error run getregional_datasets.py for 0.125x0.125 sufdata using

[Maggie Xia]

Dear all,

I am trying to subset a 0.125x0.125 surfdata (obtained from https://svn-ccsm-inputdata.cgd.ucar...fdata_0.125x0.125_78pfts_simyr2000_c151014.nc) using the command:

./getregional_datasets.pl -sw 25,73 -ne 41,105 -i inlist_1 -o outlist_1

and the error occured:


I have used this tool successfully several times, not sure why this happened. Can anyone give me a hint here? Thanks in advance.

 

[oleson]

Are you sure this file exists in that directory?

/disk/cesmlab/inputdata/lnd/clm2/surfdata_map/release-clm5.0.18/surfdata_0.125x0.125_78pfts_simyr2000_c151014.nc

In our inputdata, that file is not in the release-clm5.0.18 sub-directory, it is simply in the surfdata_map directory.

What is the domain file you are using (fatmlndfrc)?

 

[oleson]

I think the problem is in these lines in getregional_datasets.ncl:

getfilecoord_namenlen( filelatnames, dimnames, dimlens, nlat, filelatnm );
getfilecoord_namenlen( filelonnames, dimnames, dimlens, nlon, filelonnm );

The function is expecting an integer for dimlens, but is getting a long integer. This is occurring for the variable MONTHLY_LAI on the surface dataset. The reason this variable is different from the other variables in the file from the perspective of NCL is likely this:

As of version 6.0.0, this function will return type "long" if any of the individual dimension sizes or the product of the dimension sizes is greater than or equal to 2 GB.

MONTHLY_LAI is likely > 2GB for this 0.125x0.125 global surface dataset.

I think you can fix this by changing these lines as:

getfilecoord_namenlen( filelatnames, dimnames, longtoint(dimlens), nlat, filelatnm );
getfilecoord_namenlen( filelonnames, dimnames, longtoint(dimlens), nlon, filelonnm );

This should convert any longs to an integer.