后使用快捷导航没有帐号？
扫一扫，访问微社区
只需一步, 快速开始
只需一步，快速开始
查看: 586|回复: 8
wrfout数据绘高分辨率三维地形图
在线时间 小时
登录后查看更多精彩内容~
才可以下载或查看，没有帐号？
一直想要画三维的地形图，因为是搞中尺度的，要求分辨率比较高，网上下的数据画不出想要的高分辨率三维地形图，前几天在论坛里看到了这篇matlab的帖子&&&&根据他的思路弄了一下，画出来的海南地形图可以用了，就想着分享到论坛里，给遇到同样困扰的朋友参考一下。
先用grads2ascii将wrfout中的地形数据提取出来（可参见兰溪的）
'open e:\7.ctl'
'grads2ascii hgt e:\hgt\hgt.txt %8.0f 375 1'
然后读入matlab中绘图(具体程序见hgt.m)
09:54 上传
点击文件名下载附件
下载积分: 金钱 -5
959 Bytes, 下载次数: 14, 下载积分: 金钱 -5
给自己设置一个个性的签名吧~
在线时间 小时
顶一下。谢谢分享
给自己设置一个个性的签名吧~
在线时间 小时
A small hint:
&& load topo
and there are two colormaps (topomap1, topomap2) designed for topography. Land and ocean in your fancy stuff are not perfectly distinguished.
给自己设置一个个性的签名吧~
在线时间 小时
谢谢分享& &
给自己设置一个个性的签名吧~
在线时间 小时
漂亮，楼主，学习学习哈
给自己设置一个个性的签名吧~
在线时间 小时
谢谢楼主分享
给自己设置一个个性的签名吧~
在线时间 小时
这张可能看得更清楚
给自己设置一个个性的签名吧~
在线时间 小时
涛哥，你好棒
给自己设置一个个性的签名吧~
在线时间 小时
给自己设置一个个性的签名吧~后使用快捷导航没有帐号？
扫一扫，访问微社区
只需一步, 快速开始
只需一步，快速开始
查看: 1895|回复: 22
请教只用用ncl画wrfout的散度，涡度甚至湿位涡
在线时间 小时
登录后查看更多精彩内容~
才可以下载或查看，没有帐号？
求教用ncl直接处理wrfout数据画速度散度，涡度以及湿位涡？请大家给个建议吧..非常感谢
已有1人关注本帖
给自己设置一个个性的签名吧~
在线时间 小时
呵呵，帮你顶了
给自己设置一个个性的签名吧~
在线时间 小时
随缘 发表于
呵呵，帮你顶了
哈哈，谢谢~
给自己设置一个个性的签名吧~
在线时间 小时
prodigyzl 发表于
哈哈，谢谢~
官网没有脚本的例子么？
给自己设置一个个性的签名吧~
在线时间 小时
官网上有脚本的
给自己设置一个个性的签名吧~
在线时间 小时
薇风拂尘 发表于
官网上有脚本的
官网上我只看到涡度和螺旋度，散度好像没有看到。。还是我没有找到？
给自己设置一个个性的签名吧~
在线时间 小时
楼主找到MPV脚本没？
给自己设置一个个性的签名吧~
在线时间 小时
报错，没有找到表文件是啥意思，是我与.ncl程序文件相应的文件没考全，还是什么
给自己设置一个个性的签名吧~
在线时间 小时
散度的没看见
给自己设置一个个性的签名吧~
在线时间 小时
报错，没有找到表文件是啥意思，是我与.ncl程序文件相应的文件没考全，还是什么
没load上函数吧
给自己设置一个个性的签名吧~后使用快捷导航没有帐号？
扫一扫，访问微社区
只需一步, 快速开始
只需一步，快速开始
查看: 1189|回复: 13
grads画图时如何选择wrfout的输出变量
在线时间 小时
登录后查看更多精彩内容~
才可以下载或查看，没有帐号？
wrfout文件用ARWpost处理后，生成dat和ctl文件，用grads画图时，查看了下ctl文件的内容发现有100个变量，
请问我该如何找到我要的那个变量啊？
请各位大虾指点一下，谢谢啦！
给自己设置一个个性的签名吧~
在线时间 小时
自己一个个的读vars呗，反正都有说明，论坛也专门有帖子讲wrf的输出变量。
你要是嫌多就再namelist.arwpost里面用list控制一下输出不就好了。
给自己设置一个个性的签名吧~
在线时间 小时
100个变量呢，我原来不是学气象的，刚刚接触这个，
每个变量后面的说明我很多都不太明白是什么意思，
请问你在哪看到有ctl文件里变量的名称的啊
给自己设置一个个性的签名吧~
在线时间 小时
bailing1175 发表于
100个变量呢，我原来不是学气象的，刚刚接触这个，
每个变量后面的说明我很多都不太明白是什么意思，
那就是说你不会用grads？这样的话我就没什么可说的了，你要先好好看看手册，多在论坛搜索。
ctl的vars后面都是变量，看不懂说明的话那是英语水平问题，这个没人能帮的。
ps：你要点回复，我才看得到
给自己设置一个个性的签名吧~
在线时间 小时
lqouc 发表于
那就是说你不会用grads？这样的话我就没什么可说的了，你要先好好看看手册，多在论坛搜索。
ctl的vars后 ...
非常感谢！
你说的手册是哪种手册啊？
给自己设置一个个性的签名吧~
在线时间 小时
bailing1175 发表于
非常感谢！
你说的手册是哪种手册啊？
lasg的手册，全称叫GRADS使用手册，论坛资料站里面有，非常简明易懂。想多学点的话可以再找找马红云和朱禾的书。
给自己设置一个个性的签名吧~
在线时间 小时
lqouc 发表于
lasg的手册，全称叫GRADS使用手册，论坛资料站里面有，非常简明易懂。想多学点的话可以再找找马红云和朱禾 ...
知道了，可那本书上我没看到有介绍怎么选wrfout输出变量的啊，谢谢
给自己设置一个个性的签名吧~
在线时间 小时
bailing1175 发表于
知道了，可那本书上我没看到有介绍怎么选wrfout输出变量的啊，谢谢
不明白你说的什么叫怎么选，你想用哪个就用哪个啊，找到变量名用就可以了。
wrfout的变量介绍论坛有帖子可以搜一下。
给自己设置一个个性的签名吧~
在线时间 小时
ctl 的描述中，vars这行之后的都是为变量，每行一个，每个后面有英文描述，你自己翻译着找呗。
给自己设置一个个性的签名吧~
在线时间 小时
做个霸气的木头 发表于
ctl 的描述中，vars这行之后的都是为变量，每行一个，每个后面有英文描述，你自己翻译着找呗。
我英语弱，用这个方法太慢了，所以想看看有没有简单、快速一点的办法。谢谢啦！！
给自己设置一个个性的签名吧~后使用快捷导航没有帐号？
扫一扫，访问微社区
只需一步, 快速开始
只需一步，快速开始
查看: 739|回复: 8
wrfout转换为ctl和dat之后grads画阴影图所有数据都没有
在线时间 小时
自己思考时长(天):
我看过提问的智慧:
ncdump wrfout都是有数据的，我用了wrf2grads和arwpostv2 v3都进行了转换，用grads2.1家园版画图的结果都如图，操作命令就是set t 25；set gxout shaded；d rainnc。这里强调下，grads应该无问题，我画以前留存的数据可以，我后来换了大型机重新安装的arwpost和wrf2grads，画新转换的ctl和dat就不行。有无方法看dat文件如何看里面内容对不对？还有就是我是在大型机上run的数据，转换之后拷贝到本地，再用grads进行画图。我第一次发帖，已经在网上找了很多信息，有类似的问题但是都没有答案，特来求助。
aix； windows
登录后查看更多精彩内容~
才可以下载或查看，没有帐号？
grads输完命令后显示contouring：-2e+38 to 2.5e+38 interval 5e+37
arwpost转换后的ctl文件：
dset D:/wrfoutctr.dat
undef 1.e30
title&&OUTPUT FROM WRF V2.2.1 MODEL
pdef&&199 179 lcc&&36.000&&117.500&&100.000& &90.000&&60.000&&117.5.000&&
xdef&&586 linear&&104.29045& &0.
ydef&&394 linear& &26.93918& &0.
zdef& &30 linear 1 1&&
tdef& &25 linear 06Z18SEP2005& && &60MN& && &
RAINC& && && & 1&&0&&ACCUMULATED TOTAL CUMULUS PRECIPITATION (mm)
RAINNC& && && &1&&0&&ACCUMULATED TOTAL GRID SCALE PRECIPITATION (mm)
给自己设置一个个性的签名吧~
在线时间 小时
RE: wrfout转换为ctl和dat之后grads画阴影图所有数据都没有
补充下namelist.arw:
start_date = '_06:00:00',
end_date& &= '_06:00:00',
interval_seconds = 3600,
debug_level = 0,
io_form_input&&= 2,
input_root_name = '/cmb/g5/shaochl/WRFV2/WRFV2/run/wrfout_d01__06v2'
output_root_name = '/cmb/g5/shaochl/wrfoutctr'
plot = 'list'
fields = 'RAINC,RAINNC'
output_type = 'grads'
mercator_defs = .false.
! split_output = .false.
! frames_per_outfile = 2
! output_type = 'grads'
! output_type = 'v5d'
! plot = 'all'
! plot = 'list'
! plot = 'all_list'
! Below is a list of all available diagnostics
! fields = 'height,geopt,theta,tc,tk,td,td2,rh,rh2,umet,vmet,pressure,u10m,v10m,wdir,wspd,wd10,ws10,slp,mcape,mcin,lcl,lfc,cape,cin,dbz,max_dbz,clfr'
interp_method = 0,
! interp_levels = .,900.,850.,800.,750.,700.,650.,600.,550.,500.,450.,400.,350.,300.,250.,200.,150.,100.,
!extrapolate = .false.
! interp_method = 0,& &&&! 0 is model levels, -1 is nice height levels, 1 is user specified pressure/height
! interp_levels = .,900.,850.,800.,750.,700.,650.,600.,550.,500.,450.,400.,350.,300.,250.,200.,150.,100.,
! interp_levels = 0.25, 0.50, 0.75, 1.00, 2.00, 3.00, 4.00, 5.00, 6.00, 7.00, 8.00, 9.00, 10.0, 11.0, 12.0, 13.0, 14.0, 15.0, 16.0, 17.0, 18.0, 19.0, 20.0,
给自己设置一个个性的签名吧~
在线时间 小时
&成长值: 690
options byteswapped
给自己设置一个个性的签名吧~
在线时间 小时
&成长值: 7365
你可以看看里面的数据到底有没有
给自己设置一个个性的签名吧~
在线时间 小时
本帖最后由 大侠MVP 于
14:52 编辑
options byteswapped
我用的ibm的xlf编译器，ctl里本来就没有这项。而且我也试过了把这个加上，输出还是一样
给自己设置一个个性的签名吧~
在线时间 小时
你可以看看里面的数据到底有没有
dat里面的数据吗？我只能用写字板打开，里面都是非零的，但是不知道具体数字是多少，有其他方法看吗？用set gxout contour画出来的图满屏幕都是很小很小孤立的圈，用shaded输出的就是我传上来这张图
给自己设置一个个性的签名吧~
在线时间 小时
我也遇到了同样的问题，直接用ncl看wrfout的数据时有的，就是不知道为什么画出来就是一个一个小圈？
给自己设置一个个性的签名吧~
在线时间 小时
options byteswapped
感谢兰溪。。按照你的提示问题已经得到解决
给自己设置一个个性的签名吧~
在线时间 小时
看你的网格中心位置，你是河北的，还是山东的？
给自己设置一个个性的签名吧~}