Project

General

Profile

wrong grid for LM external parameters for ERA5

Added by Mostafa Hamouda 7 months ago

Dear all,

I am trying to downscale ERA5 cas files using cosmo. I run into an error using int2lm_2_clm4, I tried the suggestions posted in the forum (such as changing startlon_tot by 0.001) but it still didn’t work. I made sure that EXTPAR domain is larger than my requeted LMout. I hope you can help me.
Thank you very much.

Input namelist:


&GRID_IN
  pcontrol_fi = 30000.,
  ie_in_tot = 544,
  je_in_tot = 224,
  ke_in_tot = 98,
  ke_soil_in= 3,
  czml_soil_in = 0.035,0.175,0.64,1.945,
  lushift_in=.FALSE.,.FALSE., lvshift_in=.FALSE.,.FALSE.,
  pollat_in = 90.0,
  pollon_in = 180.0,
  startlat_in_tot = 19.25058,
  startlon_in_tot = -68.90625,
  endlat_in_tot = 81.920,
  endlon_in_tot = 83.812,
/END
&LMGRID
  ivctype = 2,
  irefatm = 2,
  ielm_tot = 300, jelm_tot = 230, kelm_tot = 40,
  pollat =  39.25, pollon = -162.0, polgam = 0, 
  dlon = 0.11, dlat= 0.11, 
  startlat_tot  = -17.00, startlon_tot  = -20.00,
  ke_soil_lm = 9,   
  czml_soil_lm = 0.005,0.025,0.07,0.16,0.34,0.7,1.42,2.86,5.74,11.5,
 /END
 &DATABASE
 /
 &DATA
  yinput_type='analysis',
  ytunit_in='d',
  ytunit_out='d',
  ylmext_form_read='ncdf',
  yinext_form_read='ncdf',
  yin_form_read='ncdf',
  ylm_form_write='ncdf',
  ie_ext = 480, je_ext = 380,
  ylmext_lfn='era5_newest.nc',
  ylmext_cat='${EXPAR_DIR}',
  yinext_lfn='cas${YDATE_INI}.nc',
  yinext_cat='${INPUT_DIR}',
  yin_cat= '${INPUT_DIR}',
  ylm_cat='${OUTPUT_DIR}',
  nprocess_ini = 131, nprocess_bd = 132, 
 /END

outputlog:


 0:  OPEN: ncdf-file:
 0:  /pf/b/b381564/work/DATA_CCLM/SPATE_4DAS/EXTPARAM/era5_newest.nc
 0:    external data file        namelist input
 0:  ie_ext              480                 480
 0:  je_ext              380                 380
 0:  lonll        -30.0000000000000               -20.0000000000000
 0:  latll        -27.0000000000000               -17.0000000000000
 0:  lonur         22.6900005340576                12.8900000000000
 0:  latur         14.6899995803833                8.19000000000000
 0:  dlon         0.110000001114943               0.110000000000000
 0:  dlat         0.109999998892832               0.110000000000000
 0:  pollat        39.2500000000000                39.2500000000000
 0:  pollon       -162.000000000000               -162.000000000000
 0:  polgam       0.000000000000000E+000          0.000000000000000E+000
 0:  wrong grid for LM external parameters
 0:  *------------------------------------------------------------*
 0:  *    PROGRAM TERMINATED BECAUSE OF ERRORS DETECTED
 0:  *              IN ROUTINE:   int2lm_org
 0:  *
 0:  *    ERROR CODE is            2
 0:  *    Error in read_nc_gdefs_ext_lm
 0:  *------------------------------------------------------------*

My EXTPAR report from WebPEP:

Model version = EXTPAR-5.2.1

pollon=-162
pollat=39.25
polgam=0
ie_tot=480
je_tot=380
startlon_tot=-30
startlat_tot=-27
dlon=0.11
dlat=0.11
oro=1
orofilter=0
sgsl=0
landuse=1
soil=1
tcl=1
aot=1
albedo=2
urban=0


Replies (3)

RE: wrong grid for LM external parameters for ERA5 - Added by Burkhardt Rockel 7 months ago

The grid you defined in LMGRID is within the area of the external data set and has the same grid mesh size, however, the nots of the LMGRID are shifted relativ to the external data set grid nots.

RE: wrong grid for LM external parameters for ERA5 - Added by Burkhardt Rockel 7 months ago

startlat_tot  = -16.99, startlon_tot  = -19.99,

may work

RE: wrong grid for LM external parameters for ERA5 - Added by Mostafa Hamouda 7 months ago

It worked!
Thank you for your help

    (1-3/3)