NCL提取wrfout的变量值

驍放假了

★ESPN★ 发表于 2014-4-10 11:00
试过了，不行呢，总有错误

你好，你可以把错误的提示信息发上来，或许大家可以帮帮你。

★ESPN★

驍放假了 发表于 2014-4-10 11:46
你好，你可以把错误的提示信息发上来，或许大家可以帮帮你。

我看了官网的那个脚本，它是画风场的，而不是提取风数据的呢，而且也没有固定经纬度设置，只有固定高度设置，不过也无法提取呢，还有很多错误的呀

老郭爱骑车

★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

virl

{:lxm_20:}{:lxm_20:}{:lxm_20:}{:lxm_20:}

★ESPN★

小__天 发表于 2014-3-4 11:13
ARWpost处理wrfout也行啊，很简单

怎么处理呢？我是提取数据，把数据提取出来放到Excel里

帅帅帅宇

目前我也是遇到了提取这个方面的问题，楼主解决了吗？

中尺度风速模拟

请问楼主解决了吗？

千裡快哉偑

请问楼主解决了吗？求指导

hahahahaha

老郭爱骑车 发表于 2014-4-12 14:16
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;本段代码用来显示指定的经纬度在计算网格中最近的点的编号和 ...

你好，我用你的代码，然后出现了如下错误，错误提示是t=T2(:,locationij(1)-1,locationij(0)-1)这一行，请问如何让解决呢？
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 527.ncl

abc258886@

读取WRF数据有没有教程？