Unable to inquire about contents of netCDF in svn881

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
thiago.pdpaula
Posts: 14
Joined: Thu Oct 03, 2013 6:44 pm
Location: Personal

Unable to inquire about contents of netCDF in svn881

#1 Unread post by thiago.pdpaula »

Hi,

I'm running a realistic application spanning 6 years (2008 - 2014) and I was able to execute the entire simulation using ROMS version svn783.

However, after upgrading to svn881, I am getting the following error:

Code: Select all

     44160 2009-01-01 00:00:00.00  2.761819E-03  1.905423E+04  1.905423E+04  6.367074E+15
                     (098,094,23)  6.296839E-03  2.155303E-04  9.631030E-02  1.162009E+00
      WRT_HIS     - wrote history     fields (Index=1,1) in record = 0000185
      WRT_RST     - wrote re-start    fields (Index=1,1) in record = 0000002
 Found Error: **   Line: 330      Source: ROMS/Modules/mod_netcdf.F, netcdf_get_dim

 NETCDF_GET_DIM - unable to inquire about contents of input NetCDF file:  /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Uwind_2009.nc
                  call from:  ROMS/Utility/inquiry.F
 Found Error: 02   Line: 462      Source: ROMS/Modules/mod_netcdf.F, netcdf_check_dim
 Found Error: 02   Line: 176      Source: ROMS/Utility/inquiry.F
 Found Error: 02   Line: 119      Source: ROMS/Utility/get_2dfld.F
 Found Error: 02   Line: 134      Source: ROMS/Nonlinear/get_data.F
 Found Error: 02   Line: 205      Source: ROMS/Nonlinear/main3d.F
 Found Error: 02   Line: 274      Source: ROMS/Drivers/nl_ocean.h
ROMS is not able to get information on the Uwind forcing file for the year 2009.

I'm using BULK_FLUXES and I have 8 unique forcing files, one for each year and for each variable (Uwind, Vwind, swrflux, lwrad_rad, Qair, Pair, rain, Tair) defined as follows on ocean.in:

Code: Select all

     NFFILES == 8                          ! number of unique forcing files

     FRCNAME == /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_lwrad_down_2008.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_lwrad_down_2009.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_lwrad_down_2010.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_lwrad_down_2011.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_lwrad_down_2012.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_lwrad_down_2013.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_lwrad_down_2014.nc \

                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_rain_2008.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_rain_2009.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_rain_2010.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_rain_2011.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_rain_2012.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_rain_2013.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_rain_2014.nc \

                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Tair_2008.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Tair_2009.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Tair_2010.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Tair_2011.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Tair_2012.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Tair_2013.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Tair_2014.nc \

                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Qair_2008.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Qair_2009.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Qair_2010.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Qair_2011.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Qair_2012.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Qair_2013.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Qair_2014.nc \

                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Pair_2008.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Pair_2009.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Pair_2010.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Pair_2011.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Pair_2012.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Pair_2013.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Pair_2014.nc \

                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_net_swrad_2008.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_net_swrad_2009.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_net_swrad_2010.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_net_swrad_2011.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_net_swrad_2012.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_net_swrad_2013.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_net_swrad_2014.nc \

                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Uwind_2008.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Uwind_2009.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Uwind_2010.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Uwind_2011.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Uwind_2012.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Uwind_2013.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Uwind_2014.nc \

                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Vwind_2008.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Vwind_2009.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Vwind_2010.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Vwind_2011.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Vwind_2012.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Vwind_2013.nc |
                /nethome/Teo_NetApp/u4mt/roms/grid_SSE_112/inputs/FRC/roms_SSE_112_CFSRR_Vwind_2014.nc 
I can inquire the mentioned netCDF file with ncdump, python, matlab and I'm quite sure every thing seems fine and in accordance to varinfo.dat. ncdump -h outputs the following:

Code: Select all

netcdf roms_SSE_112_CFSRR_Uwind_2009 {
dimensions:
        wind_time = UNLIMITED ; // (8761 currently)
        eta_rho = 222 ;
        xi_rho = 217 ;
variables:
        double wind_time(wind_time) ;
                wind_time:long_name = "time" ;
                wind_time:units = "days since 1970-01-01 00:00:00 GMT" ;
        float lon(eta_rho, xi_rho) ;
                lon:long_name = "longitude" ;
                lon:units = "degree_east" ;
                lon:standard_name = "longitude" ;
        float lat(eta_rho, xi_rho) ;
                lat:long_name = "latitude" ;
                lat:units = "degree_north" ;
                lat:standard_name = "latitude" ;
        float Uwind(wind_time, eta_rho, xi_rho) ;
                Uwind:long_name = "10-meter u-wind component" ;
                Uwind:units = "meter second-1" ;
                Uwind:time = "wind_time" ;
                Uwind:obs = "wind is given in rho-points, wind-stress in u/v points" ;
vimdiff in the routines mod_netcdf.F for versions 783 and 881 shows that there are a few differences between them, but that doesn't explain why the model can read forcing files for 2008 but not for 2009.

Did anybody have similar issue? Any ideas on how to overcome this problem?

I'm really interested in moving to the new version since it has several upgrades.

Regards.

User avatar
m.hadfield
Posts: 521
Joined: Tue Jul 01, 2003 4:12 am
Location: NIWA

Re: Unable to inquire about contents of netCDF in svn881

#2 Unread post by m.hadfield »

> Did anybody have similar issue?

Yes, a week or two ago. ROMS seemed to object to one forcing file in a sequence, the third or fourth as I recall, but the file was similar to all the others as far as I could see. I resolved to look into this at some stage, but haven't got around to it.

I have deleted the stdout log and it would take an hour or two to reproduce it, but as I recall it suggested increasing Mdims in mod_netcdf.F from its current value of 50. I haven't tried that, and I can't see that it should be necessary as none of my files have anything like 50 dimensions. But I guess this suggests a problem in the management of dimension IDs.

> Any ideas on how to overcome this problem?

Nothing beyond stepping laboriously through the code. Unfortunately, I struck the error an hour or two into a simulation, and it would need a bit of work to reproduce it with a shorter run.

Post Reply