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

asked 2021-01-23 07:23:33 +0200

updated 2021-01-23 07:26:40 +0200

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!

edit retag flag offensive close merge delete

add a comment

0

answered 2021-01-25 15:12:45 +0200

aakashsane
1

Hi Andres, I don't have year 2004 files. I have whatever I downloaded from CROCO website: https://www.croco-ocean.org/documenta... >>> 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 ; }