请教如何查看wrf模式结果某个经纬度的T2变量值

shangbn305

想请教大家，我运行完WRF后，想查看 wrfout_d04_2011-07-20_00:00:00.nc文件中某个经纬度的T2变量的值，该如何实现呢？
感觉应该简单，但是大脑一片混乱，完全不知道用什么方法，任务比较急，大家有知道的可以告诉我么，万分万分感谢。

lsd1984

送一个ncl的脚本
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;本段代码用来显示指定的经纬度在计算网格中最近的点的编号和经纬度
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
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("//wrfout_d02_2011-08-04_00:00:00.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.7854194,31.19645278,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)


end

shangbn305

lsd1984 发表于 2013-6-21 17:03

                               
登录/注册后可看大图

送一个ncl的脚本
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;本段代码用来显示指定的经纬度在计算网格中最近 ...

非常非常非常感谢！！

Forcast

如果你有输出的CTL文件的话，在GRADS里打开这个CTL文件
第1步 set t 时间
第2步 set lon 经度
第3步 set lat 纬度
第4步 d t2
结果就会显示出来了

shangbn305

Forcast 发表于 2013-6-22 09:25

                               
登录/注册后可看大图

如果你有输出的CTL文件的话，在GRADS里打开这个CTL文件
第1步 set t 时间
第2步 set lon 经度

嗯，谢谢您，恍然大悟

瓶子

lsd1984 发表于 2013-6-21 17:03

                               
登录/注册后可看大图

送一个ncl的脚本
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;本段代码用来显示指定的经纬度在计算网格中最近 ...

正需要这个，非常感谢分享，我是ncl菜鸟，嘿嘿

leavesxxm

lsd1984 发表于 2013-6-21 17:03
送一个ncl的脚本
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;本段代码用来显示指定的经纬度在计算网格中最 ...

正需要这个，真的非常感谢！

hahahahaha

lsd1984 发表于 2013-6-21 17:03
送一个ncl的脚本
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;本段代码用来显示指定的经纬度在计算网格中最 ...

你好，我用你的代码试了一下，但是有一些问题，你知道是为什么吗？

Variable: lonlat
Type: float
Total Size: 8 bytes
            2 values
Number of Dimensions: 1
Dimensions and sizes:   [lon_lat_location | 2]
Coordinates:
(0)     121.6453
(1)     31.39188
fatal:["Execute.c":7432]:performASSIGN_VAR_VAR_OP: Number of subscripts on rhs do not match
                        number of dimensions of variable, (3) Subscripts used, (2)

fatal:["Execute.c":8128]:Execute: Error occurred at or near line 14 in file test.ncl

lsd1984

hahahahaha 发表于 2016-5-31 14:21
你好，我用你的代码试了一下，但是有一些问题，你知道是为什么吗？

Variable: lonlat

你确定一下你的时次数量和经纬度坐标范围。我不知道你里边的14行是哪一行

hahahahaha

lsd1984 发表于 2016-6-2 08:23
你确定一下你的时次数量和经纬度坐标范围。我不知道你里边的14行是哪一行

是这一行 t=T2(:,locationij(1)-1,locationij(0)-1)