unable to find requested variable: river

[adundunis]

Hi everyone.

I'm trying to run my model in an estuary. Compiling was successful but I encountered a problem during the run.

Pasted below are the last parts of the print out in the terminal:

INQUIRE - unable to find requested variable: river
in file:
Data/NAKDONG_river1.nc

Elapsed CPU time (seconds):

Thread # 0 CPU: 0.319
Total: 0.319

Nonlinear model elapsed time profile, Grid: 01

Allocation and array initialization .............. 0.121 (37.9038 %)
2D/3D coupling, vertical metrics ................. 0.004 ( 1.3507 %)
Omega vertical velocity .......................... 0.002 ( 0.4797 %)
Equation of state for seawater ................... 0.010 ( 3.2412 %)
Total: 0.137 42.9754

All percentages are with respect to total time = 0.319


ROMS/TOMS - Output NetCDF summary for Grid 01:

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


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

I can see that it says that the problem is in my river discharge forcing file, specifically, variable "river" cannot be found. However, I checked and recreated that input file many times and made sure that the said variable is indeed in there.

Pasted below is the ncdisp of my river discharge file.

Source:
/Users/Dondon/Documents/NAKDONG/ROMS/Input/Forcing/Discharge/FINAL/NAKDONG_river1.nc
Format:
netcdf4_classic
Dimensions:
river_time = 52563
river = 20
s_rho = 11
Variables:
river_time
Size: 52563x1
Dimensions: river_time
Datatype: double
Attributes:
long_name = 'river runoff time of river 1'
units = 'Days since 01 Dec 2014'
river_transport
Size: 20x52563
Dimensions: river,river_time
Datatype: double
Attributes:
long_name = 'river 1 runoff vertically integrated mass transport'
units = 'meter3 second-1'
river_salt
Size: 20x11x52563
Dimensions: river,s_rho,river_time
Datatype: double
Attributes:
long_name = 'river 1 runoff salinity'
units = 'PSU'
river_temp
Size: 20x11x52563
Dimensions: river,s_rho,river_time
Datatype: double
Attributes:
long_name = 'river 1 runoff potential temperature'
units = 'Celsius'
river_flag
Size: 20x1
Dimensions: river
Datatype: double
Attributes:
long_name = 'river 1 runoff tracer flag'
units = 'Dimensionless'
river_direction
Size: 20x1
Dimensions: river
Datatype: double
Attributes:
long_name = 'river 1 runoff direction'
units = 'Dimensionless'
river_Xposition
Size: 20x1
Dimensions: river
Datatype: double
Attributes:
long_name = 'river 1 XI-position at RHO-points'
units = 'Dimensionless'
river_Eposition
Size: 20x1
Dimensions: river
Datatype: double
Attributes:
long_name = 'river 1 Eta-position at RHO-points'
units = 'Dimensionless'
river_Vshape
Size: 11x1
Dimensions: s_rho
Datatype: double
Attributes:
long_name = 'river 1 runoff mass transport vertical profile'
units = 'Dimensionless'
river
Size: 20x1
Dimensions: river
Datatype: double
Attributes:
long_name = 'river 1 runoff id'
units = 'Dimensionless'

Can I be enlightened of the possible reasons why that error was shown?

Thanks in advance!

[kate]

I'd have to watch it in the debugger. What is the last field it reports on having read before the rivers? My river reading output looks like:

Code: Select all

    GET_NGFLD   - river runoff  XI-positions at RHO-points
                   (Min =  8.50000000E+01 Max =  3.19000000E+02)
    GET_NGFLD   - river runoff ETA-positions at RHO-points
                   (Min =  1.00000000E+00 Max =  2.60000000E+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 =  4.00000000E-02)
    GET_NGFLD   - river runoff mass transport,               t = 39746 00:00:00
                   (Rec=0010650, Index=2, File: NGOA_rivers_1980_2014.nc)
                   (Tmin=      29097.0000 Tmax=      41880.0000)
                   (Min = -7.67441412E+02 Max =  4.20762235E+02)
    GET_NGFLD   - river runoff potential temperature,        t =   299 00:00:00
                   (Rec=0000299, Index=1, File: NGOA_rivers_1980_2014.nc)
                   (Tmin=          1.0000 Tmax=        365.0000)
                   (Min =  4.17200000E+00 Max =  4.17200000E+00)
    GET_NGFLD   - river runoff salinity,                     t =   299 00:00:00
                   (Rec=0000299, Index=1, File: NGOA_rivers_1980_2014.nc)
                   (Tmin=          1.0000 Tmax=        365.0000)
                   (Min =  0.00000000E+00 Max =  0.00000000E+00)

so it doesn't explicitly list reading the "river" variable.

The river variable should be an integer, however.

[adundunis]

Hi Kate!

My debugger stopped at "river runoff direction" and didn't display anything about "river runoff mass transport vertical profile" like yours did, so I suspected it had something to do with my river_Vshape variable. I checked and indeed, I made a mistake configuring it.

I corrected it and now the run reads my discharge file with no error.

I just want you to know that your response helped me a lot in solving this problem. Thanks a bunch!

Now I'll proceed trying to solving other problems in my model because it is still blowing up. I'll probably post more things on this forum in the future whenever I feel like I've exhausted the solutions I can think of.

Regards! And thanks again!