ncl处理逐日数据提取指定月份

是冉冉升起的冉

;NCL处理逐日资料，提取多年的指定月份数据，使用cd_calendar函数，提取冬季数据（12月和次年1、2月）
;没有检查提取后的数据是否正确

h=new((/3,90,11,144/),float)
f_1958=addfile("/data/hgt/hgt.1958.nc","r")
t   = f_1958->time
utc_date = cd_calendar(t, 0)
mo  = tointeger(utc_date(:,1)); 月份

begin
level=10;指定高度
do j=0,2,1
do i=1958,1960,1
    yr=i
    f=addfile("/data/hgt/hgt."+yr+".nc","r")
    De     = ind(mo.eq.12); 12月
    h_De   =  f->hgt(De,{level},{65:90},:)   ; 12月份
    yr_JF  =i+1

    f1=addfile("/data/hgt/hgt."+yr_JF+".nc","r")
    JF  = ind(mo.eq.1.or.mo.eq.2)
    h_JF   =  f1->hgt(JF,{level},{65:90},:)  ;次年1月、2月份
    h_DJF=array_append_record(h_De,h_JF,0)

end do
   h(j,:,:,:)=h_DJF(:,:,:)

end do
end