- 积分
- 15746
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-3-8
- 最后登录
- 1970-1-1
|

楼主 |
发表于 2016-7-4 09:49:48
|
显示全部楼层
谢谢版主的回复,我检查了一下,的确拼接后的数组好像没有数据,可是为什么呢?
可以麻烦您帮忙看一下吗,谢谢版主
min_lat = -30.0
max_lat =90.0
min_lon = 30.0
max_lon = 180.0
yrStrt1 = 1901
yrLast1 = 1978
yrStrt2 = 1979
yrLast2 = 2000
season = "DJF" ; choose chun seasonal
f1=addfile("./era20monofdaily.nc","r")
f2=addfile("./era.interim.monofdaily1.1.nc","r")
TIME1 = f1->time
TIME2 = f2->time
YYYY1 = cd_calendar(TIME1,-1)/100 ; entire file
iYYYY1 = ind(YYYY1.ge.yrStrt1 .and. YYYY1.le.yrLast1)
YYYY2 = cd_calendar(TIME2,-1)/100 ; entire file
iYYYY2 = ind(YYYY2.ge.yrStrt2 .and. YYYY2.le.yrLast2)
ishf1 = short2flt(f1->ishf)
ishf2 = short2flt(f2->ishf)
gt1=ishf1(iYYYY1,:,:)
gt2=ishf2(iYYYY2,:,:)
printVarSummary(gt1)
printVarSummary(gt2)
lat=f1->latitude
lon=f1->longitude
dsizes_x1=dimsizes(gt1)
ntim1=dsizes_x1(0)
nlat1=dsizes_x1(1)
nlon1=dsizes_x1(2)
dsizes_x2=dimsizes(gt2)
ntim2=dsizes_x2(0)
nlat2=dsizes_x2(1)
nlon2=dsizes_x2(2)
gt = new((/ntim1+ntim2,nlat1,nlon1/),"float")
time=fspan(0,ntim1+ntim2-1,ntim1+ntim2)
gt!0="time"
gt!1="lat"
gt!2="lon"
gt&time=time
gt&lat=gt1&latitude
gt&lon=gt2&longitude
time!0 = "time"
time@long_name = "time"
time@units = "yr"
lat!0 = "lat"
lat@long_name = "latitude"
lat@units = "degrees_north"
lon!0 = "lon"
lon@long_name = "longitude"
lon@units = "degrees_east"
printVarSummary(gt)
|
|