Large Initial Current caused by temperature mixing in shallow water regions?

Cool Findings and Plots

Moderators: arango, robertson

Post Reply
Message
Author
tommy307507
Posts: 5
Joined: Tue Jun 01, 2021 6:11 pm
Location: CUHK

Large Initial Current caused by temperature mixing in shallow water regions?

#1 Post by tommy307507 »

I am an ocean modelling amatur and doesnt have much knowledge in running ocean models.
I am running ROMS model on a 10km resolution domain with a 10 sec timestep
I observed that my free surface behaves weirdly like this at about 2 minutes and couldn't figure out why
Image
I have low Rx0 (max 0.5) and Rx1 (only large in a few grid points, less than 5 for most points)
and the large Rx1 points are different from what's shown here.
This is what I get in my temperature after 30 mins, all of the points are along shallow water area
Image

I think this is considered an "upwelling" current as cold bottom water rushes up and replaces the warm surface water and causing my free surface to go wild...
However this also produces unrealistic SSH values up to 12m and u v velocity up to 4 ms-1 ... I am not sure if this is considered normal in ROMS but this looks wrong to me
I tried plotting the bottom u v stresses and they corrsponds to my free surface as well...

Initially I thought this is due to my grid being too steep and I tried using some unrealistic bathy modificiations , then I found out the more shallow water I had ( due to the bathy modifications ) the worse it gets so I used the LP-smoothing in Pyroms instead to produce a smooth grid (although it does take overnight)
This current makes my model very unstable at the beginning and I have to use very small timesteps to avoid the vertical courant number being > than 1, it does go down eventually

Here is a log of my test run

Code: Select all

 Model Input Parameters:  ROMS/TOMS version 3.7
                          Wednesday - September 8, 2021 - 10:16:09 AM
 -----------------------------------------------------------------------------

 AAMC

 Operating system : Linux
 CPU/hardware     : x86_64
 Compiler system  : ifort
 Compiler command : /work/intel/compilers_and_libraries_2018.5.274/linux/mpi/intel64/bin/mpiifort
 Compiler flags   : -heap-arrays -fp-model precise -ip -O3  -free -free

 Input Script  :

 SVN Root URL  : https:://myroms.org/svn/src
 SVN Revision  : Unversioned directory

 Local Root    : /avm/home/dev/COAWST/v3.2/COAWST-master
 Header Dir    : /avm/home/dev/COAWST/v3.2/COAWST-master/Projects/Hongkong
 Header file   : hongkong_roms.h
 Analytical Dir: /avm/home/dev/COAWST/v3.2/COAWST-master/Projects/Hongkong

 Resolution, Grid 01: 1248x1098x035,  Parallel Nodes:  80,  Tiling: 010x008


 Physical Parameters, Grid: 01
 =============================

      86400  ntimes            Number of timesteps for 3-D equations.
     10.000  dt                Timestep size (s) for 3-D equations.
        180  ndtfast           Number of timesteps for 2-D equations between
                                 each 3D timestep.
          1  ERstr             Starting ensemble/perturbation run number.
          1  ERend             Ending ensemble/perturbation run number.
          0  nrrec             Number of restart records to read from disk.
          F  LcycleRST         Switch to recycle time-records in restart file.
       1440  nRST              Number of timesteps between the writing of data
                                 into restart fields.
          1  ninfo             Number of timesteps between print of information
                                 to standard output.
          T  ldefout           Switch to create a new output NetCDF file(s).
         10  nHIS              Number of timesteps between the writing fields
                                 into history file.
          0  nQCK              Number of timesteps between the writing fields
                                 into quicksave file.
          1  ntsAVG            Starting timestep for the accumulation of output
                                 time-averaged data.
          0  nAVG              Number of timesteps between the writing of
                                 time-averaged data into averages file.
 0.0000E+00  nl_tnu2(01)       NLM Horizontal, harmonic mixing coefficient
                                 (m2/s) for tracer 01: temp
 0.0000E+00  nl_tnu2(02)       NLM Horizontal, harmonic mixing coefficient
                                 (m2/s) for tracer 02: salt
 5.0000E+00  nl_visc2          NLM Horizontal, harmonic mixing coefficient
                                 (m2/s) for momentum.
          F  LuvSponge         Turning OFF sponge on horizontal momentum.
          F  LtracerSponge(01) Turning OFF sponge on tracer 01: temp
          F  LtracerSponge(02) Turning OFF sponge on tracer 02: salt
 1.0000E-06  Akt_bak(01)       Background vertical mixing coefficient (m2/s)
                                 for tracer 01: temp
 1.0000E-06  Akt_bak(02)       Background vertical mixing coefficient (m2/s)
                                 for tracer 02: salt
 1.0000E-05  Akv_bak           Background vertical mixing coefficient (m2/s)
                                 for momentum.
 5.0000E-06  Akk_bak           Background vertical mixing coefficient (m2/s)
                                 for turbulent energy.
 5.0000E-06  Akp_bak           Background vertical mixing coefficient (m2/s)
                                 for turbulent generic statistical field.
 0.0000E+00  rdrg              Linear bottom drag coefficient (m/s).
 3.0000E-03  rdrg2             Quadratic bottom drag coefficient.
 0.0000E+00  Zob               Bottom roughness (m).
 1.0000E+01  blk_ZQ            Height (m) of surface air humidity measurement.
 1.0000E+01  blk_ZT            Height (m) of surface air temperature measurement.
 1.0000E+01  blk_ZW            Height (m) of surface winds measurement.
          1  lmd_Jwt           Jerlov water type.
          2  Vtransform        S-coordinate transformation equation.
          4  Vstretching       S-coordinate stretching function.
 3.0000E+00  theta_s           S-coordinate surface control parameter.
 4.0000E-01  theta_b           S-coordinate bottom  control parameter.
     10.000  Tcline            S-coordinate surface/bottom layer width (m) used
                                 in vertical coordinate stretching.
   1025.000  rho0              Mean density (kg/m3) for Boussinesq approximation.
  59396.000  dstart            Time-stamp assigned to model initialization (days).
18581117.00  time_ref          Reference time for units attribute (yyyymmdd.dd)
 0.0000E+00  Tnudg(01)         Nudging/relaxation time scale (days)
                                 for tracer 01: temp
 0.0000E+00  Tnudg(02)         Nudging/relaxation time scale (days)
                                 for tracer 02: salt
 9.0000E+01  Tnudg_SSS         Nudging/relaxation time scale (days)
                                 for sea surface salinity.
 0.0000E+00  Znudg             Nudging/relaxation time scale (days)
                                 for free-surface.
 0.0000E+00  M2nudg            Nudging/relaxation time scale (days)
                                 for 2D momentum.
 0.0000E+00  M3nudg            Nudging/relaxation time scale (days)
                                 for 3D momentum.
 0.0000E+00  obcfac            Factor between passive and active
                                 open boundary conditions.
          F  VolCons(1)        NLM western  edge boundary volume conservation.
          F  VolCons(2)        NLM southern edge boundary volume conservation.
          F  VolCons(3)        NLM eastern  edge boundary volume conservation.
          F  VolCons(4)        NLM northern edge boundary volume conservation.
     14.000  T0                Background potential temperature (C) constant.
     35.000  S0                Background salinity (PSU) constant.
      1.000  gamma2            Slipperiness variable: free-slip (1.0) or
                                                      no-slip (-1.0).
          F  LuvSrc            Turning OFF momentum point Sources/Sinks.
          F  LwSrc             Turning OFF volume influx point Sources/Sinks.
          F  LtracerSrc(01)    Turning OFF point Sources/Sinks on tracer 01: temp
          F  LtracerSrc(02)    Turning OFF point Sources/Sinks on tracer 02: salt
          F  LsshCLM           Turning OFF processing of SSH climatology.
          F  Lm2CLM            Turning OFF processing of 2D momentum climatology.
          F  Lm3CLM            Turning OFF processing of 3D momentum climatology.
          F  LtracerCLM(01)    Turning OFF processing of climatology tracer 01: temp
          F  LtracerCLM(02)    Turning OFF processing of climatology tracer 02: salt
          F  LnudgeM2CLM       Turning OFF nudging of 2D momentum climatology.
          F  LnudgeM3CLM       Turning OFF nudging of 3D momentum climatology.
          F  LnudgeTCLM(01)    Turning OFF nudging of climatology tracer 01: temp
          F  LnudgeTCLM(02)    Turning OFF nudging of climatology tracer 02: salt

          T  Hout(idFsur)      Write out free-surface.
          T  Hout(idUbar)      Write out 2D U-momentum component.
          T  Hout(idVbar)      Write out 2D V-momentum component.
          T  Hout(idUvel)      Write out 3D U-momentum component.
          T  Hout(idVvel)      Write out 3D V-momentum component.
          T  Hout(idWvel)      Write out W-momentum component.
          T  Hout(idTvar)      Write out tracer 01: temp
          T  Hout(idTvar)      Write out tracer 02: salt
          T  Hout(idUsms)      Write out surface U-momentum stress.
          T  Hout(idVsms)      Write out surface V-momentum stress.
          T  Hout(idUbms)      Write out bottom U-momentum stress.
          T  Hout(idVbms)      Write out bottom V-momentum stress.
          T  Hout(idDano)      Write out density anomaly.

 Output/Input Files:

             Output Restart File:  ocean_rst_2.nc
             Output History File:  ocean_his.nc
            Output Averages File:  ocean_avg.nc
                 Input Grid File:  Projects/Hongkong/roms_grid.nc
    Input Nonlinear Initial File:  Projects/Hongkong/roms_ini.nc
           Input Forcing File 01:  Projects/Hongkong/AAMC20210601-20210630/tide_aamc.nc
           Input Forcing File 02:  Projects/Hongkong/roms_forc.nc
         Input Boundary File 01:  Projects/Hongkong/roms_bdy.nc

 Tile partition information for Grid 01:  1248x1098x0035  tiling: 010x008

     tile     Istr     Iend     Jstr     Jend     Npts

 Number of tracers:            2
        0        1      124        1      135   585900
        1      125      249        1      135   590625
        2      250      374        1      135   590625
        3      375      499        1      135   590625
        4      500      624        1      135   590625
        5      625      749        1      135   590625
        6      750      874        1      135   590625
        7      875      999        1      135   590625
        8     1000     1124        1      135   590625
        9     1125     1248        1      135   585900
       10        1      124      136      273   598920
       11      125      249      136      273   603750
       12      250      374      136      273   603750
       13      375      499      136      273   603750
       14      500      624      136      273   603750
       15      625      749      136      273   603750
       16      750      874      136      273   603750
       17      875      999      136      273   603750
       18     1000     1124      136      273   603750
       19     1125     1248      136      273   598920
       20        1      124      274      411   598920
       21      125      249      274      411   603750
       22      250      374      274      411   603750
       23      375      499      274      411   603750
       24      500      624      274      411   603750
       25      625      749      274      411   603750
       26      750      874      274      411   603750
       27      875      999      274      411   603750
       28     1000     1124      274      411   603750
       29     1125     1248      274      411   598920
       30        1      124      412      549   598920
       31      125      249      412      549   603750
       32      250      374      412      549   603750
       33      375      499      412      549   603750
       34      500      624      412      549   603750
       35      625      749      412      549   603750
       36      750      874      412      549   603750
       37      875      999      412      549   603750
       38     1000     1124      412      549   603750
       39     1125     1248      412      549   598920
       40        1      124      550      687   598920
       41      125      249      550      687   603750
       42      250      374      550      687   603750
       43      375      499      550      687   603750
       44      500      624      550      687   603750
       45      625      749      550      687   603750
       46      750      874      550      687   603750
       47      875      999      550      687   603750
       48     1000     1124      550      687   603750
       49     1125     1248      550      687   598920
       50        1      124      688      825   598920
       51      125      249      688      825   603750
       52      250      374      688      825   603750
       53      375      499      688      825   603750
       54      500      624      688      825   603750
       55      625      749      688      825   603750
       56      750      874      688      825   603750
       57      875      999      688      825   603750
       58     1000     1124      688      825   603750
       59     1125     1248      688      825   598920
       60        1      124      826      963   598920
       61      125      249      826      963   603750
       62      250      374      826      963   603750
       63      375      499      826      963   603750
       64      500      624      826      963   603750
       65      625      749      826      963   603750
       66      750      874      826      963   603750
       67      875      999      826      963   603750
       68     1000     1124      826      963   603750
       69     1125     1248      826      963   598920
       70        1      124      964     1098   585900
       71      125      249      964     1098   590625
       72      250      374      964     1098   590625
       73      375      499      964     1098   590625
       74      500      624      964     1098   590625
       75      625      749      964     1098   590625
       76      750      874      964     1098   590625
       77      875      999      964     1098   590625
       78     1000     1124      964     1098   590625
       79     1125     1248      964     1098   585900

 Tile minimum and maximum fractional coordinates for Grid 01:
   (interior points only)

     tile     Xmin     Xmax     Ymin     Ymax     grid

        0     0.50   124.50     0.50   135.50  RHO-points
        1   124.50   249.50     0.50   135.50  RHO-points
        2   249.50   374.50     0.50   135.50  RHO-points
        3   374.50   499.50     0.50   135.50  RHO-points
        4   499.50   624.50     0.50   135.50  RHO-points
        5   624.50   749.50     0.50   135.50  RHO-points
        6   749.50   874.50     0.50   135.50  RHO-points
        7   874.50   999.50     0.50   135.50  RHO-points
        8   999.50  1124.50     0.50   135.50  RHO-points
        9  1124.50  1248.50     0.50   135.50  RHO-points
       10     0.50   124.50   135.50   273.50  RHO-points
       11   124.50   249.50   135.50   273.50  RHO-points
       12   249.50   374.50   135.50   273.50  RHO-points
       13   374.50   499.50   135.50   273.50  RHO-points
       14   499.50   624.50   135.50   273.50  RHO-points
       15   624.50   749.50   135.50   273.50  RHO-points
       16   749.50   874.50   135.50   273.50  RHO-points
       17   874.50   999.50   135.50   273.50  RHO-points
       18   999.50  1124.50   135.50   273.50  RHO-points
       19  1124.50  1248.50   135.50   273.50  RHO-points
       20     0.50   124.50   273.50   411.50  RHO-points
       21   124.50   249.50   273.50   411.50  RHO-points
       22   249.50   374.50   273.50   411.50  RHO-points
       23   374.50   499.50   273.50   411.50  RHO-points
       24   499.50   624.50   273.50   411.50  RHO-points
       25   624.50   749.50   273.50   411.50  RHO-points
       26   749.50   874.50   273.50   411.50  RHO-points
       27   874.50   999.50   273.50   411.50  RHO-points
       28   999.50  1124.50   273.50   411.50  RHO-points
       29  1124.50  1248.50   273.50   411.50  RHO-points
       30     0.50   124.50   411.50   549.50  RHO-points
       31   124.50   249.50   411.50   549.50  RHO-points
       32   249.50   374.50   411.50   549.50  RHO-points
       33   374.50   499.50   411.50   549.50  RHO-points
       34   499.50   624.50   411.50   549.50  RHO-points
       35   624.50   749.50   411.50   549.50  RHO-points
       36   749.50   874.50   411.50   549.50  RHO-points
       37   874.50   999.50   411.50   549.50  RHO-points
       38   999.50  1124.50   411.50   549.50  RHO-points
       39  1124.50  1248.50   411.50   549.50  RHO-points
       40     0.50   124.50   549.50   687.50  RHO-points
       41   124.50   249.50   549.50   687.50  RHO-points
       42   249.50   374.50   549.50   687.50  RHO-points
       43   374.50   499.50   549.50   687.50  RHO-points
       44   499.50   624.50   549.50   687.50  RHO-points
       45   624.50   749.50   549.50   687.50  RHO-points
       46   749.50   874.50   549.50   687.50  RHO-points
       47   874.50   999.50   549.50   687.50  RHO-points
       48   999.50  1124.50   549.50   687.50  RHO-points
       49  1124.50  1248.50   549.50   687.50  RHO-points
       50     0.50   124.50   687.50   825.50  RHO-points
       51   124.50   249.50   687.50   825.50  RHO-points
       52   249.50   374.50   687.50   825.50  RHO-points
       53   374.50   499.50   687.50   825.50  RHO-points
       54   499.50   624.50   687.50   825.50  RHO-points
       55   624.50   749.50   687.50   825.50  RHO-points
       56   749.50   874.50   687.50   825.50  RHO-points
       57   874.50   999.50   687.50   825.50  RHO-points
       58   999.50  1124.50   687.50   825.50  RHO-points
       59  1124.50  1248.50   687.50   825.50  RHO-points
       60     0.50   124.50   825.50   963.50  RHO-points
       61   124.50   249.50   825.50   963.50  RHO-points
       62   249.50   374.50   825.50   963.50  RHO-points
       63   374.50   499.50   825.50   963.50  RHO-points
       64   499.50   624.50   825.50   963.50  RHO-points
       65   624.50   749.50   825.50   963.50  RHO-points
       66   749.50   874.50   825.50   963.50  RHO-points
       67   874.50   999.50   825.50   963.50  RHO-points
       68   999.50  1124.50   825.50   963.50  RHO-points
       69  1124.50  1248.50   825.50   963.50  RHO-points
       70     0.50   124.50   963.50  1098.50  RHO-points
       71   124.50   249.50   963.50  1098.50  RHO-points
       72   249.50   374.50   963.50  1098.50  RHO-points
       73   374.50   499.50   963.50  1098.50  RHO-points
       74   499.50   624.50   963.50  1098.50  RHO-points
       75   624.50   749.50   963.50  1098.50  RHO-points
       76   749.50   874.50   963.50  1098.50  RHO-points
       77   874.50   999.50   963.50  1098.50  RHO-points
       78   999.50  1124.50   963.50  1098.50  RHO-points
       79  1124.50  1248.50   963.50  1098.50  RHO-points

        0     1.00   124.50     0.50   135.50    U-points
        1   124.50   249.50     0.50   135.50    U-points
        2   249.50   374.50     0.50   135.50    U-points
        3   374.50   499.50     0.50   135.50    U-points
        4   499.50   624.50     0.50   135.50    U-points
        5   624.50   749.50     0.50   135.50    U-points
        6   749.50   874.50     0.50   135.50    U-points
        7   874.50   999.50     0.50   135.50    U-points
        8   999.50  1124.50     0.50   135.50    U-points
        9  1124.50  1248.00     0.50   135.50    U-points
       10     1.00   124.50   135.50   273.50    U-points
       11   124.50   249.50   135.50   273.50    U-points
       12   249.50   374.50   135.50   273.50    U-points
       13   374.50   499.50   135.50   273.50    U-points
       14   499.50   624.50   135.50   273.50    U-points
       15   624.50   749.50   135.50   273.50    U-points
       16   749.50   874.50   135.50   273.50    U-points
       17   874.50   999.50   135.50   273.50    U-points
       18   999.50  1124.50   135.50   273.50    U-points
       19  1124.50  1248.00   135.50   273.50    U-points
       20     1.00   124.50   273.50   411.50    U-points
       21   124.50   249.50   273.50   411.50    U-points
       22   249.50   374.50   273.50   411.50    U-points
       23   374.50   499.50   273.50   411.50    U-points
       24   499.50   624.50   273.50   411.50    U-points
       25   624.50   749.50   273.50   411.50    U-points
       26   749.50   874.50   273.50   411.50    U-points
       27   874.50   999.50   273.50   411.50    U-points
       28   999.50  1124.50   273.50   411.50    U-points
       29  1124.50  1248.00   273.50   411.50    U-points
       30     1.00   124.50   411.50   549.50    U-points
       31   124.50   249.50   411.50   549.50    U-points
       32   249.50   374.50   411.50   549.50    U-points
       33   374.50   499.50   411.50   549.50    U-points
       34   499.50   624.50   411.50   549.50    U-points
       35   624.50   749.50   411.50   549.50    U-points
       36   749.50   874.50   411.50   549.50    U-points
       37   874.50   999.50   411.50   549.50    U-points
       38   999.50  1124.50   411.50   549.50    U-points
       39  1124.50  1248.00   411.50   549.50    U-points
       40     1.00   124.50   549.50   687.50    U-points
       41   124.50   249.50   549.50   687.50    U-points
       42   249.50   374.50   549.50   687.50    U-points
       43   374.50   499.50   549.50   687.50    U-points
       44   499.50   624.50   549.50   687.50    U-points
       45   624.50   749.50   549.50   687.50    U-points
       46   749.50   874.50   549.50   687.50    U-points
       47   874.50   999.50   549.50   687.50    U-points
       48   999.50  1124.50   549.50   687.50    U-points
       49  1124.50  1248.00   549.50   687.50    U-points
       50     1.00   124.50   687.50   825.50    U-points
       51   124.50   249.50   687.50   825.50    U-points
       52   249.50   374.50   687.50   825.50    U-points
       53   374.50   499.50   687.50   825.50    U-points
       54   499.50   624.50   687.50   825.50    U-points
       55   624.50   749.50   687.50   825.50    U-points
       56   749.50   874.50   687.50   825.50    U-points
       57   874.50   999.50   687.50   825.50    U-points
       58   999.50  1124.50   687.50   825.50    U-points
       59  1124.50  1248.00   687.50   825.50    U-points
       60     1.00   124.50   825.50   963.50    U-points
       61   124.50   249.50   825.50   963.50    U-points
       62   249.50   374.50   825.50   963.50    U-points
       63   374.50   499.50   825.50   963.50    U-points
       64   499.50   624.50   825.50   963.50    U-points
       65   624.50   749.50   825.50   963.50    U-points
       66   749.50   874.50   825.50   963.50    U-points
       67   874.50   999.50   825.50   963.50    U-points
       68   999.50  1124.50   825.50   963.50    U-points
       69  1124.50  1248.00   825.50   963.50    U-points
       70     1.00   124.50   963.50  1098.50    U-points
       71   124.50   249.50   963.50  1098.50    U-points
       72   249.50   374.50   963.50  1098.50    U-points
       73   374.50   499.50   963.50  1098.50    U-points
       74   499.50   624.50   963.50  1098.50    U-points
       75   624.50   749.50   963.50  1098.50    U-points
       76   749.50   874.50   963.50  1098.50    U-points
       77   874.50   999.50   963.50  1098.50    U-points
       78   999.50  1124.50   963.50  1098.50    U-points
       79  1124.50  1248.00   963.50  1098.50    U-points

        0     0.50   124.50     1.00   135.50    V-points
        1   124.50   249.50     1.00   135.50    V-points
        2   249.50   374.50     1.00   135.50    V-points
        3   374.50   499.50     1.00   135.50    V-points
        4   499.50   624.50     1.00   135.50    V-points
        5   624.50   749.50     1.00   135.50    V-points
        6   749.50   874.50     1.00   135.50    V-points
        7   874.50   999.50     1.00   135.50    V-points
        8   999.50  1124.50     1.00   135.50    V-points
        9  1124.50  1248.50     1.00   135.50    V-points
       10     0.50   124.50   135.50   273.50    V-points
       11   124.50   249.50   135.50   273.50    V-points
       12   249.50   374.50   135.50   273.50    V-points
       13   374.50   499.50   135.50   273.50    V-points
       14   499.50   624.50   135.50   273.50    V-points
 Node # 41 (pid=   77568) is active.
       15   624.50   749.50   135.50   273.50    V-points
       16   749.50   874.50   135.50   273.50    V-points
       17   874.50   999.50   135.50   273.50    V-points
       18   999.50  1124.50   135.50   273.50    V-points
       19  1124.50  1248.50   135.50   273.50    V-points
       20     0.50   124.50   273.50   411.50    V-points
       21   124.50   249.50   273.50   411.50    V-points
       22   249.50   374.50   273.50   411.50    V-points
       23   374.50   499.50   273.50   411.50    V-points
 Node # 42 (pid=   77569) is active.
       24   499.50   624.50   273.50   411.50    V-points
       25   624.50   749.50   273.50   411.50    V-points
       26   749.50   874.50   273.50   411.50    V-points
       27   874.50   999.50   273.50   411.50    V-points
       28   999.50  1124.50   273.50   411.50    V-points
       29  1124.50  1248.50   273.50   411.50    V-points
       30     0.50   124.50   411.50   549.50    V-points
       31   124.50   249.50   411.50   549.50    V-points
       32   249.50   374.50   411.50   549.50    V-points
 Node # 50 (pid=   77577) is active.
       33   374.50   499.50   411.50   549.50    V-points
       34   499.50   624.50   411.50   549.50    V-points
       35   624.50   749.50   411.50   549.50    V-points
       36   749.50   874.50   411.50   549.50    V-points
       37   874.50   999.50   411.50   549.50    V-points
       38   999.50  1124.50   411.50   549.50    V-points
       39  1124.50  1248.50   411.50   549.50    V-points
       40     0.50   124.50   549.50   687.50    V-points
       41   124.50   249.50   549.50   687.50    V-points
 Node # 51 (pid=   77578) is active.
       42   249.50   374.50   549.50   687.50    V-points
       43   374.50   499.50   549.50   687.50    V-points
       44   499.50   624.50   549.50   687.50    V-points
       45   624.50   749.50   549.50   687.50    V-points
       46   749.50   874.50   549.50   687.50    V-points
       47   874.50   999.50   549.50   687.50    V-points
       48   999.50  1124.50   549.50   687.50    V-points
       49  1124.50  1248.50   549.50   687.50    V-points
       50     0.50   124.50   687.50   825.50    V-points
 Node # 55 (pid=   77582) is active.
       51   124.50   249.50   687.50   825.50    V-points
       52   249.50   374.50   687.50   825.50    V-points
       53   374.50   499.50   687.50   825.50    V-points
       54   499.50   624.50   687.50   825.50    V-points
       55   624.50   749.50   687.50   825.50    V-points
       56   749.50   874.50   687.50   825.50    V-points
       57   874.50   999.50   687.50   825.50    V-points
       58   999.50  1124.50   687.50   825.50    V-points
       59  1124.50  1248.50   687.50   825.50    V-points
 Node # 59 (pid=   77586) is active.
       60     0.50   124.50   825.50   963.50    V-points
       61   124.50   249.50   825.50   963.50    V-points
       62   249.50   374.50   825.50   963.50    V-points
       63   374.50   499.50   825.50   963.50    V-points
       64   499.50   624.50   825.50   963.50    V-points
       65   624.50   749.50   825.50   963.50    V-points
       66   749.50   874.50   825.50   963.50    V-points
       67   874.50   999.50   825.50   963.50    V-points
       68   999.50  1124.50   825.50   963.50    V-points
 Node # 69 (pid=   77596) is active.
       69  1124.50  1248.50   825.50   963.50    V-points
       70     0.50   124.50   963.50  1098.00    V-points
       71   124.50   249.50   963.50  1098.00    V-points
       72   249.50   374.50   963.50  1098.00    V-points
       73   374.50   499.50   963.50  1098.00    V-points
       74   499.50   624.50   963.50  1098.00    V-points
       75   624.50   749.50   963.50  1098.00    V-points
       76   749.50   874.50   963.50  1098.00    V-points
       77   874.50   999.50   963.50  1098.00    V-points
 Node # 74 (pid=   77601) is active.
       78   999.50  1124.50   963.50  1098.00    V-points
       79  1124.50  1248.50   963.50  1098.00    V-points

 Maximum halo size in XI and ETA directions:

               HaloSizeI(1) =     417
               HaloSizeJ(1) =     456
                TileSide(1) =     146
                TileSize(1) =   19418


 Lateral Boundary Conditions: NLM
 ============================

 Variable               Grid  West Edge    South Edge   East Edge    North Edge
 ---------              ----  ----------   ----------   ----------   ----------

 zeta                     1   Chapman Imp  Chapman Imp  Chapman Imp  Chapman Imp

 ubar                     1   Flather      Flather      Flather      Flather
 Node # 40 (pid=   77567) is active.

 vbar                     1   Flather      Flather      Flather      Flather

 u                        1   Rad + Nud    Rad + Nud    Rad + Nud    Rad + Nud

 v                        1   Rad + Nud    Rad + Nud    Rad + Nud    Rad + Nud

 temp                     1   Rad + Nud    Rad + Nud    Rad + Nud    Rad + Nud

 salt                     1   Rad + Nud    Rad + Nud    Rad + Nud    Rad + Nud

 tke                      1   Radiation    Radiation    Radiation    Radiation

 Activated C-preprocessing Options:

 HONGKONG            AAMC
 ANA_BSFLUX          Analytical kinematic bottom salinity flux.
 ANA_BTFLUX          Analytical kinematic bottom temperature flux.
 ASSUMED_SHAPE       Using assumed-shape arrays.
 ATM_PRESS           Impose atmospheric pressure onto sea surface.
 AVERAGES            Writing out time-averaged nonlinear model fields.
 BULK_FLUXES         Surface bulk fluxes parameterization.
 DJ_GRADPS           Parabolic Splines density Jacobian (Shchepetkin, 2002).
 DOUBLE_PRECISION    Double precision arithmetic.
 EMINUSP             Compute Salt Flux using E-P.
 KANTHA_CLAYSON      Kantha and Clayson stability function formulation.
 MASKING             Land/Sea masking.
 Node # 44 (pid=   77571) is active.
 MIX_ISO_TS          Mixing of tracers along isopycnal surfaces.
 MIX_S_UV            Mixing of momentum along constant S-surfaces.
 MPI                 MPI distributed-memory configuration.
 MY25_MIXING         Mellor/Yamada Level-2.5 mixing closure.
 NONLINEAR           Nonlinear Model.
 NONLIN_EOS          Nonlinear Equation of State for seawater.
 NO_LBC_ATT          Not checking NetCDF global attribute NLM_LBC during restart.
 N2S2_HORAVG         Horizontal smoothing of buoyancy and shear.
 POWER_LAW           Power-law shape time-averaging barotropic filter.
 PROFILE             Time profiling activated .
 K_C4ADVECTION       Fourth-order centered differences advection of TKE fields.
 Node # 45 (pid=   77572) is active.
 RADIATION_2D        Use tangential phase speed in radiation conditions.
 READ_WATER          Reading data at water points only.
 RI_SPLINES          Parabolic Spline Reconstruction for Richardson Number.
 !RST_SINGLE         Double precision fields in restart NetCDF file.
 SALINITY            Using salinity.
 SOLAR_SOURCE        Solar Radiation Source Term.
 SOLVE3D             Solving 3D Primitive Equations.
 TS_MPDATA           Recursive flux corrected MPDATA 3D advection of tracers.
 TS_DIF2             Harmonic mixing of tracers.
 UV_ADV              Advection of momentum.
 UV_COR              Coriolis term.
 UV_U3HADVECTION     Third-order upstream horizontal advection of 3D momentum.
 UV_C4VADVECTION     Fourth-order centered vertical advection of momentum.
 Node # 46 (pid=   77573) is active.
 UV_QDRAG            Quadratic bottom stress.
 UV_VIS2             Harmonic mixing of momentum.
 UV_SMAGORINSKY      Smagorinksy-like time-dependent viscosity coefficients.
 VAR_RHO_2D          Variable density barotropic mode.
 VISC_3DCOEF         Horizontal, time-dependent 3D viscosity coefficient.

 CHECKADJ - use caution when activating: UV_SMAGORINSKY
            REASON: stability problems, WARNING.

 
 INITIAL: Configuring and initializing forward nonlinear model ...
 *******

 Vertical S-coordinate System, Grid 01:

 level   S-coord     Cs-curve   Z   at hmin       at hc    half way     at hmax

    35   0.0000000   0.0000000        0.000       0.000       0.000       0.000
    34  -0.0285714  -0.0004918       -0.145      -0.145      -2.050      -3.815
    33  -0.0571429  -0.0019703       -0.296      -0.296      -7.641     -14.714
    32  -0.0857143  -0.0044449       -0.451      -0.451     -16.808     -32.763
    31  -0.1142857  -0.0079314       -0.611      -0.611     -29.607     -58.075
    30  -0.1428571  -0.0124518       -0.777      -0.777     -46.116     -90.808
    29  -0.1714286  -0.0180347       -0.947      -0.947     -66.438    -131.169
    28  -0.2000000  -0.0247149       -1.124      -1.124     -90.699    -179.406
    27  -0.2285714  -0.0325342       -1.306      -1.306    -119.049    -235.820
    26  -0.2571429  -0.0415407       -1.493      -1.493    -151.659    -300.756
    25  -0.2857143  -0.0517896       -1.688      -1.688    -188.729    -374.610
    24  -0.3142857  -0.0633428       -1.888      -1.888    -230.480    -457.827
    23  -0.3428571  -0.0762693       -2.096      -2.096    -277.159    -550.901
    22  -0.3714286  -0.0906450       -2.310      -2.310    -329.041    -654.378
    21  -0.4000000  -0.1065532       -2.533      -2.533    -386.422    -768.856
    20  -0.4285714  -0.1240840       -2.763      -2.763    -449.627    -894.981
    19  -0.4571429  -0.1433350       -3.002      -3.002    -519.006   -1033.454
    18  -0.4857143  -0.1644105       -3.251      -3.251    -594.934   -1185.024
    17  -0.5142857  -0.1874218       -3.509      -3.509    -677.809   -1350.489
    16  -0.5428571  -0.2124868       -3.777      -3.777    -768.056   -1530.696
    15  -0.5714286  -0.2397294       -4.056      -4.056    -866.118   -1726.535
    14  -0.6000000  -0.2692794       -4.346      -4.346    -972.461   -1938.937
    13  -0.6285714  -0.3012713       -4.649      -4.649   -1087.569   -2168.867
    12  -0.6571429  -0.3358436       -4.965      -4.965   -1211.939   -2417.319
    11  -0.6857143  -0.3731377       -5.294      -5.294   -1346.077   -2685.309
    10  -0.7142857  -0.4132962       -5.638      -5.638   -1490.496   -2973.861
     9  -0.7428571  -0.4564613       -5.997      -5.997   -1645.706   -3283.995
     8  -0.7714286  -0.5027728       -6.371      -6.371   -1812.209   -3616.714
     7  -0.8000000  -0.5523651       -6.762      -6.762   -1990.486   -3972.983
     6  -0.8285714  -0.6053646       -7.170      -7.170   -2180.994   -4353.711
     5  -0.8571429  -0.6618863       -7.595      -7.595   -2384.142   -4759.721
     4  -0.8857143  -0.7220296       -8.039      -8.039   -2600.289   -5191.727
     3  -0.9142857  -0.7858738       -8.501      -8.501   -2829.719   -5650.300
     2  -0.9428571  -0.8534736       -8.982      -8.982   -3072.628   -6135.831
     1  -0.9714286  -0.9248530       -9.481      -9.481   -3329.103   -6648.493
     0  -1.0000000  -1.0000000      -10.000     -10.000   -3599.100   -7188.199


 ndtfast, nfast =  180 241   nfast/ndtfast =  1.33889

 Centers of gravity and integrals (values must be 1, 1, approx 1/2, 1, 1):

    1.000000000000 1.011547989092 0.505773994546 1.000000000000 1.000000000000

 Power filter parameters, Fgamma, gamma =  0.28400   0.26822

 Metrics information for Grid 01:
 ===============================

 Minimum X-grid spacing, DXmin =  1.08226339E+01 km
 Maximum X-grid spacing, DXmax =  1.08234832E+01 km
 Minimum Y-grid spacing, DYmin =  1.07500622E+01 km
 Maximum Y-grid spacing, DYmax =  1.08209760E+01 km
 Minimum Z-grid spacing, DZmin =  1.45316090E-01 m
 Maximum Z-grid spacing, DZmax =  5.39706365E+02 m

 Minimum barotropic Courant Number =  7.19051938E-05
 Maximum barotropic Courant Number =  1.93400940E-03
 Maximum Coriolis   Courant Number =  1.27985656E-03


 Minimum horizontal viscosity coefficient =  8.98857253E+01 m2/s
 Maximum horizontal viscosity coefficient =  9.01852461E+01 m2/s

 NLM: GET_STATE - Read state initial conditions,             t = 59448 00:00:00
                   (Grid 01, File: roms_ini.nc, Rec=0001, Index=1)
                - free-surface
                   (Min = -5.23771516E-01 Max =  1.68767574E+00)
                - vertically integrated u-momentum component
                   (Min = -5.88693481E-01 Max =  8.93720017E-01)
                - vertically integrated v-momentum component
                   (Min = -8.99220987E-01 Max =  8.94592932E-01)
                - u-momentum component
                   (Min = -1.64256394E+00 Max =  2.86598516E+00)
                - v-momentum component
                   (Min = -2.28817463E+00 Max =  1.81433034E+00)
                - potential temperature
                   (Min = -1.89512908E+00 Max =  3.59746475E+01)
                - salinity
                   (Min =  0.00000000E+00 Max =  4.39164734E+01)
    GET_2DFLD   - surface u-wind component,                  t = 59448 00:00:00
                   (Rec=0000001, Index=1, File: roms_forc.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -2.12692261E+01 Max =  1.47474135E+01)
    GET_2DFLD   - surface v-wind component,                  t = 59448 00:00:00
                   (Rec=0000001, Index=1, File: roms_forc.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -1.65383660E+01 Max =  2.22907516E+01)
    GET_2DFLD   - surface air pressure,                      t = 59448 00:00:00
                   (Rec=0000001, Index=1, File: roms_forc.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  5.02593040E+02 Max =  1.02620879E+03)
    GET_2DFLD   - surface air temperature,                   t = 59448 00:00:00
                   (Rec=0000001, Index=1, File: roms_forc.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -1.72279855E+01 Max =  4.12276120E+01)
    GET_2DFLD   - surface air relative humidity,             t = 59448 00:00:00
                   (Rec=0000001, Index=1, File: roms_forc.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  8.67188807E-02 Max =  9.99982713E-01)
    GET_2DFLD   - solar shortwave radiation flux,            t = 59448 00:00:00
                   (Rec=0000001, Index=1, File: roms_forc.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  0.00000000E+00 Max =  2.26388896E-04)
    GET_2DFLD   - net longwave radiation flux,               t = 59448 00:00:00
                   (Rec=0000001, Index=1, File: roms_forc.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -6.33245609E-05 Max =  2.18022298E-06)
    GET_2DFLD   - rain fall rate,                            t = 59448 00:00:00
                   (Rec=0000001, Index=1, File: roms_forc.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  0.00000000E+00 Max =  4.18852894E-03)
    GET_NGFLD   - free-surface western boundary condition,   t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  0.00000000E+00 Max =  9.18808910E-01)
    GET_NGFLD   - free-surface eastern boundary condition,   t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -3.72236876E-01 Max =  9.85347354E-01)
    GET_NGFLD   - free-surface southern boundary condition,  t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  0.00000000E+00 Max =  9.30096322E-01)
    GET_NGFLD   - free-surface northern boundary condition,  t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -1.54842855E-01 Max =  0.00000000E+00)
    GET_NGFLD   - 2D u-momentum western boundary condition,  t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -1.93881100E-01 Max =  3.98690294E-01)
    GET_NGFLD   - 2D v-momentum western boundary condition,  t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -1.24536020E-01 Max =  5.05757189E-01)
    GET_NGFLD   - 2D u-momentum eastern boundary condition,  t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -2.22985159E-01 Max =  1.69643588E-01)
    GET_NGFLD   - 2D v-momentum eastern boundary condition,  t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -2.49782698E-01 Max =  2.32015508E-01)
    GET_NGFLD   - 2D u-momentum southern boundary condition, t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -1.27306205E-01 Max =  1.45040155E-01)
    GET_NGFLD   - 2D v-momentum southern boundary condition, t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -3.59970067E-01 Max =  1.85509109E-01)
    GET_NGFLD   - 2D u-momentum northern boundary condition, t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -8.55978581E-02 Max =  3.09511009E-02)
    GET_NGFLD   - 2D v-momentum northern boundary condition, t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -1.53482416E-01 Max =  5.01764663E-02)
    GET_NGFLD   - 3D u-momentum western boundary condition,  t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -8.56141626E-01 Max =  7.60081740E-01)
    GET_NGFLD   - 3D v-momentum western boundary condition,  t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -8.91905015E-01 Max =  9.14233817E-01)
    GET_NGFLD   - 3D u-momentum eastern boundary condition,  t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -1.11849057E+00 Max =  8.93760105E-01)
    GET_NGFLD   - 3D v-momentum eastern boundary condition,  t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -1.23631883E+00 Max =  5.82368326E-01)
    GET_NGFLD   - 3D u-momentum southern boundary condition, t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -5.19719120E-01 Max =  7.16236399E-01)
    GET_NGFLD   - 3D v-momentum southern boundary condition, t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -7.65520368E-01 Max =  7.96315578E-01)
    GET_NGFLD   - 3D u-momentum northern boundary condition, t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -1.61467013E-01 Max =  9.47605396E-02)
    GET_NGFLD   - 3D v-momentum northern boundary condition, t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -1.91182274E-01 Max =  9.12230464E-02)
    GET_NGFLD   - temperature western boundary condition,    t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  0.00000000E+00 Max =  3.38317888E+01)
    GET_NGFLD   - salinity western boundary condition,       t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  0.00000000E+00 Max =  4.05570047E+01)
    GET_NGFLD   - temperature eastern boundary condition,    t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  0.00000000E+00 Max =  3.00975357E+01)
    GET_NGFLD   - salinity eastern boundary condition,       t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  0.00000000E+00 Max =  3.59606162E+01)
    GET_NGFLD   - temperature southern boundary condition,   t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  0.00000000E+00 Max =  2.28045448E+01)
    GET_NGFLD   - salinity southern boundary condition,      t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  0.00000000E+00 Max =  3.60215520E+01)
    GET_NGFLD   - temperature northern boundary condition,   t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  0.00000000E+00 Max =  1.16957574E+01)
    GET_NGFLD   - salinity northern boundary condition,      t = 59448 00:00:00
                   (Rec=0000001, Index=2, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  0.00000000E+00 Max =  3.27917837E+01)

 Basin information for Grid 01:

 Maximum grid stiffness ratios:  rx0 =   5.189205E-01 (Beckmann and Haidvogel)
                                 rx1 =   1.641224E+01 (Haney)

 Initial basin volumes: TotVolume =  2.7396633652E+17 m3
                        MinVolume =  1.6906680057E+07 m3
                        MaxVolume =  6.2809235998E+10 m3
                          Max/Min =  3.7150543919E+03

 NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 00449281 - 00086400)

    GET_2DFLD   - surface u-wind component,                  t = 59448 06:00:00
                   (Rec=0000002, Index=2, File: roms_forc.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -1.62967982E+01 Max =  1.49250742E+01)
    GET_2DFLD   - surface v-wind component,                  t = 59448 06:00:00
                   (Rec=0000002, Index=2, File: roms_forc.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -1.42649994E+01 Max =  2.24804682E+01)
    GET_2DFLD   - surface air pressure,                      t = 59448 06:00:00
                   (Rec=0000002, Index=2, File: roms_forc.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  5.03315126E+02 Max =  1.02483943E+03)
    GET_2DFLD   - surface air temperature,                   t = 59448 06:00:00
                   (Rec=0000002, Index=2, File: roms_forc.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  9.44074498E-02 Max =  6.65694407E+01)
    GET_2DFLD   - surface air relative humidity,             t = 59448 06:00:00
                   (Rec=0000002, Index=2, File: roms_forc.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  3.59224866E-02 Max =  1.00000000E+00)
    GET_2DFLD   - solar shortwave radiation flux,            t = 59448 06:00:00
                   (Rec=0000002, Index=2, File: roms_forc.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  5.05144004E-07 Max =  2.26292645E-04)
    GET_2DFLD   - net longwave radiation flux,               t = 59448 06:00:00
                   (Rec=0000002, Index=2, File: roms_forc.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -8.55175232E-05 Max =  1.92112509E-06)
    GET_2DFLD   - rain fall rate,                            t = 59448 06:00:00
                   (Rec=0000002, Index=2, File: roms_forc.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  0.00000000E+00 Max =  3.67907220E-03)
    GET_NGFLD   - free-surface western boundary condition,   t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  0.00000000E+00 Max =  9.13595725E-01)
    GET_NGFLD   - free-surface eastern boundary condition,   t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -3.80341571E-01 Max =  1.00688323E+00)
    GET_NGFLD   - free-surface southern boundary condition,  t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  0.00000000E+00 Max =  9.21159827E-01)
    GET_NGFLD   - free-surface northern boundary condition,  t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -1.47869733E-01 Max =  0.00000000E+00)
    GET_NGFLD   - 2D u-momentum western boundary condition,  t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -1.94615663E-01 Max =  3.86651699E-01)
    GET_NGFLD   - 2D v-momentum western boundary condition,  t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -1.23470110E-01 Max =  5.13688864E-01)
    GET_NGFLD   - 2D u-momentum eastern boundary condition,  t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -2.25827859E-01 Max =  1.70569687E-01)
    GET_NGFLD   - 2D v-momentum eastern boundary condition,  t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -2.53845765E-01 Max =  2.28138624E-01)
    GET_NGFLD   - 2D u-momentum southern boundary condition, t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -1.28899920E-01 Max =  1.42907998E-01)
    GET_NGFLD   - 2D v-momentum southern boundary condition, t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -3.63704986E-01 Max =  1.85439025E-01)
    GET_NGFLD   - 2D u-momentum northern boundary condition, t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -7.87085736E-02 Max =  3.78793816E-02)
    GET_NGFLD   - 2D v-momentum northern boundary condition, t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -1.53680235E-01 Max =  5.56665103E-02)
    GET_NGFLD   - 3D u-momentum western boundary condition,  t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -8.73040843E-01 Max =  7.62176647E-01)
    GET_NGFLD   - 3D v-momentum western boundary condition,  t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -9.38988573E-01 Max =  8.87503651E-01)
    GET_NGFLD   - 3D u-momentum eastern boundary condition,  t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -1.25736728E+00 Max =  8.33861637E-01)
    GET_NGFLD   - 3D v-momentum eastern boundary condition,  t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -1.07589611E+00 Max =  5.47690883E-01)
    GET_NGFLD   - 3D u-momentum southern boundary condition, t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -5.72394462E-01 Max =  8.20241750E-01)
    GET_NGFLD   - 3D v-momentum southern boundary condition, t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -7.96319550E-01 Max =  7.89201078E-01)
    GET_NGFLD   - 3D u-momentum northern boundary condition, t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -1.60715406E-01 Max =  8.72553923E-02)
    GET_NGFLD   - 3D v-momentum northern boundary condition, t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min = -1.97438264E-01 Max =  8.77242500E-02)
    GET_NGFLD   - temperature western boundary condition,    t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  0.00000000E+00 Max =  3.37482743E+01)
    GET_NGFLD   - salinity western boundary condition,       t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  0.00000000E+00 Max =  4.05570047E+01)
    GET_NGFLD   - temperature eastern boundary condition,    t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  0.00000000E+00 Max =  3.02205660E+01)
    GET_NGFLD   - salinity eastern boundary condition,       t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  0.00000000E+00 Max =  3.59577041E+01)
    GET_NGFLD   - temperature southern boundary condition,   t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  0.00000000E+00 Max =  2.31599658E+01)
    GET_NGFLD   - salinity southern boundary condition,      t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  0.00000000E+00 Max =  3.60048515E+01)
    GET_NGFLD   - temperature northern boundary condition,   t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  0.00000000E+00 Max =  1.17224400E+01)
    GET_NGFLD   - salinity northern boundary condition,      t = 59448 06:00:00
                   (Rec=0000002, Index=1, File: roms_bdy.nc)
                   (Tmin=      59448.0000 Tmax=      59451.2500)
                   (Min =  0.00000000E+00 Max =  3.27933619E+01)

   STEP   Day HH:MM:SS  KINETIC_ENRG   POTEN_ENRG    TOTAL_ENRG    NET_VOLUME
          C => (i,j,k)       Cu            Cv            Cw         Max Speed

  449280 59448 00:00:00  6.330212E-03  2.194455E+04  2.194456E+04  2.747863E+17
        (0139,0370,34)  1.823549E-03  1.992562E-03  0.000000E+00  2.914745E+00
      DEF_HIS     - creating  history      file, Grid 01: ocean_his.nc
      WRT_HIS     - wrote history     fields (Index=1,1) in record = 0000001
  449281 59448 00:00:10  6.330793E-03  2.194455E+04  2.194456E+04  2.747863E+17
        (0858,0349,24)  6.361238E-04  2.791231E-04  3.430988E-03  2.914673E+00
  449282 59448 00:00:20  6.332671E-03  2.194455E+04  2.194456E+04  2.747863E+17
        (1080,0197,35)  2.285632E-04  3.417999E-04  5.214773E-01  2.914633E+00
  449283 59448 00:00:30  6.335590E-03  2.194455E+04  2.194456E+04  2.747863E+17
        (0865,0399,35)  4.783167E-05  1.043398E-05  6.417421E-01  2.914633E+00
  449284 59448 00:00:40  6.339342E-03  2.194455E+04  2.194456E+04  2.747863E+17
        (0953,0375,35)  5.513642E-05  4.775511E-05  8.049609E-01  2.914690E+00
  449285 59448 00:00:50  6.343717E-03  2.194455E+04  2.194456E+04  2.747863E+17
        (0953,0375,35)  6.610631E-05  3.569643E-05  9.157686E-01  2.914818E+00
How do I resolve this problem or is there anyway I can deal with this?

tommy307507
Posts: 5
Joined: Tue Jun 01, 2021 6:11 pm
Location: CUHK

Re: Large Initial Current caused by temperature mixing in shallow water regions?

#2 Post by tommy307507 »

Update ... I tried the "HPGE" test by setting all forcing except temp and Pair to zero, defining TS_FIXED , letting initial zeta, u ,v be 0 and found out the current will still happen, this is the velocity after a few minutes.. Does that mean this is purely due to my grid being not smooth and I should do more smoothing around those areas? But their Rx0 and Rx1 are already small, Rx0 max is 0.1 (darker area in pic) , Rx1 max is around 25 in a few grid points
Image
Image
For example ,Phillipines is around the worse I was getting with my zeta but the Rx1 is 2.57 around the land boundaries, which is not high at all according to all sources I read in the forum. Maybe my HPGE test has a few problems as I also see the effect of the temperature field on my u and v around Japan ... But over all my zeta field and u v field has the exact artifacts as described with the first post
Image

Image

Image
I am really confused by this... is this because of the temperature and salinity field?

tommy307507
Posts: 5
Joined: Tue Jun 01, 2021 6:11 pm
Location: CUHK

Re: Large Initial Current caused by temperature mixing in shallow water regions?

#3 Post by tommy307507 »

Update again , after doing some reading I realized that this is due to the initial temp and salinity fields that I inputted, the actual HPGE test requires a initially stable analytic T,S field.
So this is purely due to the density gradient and can be called the upwelling currents?
it causes some unrealistic zeta and u v numbers and bounces waves around the whole domain. is there anyway to change the initial values to fix this issue , or should I take a 'stable' version of the T,S field from running the model for a long time and then restarting it on those profiles to get a better initial condition?

tommy307507
Posts: 5
Joined: Tue Jun 01, 2021 6:11 pm
Location: CUHK

Re: Large Initial Current caused by temperature mixing in shallow water regions?

#4 Post by tommy307507 »

Final Update
I think all the trouble is because of the LP-Smoothing method that I used, although it preserved the overall features and the smothed my Rx0 and Rx1 correctly , I forgot to take in to account for the difference bewteen my finished bathy vs the bathy of the Hycom model and did some bad vertical extrapolation. Some shallow water essentially became deep ocean after my smoothing and they all corrosponds to the area that has the high gravity waves initially.... I might do some more digging in my pyroms code to make sure that the smoothed bathy isnt too far from the model input .... this explains the huge change in temperature after running the test case for 24 hours .. I guess the test case eventually did it job for me anyways

Post Reply