爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 8190|回复: 4

分享一个ncl处理wrfout,计算P,rh,T等气象元素的脚本

[复制链接]
发表于 2017-8-12 17:04:16 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

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

x
nb1 = addfile("/home/alex/test/wrfout_d01_2015-07-20_00:00:00","r")

;TK为温度,P为压强
P = nb1->P
pb = nb1->PB
P = P + pb

T = nb1->T
T = T+300
TK = wrf_tk(P,T)


;RH为相对湿度
qv = nb1->QVAPOR
RH = wrf_rh(qv,P,TK)

;位势高
PH = nb1->PH
PHB=nb1->PHB
PH=PH+PHB
z=wrf_user_unstagger(PH,PH@stagger)
z=z/9.81

fout = addfile("输出nc文件的完整路径","c")
lat = nb1->XLAT
lon = nb1->XLONG
times=nb1->Times
fout->lat = lat
fout->lon = lon
fout->time =times
fout->P=P
fout->T=T
fout->RH=RH
fout->hgt = z
密码修改失败请联系微信:mofangbao
发表于 2017-8-12 20:51:51 | 显示全部楼层
嘛,怎么说呢……
先说个“谢谢分享”吧



我本身是不用NCL的,所以后面的就当我是在鬼扯了

为什么不用  wrf_user_getvar  系列呢?

是为了提示自己 wrfout 文件里变量“T”是减掉定值300K的么?
是为了提示自己有 wrf_tk wrf_rh 这样的函数么?
(z印象中似乎也是直接会算的,但不确定了,想想还是不喷了)

我觉着知道这些函数是怎么算的,挺好;
即使不知道怎么算的,知道有这些函数,知道这些函数的输入输出,也挺好;
但正常用的话有封装好的还是用封装好的,里面不也是这么调用的么
密码修改失败请联系微信:mofangbao
回复 支持 0 反对 1

使用道具 举报

发表于 2017-8-13 10:43:57 | 显示全部楼层
现在好多气象台的人都用ncl啊
密码修改失败请联系微信:mofangbao
发表于 2017-8-18 18:45:24 | 显示全部楼层
赞一个!!!
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

发表于 2025-1-23 17:39:40 | 显示全部楼层
楼主有python处理的脚本吗?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表