ncl出错问题

水霖玲

[img]file:///C:\Users\lenovo\Documents\Tencent Files\569417977\Image\NL[A32AIB32MP@]U35ZJ{O6.jpg[/img]


请问各大神们如何解决后两条警告啊。

longlivehj

如果你用到了onedtond，那么你一维数组的元素个数大于输出的多维数组元素个数，检查各维维数相乘是否等于一维数组元素个数。

水霖玲

longlivehj 发表于 2014-4-22 14:02
如果你用到了onedtond，那么你一维数组的元素个数大于输出的多维数组元素个数，检查各维维数相乘是否等于一 ...

这个解决了，后两个警告呢

水霖玲

longlivehj 发表于 2014-4-22 14:02
如果你用到了onedtond，那么你一维数组的元素个数大于输出的多维数组元素个数，检查各维维数相乘是否等于一 ...

dim = dimsizes(mmy)  
var_3d_summer=new((/30,12/),"float")
var_3d_cli=new((/12,30/),"float")
mmt=new(30,"float")
var_2d_winter1=new(30,"float")
var_2d_winter2=new(30,"float")
var_3d_summer = onedtond(mmy,(/dim/12,12/))  ;ti qu chu ji jie
var_3d_summer!1="season"
var_3d_summer!0="time"
var_3d_cli= var_3d_summer(season|:,time|:)   
var_2d_winter1 = dim_avg_n(var_3d_cli(0:1,:),0)
var_2d_winter2 = dim_avg_n(var_3d_cli(11,:),0)     
mmt=(var_2d_winter1+var_2d_winter2)/2.    ;DJF

;************************************************
; extract time series from 3d data
;************************************************
   in = addfile("$NCARG_ROOT/lib/ncarg/data/cdf/hadlsst.nc","r")
  
   ts =in->sst(1331:1690,{-80.5:80.5},{-179.5:179.5})
   ts!0="time"       ;1870.1
   ts!1="lat"
   ts!2="lon"
   laa=dimsizes(ts&lat)
   loo=dimsizes(ts&lon)
   tt1 = ndtooned(ts)
   tt=new((/30,12,laa,loo/),"float")
   tt = onedtond(tt1,(/30,12,laa,loo/))
   var1=new((/30,laa,loo/),"float")
   var2=new((/30,laa,loo/),"float")
   var3=new((/30,laa,loo/),"float")
   var1 = dim_avg_n(tt(:,0:1,:,:),1)
   var2 = dim_avg_n(tt(:,11,:,:),1)   
   var3=(var1+var2)/2.    ;DJF

倒数第二行出错了，说是左右维度不符，可为何正数12行类似的却没出错。

longlivehj

水霖玲 发表于 2014-4-22 16:11
dim = dimsizes(mmy)  
var_3d_summer=new((/30,12/),"float")
var_3d_cli=new((/12,30/),"float")

如果是求DJF的均值，下面一个语句就可以了：  
var3 = dim_avg_n(tt(:,(/0, 1, 11/), :, :),1)

如果非要按你程序里面的逻辑，tt(:,11,:,:)是3维数组，也就是月份维没有了，dim_avg_n(tt(:,11,:,:),1)是对lat维求均值，出来是(/30, loo/)大小的数组，当然与你的var2定义不符！

水霖玲

longlivehj 发表于 2014-4-22 16:23
如果是求DJF的均值，下面一个语句就可以了：  
var3 = dim_avg_n(tt(:,(/0, 1, 11/), :, :),1)

是的。明白了。太感谢了

水霖玲

longlivehj 发表于 2014-4-22 16:23
如果是求DJF的均值，下面一个语句就可以了：  
var3 = dim_avg_n(tt(:,(/0, 1, 11/), :, :),1)

还有个问题，我设置的中心经度为180，就是画出来的图中间有一条白色的竖线，即使我设置了res@gsnAddCyclic=True也没用。这是与我数据资料只有-179.5-179.5有关吗

longlivehj

水霖玲 发表于 2014-4-22 17:39
还有个问题，我设置的中心经度为180，就是画出来的图中间有一条白色的竖线，即使我设置了res@gsnAddCycli ...

是的，有关系！