请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 67|回复: 0

NCL将wrf格点数据插值到站点,除了一两个站点都是缺测值,求求大家帮我看看!

[复制链接]

新浪微博达人勋

发表于 3 天前 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

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)

密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

快速回复 返回顶部 返回列表