- 积分
- 462
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2021-9-17
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
file_wrf=addfile(dir_wrf+"wrfout_d01_2011-03-15_06_00_00","r")
f=addfile("wrfchemi_d01_2011-03-15_06_00_00","r")
f_obs=asciiread(dir_obs+file_obs, (/32,8/), "float");(:,3)
lat2d=file_wrf->XLAT(0,:,:)
lon2d=file_wrf->XLONG(0,:,:)
lat=lat2d(:,0)
lon=lon2d(0,:)
e_ces1=f->E_CES(:,0,:,:)
cesnd=e_ces1
;;e_ces=3600*(3000*3000)*e_ces1 ;;单位换算
;copy_VarMeta(e_ces1,e_ces)
cesnd!0="lat"
cesnd!1="lon"
cesnd&lat=lat
cesnd&lon=lon
;printVarSummary(f_obs)
;站点的经纬度-------------------------
LAT=f_obs(:,0)
LON=f_obs(:,1)
;-----------------------------------------
;fo = rcm2points_Wrap (lat2d,lon2d,cesnd,LON,LAT, 0) ; => fo(ntim,nlvl,5)
fo=linint2_points_Wrap(lon,lat,cesnd, False, LON,LAT, 0)
|
|