NCL脚本求教

迷幻城堡

在下NCL小白，今天在学习NCL尝试用ncl脚本将500hpa的位势高度的月平均资料处理成年平均资料，然后画出01961-1990年气候场遇到了一点问题，请教各位大神！下面是我编辑的脚本：系统显示在脚本中红色那行左右有这样一个错误：Subscript out of range, error in subscript #0，应该是哪里下标出了问题 但是我就是找不到。对了 关于读出来的hgt500文件格式如下：dimt （这个dimt是dimt=dimsizes(hgt500) print(dimt)得到的）
Type: integer
Total Size: 12 bytes
            3 values
Number of Dimensions: 1
Dimensions and sizes:   [3]
Coordinates:
(0)     1704
(1)     91
(2)     180

脚本：

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
begin
f=addfile("C:/dream/bj/summer/data/hgt.mon.mean.nc","r")
hgt500=f->hgt(:,10,:,:)
hgt500=hgt500/10
dimt=dimsizes(hgt500)
print(dimt)
lon=f->lon
lat=f->lat
;calculate winterhgt
winterhgt500=new((/141,91,180/),float)
do i=0,90
  do j=0,179
    do it=0,1702
      if(mod(it,12).eq.11)then
       winterhgt500((it+1)/12,i,j)=(hgt500(it,i,j)+hgt500(it+1,i,j)+hgt500(it+2,i,j))/3
       end if
    end do
  end do
end do
copy_VarMeta(hgt500,winterhgt500)
;calculate ave winterhgt for 1961-1990
avehgt6190=new((/91,180/),float)
do i=0,90
  do j=0,179
    do it=90,119
      avehgt6190(i,j)=avehgt6190(i,j)+winterght500(it,i,j)
    end do
    avehgt6190(i,j)=avehgt6190(i,j)/30
  end do
end do
avehgt6190@op="time average:"+winterght500@long_name
wks=gsn_open_wks("png","C:/dream/bj/summer/slp/hgt/avehgt6191two")       ; send graphics to PNG file
res            = True                          ; plot mods desired
res@gsnPolar   = "NH"                          ; specify the hemisphere
res@gsnLeftString         = ""       ; turn off subtitles
res@gsnRightString        = ""
res@gsnCenterString       = ""
plot = gsn_csm_contour_map_polar(wks,avehgt6190,res)    ; create the plot
end

流沙

标红的代码写错了。一共142年的数据，提取141的冬季数据。数组初始化是141大小，那么索引是从0到140。而你写的是从1到141,141自然就越界了。应改为winterhgt500((it+1)/12-1,i,j)=(hgt500(it,i,j)+hgt500(it+1,i,j)+hgt500(it+2,i,j))/3
另外，楼主全是硬编码写的，既然已经dimsizes取了坐标维数，为什么你的dimt(0),dimt(1),dimt(2)来分别代替你的数字呢。

freekiller

1702+2 数组越界了，你既然知道从0开始，应该最后是1073，而不是1074

迷幻城堡

流沙 发表于 2016-8-7 09:55
标红的代码写错了。一共142年的数据，提取141的冬季数据。数组初始化是141大小，那么索引是从0到140。而你 ...

嗯呐 这个问题昨天晚上被检查出来了 nc才学用有时候就会跟matlab搞混了。大神我还想请教您一下 我这样在赋值元数据的地方还是有问题？我把月数据处理成年数据那之类长度从1704*91*180变成看141*91*180，那这样元数据改怎么复制呢？用copy_VarMeta好像只能复制属性不能复制坐标变量额  麻烦您了

迷幻城堡

freekiller 发表于 2016-8-7 10:13
1702+2 数组越界了，你既然知道从0开始，应该最后是1073，而不是1074

嗯呐 这个问题昨天晚上被检查出来了 nc才学用有时候就会跟matlab搞混了。大神我还想请教您一下 我这样在赋值元数据的地方还是有问题？我把月数据处理成年数据那之类长度从1704*91*180变成看141*91*180，那这样元数据改怎么复制呢？用copy_VarMeta好像只能复制属性不能复制坐标变量额

流沙

迷幻城堡 发表于 2016-8-7 10:31
嗯呐 这个问题昨天晚上被检查出来了 nc才学用有时候就会跟matlab搞混了。大神我还想请教您一下 我这样在 ...

坐标元数据需要自己设置。因为时间维已经改了。
如 a&lat = lat等

dakangshuji

同为ncl小白一个