提取某纬度的位温值出现问题

kanyu · 发表于 2015-3-3 17:13:10

您需要登录才可以下载或查看，没有帐号？立即注册

x

需要提取某一个纬度的位温值，但提取过程中出现问题，代码及问题如下
代码：
lat = f->XLAT
ind= ind(lat .eq. 31.2052)
t    = f->T(0,:,ind,:)

问题：

longlivehj · 发表于 2015-3-3 17:35:02

lat是多维的，不能用ind去找。楼主这种情况，可以直接用下面的语句：
t = f->T(0,:,{31.2052},:)

kanyu · 发表于 2015-3-3 19:15:04

longlivehj 发表于 2015-3-3 17:35
lat是多维的，不能用ind去找。楼主这种情况，可以直接用下面的语句：
t = f->T(0,:,{31.2052},:)

我这样修改之后出现这样的问题：
fatal:Dimension (south_north) of file (wrfout_d02_2013-09-13_08) does not have an associated coordinate variable

longlivehj · 发表于 2015-3-3 21:11:20

kanyu 发表于 2015-3-3 19:15
我这样修改之后出现这样的问题：
fatal:Dimension (south_north) of file (wrfout_d02_2013-09-13_08) d ...

不好意思，是我弄错了。
可以考虑先用rcm2rgrid转换一下。

longlivehj · 发表于 2015-3-3 21:11:20

kanyu 发表于 2015-3-3 19:15
我这样修改之后出现这样的问题：
fatal:Dimension (south_north) of file (wrfout_d02_2013-09-13_08) d ...

不好意思，是我弄错了。
可以考虑先用rcm2rgrid转换一下。

超台十八 · 发表于 2015-3-4 06:40:10

kanyu 发表于 2015-3-3 19:15
我这样修改之后出现这样的问题：
fatal:Dimension (south_north) of file (wrfout_d02_2013-09-13_08) d ...

资料不大的话，我整个先读出来 t = f->T

再定義 t 的纬度。

kanyu · 发表于 2015-3-4 09:11:51

超台十八发表于 2015-3-4 06:40
资料不大的话，我整个先读出来 t = f->T

再定義 t 的纬度。

可以具体说一下吗？

kanyu · 发表于 2015-3-4 09:12:36

longlivehj 发表于 2015-3-3 21:11
不好意思，是我弄错了。
可以考虑先用rcm2rgrid转换一下。

是将相关变量都转换一下还是只转换纬度的

longlivehj · 发表于 2015-3-4 09:30:13

kanyu 发表于 2015-3-4 09:12
是将相关变量都转换一下还是只转换纬度的

wrf坐标是curvilinear的，就是同一列上格点的经度是不一样的，所以没有办法直接为变量设置坐标变量。用rcm2rgrid_Wrap将变量插值后，就拥有rectilinear的坐标变量，取某一纬度就非常方便了。
当然，也可以取一定范围的纬度，官网有例子的。

kanyu · 发表于 2015-3-4 10:01:43

longlivehj 发表于 2015-3-4 09:30
wrf坐标是curvilinear的，就是同一列上格点的经度是不一样的，所以没有办法直接为变量设置坐标变量。用rc ...

我提取这个某一纬度的变量之后，需要画出剖面图，纵坐标是气压，应该怎么做呢？直接用gsn_csm_pres_hgt这个函数可以吗？

		自动登录	找回密码
密码			立即注册

[作图] 提取某纬度的位温值出现问题