- 积分
- 634
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-10-27
- 最后登录
- 1970-1-1
|
发表于 2020-7-2 23:36:06
|
显示全部楼层
轨道数据网格化可以参考这个例子 :http://www.ncl.ucar.edu/Applications/Scripts/hdf5eos_3a.ncl
;---Set up options for regridding
Opt = True
Opt@SrcGridLat = lat ; unstructured
Opt@SrcGridLon = lon ;
Opt@InterpMethod = "bilinear"
;---If you don't set these two, the regridding will be VERY slow
if (DstGridType.eq."1x1") then
Opt@DstURCorner = (/ 89.5d, 179.5d/) ; for 2x2
end if
if (DstGridType.eq."2x2") then
Opt@DstURCorner = (/ 89d, 179d/) ; for 2x2
end if
if (DstGridType.eq."5x5") then
Opt@DstURCorner = (/ 87.5d, 177.5d/)
end if
Opt@DstLLCorner = -Opt@DstURCorner
Opt@DstGridType = DstGridType ; destination grid
Opt@SrcRegional = False
Opt@DstRegional = False
Opt@RemoveSrcFile = True ; remove SCRIP grid destination files
Opt@RemoveDstFile = True ; remove SCRIP grid destination files
Opt@ForceOverwrite = True
Opt@Debug = True
;---Extract desired level for regridding
xp = x(:,{plvl})
stat_xp = stat_dispersion(xp, opt ) ; ALWAYS look at satellite data
; Here, just 'xp'
Opt@SrcGridMask = where(ismissing(lat) .or. ismissing(xp) .or. \
abs(lat).gt.90 .or. abs(lon).gt.180, 0,1)
Opt@WgtFileName = "MLS-Aura_L2GP-"+app+"."+Opt@DstGridType+".nc"
xp_regrid = ESMF_regrid(xp,Opt)
把每天的轨道数据转为格点数据,再求格点数据的月平均。 |
|