- 积分
- 184
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-5-14
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
begin
year=ispan(1998,2011,1)
month=(/"01","02","03","04","05","06","07","08","09","10","11","12"/)
n=0
out=addfile("/out.hdf","c")
do i=0,13
do j=0,11
fnl=addfile("/PF.RPF."+year(i)+month(j)+"_LEVEL.EastAsia.HDF","r")
;lat= fnl->LAT(1:536767) ;数据太大
;lon= fnl->LON(1:536767)
;flashcount=fnl->FLASHCOUNT(1:536767)
dimNames=(/"lat","lon","flashcount"/)
dimSizes=(/-1,-1,-1/)
do k=0,536766
if (fnl->LAT(k).gt.0.and.fnl->LAT(k).lt.34)then
if (fnl->LON(k).gt.80.and.fnl->LON(k).lt.144)then
count&lat(n)=fnl->LAT(k)
count&lon(n)=fnl->LON(k)
count&flashcount(n)=flashcount(k)
n=n+1
end if
end if
end do
end do
end do
out->lat=(/lat/)
out->lon=(/lon/)
out->flashcount=(/flashcount/)
end
错误fatal:Could not create (/out.hdf)
fatal:["Execute.c":8575]:Execute: Error occurred at or near line 25 in file flashcount.ncl
25行 count&lat(n)=fnl->LAT(k)
|
|