- 积分
- 54
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-4-27
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
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
|
|