请教大家ncl如何写netcdf数据

天边的那边

我想把数据写成nc格式的，不知道怎么写，向大家请教。这个是参照ncl手册写，出错了，不知道问题在哪begin
f=addfile("wrfinput_d01","r")
xlon =f->XLONG(0,:,:)
xlat = f->XLAT(0,:,:)

ff=asciiread("SEVP_CLI_CHN_PRE_DAY_GRID_0.25-20140826.asc",-1,"float")
x2=onedtond(ff(6:),(/144,256/))
x2@_FillValue = -999.00
nx=256
ny=144
rlat=new(ny,"float")
rlon=new(nx,"float")

rlon(0)=72.125
rlat(0)=18.125

do i=0,nx-1
rlon(i)=72.125+i*0.25
end do

do j=0,ny-1
rlat(j)=18.125+j*0.25
end do

x2!0="rlat"
x2!1="rlon"
x2&rlat=rlat(:)
x2&rlon=rlon(:)

wrf_tem=rgrid2rcm(x2&rlat,x2&rlon,x2,lat,lon,0)

rain=addfile("raintest.nc","c")
fileAtt = True
fileAtt@title = "raintogrid"
fileAtt@Conventions = "None"
fileAtt@creation_date = systemfunc("date")
fileattdef(rain,fileAtt)

dimNames = (/"time","lat","lon"/)
dimSizes = (/-1,ny,nx/)
dimUnlim = (/True,False,False/)

filedimdef(rain,dimNames,dimSizes,dimUnlim)
filevardef(rain,"time",typeof(time),getvardims(time))
filevardef(rain,"lat",typeof(lat),"lat")
filevardef(rain,"lon","float","lon")
filevardef(rain,"TMP",typeof(T),getvardims(T))

filevarattdef(rain,"time",time)
filevarattdef(rain,"lat",xlat)
filevarattdef(rain,"lon",xlon)
filevarattdef(rain,"TMP",T)
setfileoption(rain,"SuppressDefineMode",False)

rain->time = (/time/)
rain->wrf_tem = wrf_tem
rain->xlat = xlat
rain->xlon = xlon
print(wrf_tem)
end

这是ncl脚本，load部分没有给出

天边的那边

怎么没有人帮我看问题出在哪里呀

longlivehj

天边的那边 发表于 2014-9-9 11:41
怎么没有人帮我看问题出在哪里呀

ncl程序出错，命令行提示的错误信息对查错是非常重要的。很可惜，楼主用“出错了”三个字概括了。这样的提问，是很难得到有效帮助的。

简单看了一下代码片段，就发现了一个问题。赋值时引用的变量，之前根本就没有定义！
filevardef(rain,"lat",typeof(lat),"lat")
filevardef(rain,"lon","float","lon")
filevardef(rain,"TMP",typeof(T),getvardims(T))

rain->wrf_tem = wrf_tem
rain->xlat = xlat
rain->xlon = xlon

天边的那边

longlivehj 发表于 2014-9-9 12:01
ncl程序出错，命令行提示的错误信息对查错是非常重要的。很可惜，楼主用“出错了”三个字概括了。这样的 ...

之前是定义了的，没有把脚本写全

天边的那边

fatal:Could not create (raintest.nc)
fatal:["Execute.c":8128]:Execute: Error occurred at or near line 45 in file test.ncl
这个是错误信息，第45行内容fileattdef(rain,fileAtt)

TMP的确是没有定义，我看手册中没有明确给出TMP代表什么

innocent

楼主，写nc文件的官网脚本在哪？没找到，麻烦楼主发一下网址，正好要用~