- 积分
- 165
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2018-9-13
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
我生成文件,用HDFEXP不能读取是怎样?希望大家给建议。代码如下
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl"
outfile = "/output/pr/pr_regrid.nc"
;diri = "/input/data/pr/"
fil = "/input/data/pr/pr_GPCP-SG_L3_v2.3_197901-201710.nc"
f = addfile(fil, "r")
tim = f->time
lat = f->lat
lon = f->lon
n1 = dimsizes(tim)
n2 = dimsizes(lat)
n3 = dimsizes(lon)
;printVarSummary(pr1)
pr = new((/n1,n2,n3/), float)
do i = 0, n1-1
pr(i,:,:) = f->pr(i,:,:)
; write loop content
end do
Opt = True
Opt@SrcFileName = "src_SCRIP.nc"
Opt@DstFileName = "dst_SCRIP.nc"
Opt@WgtFileName = "NCEP_2_Rect.nc"
Opt@ForceOverwrite = True
Opt@SrcGridMask = where(.not.ismissing(pr),1,0) ; use SrcMask2D in NCL V6.2.1 and earlier
Opt@DstGridType = "1x1"
;Opt@DstLLCorner = (/-89.5.d, 89.d/)
;Opt@DstURCorner = (/-180.d, 180.d/)
pr_regrid = new((/n1,181,361/), float)
do i = 0, n1-1
pr_regrid(i,:,:) =ESMF_regrid(pr(i,:,:),Opt)
; write loop content
end do
copy_VarAtts(pr, pr_regrid)
printVarSummary(pr_regrid)
;printVarSummary(src_SCRIP.nc)
;............................................
;write .nc
;............................................
fout = addfile(outfile,"c")
filedimdef(fout,"time",-1,True) ;-- make time and UNLIMITED dimension
fout->pr = pr_regrid
|
|