IO error in reading netcdf file

[sankaras]

Hello,

I getting the following error, while running ROMS for an estuary application, with tides and river flow. Any clues about where the problem is?. I generated the grid using seagrid and edited the grid using editmask.


Thanks,

Sankar


OPENCDF - inconsistent size of dimension: xi_psi 59 61

ROMS/TOMS - IO error ................ exit_flag: 4


Elapsed CPU time (seconds):

Thread # 0 CPU: 0.078
Total: 0.078

Nonlinear model elapsed time profile:

Total: 0.000 0.0000

All percentages are with respect to total time = 0.078

ROMS/TOMS - Output NetCDF summary for Grid 01:

ERROR: OPENCDF - Can not open NetCDF file.
REASON: No error


My cpp options were:

ESTUARY_TEST Suspended Sediment Test in an Estuary
ANA_BSFLUX Analytical kinematic bottom salinity flux.
ANA_BTFLUX Analytical kinematic bottom temperature flux.
ANA_SSFLUX Analytical kinematic surface salinity flux.
ANA_STFLUX Analytical kinematic surface temperature flux.
AVERAGES Writing out time-averaged fields.
AVERAGES_AKS Writing out time-averaged vertical S-diffusion.
CURVGRID Orthogonal curvilinear grid.
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.
NONLINEAR Nonlinear Model.
!NONLIN_EOS Linear Equation of State for seawater.
NORTHERN_WALL Wall boundary at Northern edge.
N2S2_HORAVG Horizontal smoothing of buoyancy and shear.
POWER_LAW Power-law shape time-averaging barotropic filter.
PRSGRD31 Standard density Jacobian formulation (Song, 1998).
PROFILE Time profiling activated .
K_GSCHEME Third-order upstream bias advection of TKE fields.
RHO_SURF Include difference between rho0 and surface density.
!RST_SINGLE Double precision fields in restart NetCDF file.
SALINITY Using salinity.
SOLVE3D Solving 3D Primitive Equations.
SOUTH_FSCLAMPED Southern edge, free-surface, Clamped condition.
SOUTH_M2REDUCED Southern edge, 2D momentum, reduced-physics condition.
SOUTH_M3RADIATION Southern edge, 3D momentum, radiation condition.
SOUTH_TCLAMPED Southern edge, tracers, Clamped condition.
SPLINES Conservative parabolic spline reconstruction.
TS_U3HADVECTION Third-order upstream bias horizontal advection of tracers.
TS_C4VADVECTION Fourth-order centered vertical advection of tracers.
UV_ADV Advection of momentum.
UV_U3HADVECTION Third-order upstream bias advection of momentum.
UV_LOGDRAG Logarithmic bottom stress.
VAR_RHO_2D Variable density barotropic mode.
WESTERN_WALL Wall boundary at Western edge.

[jcwarner]

the error is inconsistent size of dimensions.
So your grid has a set of dimensions, the input file has a set of dimensions, the forcing files have grid dimensions, etc etc.
They all need to be consistent.
One needs to be careful. In the x-direction, the grid dimensions for the number of rho points goes from 0 to L.
LP=L+1, so 0:L = LP.
The input file wants Lm , and Lm = L-1.
Similar in y.
There have been other posts in this for you to search.

This is the kind of info that needs to be added to the wiki, so that we only have to do this once. Kate or Hernan had a nice set of grid patterns that could be put on
https://www.myroms.org/wiki/index.php?t ... Generation

[sankaras]

John,

Yes, I had the wrong dimensions in the input file. I should have carefully read the glossary at the end of the input file.

Lm Number of INTERIOR grid RHO-points in the XI-direction for
! each nested grid, [1:Ngrids]. If using NetCDF files as
! input, Lm=xi_rho-2 where "xi_rho" is the NetCDF file
! dimension of RHO-points. Recall that all RHO-point
! variables have a computational I-range of [0:Lm+1].
!
! Mm Number of INTERIOR grid RHO-points in the ETA-direction for
! each nested grid, [1:Ngrids]. If using NetCDF files as
! input, Mm=eta_rho-2 where "eta_rho" is the NetCDF file
! dimension of RHO-points. Recall that all RHO-point
! variables have a computational J-range of [0:Mm+1].


Thanks very much for your help.

Sankar