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 ?
initial-nest
Re: initial-nest
Your answer will make me sure about what i have done.
Cheers
jafar
Cheers
jafar
Re: initial-nest
Maybe you should tell us what you've tried and what works for you.
-
arango
- Site Admin
- Posts: 1367
- Joined: Wed Feb 26, 2003 4:41 pm
- Location: DMCS, Rutgers University
- Contact:
Re: initial-nest
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:
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.
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
Re: initial-nest
Dear Arango
many many thanks. i found my answer with your nice post.
cheers
jafar
many many thanks. i found my answer with your nice post.
cheers
jafar
Re: initial-nest
Base on your post i have set
but it gives errors which show that model is finding initial data for river in refined grid 
and when i set
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
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 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.ncCode: 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 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 Re: initial-nest
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
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(
):
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 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 Re: initial-nest
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.
Re: initial-nest
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.
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
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu
Re: initial-nest
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
..
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
..
Re: initial-nest
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
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
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
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
Re: initial-nest
To clarify, since this was not interpreted correctly:Wilkin reply said that whenever UV_psource is defined, both grids need river as force file.
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.
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:Therefor child grid needs river data even if it does not have any river cell.
! 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
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu
Re: initial-nest
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 :
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
Many thanks again for your patient
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.sorry if i am too confused about this concept
Many thanks again for your patient
Re: initial-nest
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.
Or you might look into what the ROMS-AGRIF people support.