请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 27794|回复: 37

[分享资料] 整层水汽通量及散度的计算

[复制链接]

新浪微博达人勋

发表于 2013-9-9 12:10:57 | 显示全部楼层 |阅读模式

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

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

x
以下分别计算出各个层次的qu,qv,hd,其中所用的数据是1948年1月到2013年1月的月平均年数据,即共781月的数据
'reinit'
'sdfopen e:\1year\air.mon.mean.nc'
'sdfopen e:\1year\rhum.mon.mean.nc'
'sdfopen e:\1year\uwnd.mon.mean.nc'
'sdfopen e:\1year\vwnd.mon.mean.nc'
'set fwrite e:\1year\1flux.dat'
'set gxout fwrite'
nt=1
'set t 'nt
while(nt<=781)
nz=1
while(nz<=8)
'set z 'nz
'define p=lev'
'define es=6.112*exp(17.67*air/(air+273.15-29.65))'
'define qs=0.622*es/(p-0.378*es)'  
'define q=rhum.2*qs/100'
'define qu=(uwnd.3*q/9.8)'
'd qu'
nz=nz+1
endwhile
nt=nt+1
endwhile
nt=1
'set t 'nt
while(nt<=781)
nz=1
while(nz<=8)
'set z 'nz
'define p=lev'
'define es=6.112*exp(17.67*air/(air+273.15-29.65))'
'define qs=0.622*es/(p-0.378*es)'
'define q=rhum.2*qs/100'
'define qv=(vwnd.4*q/9.8)'
'd qv'
nz=nz+1
endwhile
nt=nt+1
endwhile
nt=1
'set t 'nt
while(nt<=781)
nz=1
while(nz<=8)
'set z 'nz
'define p=lev'
'define es=6.112*exp(17.67*air/(air+273.15-29.65))'
'define qs=0.622*es/(p-0.378*es)'
'define q=rhum.2*qs/100'
'define qu=(uwnd.3*q/9.8)'
'define qv=(vwnd.4*q/9.8)'
'define hd=hdivg(qu,qv)'
'd hd'
nz=nz+1
endwhile
nt=nt+1
endwhile
'disable fwrite'
'reinit'
;
以下是进行整层计算前的ctl文件
dset e:/1year/1flux.dat
title vertically integrated moisture flux
UNDEF -9.99E33
xdef 144 LINEAR   0 2.5
ydef 73 LINEAR   -90 2.5
zdef 8 LEVELS 1000 925 850 700 600 500 400 300
tdef 781 LINEAR JAN1948 1mo
vars 3
qu 1 99 u_direction
qv 1 99 v_direction
hd 1 99 horizontal divergence
endvars
以下为整层计算
'reinit'
'open e:\1year\1fluxdata.ctl'
'sdfopen e:\1year\pres.sfc.mon.mean.nc'
'set fwrite e:\1year\flux_all.dat'
'set gxout fwrite'
'set z 1'
'define quall=vint(pres.2,qu,300)*9.8'
'define qvall=vint(pres.2,qv,300)*9.8'
'd quall'
'd qvall'
'define hdall=hdivg(quall,qvall)'
'd hdall'
'disable fwrite'
'reinit'
;
以下为整层文件所写的ctl
dset e:/1year/flux_all.dat
title vertically integrated moisture
UNDEF -9.99E33
xdef 144 LINEAR   0 2.5
ydef 73 LINEAR   -90 2.5
zdef 8 LEVELS 1000 925 850 700 600 500 400 300
tdef 781 LINEAR JAN1948 1mo
vars 3
quall 1 99 u_direction
qvall 1 99 v_direction
hdall 1 99 horizontal divergence
endvars
玛法麻烦大家帮我看下上面的程序哪有问题啊?当我最后单独绘制quall,qvall, hdall是发现整个图形一个颜色,期待高手指点啊,在线等待中。。。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-9-9 13:10:30 | 显示全部楼层
顶一下,希望大家帮我看下我的错误啊
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-9-9 13:11:12 | 显示全部楼层
自顶一下,希望大家帮我看下我的错误啊
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-9-9 14:17:43 | 显示全部楼层
。。。。。。。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-9-9 19:23:05 | 显示全部楼层
帮顶
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2013-9-10 10:03:58 | 显示全部楼层
为何没人关注我的问题啊,太长了吗?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-10-18 12:51:30 | 显示全部楼层
老哥,时间应该是最外层循环,你怎么每个量都在循环呀!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2013-10-19 14:55:49 | 显示全部楼层
前面提取qu,qv的那些数据你检验过可以出图不没?还有你最后画图的时候的gs文件贴出来一下赛
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-6-3 22:43:00 | 显示全部楼层
我也在学习,看不懂你这个
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2014-6-21 12:45:45 | 显示全部楼层
帮顶!!!求大神解释
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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