CLM3.5 runtime error at reading of the forcing data

[asanuma@suiri_tsukuba_ac_jp]

Hi,

I am trying to do a single point simulation with CLM3.5, but encountered a runtime error. The
last portion of the log file says,

==BEGIN===
Successfully made atm -> srf interpolation
Successfully initialized area-averaging interpolation

Successfully set up atmospheric grid
GPTLstop: thread 0 timer for init had not been started.
(GETFIL): attempting to find local file 2003-01.nc
(GETFIL): using /home3/blm/clm/atmfcng/clm3.5//2003-01.nc

------------------------------------------------------------
nstep= 1 date= 20030101 sec= 1800
ATMDRV: attempting to read data from /home3/blm/clm/atmfcng/clm3.5//2003-01.nc
------------------------------------------------------------

MCT::m_AttrVect::init_: allocate() error, stat =179
000.MCT(MPEU)::die.: from MCT::m_AttrVect::init_()
MPI_Abort: error code = 2
==END===

It seems the coupler exits at reading of the forcing data, for which I am using those come with the source files.

The namelist to CLM is as follows:
==BEGIN===
&clm_inparm
caseid = '$case'
ctitle = '${case}: spun -up initialization'
finidat = ''
fsurdat = 'surfdata_0001x0001.nc'
fatmgrid = 'griddata_0001x0001.nc'
fatmlndfrc = 'fracdata_0001x0001.nc'
fpftcon = '$CSMDATA/pftdata/pft-physiology.c070207'
offline_atmdir = '$BLM/clm/atmfcng/clm3.5/'
frivinp_rtm = '$CSMDATA/rtmdata/rdirc.05.061026'
nsrest = 0
dtime = 1800
start_ymd = 20030101
start_tod = 0
stop_ymd = 20040101
stop_tod = 0
irad = -1
wrtdia = .true.
mss_irt = 0
hist_dov2xy = .true.
hist_nhtfrq = 0
hist_mfilt = 36
hist_crtinic = 'MONTHLY'
rpntpath = './lnd.$case.rpointer'
/
&prof_inparm
/
==END====

The griddata "griddata_0001x0001.nc" was constructed by using tools/griddata/griddata with the namelist, which is 1 degree grid box, as follows

==BEGIN===
&clmexp
mksrf_fnavyoro="/home/blm/clm/clm3.5in/lnd/clm2/rawdata/mksrf_navyoro_20min.c010129.nc"
mksrf_lsmlon = 1
mksrf_lsmlat = 1
mksrf_edgew = 108.5
mksrf_edgee = 109.5
mksrf_edges = 46.5
mksrf_edgen = 47.5
/
==END===

Howver, to my surprise, a 2 degree grid box constructed with
====
mksrf_edgew = 107.5
mksrf_edgee = 109.5
mksrf_edges = 45.5
mksrf_edgen = 47.5
/
====
gives a successfull run.

Therefore the problem seems to be linked with the resolution.:confused:

Would you give me a suggestion?

 

[asanuma@suiri_tsukuba_ac_jp]

After further investivation, I found that this runtime error ocurrs at interpa2s() in atmdrvMod.F90, when it calls mct_Smat_AvMult(), and the occurrence depends on the geographical location of the atm and land grids, and their horizontal interval.

So far, I have resolved this by changing (interpolating) the resolution of atm forcing data.

Probably, ver 3.5 does not fully support single column mode.. :confused:

 

[chen1718@umn_edu]

Hi, asanuma:

I'm curious about how did you change the resolution of the forcing data? Did you use another forcing data instead or did you modify the testing data that released with the model?

Thanks,

Ming

 

[asanuma@suiri_tsukuba_ac_jp]

In my case, the computation is for a single point using an surface observations as the forcing. Therefore, I can put the same data at any resolution in space.

If you are looking for a tool to interpolate global data set, CDO (Climate Data Operator) developped at Max Plank is works well.

 

[chen1718@umn_edu]

Thank you for your information! CDO is really useful for me.

Ming

In my case, the computation is for a single point using an surface observations as the forcing. Therefore, I can put the same data at any resolution in space.

If you are looking for a tool to interpolate global data set, CDO (Climate Data Operator) developped at Max Plank is works well.