CLM4.0创建地表数据报错，求解？

奥兰多 · 发表于 2013-7-25 14:52:31

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

各位大侠，本人做单点，我是按照Tools及分目录下的README一步一步做的，先是生成网格数据，然后放入inputdata下的griddata里面，再在mksurfdata里面修改mksurfdata.singlept为：
&clmexp
mksrf_fgrid       = '/snfs01/chenxl/inputdata/lnd/clm2/griddata/griddata_1x1_130723.nc'
mksrf_fsoitex    = '/snfs01/chenxl/inputdata/lnd/clm2/rawdata/mksrf_soitex.10level.c010119.nc'
mksrf_forganic    = '/snfs01/chenxl/inputdata/lnd/clm2/rawdata/mksrf_organic.10level.0.5deg.081112.nc'
mksrf_flanwat    = '/snfs01/chenxl/inputdata/lnd/clm2/rawdata/mksrf_lanwat.050425.nc'
mksrf_fmax       = '/snfs01/chenxl/inputdata/lnd/clm2/rawdata/mksrf_fmax.070406.nc'
mksrf_fglacier    = '/snfs01/chenxl/inputdata/lnd/clm2/rawdata/mksrf_glacier.060929.nc'
mksrf_ftopo       = '/snfs01/chenxl/inputdata/lnd/clm2/rawdata/mksrf_topo.10min.c080912.nc'
mksrf_ffrac       = '/snfs01/chenxl/inputdata/lnd/clm2/griddata/fracdata_1x1_130723.nc'
mksrf_fvocef    = '/snfs01/chenxl/inputdata/lnd/clm2/rawdata/mksrf_vocef.c060502.nc'
mksrf_furban    = '/snfs01/chenxl/inputdata/lnd/clm2/rawdata/mksrf_urban_3den_0.5x0.5_simyr2000.c090223_v1.nc'
mksrf_gridtype    = 'regional'
mksrf_fvegtyp    = '/snfs01/chenxl/inputdata/lnd/clm2/rawdata/pftlandusedyn.0.5x0.5.simyr1850-2005.c090630/mksrf_landuse_rc2000_c090630.nc'
mksrf_fsoicol    = '/snfs01/chenxl/inputdata/lnd/clm2/rawdata/pftlandusedyn.0.5x0.5.simyr1850-2005.c090630/mksrf_soilcol_global_c090324.nc'
mksrf_flai       = '/snfs01/chenxl/inputdata/lnd/clm2/rawdata/pftlandusedyn.0.5x0.5.simyr1850-2005.c090630/mksrf_lai_global_c090506.nc'
mksrf_firrig    = ' '
mksrf_fdynuse    = 'pftdyn_simyr2000.txt'
outnc_double    = .true.

执行./mksurfdata < mksurfdata.singlept, 运行到这步就出问题了，
Attempting to make elevation .....

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          1  x          1

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          1       2160

关键是最后一句，求各位大侠解答原因？

xo2001 · 发表于 2013-7-25 14:59:30

你把grid的也贴出来看看~

奥兰多 · 发表于 2013-7-25 15:04:09

xo2001 发表于 2013-7-25 14:59

登录/注册后可看大图

你把grid的也贴出来看看~

netcdf griddata_1x1_130723 {
dimensions:
lsmlon = 1 ;
lsmlat = 1 ;
nchar = 128 ;
variables:
int NUMLON(lsmlat) ;
NUMLON:long_name = "number of grid cells at each latitude" ;
NUMLON:units = "unitless" ;
double LONGXY(lsmlat, lsmlon) ;
LONGXY:long_name = "longitude" ;
LONGXY:units = "degrees east" ;
double LATIXY(lsmlat, lsmlon) ;
LATIXY:long_name = "latitude" ;
LATIXY:units = "degrees north" ;
double EDGEN ;
EDGEN:long_name = "northern edge of surface grid" ;
EDGEN:units = "degrees north" ;
double EDGEE ;
EDGEE:long_name = "eastern edge of surface grid" ;
EDGEE:units = "degrees east" ;
double EDGES ;
EDGES:long_name = "southern edge of surface grid" ;
EDGES:units = "degrees north" ;
double EDGEW ;
EDGEW:long_name = "western edge of surface grid" ;
EDGEW:units = "degrees east" ;
double LATN(lsmlat, lsmlon) ;
LATN:long_name = "latitude of north edge" ;
LATN:units = "degrees north" ;
double LONE(lsmlat, lsmlon) ;
LONE:long_name = "longitude of east edge" ;
LONE:units = "degrees east" ;
double LATS(lsmlat, lsmlon) ;
LATS:long_name = "latitude of south edge" ;
LATS:units = "degrees north" ;
double LONW(lsmlat, lsmlon) ;
LONW:long_name = "longitude of west edge" ;
LONW:units = "degrees east" ;
double AREA(lsmlat, lsmlon) ;
AREA:long_name = "area" ;
AREA:units = "km^2" ;

// global attributes:
:Conventions = "NCAR-CSM" ;
:History_Log = "created on: 07-23-13 11:38:46" ;
:Logname = "chenxl" ;
:Host =  ;
:fcamfile =  ;
:fccsmdom =  ;
:fcamtopo =  ;
:fclmgrid =  ;
:fnavyoro = "/snfs01/chenxl/inputdata/lnd/clm2/rawdata/mksrf_navyoro_20min.c010129.nc" ;
:frawtopo =  ;
:File_used = "/snfs01/chenxl/inputdata/lnd/clm2/rawdata/mksrf_navyoro_20min.c010129.nc" ;
:Source = "Community Land Model: CLM3" ;
:Version = "$HeadURL: https://svn-ccsm-models.cgd.ucar ... a/creategridMod.F90 $" ;
:Revision_Id = "$Id: creategridMod.F90 28715 2011-05-27 21:18:44Z erik $" ;
data:

NUMLON = 1 ;

LONGXY =
  91.4 ;

LATIXY =
  31.51 ;

EDGEN = 32.51 ;

EDGEE = 92.4 ;

EDGES = 30.51 ;

EDGEW = 90.4 ;

LATN =
  32.51 ;

LONE =
  92.4 ;

LATS =
  30.51 ;

LONW =
  90.4 ;

AREA =
  42165.4995477498 ;

奥兰多 · 发表于 2013-7-25 15:45:13

xo2001 发表于 2013-7-25 14:59

登录/注册后可看大图

你把grid的也贴出来看看~

高手帮我找见问题了，这块的mksrf_ffrac = '/snfs01/chenxl/inputdata/lnd/clm2/griddata/fracdata_1x1_130723.nc'
不能用第一步生成的，只能用例子自带的，inputdata/lnd/clm2/griddata/fracdata_10min_USGS_071205.nc，运行就好了。不过不知道原因是什么。

andrewsoong · 发表于 2013-7-25 16:13:19

奥兰多发表于 2013-7-25 15:45

登录/注册后可看大图

高手帮我找见问题了，这块的mksrf_ffrac = '/snfs01/chenxl/inputdata/lnd/clm2/griddata/fracdat ...

只能用自带的啊，如果是想用自己的做的更详细的更好的数据那可怎么办？你的不是？

奥兰多 · 发表于 2013-7-25 16:43:21

andrewsoong 发表于 2013-7-25 16:13

登录/注册后可看大图

只能用自带的啊，如果是想用自己的做的更详细的更好的数据那可怎么办？你的不是？

恩，我做的是单点的，不能用我第一步生成的fracdata，用了就报错，换成模式自身的例子，就好了，不过等生成了，地表数据可以自己修改，你可以修改的具体点。

andrewsoong · 发表于 2013-7-25 16:45:56

奥兰多发表于 2013-7-25 16:43

登录/注册后可看大图

恩，我做的是单点的，不能用我第一步生成的fracdata，用了就报错，换成模式自身的例子，就好了，不过等生 ...

我还设想着要使用自己做的数据来做区域的，看来是不行啊，自己一个个修改，多麻烦~~~

奥兰多 · 发表于 2013-7-25 17:42:13

andrewsoong 发表于 2013-7-25 16:45

登录/注册后可看大图

我还设想着要使用自己做的数据来做区域的，看来是不行啊，自己一个个修改，多麻烦~~~

这可能是一个bug，不过要是能解决就好了。估计是咱生成的缺了什么东西吧。有时间可以试着找找看，说不定能解决这个问题。

xo2001 · 发表于 2013-7-25 21:37:29

肯定是你们做的有问题了，我做的很OK啊...

奥兰多 · 发表于 2013-7-26 09:38:07

xo2001 发表于 2013-7-25 21:37

登录/注册后可看大图

肯定是你们做的有问题了，我做的很OK啊...

恩，那就好，我只是试了单点的，区域的没试呢。

		自动登录	找回密码
密码			立即注册