initial-nest

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
jafar1979
Posts: 37
Joined: Sat Feb 18, 2012 6:51 pm
Location: inioas

initial-nest

#1 Unread post by jafar1979 »

1-please tell me that if initial condition provide by netcdf file, parent and child need their own initial files, separately?
2-with regard that in NFFILES just number of parent's force files is set: should number of force files for parent and child be same ?

jafar1979
Posts: 37
Joined: Sat Feb 18, 2012 6:51 pm
Location: inioas

Re: initial-nest

#2 Unread post by jafar1979 »

Your answer will make me sure about what i have done.
Cheers
jafar

User avatar
kate
Posts: 4088
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: initial-nest

#3 Unread post by kate »

Maybe you should tell us what you've tried and what works for you.

User avatar
arango
Site Admin
Posts: 1347
Joined: Wed Feb 26, 2003 4:41 pm
Location: DMCS, Rutgers University
Contact:

Re: initial-nest

#4 Unread post by arango »

The answer for both questions is that you need separate files. Nesting is complex and users need to explore first with the larger grid before trying to run with nested grid. It requires expertise to set the nested grids correctly. Users need lots of experience with ROMS, ocean modeling, and the ocean circulation in the area to be modeled.

In one of my refined-grid applications, I have:

Code: Select all

     GRDNAME == ../Data/grd_SCS7km_jl3f.nc \
                ../Data/grd_SCS2km_jw2.nc
     ININAME == ../Data/scs7km-run201-2006-10-03.nc \
                ../Data/scs2km_ini_run201_10-Mar-2006.nc

...

    NFFILES == 6 5

     FRCNAME == ../Data/sustr-scs-era-interim-2006-v3e.nc \
                ../Data/svstr-scs-era-interim-2006-v3e.nc \
                ../Data/shflux-scs-era-interim-2006-v3e.nc \
                ../Data/swflux-scs-era-interim-2006-v3e.nc \
                ../Data/swrad-scs-era-interim-2006-v3e.nc \
                ../Data/scs7km_tide_9harmonics_20060101v2.nc \
                ../Data/sustr-scs-era-interim-2006-v3e.nc \
                ../Data/svstr-scs-era-interim-2006-v3e.nc \
                ../Data/shflux-scs-era-interim-2006-v3e.nc \
                ../Data/swflux-scs-era-interim-2006-v3e.nc \
                ../Data/swrad-scs-era-interim-2006-v3e.nc

! Output NetCDF file names, [1:Ngrids].

     RSTNAME == scs_rst.nc   \
                luzon_rst.nc
     HISNAME == scs_his.nc   \
                luzon_his.nc
     AVGNAME == scs_avg.nc   \
                luzon_avg.nc
     DIANAME == scs_dia.nc   \
                luzon_dia.nc
Notice that the tide harmonics (file 6) is only used in the first (coarse) grid. Files 7 to 11 are the forcing files for the finer grid.

jafar1979
Posts: 37
Joined: Sat Feb 18, 2012 6:51 pm
Location: inioas

Re: initial-nest

#5 Unread post by jafar1979 »

Dear Arango
many many thanks. i found my answer with your nice post.
cheers
jafar

jafar1979
Posts: 37
Joined: Sat Feb 18, 2012 6:51 pm
Location: inioas

Re: initial-nest

#6 Unread post by jafar1979 »

Base on your post i have set

Code: Select all

LtracerSrc == T T F F       !there are river cells just in parent grid
NFFILES ==5 4                    ! number of forcing files
     FRCNAME == /home/jafar/data/Forcing/CS_shflux2_era_06.nc \
                /home/jafar/data/Forcing/wind2stress/CS_wind2stress_06.nc \
                 /home/jafar/data/Forcing/CS_swflux2_era_06.nc \
                 /home/jafar/data/Forcing/CS_swrad2_era_06.nc \
                 /home/jafar/data/river/rivers_force_06.nc \
                      /home/jafar/data/Forcing/CS_shflux2_era_06.nc \
                      /home/jafar/data/Forcing/wind2stress/CS_wind2stress_06.nc \
                      /home/jafar/data/Forcing/CS_swflux2_era_06.nc \
                      /home/jafar/data/Forcing/CS_swrad2_era_06.nc 
but it gives errors which show that model is finding initial data for river in refined grid :shock:

Code: Select all

 NLM: GET_STATE - Read state initial conditions,             t =   303 24:00:00
                   (File: CS_initial_Coh_01nov06_ref1.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.24582748E-06 Max =  1.90413388E+01)
                - salinity
                   (Min =  0.00000000E+00 Max =  1.09521994E+01)
 GET_NGFLD   - unable to find requested variable: river_Xposition
               in file: /home/jafar/data/Forcing/CS_swflux2_era_06.nc
 GET_NGFLD   - unable to find requested variable: river_Xposition
               in file: /home/jafar/data/Forcing/CS_swflux2_era_06.nc
 GET_NGFLD   - unable to find requested variable: river_Xposition
               in file: /home/jafar/data/Forcing/CS_swflux2_era_06.nc
 Elapsed CPU time (seconds):
 GET_NGFLD   - unable to find requested variable: river_Xposition
               in file: /home/jafar/data/Forcing/CS_swflux2_era_06.nc
and when i set

Code: Select all

NFFILES ==5 5                    ! number of forcing files
     FRCNAME == /home/jafar/data/Forcing/CS_shflux2_era_06.nc \
                /home/jafar/data/Forcing/wind2stress/CS_wind2stress_06.nc \
                 /home/jafar/data/Forcing/CS_swflux2_era_06.nc \
                 /home/jafar/data/Forcing/CS_swrad2_era_06.nc \
                 /home/jafar/data/river/rivers_force_06.nc \
                      /home/jafar/data/Forcing/CS_shflux2_era_06.nc \
                      /home/jafar/data/Forcing/wind2stress/CS_wind2stress_06.nc \
                      /home/jafar/data/Forcing/CS_swflux2_era_06.nc \
                      /home/jafar/data/Forcing/CS_swrad2_era_06.nc \
                     /home/jafar/data/river/rivers_force_06.nc 
model run with out problem :!: why it need river data for child grid?does because of UV_PSOURCE and TS_PSOURCE options, it need river's initial for both grids? Now, Would you please tell me what can i do?
P.S

Code: Select all

Coupled Input File name = coupling_rip_current.in
 Model Coupling: 
       Ocean Model MPI nodes: 000 - 003

       Waves Model MPI nodes: 004 - 007
..

..
 Resolution, Grid 01: 0194x0344x032,  Parallel Nodes:   4,  Tiling: 002x002
 Physical Parameters, Grid: 01
 =============================

      21600  ntimes          Number of timesteps for 3-D equations.
    120.000  dt              Timestep size (s) for 3-D equations.
         30  ndtfast         Number of timesteps for 2-D equations between
                               each 3D timestep.
          2  Numgrids        Number of refined grids.
    960.000  TI_WAV_OCN      Time interval (s) between coupling WAV-OCN models.
          8  nOCN_WAV        Number of OCN timesteps between coupling to WAV.
          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.
       9600  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).
        360  nHIS            Number of timesteps between the writing fields
                               into history file.
          1  ntsAVG          Starting timestep for the accumulation of output
                               time-averaged data.
       9600  nAVG            Number of timesteps between the writing of
                               time-averaged data into averages file.
          1  ntsDIA          Starting timestep for the accumulation of output
                               time-averaged diagnostics data.
      28800  nDIA            Number of timesteps between the writing of
                               time-averaged data into diagnostics file.
 1.0000E+02  nl_tnu2(01)     NLM Horizontal, harmonic mixing coefficient
                               (m2/s) for tracer 01: temp
 1.0000E+02  nl_tnu2(02)     NLM Horizontal, harmonic mixing coefficient
                               (m2/s) for tracer 02: salt
 1.0000E+01  nl_visc2        NLM Horizontal, harmonic mixing coefficient
                               (m2/s) for momentum.
 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
 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.
      3.000  gls_p           GLS stability exponent.
      1.500  gls_m           GLS turbulent kinetic energy exponent.
     -1.000  gls_n           GLS turbulent length scale exponent.
 7.6000E-06  gls_Kmin        GLS minimum value of turbulent kinetic energy.
 1.0000E-12  gls_Pmin        GLS minimum value of dissipation.
 5.4770E-01  gls_cmu0        GLS stability coefficient.
 1.4400E+00  gls_c1          GLS shear production coefficient.
 1.9200E+00  gls_c2          GLS dissipation coefficient.
-4.0000E-01  gls_c3m         GLS stable buoyancy production coefficient.
 1.0000E+00  gls_c3p         GLS unstable buoyancy production coefficient.
 1.0000E+00  gls_sigk        GLS constant Schmidt number for TKE.
 1.3000E+00  gls_sigp        GLS constant Schmidt number for PSI.
   1400.000  charnok_alpha   Charnok factor for Zos calculation.
      0.500  zos_hsig_alpha  Factor for Zos calculation using Hsig(Awave).
      0.020  sz_alpha        Factor for Wave dissipation surface tke flux .
    100.000  crgban_cw       Factor for Craig/Banner surface tke flux.
      0.500  wec_alpha       Factor for WEC roller/breaking dissipation.
 2.0000E-04  rdrg            Linear bottom drag coefficient (m/s).
 3.3000E-03  rdrg2           Quadratic bottom drag coefficient.
 4.0000E-03  Zob             Bottom roughness (m).
 4.0000E-03  Zos             Surface roughness (m).
          1  lmd_Jwt         Jerlov water type.
          2  Vtransform      S-coordinate transformation equation.
          4  Vstretching     S-coordinate stretching function.
 8.0000E+00  theta_s         S-coordinate surface control parameter.
 4.0000E-01  theta_b         S-coordinate bottom  control parameter.
      6.000  Tcline          S-coordinate surface/bottom layer width (m) used
                               in vertical coordinate stretching.
   1025.000  rho0            Mean density (kg/m3) for Boussinesq approximation.
    304.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  LtracerSrc(01)  Processing point sources/Sink on tracer 01: temp
          T  LtracerSrc(02)  Processing point sources/Sink on 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
          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(idWamp)    Write out wave height.
          T  Hout(idWlen)    Write out mean wavelength.
          T  Hout(idWdir)    Write out wave direction.
          T  Hout(idTsur)    Write out surface net heat flux.
          T  Hout(idTsur)    Write out surface net salt flux.
          T  Hout(idSrad)    Write out shortwave radiation flux.

 Output/Input Files:

             Output Restart File:  ocean_rip_current_rst.nc
             Output History File:  ocean_rip_current_his.nc
            Output Averages File:  ocean_rip_current_avg.nc
         Output Diagnostics File:  ocean_rip_current_dia.nc
        Physical parameters File:  ocean_rip_current.in
                 Input Grid File:  /home/jafar/data/grid/caspian_grd_4.nc
    Input Nonlinear Initial File:  /home/jafar/data/clim-init/CS_initial_Coh_01nov06.nc
           Input Forcing File 01:  /home/jafar/data/Forcing/CS_shflux2_era_06.nc
           Input Forcing File 02:  /home/jafar/data/Forcing/wind2stress/CS_wind2stress_06.nc
           Input Forcing File 03:  /home/jafar/data/Forcing/CS_swflux2_era_06.nc
           Input Forcing File 04:  /home/jafar/data/Forcing/CS_swrad2_era_06.nc
           Input Forcing File 05:  /home/jafar/data/river/rivers_force_06.nc

 Resolution, Grid 02: 0100x0225x032,  Parallel Nodes:   4,  Tiling: 002x002


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

     108000  ntimes          Number of timesteps for 3-D equations.
     24.000  dt              Timestep size (s) for 3-D equations.
         30  ndtfast         Number of timesteps for 2-D equations between
                               each 3D timestep.
          2  Numgrids        Number of refined grids.
    960.000  TI_WAV_OCN      Time interval (s) between coupling WAV-OCN models.
         40  nOCN_WAV        Number of OCN timesteps between coupling to WAV.
          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.
       9600  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).
       1800  nHIS            Number of timesteps between the writing fields
                               into history file.
          1  ntsAVG          Starting timestep for the accumulation of output
                               time-averaged data.
       9600  nAVG            Number of timesteps between the writing of
                               time-averaged data into averages file.
          1  ntsDIA          Starting timestep for the accumulation of output
                               time-averaged diagnostics data.
      28800  nDIA            Number of timesteps between the writing of
                               time-averaged data into diagnostics file.
 1.0000E+02  nl_tnu2(01)     NLM Horizontal, harmonic mixing coefficient
                               (m2/s) for tracer 01: temp
 1.0000E+02  nl_tnu2(02)     NLM Horizontal, harmonic mixing coefficient
                               (m2/s) for tracer 02: salt
 1.0000E+01  nl_visc2        NLM Horizontal, harmonic mixing coefficient
                               (m2/s) for momentum.
 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
 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.
      3.000  gls_p           GLS stability exponent.
      1.500  gls_m           GLS turbulent kinetic energy exponent.
     -1.000  gls_n           GLS turbulent length scale exponent.
 7.6000E-06  gls_Kmin        GLS minimum value of turbulent kinetic energy.
 1.0000E-12  gls_Pmin        GLS minimum value of dissipation.
 5.4770E-01  gls_cmu0        GLS stability coefficient.
 1.4400E+00  gls_c1          GLS shear production coefficient.
 1.9200E+00  gls_c2          GLS dissipation coefficient.
-4.0000E-01  gls_c3m         GLS stable buoyancy production coefficient.
 1.0000E+00  gls_c3p         GLS unstable buoyancy production coefficient.
 1.0000E+00  gls_sigk        GLS constant Schmidt number for TKE.
 1.3000E+00  gls_sigp        GLS constant Schmidt number for PSI.
   1400.000  charnok_alpha   Charnok factor for Zos calculation.
      0.500  zos_hsig_alpha  Factor for Zos calculation using Hsig(Awave).
      0.020  sz_alpha        Factor for Wave dissipation surface tke flux .
    100.000  crgban_cw       Factor for Craig/Banner surface tke flux.
      0.500  wec_alpha       Factor for WEC roller/breaking dissipation.
 2.0000E-04  rdrg            Linear bottom drag coefficient (m/s).
 3.3000E-03  rdrg2           Quadratic bottom drag coefficient.
 4.0000E-03  Zob             Bottom roughness (m).
 4.0000E-03  Zos             Surface roughness (m).
          1  lmd_Jwt         Jerlov water type.
          2  Vtransform      S-coordinate transformation equation.
          4  Vstretching     S-coordinate stretching function.
 8.0000E+00  theta_s         S-coordinate surface control parameter.
 4.0000E-01  theta_b         S-coordinate bottom  control parameter.
      6.000  Tcline          S-coordinate surface/bottom layer width (m) used
                               in vertical coordinate stretching.
   1025.000  rho0            Mean density (kg/m3) for Boussinesq approximation.
    304.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).
          F  LtracerSrc(01)  Processing point sources/Sink on tracer 01: temp
          F  LtracerSrc(02)  Processing point sources/Sink on 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
          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(idWamp)    Write out wave height.
          T  Hout(idWlen)    Write out mean wavelength.
          T  Hout(idWdir)    Write out wave direction.
          T  Hout(idTsur)    Write out surface net heat flux.
          T  Hout(idTsur)    Write out surface net salt flux.
          T  Hout(idSrad)    Write out shortwave radiation flux.

 Output/Input Files:

             Output Restart File:  ocean_rip_current_rst_ref1.nc
             Output History File:  ocean_rip_current_his_ref1.nc
            Output Averages File:  ocean_rip_current_avg_ref1.nc
         Output Diagnostics File:  ocean_rip_current_dia_ref1.nc
        Physical parameters File:  ocean_rip_current.in
                 Input Grid File:  /home/jafar/data/grid/caspian_grd_ref1.nc
    Input Nonlinear Initial File:  /home/jafar/data/clim-init/CS_initial_Coh_01nov06_ref1.nc
           Input Forcing File 01:  /home/jafar/data/Forcing/CS_shflux2_era_06.nc
           Input Forcing File 02:  /home/jafar/data/Forcing/wind2stress/CS_wind2stress_06.nc
           Input Forcing File 03:  /home/jafar/data/Forcing/CS_swflux2_era_06.nc
           Input Forcing File 04:  /home/jafar/data/Forcing/CS_swrad2_era_06.nc

..
..

..
 Activated C-preprocessing Options:

 RIPCURRENT          Rip Current from Haas and Warner
 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.
 DIAGNOSTICS_TS      Computing and writing tracer diagnostic terms.
 DIAGNOSTICS_UV      Computing and writing momentum diagnostic terms.
 DJ_GRADPS           Parabolic Splines density Jacobian (Shchepetkin, 2002).
 DOUBLE_PRECISION    Double precision arithmetic.
 EASTERN_WALL        Wall boundary at Eastern edge.
 GLS_MIXING          Generic Length-Scale turbulence closure.
 KANTHA_CLAYSON      Kantha and Clayson stability function formulation.
 MASKING             Land/Sea masking.
 MCT_LIB             Using Model Coupling Toolkit library.
 MIX_S_TS            Mixing of tracers along constant S-surfaces.
 MIX_S_UV            Mixing of momentum along constant S-surfaces.
 MPI                 MPI distributed-memory configuration.
 WEC_VF              Vortex Force wave current interaction- UNDER DEVELOPMENT!!!!.
 WDISS_WAVEMOD       Wave energy dissipation acquired from coupled wave model.
 ROLLER_RENIERS      Wave energy roller based on Reniers 2004.
 NONLINEAR           Nonlinear Model.
 NONLIN_EOS          Nonlinear Equation of State for seawater.
 NORTHERN_WALL       Wall boundary at Northern edge.
 N2S2_HORAVG         Horizontal smoothing of buoyancy and shear.
 OUT_DOUBLE          Double precision output fields in NetCDF files.
 PERFECT_RESTART     Processing perfect restart variables.
 POWER_LAW           Power-law shape time-averaging barotropic filter.
 PROFILE             Time profiling activated .
 K_GSCHEME           Third-order upstream advection of TKE fields.
 REFINED_GRID        Refined grids option selected.
 !RST_SINGLE         Double precision fields in restart NetCDF file.
 SALINITY            Using salinity.
 SOLAR_SOURCE        Solar Radiation Source Term.
 SOLVE3D             Solving 3D Primitive Equations.
 SOUTHERN_WALL       Wall boundary at Southern edge.
 SPLINES             Conservative parabolic spline reconstruction.
 SWAN_COUPLING       SWAN model coupling.
 THREE_GHOST         Using three Ghost Points in halo regions.
 TS_U3HADVECTION     Third-order upstream horizontal advection of tracers.
 TS_C4VADVECTION     Fourth-order centered vertical advection of tracers.
 TS_DIF2             Harmonic mixing of tracers.
 TS_PSOURCE          Tracers point sources and sinks.
 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_KIRBY            Compute uwave and vwave Kirby avg velocities.
 UV_QDRAG            Quadratic bottom stress.
 UV_PSOURCE          Mass point sources and sinks.
 UV_VIS2             Harmonic mixing of momentum.
 VAR_RHO_2D          Variable density barotropic mode.
 WAVES_OCEAN         Two-way wave-ocean models coupling.
 WESTERN_WALL        Wall boundary at Western edge.
 EAST_M2SGRADIENT    Eastern edge, 2D stokes, gradient condition.
 WEST_M2SGRADIENT    Western edge, 2D stokes, gradient condition.
 NORTH_M2SGRADIENT   Northern edge, 2D stokes, gradientcondition.
 SOUTH_M2SGRADIENT   Southern edge, 2D stokes, gradient condition.
 EAST_M3SGRADIENT    Eastern edge, 3D stokes, gradient condition.
 WEST_M3SGRADIENT    Western edge, 3D stokes, gradient condition.
 NORTH_M3SGRADIENT   Northern edge, 3D stokes, gradient condition.
 SOUTH_M3SGRADIENT   Southern edge, 3D stokes, gradient condition.
MCT::m_MCTWorld::initm_MCTERROR:  MCTWorld has already been initialized...Continuing
MCT::m_MCTWorld::initm_MCTERROR:  MCTWorld has already been initialized...Continuing
MCT::m_MCTWorld::initm_MCTERROR:  MCTWorld has already been initialized...Continuing
MCT::m_MCTWorld::initm_MCTERROR:  MCTWorld has already been initialized...Continuing

SWAN grid   2 is preparing computation

..

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

..
 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
 Minimum X-grid spacing, DXmin =  3.33514538E+00 km
 Maximum X-grid spacing, DXmax =  3.98831267E+00 km
 Minimum Y-grid spacing, DYmin =  3.33800882E+00 km
 Maximum Y-grid spacing, DYmax =  3.98739982E+00 km
 Minimum Z-grid spacing, DZmin =  1.01057987E-01 m
 Maximum Z-grid spacing, DZmax =  1.83012378E+02 m

 Minimum barotropic Courant Number =  1.18722411E-02
 Maximum barotropic Courant Number =  1.44212673E-01
 Maximum Coriolis   Courant Number =  1.27801807E-02


 NLM: GET_STATE - Read state initial conditions,             t =   303 24:00:00
                   (File: CS_initial_Coh_01nov06.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.24582748E-06 Max =  2.29034667E+01)
                - salinity
                   (Min =  0.00000000E+00 Max =  1.29654990E+01)
    GET_NGFLD   - river runoff  XI-positions at RHO-points
                   (Min =  1.80000000E+01 Max =  1.11000000E+02)
    GET_NGFLD   - river runoff ETA-positions at RHO-points
                   (Min =  8.50000000E+01 Max =  3.31000000E+02)
    GET_NGFLD   - river runoff direction
                   (Min =  0.00000000E+00 Max =  1.00000000E+00)
    GET_NGFLD   - river runoff mass transport vertical profile
                   (Min =  0.00000000E+00 Max =  7.63425333E-02)
    GET_NGFLD   - river runoff mass transport,               t =   289 03:45:00
                   (Rec=0010, Index=2, File: rivers_force_06.nc)
                   (Tmin=         15.2188 Tmax=        350.0313)
                   (Min =  1.95000000E+01 Max =  1.34300000E+02)
    GET_NGFLD   - river runoff potential temperature,        t =   289 03:45:00
                   (Rec=0010, Index=2, File: rivers_force_06.nc)
                   (Tmin=         15.2188 Tmax=        350.0313)
                   (Min =  2.00000000E+00 Max =  1.78000000E+01)
    GET_NGFLD   - river runoff salinity,                     t =   289 03:45:00
                   (Rec=0010, Index=2, File: rivers_force_06.nc)
                   (Tmin=         15.2188 Tmax=        350.0313)
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_2DFLD   - surface u-momentum stress,                 t =   303 24:00:00
                   (Rec=2432, Index=1, File: CS_wind2stress_06.nc)
                   (Tmin=          0.1250 Tmax=        365.0000)
                   (Min = -8.32389989E-06 Max =  7.00611541E-06)
    GET_2DFLD   - surface v-momentum stress,                 t =   303 24:00:00
                   (Rec=2432, Index=1, File: CS_wind2stress_06.nc)
                   (Tmin=          0.1250 Tmax=        365.0000)
                   (Min = -2.51533043E-06 Max =  9.34370765E-06)
    GET_2DFLD   - solar shortwave radiation flux,            t =   303 22:30:00
                   (Rec=2432, Index=1, File: CS_swrad2_era_06.nc)
                   (Tmin=          0.0625 Tmax=        364.9375)
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)
    GET_2DFLD   - surface net heat flux,                     t =   303 22:30:00
                   (Rec=2432, Index=1, File: CS_shflux2_era_06.nc)
                   (Tmin=          0.0625 Tmax=        364.9375)
                   (Min = -3.27425753E-05 Max = -1.92829103E-06)
    GET_2DFLD   - surface net freswater flux, (E-P),         t =   303 22:30:00
                   (Rec=2432, Index=1, File: CS_swflux2_era_06.nc)
                   (Tmin=          0.0625 Tmax=        364.9375)
                   (Min = -1.48378471E-07 Max =  2.19167402E-08)

 Maximum grid stiffness ratios:  rx0 =   3.336777E-01 (Beckmann and Haidvogel)
                                 rx1 =   3.586496E+00 (Haney)


 Initial basin volumes: TotVolume =  7.6968928605E+13 m3
                        MinVolume =  1.1438083862E+06 m3
                        MaxVolume =  2.6995814696E+09 m3
                          Max/Min =  2.3601693274E+03

 == SWAN grid  1 sent wave data to ROMS
 ** ROMS grid  1 recv data from SWAN
 ** ROMS grid  1 sent data to SWAN


 == SWAN grid  1 recvd data from ROMS

 INITIAL: Configuring and initializing forward nonlinear model ...
..
 Power filter parameters, Fgamma, gamma =  0.28400   0.18933

 Minimum X-grid spacing, DXmin =  7.78033760E-01 km
 Maximum X-grid spacing, DXmax =  7.95084769E-01 km
 Minimum Y-grid spacing, DYmin =  7.78294027E-01 km
 Maximum Y-grid spacing, DYmax =  7.95194534E-01 km
 Minimum Z-grid spacing, DZmin =  1.01057987E-01 m
 Maximum Z-grid spacing, DZmax =  1.57008050E+02 m

 Minimum barotropic Courant Number =  1.65909123E-02
 Maximum barotropic Courant Number =  1.31908486E-01
 Maximum Coriolis   Courant Number =  2.15046875E-03


 NLM: GET_STATE - Read state initial conditions,             t =   303 24:00:00
                   (File: CS_initial_Coh_01nov06_ref1.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.24582748E-06 Max =  1.90413388E+01)
                - salinity
                   (Min =  0.00000000E+00 Max =  1.09521994E+01)

 GET_NGFLD   - unable to find requested variable: river_Xposition
               in file: /home/jafar/data/Forcing/CS_swflux2_era_06.nc

 GET_NGFLD   - unable to find requested variable: river_Xposition
               in file: /home/jafar/data/Forcing/CS_swflux2_era_06.nc

 GET_NGFLD   - unable to find requested variable: river_Xposition
               in file: /home/jafar/data/Forcing/CS_swflux2_era_06.nc


 Elapsed CPU time (seconds):

 GET_NGFLD   - unable to find requested variable: river_Xposition
               in file: /home/jafar/data/Forcing/CS_swflux2_era_06.nc
 ..
 ROMS/TOMS - Output NetCDF summary for Grid 01:

 ROMS/TOMS - Output NetCDF summary for Grid 02:

 ROMS/TOMS - Input error ............. exit_flag:   2


 ERROR: Abnormal termination: NetCDF INPUT.
 REASON: No error 

jafar1979
Posts: 37
Joined: Sat Feb 18, 2012 6:51 pm
Location: inioas

Re: initial-nest

#7 Unread post by jafar1979 »

Dear Arango
sorry for same post. I cant find what I can do for above question. I hope to receive your idea about 2 follow questions
1-is there any way that refined grid does not need river data while coarse grid need? Using

Code: Select all

GRDNAME == ../data/grd_3km.nc \
              ../data/grd_ref1.nc
ININAME == ../data/CS_initial_coarse.nc \
           ../data/CS_initail_fine.nc  
NFFILES ==5 4
FRCNAME == ./data/Forcing/CS_shflux2_era_06.nc \
           ./ data/Forcing/wind2stress/CS_wind2stress_06.nc \
           ./ data/Forcing/CS_swflux2_era_06.nc \
           ./data/Forcing/CS_swrad2_era_06.nc \
           ./data/river/rivers_force_06.nc \
           ./data/Forcing/CS_shflux2_era_06.nc \
           ./data/Forcing/wind2stress/CS_wind2stress_06.nc \
           ./data/Forcing/CS_swflux2_era_06.nc \
           ./data/Forcing/CS_swrad2_era_06.nc 
cant eliminate needing of river data for fine grid.
2-Or is this guess true :when river_Xposition and river_Yposition (river's cell positions) are defined in parent grid, using of such river netcdf file in refined grid makes no problem, Therefor there is not problem to use( :?: ):

Code: Select all

LtracerSrc == T T F F
NFFILES ==5 5
     FRCNAME ==./data/Forcing/CS_shflux2_era_06.nc \
               ./ data/Forcing/wind2stress/CS_wind2stress_06.nc \
               ./ data/Forcing/CS_swflux2_era_06.nc \
               ./data/Forcing/CS_swrad2_era_06.nc \
               ./data/river/rivers_force_06.nc \
               ./data/Forcing/CS_shflux2_era_06.nc \
               ./data/Forcing/wind2stress/CS_wind2stress_06.nc \
               ./data/Forcing/CS_swflux2_era_06.nc \
               ./data/Forcing/CS_swrad2_era_06.nc \
               ./data/river/rivers_force_06.nc 

Tomasz
Posts: 23
Joined: Tue Oct 07, 2008 11:27 am
Location: Marine Institute, Ireland

Re: initial-nest

#8 Unread post by Tomasz »

Based on my tests with the refined grids you need to specify rivers in both parent and child grids. That is, the model runs when I set the river in the child grid only, but the parent solution is without the river plume. My impression was that the child's solution is interpolated onto the parent's solution, but apparently this is not the case.

User avatar
wilkin
Posts: 875
Joined: Mon Apr 28, 2003 5:44 pm
Location: Rutgers University
Contact:

Re: initial-nest

#9 Unread post by wilkin »

According to your standard output you have set

UV_PSOURCE Mass point sources and sinks

which is a deprecated option when using point sources with nesting in the myroms.org code.

This option had to be removed because the forcing file handling logic would break if the parent grid had rivers but the child did not (because both grids run the same executable binary).

The option was replaced with logical flags in ocean.in to switch whether a grid had sources, or not. Two options are provided depending on how the source is defined, LuvSrc or LwSrc. One or other has to be set true in ocean.in. They would be reported to output netcdf.

int LuvSrc ;
LuvSrc:long_name = "momentum point sources and sink activation switch" ;
LuvSrc:flag_values = 0, 1 ;
LuvSrc:flag_meanings = ".FALSE. .TRUE." ;
int LwSrc ;
LwSrc:long_name = "mass point sources and sink activation switch" ;
LwSrc:flag_values = 0, 1 ;
LwSrc:flag_meanings = ".FALSE. .TRUE." ;

The LtracerSrc flag you have set to TRUE only controls which tracers inflow in the rivers. The other flags are needed to activate the actual transport.
John Wilkin: DMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu

jafar1979
Posts: 37
Joined: Sat Feb 18, 2012 6:51 pm
Location: inioas

Re: initial-nest

#10 Unread post by jafar1979 »

Then, Is bellow set up correct when we have river just for parent grid:
LuvSrc == T F ! horizontal momentum transport *
LwSrc == F ! volume vertical influx *
LtracerSrc == T T F F
..
..
SSFNAME == /home/jafar/data/river/rivers_force_06.nc
..

jafar1979
Posts: 37
Joined: Sat Feb 18, 2012 6:51 pm
Location: inioas

Re: initial-nest

#11 Unread post by jafar1979 »

Wilkin reply said that whenever UV_psource is defined, both grids need river as force file. Therefor child grid needs river data evenif it does not have any river cell. So base on ROMS structure which way is more correct when river.nc file just contains river cell's positions and data (salt, temp, runoff ,etc) for parent grid:
1- using same river.nc file (with same data) for parent and child grid, simultaneously :?: will it effect on result :?:
2- Or setting all data in river.nc file to zero (as river_zero.nc ) and then use it for child grid while river.nc applied for parent grid :?:
Or not tehre is other way?
Sorry, those questions show that i don't completely know what will do in ROMS :oops: Also i cant use ROMS new version because several refined grids used in my test. There fore i become so appropriate to receive your experience to understand what is correct in ROMS.
Thanks for your attention and time and sorry for such basic questions
cheers
jafar

User avatar
wilkin
Posts: 875
Joined: Mon Apr 28, 2003 5:44 pm
Location: Rutgers University
Contact:

Re: initial-nest

#12 Unread post by wilkin »

Wilkin reply said that whenever UV_psource is defined, both grids need river as force file.
To clarify, since this was not interpreted correctly:

The #define UV_PSOURCE option has gone away in the updated code so if it is defined it does nothing. Its function has been replaced by logical flags in ocean.in.
Therefor child grid needs river data even if it does not have any river cell.
No, the child grid does not need river data if the flag LuvSrc is False. But there should be placeholder to the nonexistent child river file in the list of inputs under SSFNAME, something like:

! Input Sources/Sinks forcing (like river runoff) file name.

SSFNAME == ./Data/espresso_river_20030101_now.nc \
/dev/tmp
John Wilkin: DMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu

jafar1979
Posts: 37
Joined: Sat Feb 18, 2012 6:51 pm
Location: inioas

Re: initial-nest

#13 Unread post by jafar1979 »

Dear Wilkin
at first thanks for your reply. Your answer is completely correct for updated code that it has flags in ocean.in file. I have checked it when i had just one nest grid. But kate told :

Code: Select all

updated code currently only supports one parent with one child grid. No grandchildren as yet.
With regard that there are several child grids in my simulation, Therefore it forced me to not using updated code. So i become really appreciate if you guide me to understand what is the best way to do for child's river file which old code need it( while it need uv_psourceis defined in header file).
sorry if i am too confused about this concept :oops:
Many thanks again for your patient

User avatar
kate
Posts: 4088
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: initial-nest

#14 Unread post by kate »

Sorry, I didn't mean to imply that it used to support grandchildren grids. It is getting better with each update and grandchildren grids are coming.

Or you might look into what the ROMS-AGRIF people support.

Post Reply