Benguela test case with CFSR style bulk forcing gives temporal inconsistency error

inputs-outputs
1

Hi All,

I am trying out the Benguela test case with CFSR formatted bulk forcing netcdf files. The test case and forcings are downloaded from CROCO’s website and I have not modified any netcdf files.

I am getting a temporal consistency error and ERROR: Abnormal termination: netCDF INPUT
I am guessing it is an easy fix but I am unable to do it.

Can someone please tell me how to fix it?

Here is the output log from the failed simulation, more details are below:

 Online forcing: first forcing, year 2005, month  1.
 Online forcing: last forcing year 2005, month 12.
 Online forcing: datasets in ../cfsrdata_benguela/ with  4 records per day.

Activated C-preprocessing Options:

          REGIONAL
          BENGUELA_LR
          OBC_EAST
          OBC_WEST
          OBC_NORTH
          OBC_SOUTH
          CURVGRID
          SPHERICAL
          MASKING
          NEW_S_COORD
          SOLVE3D
          UV_COR
          UV_ADV
          SALINITY
          NONLIN_EOS
          UV_HADV_WENO5
          UV_VIS2
          UV_VADV_SPLINES
          TS_HADV_UP5
          TS_DIF2
          TS_MIX_S
          TS_VADV_WENO5
          GLS_MIXING
          GLS_KEPSILON
          BULK_FLUX
          BULK_FAIRALL
          BULK_LW
          BULK_EP
          BULK_SMFLUX
          ONLINE
          FRC_BRY
          Z_FRC_BRY
          M2_FRC_BRY
          M3_FRC_BRY
          T_FRC_BRY
          ANA_BSFLUX
          ANA_BTFLUX
          OBC_M2CHARACT
          OBC_M3ORLANSKI
          OBC_TORLANSKI
          SPLIT_EOS
          MPI_COMM_WORLD
          M2FILTER_POWER
          HZR
          VAR_RHO_2D
          UV_MIX_S
          NTRA_T3DMIX
          CANUTO_A
          BULK_SM_UPDATE
          CUBIC_INTERP
          NF_CLOBBER


 NUMBER OF THREADS:  1 BLOCKING:  1 x  1.

 Spherical grid detected.

 Time splitting: ndtfast = 60    nfast = 82

 Maximum grid stiffness ratios:   rx0 =  0.239536297927541      rx1 =
   7.93974854937451     
      GET_INITIAL -- Processing data for time =   0.000     record =   1
      GET_BRY    -- Read all boundary data         for time =    345.0    
      GET_BRY    -- Read all boundary data         for time =    15.00    
 Open Meteo file :
 /glade/work/asane/croco_forcfiles/back_upcfsrdata_benguela/Temperature_height_a
 bove_ground_Y2005M1.nc
      ONLINE_BULK -- for temporal consistency open CFSR file of year 2004 month 12 for Temperature_height_above_ground
 Open Meteo file :
 ../cfsrdata_benguela/Temperature_height_a
 bove_ground_Y2004M12.nc

 ONLINE_GET_BULK - ERROR: unable to open forcing NetCDF file: ../cfsrdata_benguela/Temperature_height_above_ground_Y2004M12.nc

 MAIN - number of records written into history  file(s):    0
        number of records written into restart  file(s):    0


 ERROR: Abnormal termination: netCDF INPUT

The croco.in file has the following at the beginning section:
title:

    BENGUELA TEST MODEL

time_stepping: NTIMES dt[sec] NDTFAST NINFO

            720      100      60      1

!run_start_date:
!01/04/2014 00:00:00
!run_end_date:
!01/01/2016 00:00:00
!output_time_steps: DT_HIS(H), DT_AVG(H), DT_RST(H)
1 6 12
grid: filename
…
and the last section is this:

online: byear bmonth recordsperday byearend bmonthend / data path
2005 1 4 2005 12
…/cfsrdata_benguela/
…

I have commented out run start date and end date because I am not using #define USE_CALENDAR.
And the time variables in the croco_ini.nc file has value 0.
The folder path …/cfsrdata_benguela/ contains all the files such as temperature, humidity, etc.

How to fix this netcdf input error?

Thanks!

2

Hi,

What do you get if you type

ncdump -h Temperature_height_a bove_ground_Y2004M12.nc ?

What is the size of the files in that directory, ./cfsrdata_benguela/ ?

3

Hi Andres,
I don’t have year 2004 files. I have whatever I downloaded from CROCO website: https://www.croco-ocean.org/documentation/ >>> Tutorial configurations >>> all tutos.

I have 2005 file: Temperature_height_above_ground_Y2005M4.nc which has size of 25 mb. All other files also have similar size: 23 to 25 mb.

And the ncdump looks like this:

ncdump -h Temperature_height_above_ground_Y2005M12.nc

ncdump -h Temperature_height_above_ground_Y2005M12.nc
netcdf Temperature_height_above_ground_Y2005M12 {
dimensions:
time = 124 ;
lat = 181 ;
lon = 291 ;
variables:
double time(time) ;
time:long_name = “Time” ;
time:units = “days since 1-Jan-2000 00:00:00” ;
float Temperature_height_above_ground(time, lat, lon) ;
float lat(lat) ;
lat:long_name = “latitude of RHO-points” ;
lat:units = “degree_north” ;
float lon(lon) ;
lon:long_name = “longitude of RHO-points” ;
lon:units = “degree_east” ;

// global attributes:
:history = “Thu Sep 26 13:42:42 2019: ncatted -O -a units,time,m,c,days since 1-Jan-2000 00:00:00 Temperature_height_above_ground_Y2005M12.nc\n”,
“Thu Sep 26 13:42:42 2019: ncap2 -O -s time = time - 36524.25 Temperature_height_above_ground_Y2005M12.nc Temperature_height_above_ground_Y2005M12.nc\n”,
“Wed Sep 25 09:22:17 2019: ncks -F -d time,1337,1460 -d lat,20,200 -d lon,490,780 Temperature_height_above_ground_Y2005.nc Temperature_height_above_ground_Y2005M12.nc” ;
:NCO = “4.7.1” ;
:nco_openmp_thread_number = 1 ;
}