Lake Jersey Nesting test case

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
Honolulu_Wang
Posts: 17
Joined: Mon Sep 08, 2014 4:17 pm
Location: Hohai University

Lake Jersey Nesting test case

#1 Unread post by Honolulu_Wang »

Hi I am testing the lake Jersey nesting case. However, results did show same as in the.
It is a constant temperature all the time.It might due to surface wind setting. May I know how to fix that by the analytical setting? Here is the log.file and also the my results and roms website results.
Thanks.

Code: Select all

--------------------------------------------------------------------------------
 Model Input Parameters:  ROMS/TOMS version 3.7  
                          Saturday - April 13, 2019 -  6:01:01 PM
 --------------------------------------------------------------------------------

 Lake Jersey nesting test case , 4 grids: A, C, D, and E

 Operating system : Linux
 CPU/hardware     : x86_64
 Compiler system  : ifort
 Compiler command : /data/users/htimy/mpichinst/bin/mpif90
 Compiler flags   : -fp-model precise -ip -O3 
 MPI Communicator : 1140850688  PET size = 36

 Input Script  : ocean_lake_jersey_nested_4g_acde.in

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

 Local Root    : /data/users/htimy/new_jeysey_lake/trunk
 Header Dir    : /data/users/htimy/new_jeysey_lake/trunk
 Header file   : lake_jersey.h
 Analytical Dir: /data/users/htimy/new_jeysey_lake/trunk

 Resolution, Grid 01: 0100x0080x008,  Parallel Nodes:  36,  Tiling: 006x006


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

       1920  ntimes            Number of timesteps for 3-D equations.
     90.000  dt                Timestep size (s) for 3-D equations.
         30  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.
          T  LcycleRST         Switch to recycle time-records in restart file.
         20  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).
         20  nHIS              Number of timesteps between the writing fields
                                 into history file.
          0  nQCK              Number of timesteps between the writing fields
                                 into quicksave file.
          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
 5.0000E-06  Akt_bak(01)       Background vertical mixing coefficient (m2/s)
                                 for tracer 01: temp
 5.0000E-06  Akt_bak(02)       Background vertical mixing coefficient (m2/s)
                                 for tracer 02: salt
 5.0000E-05  Akv_bak           Background vertical mixing coefficient (m2/s)
                                 for momentum.
 3.0000E-04  rdrg              Linear bottom drag coefficient (m/s).
 3.0000E-03  rdrg2             Quadratic bottom drag coefficient.
 1.5000E-02  Zob               Bottom roughness (m).
          1  Vtransform        S-coordinate transformation equation.
          1  Vstretching       S-coordinate stretching function.
 1.0000E+00  theta_s           S-coordinate surface control parameter.
 1.0000E+00  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.
      0.000  dstart            Time-stamp assigned to model initialization (days).
       0.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
 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.
     10.000  T0                Background potential temperature (C) constant.
     30.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(idOvel)      Write out omega vertical velocity.
          T  Hout(idTvar)      Write out tracer 01: temp
          T  Hout(idTvar)      Write out tracer 02: salt

 Output/Input Files:

             Output Restart File:  ../output/lake_jersey_rst_a.nc
             Output History File:  ../output/lake_jersey_his_a.nc
                 Input Grid File:  ../input/lake_jersey_grd_a.nc
  Nesting grid connectivity File:  ../input/lake_jersey_ngc_4g_acde.nc
    Input Nonlinear Initial File:  ../input/lake_jersey_ini_a.nc

 Resolution, Grid 02: 0120x0130x008,  Parallel Nodes:  36,  Tiling: 006x006


 Physical Parameters, Grid: 02
 =============================

       9600  ntimes            Number of timesteps for 3-D equations.
     18.000  dt                Timestep size (s) for 3-D equations.
         30  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.
          T  LcycleRST         Switch to recycle time-records in restart file.
        100  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).
        100  nHIS              Number of timesteps between the writing fields
                                 into history file.
          0  nQCK              Number of timesteps between the writing fields
                                 into quicksave file.
          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
 5.0000E-06  Akt_bak(01)       Background vertical mixing coefficient (m2/s)
                                 for tracer 01: temp
 5.0000E-06  Akt_bak(02)       Background vertical mixing coefficient (m2/s)
                                 for tracer 02: salt
 5.0000E-05  Akv_bak           Background vertical mixing coefficient (m2/s)
                                 for momentum.
 3.0000E-04  rdrg              Linear bottom drag coefficient (m/s).
 3.0000E-03  rdrg2             Quadratic bottom drag coefficient.
 1.5000E-02  Zob               Bottom roughness (m).
          1  Vtransform        S-coordinate transformation equation.
          1  Vstretching       S-coordinate stretching function.
 1.0000E+00  theta_s           S-coordinate surface control parameter.
 1.0000E+00  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.
      0.000  dstart            Time-stamp assigned to model initialization (days).
       0.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
 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.
     10.000  T0                Background potential temperature (C) constant.
     30.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(idOvel)      Write out omega vertical velocity.
          T  Hout(idTvar)      Write out tracer 01: temp
          T  Hout(idTvar)      Write out tracer 02: salt

 Output/Input Files:

             Output Restart File:  ../output/lake_jersey_rst_c.nc
             Output History File:  ../output/lake_jersey_his_c.nc
                 Input Grid File:  ../input/lake_jersey_grd_c.nc
  Nesting grid connectivity File:  ../input/lake_jersey_ngc_4g_acde.nc
    Input Nonlinear Initial File:  ../input/lake_jersey_ini_c.nc

 Resolution, Grid 03: 0114x0060x008,  Parallel Nodes:  36,  Tiling: 006x006


 Physical Parameters, Grid: 03
 =============================

       5760  ntimes            Number of timesteps for 3-D equations.
     30.000  dt                Timestep size (s) for 3-D equations.
         30  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.
          T  LcycleRST         Switch to recycle time-records in restart file.
         60  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).
         60  nHIS              Number of timesteps between the writing fields
                                 into history file.
          0  nQCK              Number of timesteps between the writing fields
                                 into quicksave file.
          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
 5.0000E-06  Akt_bak(01)       Background vertical mixing coefficient (m2/s)
                                 for tracer 01: temp
 5.0000E-06  Akt_bak(02)       Background vertical mixing coefficient (m2/s)
                                 for tracer 02: salt
 5.0000E-05  Akv_bak           Background vertical mixing coefficient (m2/s)
                                 for momentum.
 3.0000E-04  rdrg              Linear bottom drag coefficient (m/s).
 3.0000E-03  rdrg2             Quadratic bottom drag coefficient.
 1.5000E-02  Zob               Bottom roughness (m).
          1  Vtransform        S-coordinate transformation equation.
          1  Vstretching       S-coordinate stretching function.
 1.0000E+00  theta_s           S-coordinate surface control parameter.
 1.0000E+00  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.
      0.000  dstart            Time-stamp assigned to model initialization (days).
       0.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
 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.
     10.000  T0                Background potential temperature (C) constant.
     30.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(idOvel)      Write out omega vertical velocity.
          T  Hout(idTvar)      Write out tracer 01: temp
          T  Hout(idTvar)      Write out tracer 02: salt

 Output/Input Files:

             Output Restart File:  ../output/lake_jersey_rst_d.nc
             Output History File:  ../output/lake_jersey_his_d.nc
                 Input Grid File:  ../input/lake_jersey_grd_d.nc
  Nesting grid connectivity File:  ../input/lake_jersey_ngc_4g_acde.nc
    Input Nonlinear Initial File:  ../input/lake_jersey_ini_d.nc

 Resolution, Grid 04: 0132x0060x008,  Parallel Nodes:  36,  Tiling: 006x006


 Physical Parameters, Grid: 04
 =============================

      17280  ntimes            Number of timesteps for 3-D equations.
     10.000  dt                Timestep size (s) for 3-D equations.
         30  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.
          T  LcycleRST         Switch to recycle time-records in restart file.
        180  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).
        180  nHIS              Number of timesteps between the writing fields
                                 into history file.
          0  nQCK              Number of timesteps between the writing fields
                                 into quicksave file.
          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
 5.0000E-06  Akt_bak(01)       Background vertical mixing coefficient (m2/s)
                                 for tracer 01: temp
 5.0000E-06  Akt_bak(02)       Background vertical mixing coefficient (m2/s)
                                 for tracer 02: salt
 5.0000E-05  Akv_bak           Background vertical mixing coefficient (m2/s)
                                 for momentum.
 3.0000E-04  rdrg              Linear bottom drag coefficient (m/s).
 3.0000E-03  rdrg2             Quadratic bottom drag coefficient.
 1.5000E-02  Zob               Bottom roughness (m).
          1  Vtransform        S-coordinate transformation equation.
          1  Vstretching       S-coordinate stretching function.
 1.0000E+00  theta_s           S-coordinate surface control parameter.
 1.0000E+00  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.
      0.000  dstart            Time-stamp assigned to model initialization (days).
       0.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
 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.
     10.000  T0                Background potential temperature (C) constant.
     30.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(idOvel)      Write out omega vertical velocity.
          T  Hout(idTvar)      Write out tracer 01: temp
          T  Hout(idTvar)      Write out tracer 02: salt

 Output/Input Files:

             Output Restart File:  ../output/lake_jersey_rst_e.nc
             Output History File:  ../output/lake_jersey_his_e.nc
                 Input Grid File:  ../input/lake_jersey_grd_e.nc
  Nesting grid connectivity File:  ../input/lake_jersey_ngc_4g_acde.nc
    Input Nonlinear Initial File:  ../input/lake_jersey_ini_e.nc

 



 Maximum halo size in XI and ETA directions:

               HaloSizeI(3) =     102
               HaloSizeJ(3) =      75
                TileSide(3) =      28
                TileSize(3) =     532


 Maximum halo size in XI and ETA directions:

               HaloSizeI(4) =     111
               HaloSizeJ(4) =      75
                TileSide(4) =      31
                TileSize(4) =     589


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

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

 zeta                     1   Closed       Closed       Closed       Closed
                          2   Nested       Nested       Nested       Nested
                          3   Nested       Nested       Nested       Nested
                          4   Nested       Nested       Nested       Nested

 ubar                     1   Closed       Closed       Closed       Closed
                          2   Nested       Nested       Nested       Nested
                          3   Nested       Nested       Nested       Nested
                          4   Nested       Nested       Nested       Nested

 vbar                     1   Closed       Closed       Closed       Closed
                          2   Nested       Nested       Nested       Nested
                          3   Nested       Nested       Nested       Nested
                          4   Nested       Nested       Nested       Nested

 u                        1   Closed       Closed       Closed       Closed
                          2   Nested       Nested       Nested       Nested
                          3   Nested       Nested       Nested       Nested
                          4   Nested       Nested       Nested       Nested

 v                        1   Closed       Closed       Closed       Closed
                          2   Nested       Nested       Nested       Nested
                          3   Nested       Nested       Nested       Nested
                          4   Nested       Nested       Nested       Nested

 temp                     1   Closed       Closed       Closed       Closed
                          2   Nested       Nested       Nested       Nested
                          3   Nested       Nested       Nested       Nested
                          4   Nested       Nested       Nested       Nested

 salt                     1   Closed       Closed       Closed       Closed
                          2   Nested       Nested       Nested       Nested
                          3   Nested       Nested       Nested       Nested
                          4   Nested       Nested       Nested       Nested

 Activated C-preprocessing Options:

 LAKE_JERSEY         Lake Jersey nesting test case , 4 grids: A, C, D, and E
 ANA_BSFLUX          Analytical kinematic bottom salinity flux.
 ANA_BTFLUX          Analytical kinematic bottom temperature flux.
 ANA_SMFLUX          Analytical kinematic surface momentum flux.
 ANA_SSFLUX          Analytical kinematic surface salinity flux.
 ANA_STFLUX          Analytical kinematic surface temperature flux.
 !ASSEMBLE_ALL...    Using mpi_isend/mpi_recv in mp_assemble routine.
 ASSUMED_SHAPE       Using assumed-shape arrays.
 !BOUNDARY_ALLGATHER Using mpi_allreduce in mp_boundary routine.
 !COLLECT_ALL...     Using mpi_isend/mpi_recv in mp_collect routine.
 DJ_GRADPS           Parabolic Splines density Jacobian (Shchepetkin, 2002).
 DOUBLE_PRECISION    Double precision arithmetic.
 MASKING             Land/Sea masking.
 MPI                 MPI distributed-memory configuration.
 NESTING             Nesting grids: Composite and Refinement.
 NONLINEAR           Nonlinear Model.
 NONLIN_EOS          Nonlinear Equation of State for seawater.
 !ONE_WAY            Two-way nesting in refinement grids.
 POWER_LAW           Power-law shape time-averaging barotropic filter.
 PROFILE             Time profiling activated .
 REDUCE_ALLGATHER    Using mpi_allgather in mp_reduce routine.
 !RST_SINGLE         Double precision fields in restart NetCDF file.
 SALINITY            Using salinity.
 SOLVE3D             Solving 3D Primitive Equations.
 TS_C4HADVECTION     Fourth-order centered horizontal advection of tracers.
 TS_C4VADVECTION     Fourth-order centered vertical advection 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.
 UV_QDRAG            Quadratic bottom stress.
 VAR_RHO_2D          Variable density barotropic mode.

 Process Information:


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

    GET_GRID    - Coriolis parameter at RHO-points: f
                   (Grid = 01, File: ../input/lake_jersey_grd_a.nc)
                   (Min =  1.00000000E-04 Max =  1.00000000E-04)

    
 Vertical S-coordinate System, Grid 01:

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

     8   0.0000000   0.0000000        0.000       0.000       0.000       0.000
     7  -0.1250000  -0.1122656       -2.148      -1.250      -2.148      -2.148
     6  -0.2500000  -0.2350037       -4.380      -2.500      -4.380      -4.380
     5  -0.3750000  -0.3654529       -6.674      -3.750      -6.674      -6.674
     4  -0.5000000  -0.5000000       -9.000      -5.000      -9.000      -9.000
     3  -0.6250000  -0.6345471      -11.326      -6.250     -11.326     -11.326
     2  -0.7500000  -0.7649963      -13.620      -7.500     -13.620     -13.620
     1  -0.8750000  -0.8877344      -15.852      -8.750     -15.852     -15.852
     0  -1.0000000  -1.0000000      -18.000     -10.000     -18.000     -18.000

 Vertical S-coordinate System, Grid 02:

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

     8   0.0000000   0.0000000        0.000       0.000       0.000       0.000
     7  -0.1250000  -0.1122656       -2.148      -1.250      -2.148      -2.148
     6  -0.2500000  -0.2350037       -4.380      -2.500      -4.380      -4.380
     5  -0.3750000  -0.3654529       -6.674      -3.750      -6.674      -6.674
     4  -0.5000000  -0.5000000       -9.000      -5.000      -9.000      -9.000
     3  -0.6250000  -0.6345471      -11.326      -6.250     -11.326     -11.326
     2  -0.7500000  -0.7649963      -13.620      -7.500     -13.620     -13.620
     1  -0.8750000  -0.8877344      -15.852      -8.750     -15.852     -15.852
     0  -1.0000000  -1.0000000      -18.000     -10.000     -18.000     -18.000

 Vertical S-coordinate System, Grid 03:

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

     8   0.0000000   0.0000000        0.000       0.000       0.000       0.000
     7  -0.1250000  -0.1122656       -2.148      -1.250      -2.148      -2.148
     6  -0.2500000  -0.2350037       -4.380      -2.500      -4.380      -4.380
     5  -0.3750000  -0.3654529       -6.674      -3.750      -6.674      -6.674
     4  -0.5000000  -0.5000000       -9.000      -5.000      -9.000      -9.000
     3  -0.6250000  -0.6345471      -11.326      -6.250     -11.326     -11.326
     2  -0.7500000  -0.7649963      -13.620      -7.500     -13.620     -13.620
     1  -0.8750000  -0.8877344      -15.852      -8.750     -15.852     -15.852
     0  -1.0000000  -1.0000000      -18.000     -10.000     -18.000     -18.000

 Vertical S-coordinate System, Grid 04:

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

     8   0.0000000   0.0000000        0.000       0.000       0.000       0.000
     7  -0.1250000  -0.1122656       -2.148      -1.250      -2.148      -2.148
     6  -0.2500000  -0.2350037       -4.380      -2.500      -4.380      -4.380
     5  -0.3750000  -0.3654529       -6.674      -3.750      -6.674      -6.674
     4  -0.5000000  -0.5000000       -9.000      -5.000      -9.000      -9.000
     3  -0.6250000  -0.6345471      -11.326      -6.250     -11.326     -11.326
     2  -0.7500000  -0.7649963      -13.620      -7.500     -13.620     -13.620
     1  -0.8750000  -0.8877344      -15.852      -8.750     -15.852     -15.852
     0  -1.0000000  -1.0000000      -18.000     -10.000     -18.000     -18.000

 
 Power filter parameters, Fgamma, gamma =  0.28400   0.18933

 Time Splitting Weights for Grid 02:    ndtfast =  30    nfast =  42
 ==================================

    Primary            Secondary            Accumulated to Current Step

  

 ndtfast, nfast =   30  42   nfast/ndtfast =  1.40000

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

    1.000000000000 1.047601458608 0.523800729304 1.000000000000 1.000000000000

 Power filter parameters, Fgamma, gamma =  0.28400   0.18933

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

 Minimum X-grid spacing, DXmin =  5.00000000E-01 km
 Maximum X-grid spacing, DXmax =  5.00000000E-01 km
 Minimum Y-grid spacing, DYmin =  5.00000000E-01 km
 Maximum Y-grid spacing, DYmax =  5.00000000E-01 km
 Minimum Z-grid spacing, DZmin =  2.14812514E+00 m
 Maximum Z-grid spacing, DZmax =  2.32637641E+00 m

 Minimum barotropic Courant Number =  1.12755310E-01
 Maximum barotropic Courant Number =  1.12755310E-01
 Maximum Coriolis   Courant Number =  9.00000000E-03


 Metrics information for Grid 02:
 ===============================

 Minimum X-grid spacing, DXmin =  1.00000000E-01 km
 Maximum X-grid spacing, DXmax =  1.00000000E-01 km
 Minimum Y-grid spacing, DYmin =  1.00000000E-01 km
 Maximum Y-grid spacing, DYmax =  1.00000000E-01 km
 Minimum Z-grid spacing, DZmin =  2.14812514E+00 m
 Maximum Z-grid spacing, DZmax =  2.32637641E+00 m

 Minimum barotropic Courant Number =  1.12755310E-01
 Maximum barotropic Courant Number =  1.12755310E-01
 Maximum Coriolis   Courant Number =  1.80000000E-03


 Metrics information for Grid 03:
 ===============================

 Minimum X-grid spacing, DXmin =  1.66666667E-01 km
 Maximum X-grid spacing, DXmax =  1.66666667E-01 km
 Minimum Y-grid spacing, DYmin =  1.66666667E-01 km
 Maximum Y-grid spacing, DYmax =  1.66666667E-01 km
 Minimum Z-grid spacing, DZmin =  2.14812514E+00 m
 Maximum Z-grid spacing, DZmax =  2.32637641E+00 m

 Minimum barotropic Courant Number =  1.12755310E-01
 Maximum barotropic Courant Number =  1.12755310E-01
 Maximum Coriolis   Courant Number =  3.00000000E-03


 Metrics information for Grid 04:
 ===============================

 Minimum X-grid spacing, DXmin =  5.55555556E-02 km
 Maximum X-grid spacing, DXmax =  5.55555556E-02 km
 Minimum Y-grid spacing, DYmin =  5.55555556E-02 km
 Maximum Y-grid spacing, DYmax =  5.55555556E-02 km
 Minimum Z-grid spacing, DZmin =  2.14812514E+00 m
 Maximum Z-grid spacing, DZmax =  2.32637641E+00 m

 Minimum barotropic Courant Number =  1.12755310E-01
 Maximum barotropic Courant Number =  1.12755310E-01
 Maximum Coriolis   Courant Number =  1.00000000E-03


 Refined Nested Grid(s) Information: 
 ==================================

 Refined    Donor   Refined   Timestep   Refined
  Grid      Grid    Scale      Ratio   Timesteps

    02        01       05      5.00000      05
    03        01       03      3.00000      03
    04        03       03      3.00000      03

 WARNING: Usually the number of Refined Timesteps must be the same
          as the Refined Scale for numerical stability.


 NLM: GET_STATE - Reading state initial conditions,                       0001-01-01 00:00:00.00
                   (Grid 01, t = 0.0000, File: lake_jersey_ini_a.nc, Rec=0001, Index=1)
                - free-surface
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - vertically integrated u-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - vertically integrated v-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - u-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - v-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - potential temperature
                   (Min =  1.05902920E+01 Max =  1.94097080E+01)
                - salinity
                   (Min =  3.32000313E+01 Max =  3.67999992E+01)

 NLM: GET_STATE - Reading state initial conditions,                       0001-01-01 00:00:00.00
                   (Grid 02, t = 0.0000, File: lake_jersey_ini_c.nc, Rec=0001, Index=1)
                - free-surface
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - vertically integrated u-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - vertically integrated v-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - u-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - v-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - potential temperature
                   (Min =  1.05902920E+01 Max =  1.94097080E+01)
                - salinity
                   (Min =  3.32000313E+01 Max =  3.67999992E+01)

 NLM: GET_STATE - Reading state initial conditions,                       0001-01-01 00:00:00.00
                   (Grid 03, t = 0.0000, File: lake_jersey_ini_d.nc, Rec=0001, Index=1)
                - free-surface
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - vertically integrated u-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - vertically integrated v-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - u-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - v-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - potential temperature
                   (Min =  1.05902920E+01 Max =  1.94097080E+01)
                - salinity
                   (Min =  3.32000310E+01 Max =  3.68000000E+01)

 NLM: GET_STATE - Reading state initial conditions,                       0001-01-01 00:00:00.00
                   (Grid 04, t = 0.0000, File: lake_jersey_ini_e.nc, Rec=0001, Index=1)
                - free-surface
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - vertically integrated u-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - vertically integrated v-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - u-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - v-momentum component
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
                - potential temperature
                   (Min =  1.05902920E+01 Max =  1.94097080E+01)
                - salinity
                   (Min =  3.32000313E+01 Max =  3.67999992E+01)

 Basin information for Grid 01:

 Maximum grid stiffness ratios:  rx0 =   0.000000E+00 (Beckmann and Haidvogel)
                                 rx1 =   0.000000E+00 (Haney)

 Initial domain volumes:  TotVolume =  3.2467500000E+10 m3
                         MinCellVol =  5.3703128450E+05 m3
                         MaxCellVol =  5.8159410269E+05 m3
                            Max/Min =  1.0829799296E+00

 Basin information for Grid 02:

 Maximum grid stiffness ratios:  rx0 =   2.960595E-16 (Beckmann and Haidvogel)
                                 rx1 =   4.961602E-15 (Haney)

 Initial domain volumes:  TotVolume =  2.8080000000E+09 m3
                         MinCellVol =  2.1481251380E+04 m3
                         MaxCellVol =  2.3263764108E+04 m3
                            Max/Min =  1.0829799296E+00

 Basin information for Grid 03:

 Maximum grid stiffness ratios:  rx0 =   2.960595E-16 (Beckmann and Haidvogel)
                                 rx1 =   4.961602E-15 (Haney)

 Initial domain volumes:  TotVolume =  3.4200000000E+09 m3
                         MinCellVol =  5.9670142722E+04 m3
                         MaxCellVol =  6.4621566966E+04 m3
                            Max/Min =  1.0829799296E+00

 Basin information for Grid 04:

 Maximum grid stiffness ratios:  rx0 =   2.960595E-16 (Beckmann and Haidvogel)
                                 rx1 =   4.961602E-15 (Haney)

 Initial domain volumes:  TotVolume =  4.4000000000E+08 m3
                         MinCellVol =  6.6300158580E+03 m3
                         MaxCellVol =  7.1801741073E+03 m3
                            Max/Min =  1.0829799296E+00

 NL ROMS/TOMS: started time-stepping: (Grid: 01 TimeSteps: 000000000001 - 000000001920)
 NL ROMS/TOMS: started time-stepping: (Grid: 02 TimeSteps: 000000000001 - 000000009600)
 NL ROMS/TOMS: started time-stepping: (Grid: 03 TimeSteps: 000000000001 - 000000005760)
 NL ROMS/TOMS: started time-stepping: (Grid: 04 TimeSteps: 000000000001 - 000000017280)


 TIME-STEP YYYY-MM-DD hh:mm:ss.ss  KINETIC_ENRG   POTEN_ENRG    TOTAL_ENRG    NET_VOLUME  Grid
                     C => (i,j,k)       Cu            Cv            Cw         Max Speed

         0 0001-01-01 00:00:00.00  0.000000E+00  8.810744E+01  8.810744E+01  3.600000E+10  01
                       (000,00,0)  0.000000E+00  0.000000E+00  0.000000E+00  0.000000E+00
      DEF_HIS     - creating  history      file, Grid 01: ../output/lake_jersey_his_a.nc
      WRT_HIS     - wrote history     fields (Index=1,1) in record = 0000001               01



 Analytical header files used:

     ROMS/Functionals/ana_btflux.h
     ROMS/Functionals/ana_smflux.h
     ROMS/Functionals/ana_stflux.h

 ROMS/TOMS: DONE... Saturday - April 13, 2019 -  6:40:29 PM

Thanks.
Attachments
roms website:<br />https://www.myroms.org/wiki/Lake_Jersey_Refinement_Example_AC
roms website:
https://www.myroms.org/wiki/Lake_Jersey_Refinement_Example_AC
roms_website.png (177.33 KiB) Viewed 1324 times
my roms run
my roms run
myrun_results.png (18.1 KiB) Viewed 1324 times

Post Reply