Error in mksurfdata for regional surfaces

[m_migliavacca1@campus_unimib_it]

Dear all,

I am running mksurfdata for generating regional surfaces for CLM.
I have already generated the griddata.nc, fracdata.nc and the domain using mkgriddata and mkdatadomain tools.
The new grids have 120x101 elements (resolution 0.25x0.25 regular grid)
When I run mksurfdata.pl I get this error after a while:

Attempting to make elevation .....
(GETFIL): attempting to find local file mksrf_topo.10min.c080912.nc
(GETFIL): using /gpfs/scratch/userjrc/jrcccu10/CLM/Inputdata/lnd/clm2/rawdata/mksrf_topo.10min.c080912.nc
(GETFIL): attempting to find local file mksrf_topo.10min.c080912.nc
(GETFIL): using /gpfs/scratch/userjrc/jrcccu10/CLM/Inputdata/lnd/clm2/rawdata/mksrf_topo.10min.c080912.nc
read_domain read lsmlon and lsmlat dims
read_domain initialized domain
read_domain read LONGXY and LATIXY fields
read_domain check NUMLON for regular grid
read_domain compute lat[ns],lon[we]
celledge, using celledge_global_new
read_domain compute cellarea
cellarea, using cellarea_global
AREAINI warning: conservation check not valid for
input grid of 2160 x 1080
output grid of 120 x 101
(GETFIL): attempting to find local file fracdata_0101x0120.nc
(GETFIL): using /gpfs/scratch/userjrc/jrcccu10/CLM_Kloster/models/lnd/clm/tools/mkgriddata/fracdata_0101x0120.nc
(GETFIL): attempting to find local file fracdata_0101x0120.nc
(GETFIL): using /gpfs/scratch/userjrc/jrcccu10/CLM_Kloster/models/lnd/clm/tools/mkgriddata/fracdata_0101x0120.nc
read_domain read lsmlon and lsmlat dims
read_domain initialized domain
read_domain read LONGXY and LATIXY fields
read_domain check NUMLON for regular grid
read_domain ERROR not regular grid, stop 120 2160

Please note that fracdata_101x0120.nc has a regular grid (0.25) and the dimensions are correct (120x101). Here my namelist:

&clmexp
mksrf_fgrid = '/gpfs/scratch/userjrc/jrcccu10/CLM_Kloster/models/lnd/clm/tools/mkgriddata/griddata_0101x0120.nc'
mksrf_fsoitex = '/gpfs/scratch/userjrc/jrcccu10/CLM/Inputdata/lnd/clm2/rawdata/mksrf_soitex.10level.c010119.nc'
mksrf_forganic = '/gpfs/scratch/userjrc/jrcccu10/CLM/Inputdata/lnd/clm2/rawdata/mksrf_organic.10level.0.5deg.081112.nc'
mksrf_flanwat = '/gpfs/scratch/userjrc/jrcccu10/CLM/Inputdata/lnd/clm2/rawdata/mksrf_lanwat.050425.nc'
mksrf_fmax = '/gpfs/scratch/userjrc/jrcccu10/CLM/Inputdata/lnd/clm2/rawdata/mksrf_fmax.070406.nc'
mksrf_fglacier = '/gpfs/scratch/userjrc/jrcccu10/CLM/Inputdata/lnd/clm2/rawdata/mksrf_glacier.060929.nc'
mksrf_fvocef = '/gpfs/scratch/userjrc/jrcccu10/CLM/Inputdata/lnd/clm2/rawdata/mksrf_vocef.c060502.nc'
mksrf_ffrac = '/gpfs/scratch/userjrc/jrcccu10/CLM_Kloster/models/lnd/clm/tools/mkgriddata/fracdata_0101x0120.nc'
mksrf_ftopo = '/gpfs/scratch/userjrc/jrcccu10/CLM/Inputdata/lnd/clm2/rawdata/mksrf_topo.10min.c080912.nc'
mksrf_furban = '/gpfs/scratch/userjrc/jrcccu10/CLM/Inputdata/lnd/clm2/rawdata/mksrf_urban_3den_0.5x0.5_simyr2000.c090223_v1.nc'
mksrf_gridtype = 'regional'
mksrf_fvegtyp =
'/gpfs/scratch/userjrc/jrcccu10/CLM/Inputdata/lnd/clm2/rawdata/pftlandusedyn.0.5x0.5.simyr1850-2005.c090630/mksrf_landuse_rc2000_c090630.nc'
mksrf_fsoicol =
'/gpfs/scratch/userjrc/jrcccu10/CLM/Inputdata/lnd/clm2/rawdata/pftlandusedyn.0.5x0.5.simyr1850-2005.c090630/mksrf_soilcol_global_c090324.nc'
mksrf_flai =
'/gpfs/scratch/userjrc/jrcccu10/CLM/Inputdata/lnd/clm2/rawdata/pftlandusedyn.0.5x0.5.simyr1850-2005.c090630/mksrf_lai_global_c090506.nc'
mksrf_firrig = ' '
mksrf_fdynuse = 'pftdyn_hist_simyr2000-2000.txt'
all_urban = .false.
outnc_double = .true.
/

Do you have any idea about how to fix it?
thank you very much
mirco

 

[slevis]

Sorry, I don't know how to fix this. However, I have a suggestion (maybe you've already tried it). Follow the example in the User's Guide for building a regional dataset using mksurfdata for the Amazon region. If that works, then look for differences between that and your failed attempt.

Sam

 

[xujunhan@gmail_com]

Dear micro,

I always use:

mksrf_ffrac = '../mkgriddata/fracdata_10min_USGS_071205.nc'

I found that if I use my fraction data, mksurfdata will not work.

Xujun

 

[m_migliavacca1@campus_unimib_it]

Dear Xujun,
thank you very much, as you suggested it seems that mksurfdata doesn't work with my own fraction data. At least for the extraxtion of regional surfaces.
Using:
mksrf_ffrac = '../mkgriddata/fracdata_10min_USGS_071205.nc'
it works,
mirco