关于ncl输出 帮帮忙啊 大神们

shushamei

我用 tc_plane = wrf_user_intrp3d(tc,p,"h",pressure,0.,False)得到tc_plane后想把tc_plane输出到nc中

现在我用一下方法建变量T却总是报错
  filevardef(fout, "T",typeof(tc_plane),getvardims(tc_plane))

错误信息：
fatal:FileAddVar: Dimension (south_north) is not currently defined, can't add variable
fatal:Execute: Error occurred at or near line 113 in file wrf_PressureLevel1.ncl

hzwjy

干嘛要用这么复杂的方法创建NC文件，直接addfile，fi->var = tc不就行了么。

你这个报错是因为filevardef之前没有定义维度信息，
完整的例子在这里 http://www.ncl.ucar.edu/Applications/method_2.shtml
主要步骤
0     fout  = addfile ("foo.nc", "c")
1     setfileoption(fout,"DefineMode",True)
2     filedimdef(fout,dimNames,dimSizes,dimUnlim)
3     filevardef(fout, "time" ,typeof(time),getvardims(time))
       filevardef......

4     setfileoption(fout,"DefineMode",False)
5     fout->time   = (/time/)  

shushamei

hzwjy 发表于 2012-9-20 19:30

                               
登录/注册后可看大图

干嘛要用这么复杂的方法创建NC文件，直接addfile，fi->var = tc不就行了么。

你这个报错是因为 ...

你说这个我看了 但是还是调试不通 我想按照一定的格式输出出来

south_north=tc_plane!0
          west_east=tc_plane!1
           dimNames = (/south_north, west_east/)  
          dimSizes = (/tc_plane@south_north, tc_plane@west_east /)
           dimUnlim = (/False, False/)   
         filedimdef(fout,dimNames,dimSizes,dimUnlim)
         filevardef(fout, "south_north" ,typeof(tc_plane@south_north),getvardims(tc_plane@south_north))
       filevardef(fout, "west_east"  ,typeof(tc_plane@west_east),getvardims(tc_plane@west_east) )
       filevardef(fout, "T",typeof(tc_plane),getvardims(tc_plane))
          
           filevarattdef(fout,"south_north" ,south_north)
           filevarattdef(fout,"west_east" ,west_east)
           filevarattdef(fout,"T" ,tc_plane)


其中这一句老是过不去

imaJUNation

感谢分享~~