Skip navigation
NASA Logo, National Aeronautics and Space Administration
Modeling Guru Banner
90 Views 2 Replies Latest reply: Mar 15, 2021 9:50 PM by Bhanu Magotra RSS
Currently Being Moderated

Mar 15, 2021 5:25 AM

Error while running LIS with "genMetforc" plugin

Hi,

I am getting the following error while running LIS with "genMetforc" plugin :

 

"Error in nf90_inq_varid in metForcGen_Variables_read Stopping"

 

I checked  metForcGen_VariablesMod.F90 and I can see that the error occurs from line 347. But, I cannot figure out the reason for this.


My input netcdf file looks like below :

 

$ ncdump -h LDT_HIST_198001010000.d01.nc

netcdf LDT_HIST_198001010000.d01 {

dimensions:

        time = 1 ;

        north_south = 501 ;

        east_west = 751 ;

variables:

        float Rainf(time, north_south, east_west) ;

                Rainf:_FillValue = -9999.9f ;

                Rainf:units = "kg m-2 s-1" ;

                Rainf:coordinates = "lon lat" ;

        float LWdown(time, north_south, east_west) ;

                LWdown:_FillValue = -9999.9f ;

                LWdown:units = "W m-2" ;

                LWdown:coordinates = "lon lat" ;

        float SWdown(time, north_south, east_west) ;

                SWdown:_FillValue = -9999.9f ;

                SWdown:units = "W m-2" ;

                SWdown:coordinates = "lon lat" ;

        float Psurf(time, north_south, east_west) ;

                Psurf:_FillValue = -9999.9f ;

                Psurf:units = "Pa" ;

                Psurf:coordinates = "lon lat" ;

        float Tair(time, north_south, east_west) ;

                Tair:_FillValue = -9999.9f ;

                Tair:units = "K" ;

                Tair:coordinates = "lon lat" ;

        float Wind_E(time, north_south, east_west) ;

                Wind_E:_FillValue = -9999.9f ;

                Wind_E:units = "ms-1" ;

                Wind_E:coordinates = "lon lat" ;

        float Wind_N(time, north_south, east_west) ;

                Wind_N:_FillValue = -9999.9f ;

                Wind_N:units = "ms-1" ;

                Wind_N:coordinates = "lon lat" ;

        double lon(north_south, east_west) ;

                lon:_FillValue = NaN ;

        double lat(north_south, east_west) ;

                lat:_FillValue = -999.9 ;

                lat:units = "W/m2" ;

        int64 time(time) ;

                time:time_increment = "3600" ;

                time:units = "days since 1980-01-01 00:00:00" ;

                time:calendar = "proleptic_gregorian" ;

 

 

// global attributes:

                :description = "The Indian Monsoon Data Assimilation and Analysis reanalysis (IMDAA), Hourly forcing over Indian subcontinent" ;

                :missing_value = -9999. ;

                :MAP_PROJECTION = "EQUIDISTANT CYLINDRICAL" ;

                :SOUTH_WEST_CORNER_LAT = -15. ;

                :SOUTH_WEST_CORNER_LON = 30. ;

                :NORTH_EAST_CORNER_LAT = 45. ;

                :NORTH_EAST_CORNER_LON = 120. ;

                :DX = 0.12 ;

                :DY = 0.12 ;

                :zenith_interp = "true,false," ;

 

I am also attaching the log file for reference.

Attachments:
  • Hi Bhanu,

     

    Was this NetCDF file generated by LDT? There are some differences that indicate it was not generated by LDT.

     

    You are seeing the error because genMetForc reader expects a file generated by LDT with a particular format. For example, it expects variable names that follow the pattern *_inst, *_tavg, or *_acc. You can see this in the calls to  nf90_inq_varid in metForcGen_VariablesMod.F90 (line 348):

     

    ios = nf90_inq_varid( nid, "Tair_"//forcopts%stat_airtmp, varid )

     

     

    where the value of forcopts%stat_airtmp is appended to Tair_ to form, for example, Tair_inst. Your variables lack these suffixes and therefore cannot be found.

     

    I've attached a copy of an LDT-generated forcing file that you might be able to use as a template for formatting your input files for use with that particular reader.

     

    Best,

    Brendan

     

    (thanks to Kristi for help on this answer)

    Attachments:

More Like This

  • Retrieving data ...

Bookmarked By (0)

Legend

  • Correct Answers - 4 points
  • Helpful Answers - 2 points
USAGov logo NASA Logo - nasa.gov