日志
NCL插值，格点变多，精细化
|
begin
;这边定义的glon和glat是你想要插值的区域的经纬度以及格点数
glon=fspan(73,106,111);定义一个浮点型数据，从90-112，分为221个格点
glat=fspan(24,43,61)
glon!0="lon"
glon@long_name="lon"
glon@units="degrees_east"
glon&lon=glon
glat!0="lat"
glat@long_name="lat"
glat@units="degrees_north"
glat@lat=glat
; Files1=systemfunc("ls evspsbl_*197901-201412.nc")
; f=addfiles(Files1,"r")
f=addfile("evspsbl_ACCESS-CM2_historical_197901-201412.nc","r")
; lat1=f[0]->lat
; lon1=f[0]->lon
; evp0=f[0]->evspsbl
；经纬度一定要提取出来
lat1=f->lat
lon1=f->lon
evp0=f->evspsbl
e0=dim_avg_n_Wrap(evp0,0)
evpin=linint2(lon1,lat1,e0,False,glon,glat,0)
evpin!0="lat"
lat@long_name="lon2d"
lat@units="degrees_east"
evpin!1="lon"
a=addfile("ACCESS-CM2.historical.197901-201412.nc","r")
pr=a->pr
lat2=a->lat
lon2=a->lon
p=dim_avg_n_Wrap(pr,0)
pin=linint2(lon2,lat2,p,False,glon,glat,0)
；定义pin维度
pin!0="lat"
pin!1="lon"
pe=pin-evpin
pe=60*60*24*pe
；如果有经纬度格点数一样可以直接复制，不一样要设置lat@units
copy_VarMeta(p,pe)