unable to open existing NetCDF file: call from: inquire.

[wangzc]

hello:
I am trying to simulate tide by 2D-ROMS.
I only retain sustr and sms_time in the frc.nc. At the same time, tide is added to the frc.nc. ANA_INITIAL and ANA_BMFLUX is defined.Compiling is ok, but something is wrong when running .
Can anyone help me?

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


Minimum X-grid spacing, DXmin = 2.75248895E+00 km
Maximum X-grid spacing, DXmax = 3.20812877E+00 km
Minimum Y-grid spacing, DYmin = 2.75409673E+00 km
Maximum Y-grid spacing, DYmax = 3.20766198E+00 km

Minimum barotropic Courant Number = 1.85253133E-01
Maximum barotropic Courant Number = 9.71708342E-01
Maximum Coriolis Courant Number = 5.84041068E-03

GET_NGFLD - tidal period
(Min = 4.30819210E+04 Max = 2.38071599E+06)
GET_2DFLD - tidal elevation amplitude
(Min = 3.26393686E-05 Max = 6.68537894E+00)
GET_2DFLD - tidal elevation phase angle
(Min = 2.41898787E-06 Max = 6.28318020E+00)
GET_2DFLD - tidal current inclination angle
(Min = 4.97589188E-07 Max = 6.28318253E+00)
GET_2DFLD - tidal current phase angle
(Min = 2.59098604E-06 Max = 6.28318291E+00)
GET_2DFLD - maximum tidal current, ellipse major axis
(Min = 1.27754016E-06 Max = 1.42023097E+00)
GET_2DFLD - minimum tidal current, ellipse minor axis
(Min = -5.23102760E-01 Max = 4.38586697E-01)
GET_2DFLD - surface u-momentum stress, t = 345 00:00:00
(Rec=0000012, Index=1, File: barotropic_tide_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface v-momentum stress, t = 345 00:00:00
(Rec=0000012, Index=1, File: barotropic_tide_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)

INQUIRE - unable to assign file counter, Fcount = 0
while processing structure:
and variable; zeta_east

NETCDF_OPEN - unable to open existing NetCDF file:


call from: inquire.FQUIRE - unable to assign file counter, Fcount = 0
while processing structure:
and variable; zeta_east
NETCDF_OPEN - unable to open existing NetCDF file:


call from: inquire.F
INITIAL: Configuring and initializing forward nonlinear model ...


Minimum X-grid spacing, DXmin = 2.75248895E+00 km
Maximum X-grid spacing, DXmax = 3.20812877E+00 km
Minimum Y-grid spacing, DYmin = 2.75409673E+00 km
Maximum Y-grid spacing, DYmax = 3.20766198E+00 km

Minimum barotropic Courant Number = 1.85253133E-01
Maximum barotropic Courant Number = 9.71708342E-01
Maximum Coriolis Courant Number = 5.84041068E-03

GET_NGFLD - tidal period
(Min = 4.30819210E+04 Max = 2.38071599E+06)
GET_2DFLD - tidal elevation amplitude
(Min = 3.26393686E-05 Max = 6.68537894E+00)
GET_2DFLD - tidal elevation phase angle
(Min = 2.41898787E-06 Max = 6.28318020E+00)
GET_2DFLD - tidal current inclination angle
(Min = 4.97589188E-07 Max = 6.28318253E+00)
GET_2DFLD - tidal current phase angle
(Min = 2.59098604E-06 Max = 6.28318291E+00)
GET_2DFLD - maximum tidal current, ellipse major axis
(Min = 1.27754016E-06 Max = 1.42023097E+00)
GET_2DFLD - minimum tidal current, ellipse minor axis
(Min = -5.23102760E-01 Max = 4.38586697E-01)
GET_2DFLD - surface u-momentum stress, t = 345 00:00:00
(Rec=0000012, Index=1, File: barotropic_tide_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)
GET_2DFLD - surface v-momentum stress, t = 345 00:00:00
(Rec=0000012, Index=1, File: barotropic_tide_frc.nc)
(Tmin= 15.0000 Tmax= 345.0000)
(Min = 0.00000000E+00 Max = 0.00000000E+00)

INQUIRE - unable to assign file counter, Fcount = 0
while processing structure:
and variable; zeta_east

NETCDF_OPEN - unable to open existing NetCDF file:


call from: inquire.FQUIRE

[kate]

It is having trouble with the boundary conditions. You don't say which boundary options you are using. For tides you probably want Chapman/Flather for the barotropic mode. You then might want to do something like ADD_FSOBC and ANA_FSOBC to add the tidal elevation to an analytic zero elevation. You can look at globaldefs.h to see what the logic is asking for to turn off the searching for boundary files.

[wangzc]

Hi, Kate. Thank you very much!
I took the your suggestion and it did work.I think this is because ROMS will not read LBC values(ubar,vbar,zeta) from *_frc.nc. Either *_bry.nc or ana_fsobc(or ana_m2obc) is needed to ROMS, right? This is also why Lm(or Mn) in the *.in equals Lm(or Mn)-2 in the *_frc.nc made by roms_tools.