NCL提取wrfout的变量值

★ESPN★

怎么用NCL提取出wrfout里的各高度（米）的风速u，v值呢？？求大神赐教

老郭爱骑车

★ESPN★ 发表于 2014-4-10 11:38
这个语句是用在哪里的呢？

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;本段代码用来显示指定的经纬度在计算网格中最近的点的编号和经纬度
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
;上边是写ascii文件需要导入的脚本文件
begin
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; 输入文件
; 读入一个nc文件。如下：（注意：需要补充好文件类型，否则不识别）
  a = addfile("3.nc","r")
;读取数据
T2=wrf_user_getvar(a,"T2",-1)        ;-1取全部数据        ;slp是sea level pressure海平面气压
times = wrf_user_getvar(a,"times",-1)
xlat=wrf_user_getvar(a,"XLAT",0)
xlon=wrf_user_getvar(a,"XLONG",0)

locationij=wrf_user_ll_to_ij(a,121.95,29.2,True);                ;此命令返回给定的（纬度，经度）对应的最近的节点坐标:网格意义上的（i,j）
;print(locationij)
lonlat=wrf_user_ij_to_ll(a,locationij(0),locationij(1),True);                ;此命令返回给定的节点坐标对应的（纬度，经度）
print(lonlat)
;xx=xlat(locationij(1)-1,locationij(0)-1)                                ;注意ij坐标到编号的转换，需要-1
;print(xx)
t=T2(:,locationij(1)-1,locationij(0)-1)
print(t)
asciiwrite("jwd.txt",t)
end

andrewsoong

官网有么？
http://www.mmm.ucar.edu/wrf/OnLi ... CL/NCL_examples.htm

BDRUSH

andrewsoong 发表于 2014-3-3 13:46

                               
登录/注册后可看大图

官网有么？
http://www.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/NCL_examples.htm

好像没看到

★ESPN★

BDRUSH 发表于 2014-3-3 14:04

                               
登录/注册后可看大图

好像没看到

是的，没看到呢

小__天

ARWpost处理wrfout也行啊，很简单

驍放假了

http://www.mmm.ucar.edu/wrf/OnLi ... TERP/wrf_Height.htm
其实版主给的那个连接里面有个“Height Level”就是了，大概是插值到高度上吧~

★ESPN★

驍放假了 发表于 2014-3-4 16:34

                               
登录/注册后可看大图

http://www.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/Examples/LEVELS_INTERP/wrf_Height.htm
其实版 ...

恩，那然后怎么输出提取的东西呢？

驍放假了

★ESPN★ 发表于 2014-3-4 17:10

                               
登录/注册后可看大图

恩，那然后怎么输出提取的东西呢？

参考里面画图的脚本就可以了吧

★ESPN★

驍放假了 发表于 2014-3-4 19:30

                               
登录/注册后可看大图

参考里面画图的脚本就可以了吧

好的，请问哪个语句是提取的语句呢？我不是要画图，而是提取数据出来

驍放假了

★ESPN★ 发表于 2014-3-6 09:09

                               
登录/注册后可看大图

好的，请问哪个语句是提取的语句呢？我不是要画图，而是提取数据出来

你看看里面，大概是先用wrf_user_getvar函数提取数据，再用wrf_user_intrp3d函数插值吧