【求助】如何用NCL写nc文件。【已解决】

一人。

各位朋友，本人这两天研究官网如何写nc文件的东西，我要读取wrfinput_d02里的INDEX变量，安照官网做了，但是出现问题了。脚本如下。

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
begin

   fi   = addfile ("wrfinput_d02", "r")
   lat  = fi->XLAT
   time = fi->Times
   lon  = fi->XLONG
   in   = fi->LU_INDEX(:,:,:)

   nlat = dimsizes(lat)
   nlon = dimsizes(lon)
   ntim = dimsizes(time)


   f = addfile("out.nc","c")

   dimNames = (/"time","lat","lon"/)
   dimSizes = (/-1,nlat,nlon/)
   dimUnlim = (/True,False,False/)
   filedimdef(f,dimNames,dimSizes,dimUnlim)



   filevardef(f,"time",typeof(time),getvardims(time))
   filevardef(f,"lat" ,typeof(lat) ,"lat")
   filevardef(f,"lon" ,typeof(lon) ,"lon")
   filevardef(f,"in" ,typeof(in) , getvardims(in) )

   filevarattdef(f,"time",time)
   filevarattdef(f,"lat" ,lat)
   filevarattdef(f,"lon" ,lon)
   filevarattdef(f,"in",in)

   f->time    =  (/time/)
   f->lat     =  (/lat/)
   f->lon     =  (/lon/)
   f->in      =  (/in/)



end

报错如下：
fatal:Could not create (out.nc)
fatal:_NclBuildArray: each element of a literal array must have the same dimension sizes, at least one item doesn't
fatal:["Execute.c":8128]:Execute: Error occurred at or near line 18 in file 1.ncl


求拯救。@longlivehj

andrewsoong

system("/bin/rm -f out.nc")
fout= addfile("out.nc","c")
setfileoption(fout,"DefineMode",True)

一人。

andrewsoong 发表于 2015-12-2 15:57
system("/bin/rm -f out.nc")
fout= addfile("out.nc","c")
setfileoption(fout,"DefineMode",True)

还是不对。

longlivehj

wrf网格的原因，不能这么写。

andrewsoong

longlivehj 发表于 2015-12-2 16:38
wrf网格的原因，不能这么写。

荒川C网格，该怎么写哇，其实我也不知道了。。。。。

一人。

longlivehj 发表于 2015-12-2 16:38
wrf网格的原因，不能这么写。

嗯，是的，谢谢你。

longlivehj

andrewsoong 发表于 2015-12-2 16:48
荒川C网格，该怎么写哇，其实我也不知道了。。。。。

嘿嘿，要么就按wrfout的形式去写，借助west_east之类的维；要么就先插值成rectilinear网格，再写入nc文件。

1374203670

楼主会写了么 从nc文件中按照要求提取出来的数据如何重新写到另外一个nc文件中呢

一人。

1374203670 发表于 2016-3-15 21:25
楼主会写了么 从nc文件中按照要求提取出来的数据如何重新写到另外一个nc文件中呢

你看下帖子的其他人的回复，就用我那个脚本，提取wrfout需要的格式。最近忙，你自己研究一下。

在水一方，

我也想请教一下