LISF 7.5 error: nf90_get_var failed for ps in read_era5

[Mahza2303]

Hello everyone,

I am running LISF Version 7.5 with ERA5 forcing. I downloaded the dataset “Near surface meteorological variables from 1979 to 2019 derived from bias-corrected reanalysis” and stitched the forcing variables into monthly files following the instructions from a LISF discussion thread. LDT completes successfully and generates the inputs, but the problem appears when running LIS.

Here is the ncdump -h output for one of the stitched files (FORCING_201901.nc):

netcdf FORCING_201901 {
dimensions:
time = 744 ;
Number_of_points = 340819 ;
variables:
float Wind(time, Number_of_points) ;
float LAT(Number_of_points) ;
float LON(Number_of_points) ;
float ZS(Number_of_points) ;
float UREF(Number_of_points) ;
float ZREF(Number_of_points) ;
float Tair(time, Number_of_points) ;
float Qair(time, Number_of_points) ;
float PSurf(time, Number_of_points) ;
float LWdown(time, Number_of_points) ;
float Snowf(time, Number_of_points) ;
float Rainf(time, Number_of_points) ;
float DIR_SWdown(time, Number_of_points) ;
}

When I run LIS, the following error appears in the log:

[INFO] Reading ERA5 file (bookend, 1 -/home/bmahzabi/Desktop/1MO_RUN/ERA5/201901/FORCING_201901.nc)
[ERR] nf90_get_var failed for ps in read_era5 Stopping.
[ERR] endrun is being called

I am also attaching my one of the lislog file, lis.config file and read_era5.f90 file below.

lis_config.txt
lislog.txt
read_era5.txt

Thanks in advance for your help!

[dmocko]

Hi @Mahza2303

Sorry for the delay. There are LIS team members who are furloughed due to the shutdown as well as team members who
were traveling.

Your ERA5 netCDF files are close, but aren't exactly what is expected by LIS.

netcdf FORCING_201901 {
dimensions:
        time = 745 ;
        Number_of_points = 340819 ;
variables:
        float Wind(time, Number_of_points) ;
        float LAT(Number_of_points) ;
        float LON(Number_of_points) ;
        float ZS(Number_of_points) ;
        float UREF(Number_of_points) ;
        float ZREF(Number_of_points) ;
        float Tair(time, Number_of_points) ;
        float Qair(time, Number_of_points) ;
        float PSurf(time, Number_of_points) ;
        float LWdown(time, Number_of_points) ;
        float Snowf(time, Number_of_points) ;
        float Rainf(time, Number_of_points) ;
        float DIR_SWdown(time, Number_of_points) ;

....where time is from 0 to 744, with "0" representing 00Z on 01 Jan 2019.

We are working on updating the ERA5 reader in the coming months, to more directly handle data from ECMWF.