立即注册 登录
气象家园 返回首页

山色in的个人空间 http://bbs.06climate.com/?134279 [收藏] [复制] [分享] [RSS]

日志

NCL插值,格点变多,精细化

已有 170 次阅读2023-5-23 23:54 | 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)

评论 (0 个评论)

facelist doodle 涂鸦板

您需要登录后才可以评论 登录 | 立即注册

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

返回顶部