ROMS kinetic energy and other problems

Report or discuss software problems and other woes

Moderators: arango, robertson

Post Reply
Message
Author
scenty
Posts: 7
Joined: Wed Mar 16, 2011 2:38 pm
Location: Xiamen, China

ROMS kinetic energy and other problems

#1 Unread post by scenty »

As a new user of ROMS, I ran a model with the entire Pacific domain and single climatological forcing input file. A part of the output file is shown below:

Code: Select all

Resolution, Grid 01: 0198x0198x032,  Parallel Nodes:  16,  Tiling: 004x004


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

      86400  ntimes          Number of timesteps for 3-D equations.
   1200.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.
         72  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).
         72  nHIS            Number of timesteps between the writing fields
                               into history file.
       2160  ndefHIS         Number of timesteps between creation of new
                               history files.
          1  ntsAVG          Starting timestep for the accumulation of output
                               time-averaged data.
         72  nAVG            Number of timesteps between the writing of
                               time-averaged data into averages file.
       2160  ndefAVG         Number of timesteps between creation of new
                               time-averaged file.
 1.2000E+02  nl_tnu2(01)     NLM Horizontal, harmonic mixing coefficient
                               (m2/s) for tracer 01: temp
 1.2000E+02  nl_tnu2(02)     NLM Horizontal, harmonic mixing coefficient
                               (m2/s) for tracer 02: salt
 6.0000E+02  nl_visc2        NLM Horizontal, harmonic mixing coefficient
                               (m2/s) for momentum.
 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.
 3.0000E-04  rdrg            Linear bottom drag coefficient (m/s).
 3.0000E-03  rdrg2           Quadratic bottom drag coefficient.
 1.0000E-02  Zob             Bottom roughness (m).
          1  lmd_Jwt         Jerlov water type.
          2  Vtransform      S-coordinate transformation equation.
          1  Vstretching     S-coordinate stretching function.
 2.0000E+00  theta_s         S-coordinate surface control parameter.
 0.0000E+00  theta_b         S-coordinate bottom  control parameter.
      0.200  Tcline          S-coordinate surface/bottom layer width (m) used
                               in vertical coordinate stretching.
   1025.000  rho0            Mean density (kg/m3) for Boussinesq approximation.
      1.000  dstart          Time-stamp assigned to model initialization (days).
20060101.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.
     10.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).
          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  Aout(idFsur)    Write out averaged free-surface.
          T  Aout(idUbar)    Write out averaged 2D U-momentum component.
          T  Aout(idVbar)    Write out averaged 2D V-momentum component.
          T  Aout(idUvel)    Write out averaged 3D U-momentum component.
          T  Aout(idVvel)    Write out averaged 3D V-momentum component.
          T  Aout(idWvel)    Write out averaged W-momentum component.
          T  Aout(idOvel)    Write out averaged omega vertical velocity.
          T  Aout(idTvar)    Write out averaged tracer 01: temp
          T  Aout(idTvar)    Write out averaged tracer 02: salt

 Output/Input Files:

             Output Restart File:  /export6/Luwf/Try2/ocean_rst.nc
        Prefix for History Files:  /export6/Luwf/Try2/ocean_his
       Prefix for Averages Files:  /export6/Luwf/Try2/ocean_avg
                 Input Grid File:  /export/home/luwf/Roms_tools/Run/roms_ncfiles/roms_grd.nc
    Input Nonlinear Initial File:  /export/home/luwf/Roms_tools/Run/roms_ncfiles/roms_ini.nc
           Input Forcing File 01:  /export/home/luwf/Roms_tools/Run/roms_ncfiles/roms_frc.nc
             Input Boundary File:  /export/home/luwf/Roms_tools/Run/roms_ncfiles/roms_bry.nc

 Tile partition information for Grid 01:  0198x0198x0032  tiling: 004x004

     tile     Istr     Iend     Jstr     Jend     Npts

        0        1       49        1       49    76832
        1       50       99        1       49    78400
        2      100      149        1       49    78400
        3      150      198        1       49    76832
        4        1       49       50       99    78400
        5       50       99       50       99    80000
        6      100      149       50       99    80000
        7      150      198       50       99    78400
        8        1       49      100      149    78400
        9       50       99      100      149    80000
       10      100      149      100      149    80000
       11      150      198      100      149    78400
       12        1       49      150      198    76832
       13       50       99      150      198    78400
       14      100      149      150      198    78400
       15      150      198      150      198    76832

 Tile minimum and maximum fractional grid coordinates:
   (interior points only)

     tile     Xmin     Xmax     Ymin     Ymax     grid

        0     0.50    49.50     0.50    49.50  RHO-points
        1    49.50    99.50     0.50    49.50  RHO-points
        2    99.50   149.50     0.50    49.50  RHO-points
        3   149.50   198.50     0.50    49.50  RHO-points
        4     0.50    49.50    49.50    99.50  RHO-points
        5    49.50    99.50    49.50    99.50  RHO-points
        6    99.50   149.50    49.50    99.50  RHO-points
        7   149.50   198.50    49.50    99.50  RHO-points
        8     0.50    49.50    99.50   149.50  RHO-points
        9    49.50    99.50    99.50   149.50  RHO-points
       10    99.50   149.50    99.50   149.50  RHO-points
       11   149.50   198.50    99.50   149.50  RHO-points
       12     0.50    49.50   149.50   198.50  RHO-points
       13    49.50    99.50   149.50   198.50  RHO-points
       14    99.50   149.50   149.50   198.50  RHO-points
       15   149.50   198.50   149.50   198.50  RHO-points

        0     1.00    49.50     0.50    49.50    U-points
        1    49.50    99.50     0.50    49.50    U-points
        2    99.50   149.50     0.50    49.50    U-points
        3   149.50   198.00     0.50    49.50    U-points
        4     1.00    49.50    49.50    99.50    U-points
        5    49.50    99.50    49.50    99.50    U-points
        6    99.50   149.50    49.50    99.50    U-points
        7   149.50   198.00    49.50    99.50    U-points
        8     1.00    49.50    99.50   149.50    U-points
        9    49.50    99.50    99.50   149.50    U-points
       10    99.50   149.50    99.50   149.50    U-points
       11   149.50   198.00    99.50   149.50    U-points
       12     1.00    49.50   149.50   198.50    U-points
       13    49.50    99.50   149.50   198.50    U-points
       14    99.50   149.50   149.50   198.50    U-points
       15   149.50   198.00   149.50   198.50    U-points

        0     0.50    49.50     1.00    49.50    V-points
        1    49.50    99.50     1.00    49.50    V-points
        2    99.50   149.50     1.00    49.50    V-points
        3   149.50   198.50     1.00    49.50    V-points
        4     0.50    49.50    49.50    99.50    V-points
        5    49.50    99.50    49.50    99.50    V-points
        6    99.50   149.50    49.50    99.50    V-points
        7   149.50   198.50    49.50    99.50    V-points
        8     0.50    49.50    99.50   149.50    V-points
        9    49.50    99.50    99.50   149.50    V-points
       10    99.50   149.50    99.50   149.50    V-points
       11   149.50   198.50    99.50   149.50    V-points
       12     0.50    49.50   149.50   198.00    V-points
       13    49.50    99.50   149.50   198.00    V-points
       14    99.50   149.50   149.50   198.00    V-points
       15   149.50   198.50   149.50   198.00    V-points

 Maximum halo size in XI and ETA directions:

               HaloSizeI(1) =     124
               HaloSizeJ(1) =     124
                TileSide(1) =      56
                TileSize(1) =    3136


 Activated C-preprocessing Options:

 Tsunami             Tsunami
 ANA_BSFLUX          Analytical kinematic bottom salinity flux.
 ANA_BTFLUX          Analytical kinematic bottom temperature flux.
 ASSUMED_SHAPE       Using assumed-shape arrays.
 AVERAGES            Writing out time-averaged fields.
 CURVGRID            Orthogonal curvilinear grid.
 DIFF_GRID           Horizontal diffusion coefficient scaled by grid size.
 DJ_GRADPS           Parabolic Splines density Jacobian (Shchepetkin, 2002).
 DOUBLE_PRECISION    Double precision arithmetic.
 EASTERN_WALL        Wall boundary at Eastern edge.
 LMD_BKPP            KPP bottom boundary layer mixing.
 LMD_CONVEC          LMD convective mixing due to shear instability.
 LMD_DDMIX           LMD double-diffusive mixing.
 LMD_MIXING          Large/McWilliams/Doney interior mixing.
 LMD_NONLOCAL        LMD convective nonlocal transport.
 LMD_RIMIX           LMD diffusivity due to shear instability.
 LMD_SKPP            KPP surface boundary layer mixing.
 MASKING             Land/Sea masking.
 MIX_GEO_TS          Mixing of tracers along geopotential surfaces.
 MIX_S_UV            Mixing of momentum along constant S-surfaces.
 MPI                 MPI distributed-memory configuration.
 NONLINEAR           Nonlinear Model.
 NONLIN_EOS          Nonlinear Equation of State for seawater.
 NORTHERN_WALL       Wall boundary at Northern edge.
 POWER_LAW           Power-law shape time-averaging barotropic filter.
 PROFILE             Time profiling activated .
 QCORRECTION         Surface net heat flux correction.
 RST_SINGLE          Single precision fields in restart NetCDF file.
 SALINITY            Using salinity.
 SOLAR_SOURCE        Solar Radiation Source Term.
 SOLVE3D             Solving 3D Primitive Equations.
 SOUTH_FSRADIATION   Southern edge, free-surface, radiation condition.
 SOUTH_M2FLATHER     Southern edge, 2D momentum, Flather condition.
 SOUTH_M3CLAMPED     Southern edge, 3D momentum, Clamped condition.
 SOUTH_TCLAMPED      Southern edge, tracers, Clamped condition.
 SOUTH_VOLCONS       Southern edge, enforce mass conservation.
 SPLINES             Conservative parabolic spline reconstruction.
 SPHERICAL           Spherical grid configuration.
 SPONGE              Enhanced horizontal mixing in the sponge areas.
 TS_C4HADVECTION     Fourth-order centered horizontal advection of tracers.
 TS_C4VADVECTION     Fourth-order centered vertical advection of tracers.
 TS_DIF2             Harmonic mixing of tracers.
 UV_ADV              Advection of momentum.
 UV_COR              Coriolis term.
 UV_C4ADVECTION      Fourth-order centered differences advection of momentum.
 UV_LOGDRAG          Logarithmic bottom stress.
 UV_VIS2             Harmonic mixing of momentum.
 VAR_RHO_2D          Variable density barotropic mode.
 VISC_GRID           Horizontal viscosity coefficient scaled by grid size.
 VISC_3DCOEF         Horizontal, time-dependent 3D viscosity coefficient.
 WESTERN_WALL        Wall boundary at Western edge.
After one-year-running, the kinetic energy shows a trend to be less and less.
Kin.jpg
Kin.jpg (25.43 KiB) Viewed 6725 times
And the flow field is strange. The Kuroshio is not visible.
heichao.jpg
heichao.jpg (207.74 KiB) Viewed 6725 times
So as the tropical region:
chidao.jpg
chidao.jpg (208.09 KiB) Viewed 6725 times
Any suggestion or help is appreciated!

User avatar
donguiyo
Posts: 29
Joined: Thu Mar 11, 2004 5:17 pm
Location: BOEM

Re: ROMS kinetic energy and other problems

#2 Unread post by donguiyo »

I found an old hard drive where I **might** have the original cpp.h file. I will post/send it later tonight. Your cpp.h has some non-standard choices. Are you using tidal forcing?

User avatar
donguiyo
Posts: 29
Joined: Thu Mar 11, 2004 5:17 pm
Location: BOEM

Re: ROMS kinetic energy and other problems

#3 Unread post by donguiyo »

Below are the CPP options that I used for the 2004 Global Roms, but 1st let me suggest a perhaps easier alternative for your southern boundary. Why not extending your southern bdry to Antarctica and using EW_PERIODIC in the east-west direction for the ACC? that has been sucessfully done in the past with other models. Or with that same configuration you could also define nudging conditions for tracers (Temp and salinity) and momentum (incoming) and radiation conditions for outgoing information including sea level (no nudging for sea level though). Also, see my note below about TS_DIF2 and UV_VIS2.....with all these it should run fine.

#define GLOBAL
/*
** Options for WORLD OCEAN - resolution from 0.4 (Eq.)
** to 0.8 degrees - enhanced in +/- 35 and tropics
*/
/* output options */
#undef DIAGNOSTICS_UV /* define if writing out momentum diagnostics */
#undef DIAGNOSTICS_TS /* define if writing out tracer diagnostics */
#undef AVERAGES

/* advection and coriolis */
#define UV_ADV
#define UV_COR
#define TS_U3HADVECTION
#define TS_C4VADVECTION

/* pressure gradient */
#define DJ_GRADPS /* Splines density Jacobian (Shchepetkin, 2000) */

/* diffusion and viscosities schemes */ ** TRY TO UNDEF UV_VIS2 and TS_DIF2 and see if your solutions remain finite
#define UV_VIS2
#define VISC_GRID /* viscosity coefficient scaled by grid size */
#define MIX_S_UV

#define TS_DIF2
#define DIFF_GRID /* Diffusivity coefficient scaled by grid size */
#define MIX_GEO_TS

#define LMD_MIXING
#ifdef LMD_MIXING
# define LMD_RIMIX
# define LMD_CONVEC
# define LMD_NONLOCAL
# define LMD_SKPP
#endif
#define UV_QDRAG

/* equation of state - relaxation - salinity - coordinates - etc */ [ ** BELOW DEFINE SPHERICAL **]
#define NONLIN_EOS
#define SOLAR_SOURCE
#define SALINITY
#define CURVGRID
#define SPHERICAL
#define MASKING
#define SOLVE3D
#define SPLINES

/* boundary conditions */
#define EW_PERIODIC
#define NORTHERN_WALL /* Northern edge, closed wall condition */
#undef WESTERN_WALL /* Northern edge, closed wall condition */
#undef EASTERN_WALL /* Northern edge, closed wall condition */
#define SOUTHERN_WALL /* Southern edge, closed wall condition */

/* Climats processing and relaxations */
#define QCORRECTION
#define SRELAXATION
#undef TCLIMATOLOGY
#define ANA_BSFLUX
#define ANA_BTFLUX
/* END OF GLOBAL */

Guillermo.

turuncu
Posts: 128
Joined: Tue Feb 01, 2005 8:21 pm
Location: Istanbul Technical University (ITU)
Contact:

Re: ROMS kinetic energy and other problems

#4 Unread post by turuncu »

Hi,

I just wonder about the grid that you used for the global simulation. I try to create same kind of grid (periodic in East-West direction) to simulate the small belt around equator but i have some difficulties to setup EW periodic grid. If it is possible could you share your grid with me.

Best regards,

--ufuk

Post Reply