ARWpost运行完毕，却不知如何提取slp，求高手赐教！

敏芬

跟WRF打交道已有一年，最近才想起要用slp数据，但是只知道生成的dat文件，却不知如何提取？如果能够赐教，感激不尽～

夏朗的芒果

wrfout里面是没有slp这个变量的。
可以用ncl或arwpost插值出来。

敏芬

夏朗的芒果 发表于 2014-12-2 12:37
wrfout里面是没有slp这个变量的。
可以用ncl或arwpost插值出来。

首先谢谢你的回答。我已经用ARWpost处理过了，生成了后缀为.dat的文件，不知后面怎么处理了，是需要转化成nc读取还是直接读取？

夏朗的芒果

敏芬 发表于 2014-12-2 12:52
首先谢谢你的回答。我已经用ARWpost处理过了，生成了后缀为.dat的文件，不知后面怎么处理了，是需要转化 ...

转化后生成两个文件，一个是.ctl，还是一个是.dat。你看看生成的ctl文件变量表里面有没有slp？
如果没有的话，在namelist.ARWpost修改添加进去。
如果有的话，就可以直接用grads 读取画图了。

敏芬

夏朗的芒果 发表于 2014-12-2 13:10
转化后生成两个文件，一个是.ctl，还是一个是.dat。你看看生成的ctl文件变量表里面有没有slp？
如果没有 ...

有。不过我想把它读取出来，格式是经度 纬度 slp 三列的。

夏朗的芒果

敏芬 发表于 2014-12-2 14:39
有。不过我想把它读取出来，格式是经度 纬度 slp 三列的。

不好意思，现在才弄懂你的意思。
你可以现在namelist.arwpost中只指定生成经度，纬度和slp。

然后利用fortran读取就可以了啊。fortran读取二进制文件论坛里很多的。

敏芬

夏朗的芒果 发表于 2014-12-2 15:48
不好意思，现在才弄懂你的意思。
你可以现在namelist.arwpost中只指定生成经度，纬度和slp。

嗯嗯，我去试试，谢谢先啦

qhhzxcmm

夏朗的芒果 发表于 2014-12-2 12:37
wrfout里面是没有slp这个变量的。
可以用ncl或arwpost插值出来。

能不能说的详细点，我用ncdump看了一下，确实没有slp这个变量，请问要怎么提取

qhhzxcmm

夏朗的芒果 发表于 2014-12-2 12:37
wrfout里面是没有slp这个变量的。
可以用ncl或arwpost插值出来。

麻烦问一下，已经生成的wrfout里面没有slp,vor这个变量，看了namelist.ARWpost,但是不知道具体怎么插值出来，

deneb1979

ncl可以直接计算ARW里面的slp:
slp = wrf_user_getvar(a,"slp",time)  ; slp