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

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 8286|回复: 1

[分享资料] 想弄一个水汽通量的时间剖面,结果居然是常数?

[复制链接]

新浪微博达人勋

发表于 2022-4-20 17:37:58 | 显示全部楼层 |阅读模式

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

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

x
'reinit'
'sdfopen L:/DZQDL/nc/erainterim/u_wind.nc'
'sdfopen L:/DZQDL/nc/erainterim/v_wind.nc'
'sdfopen L:/DZQDL/nc/erainterim/specific_humidity2.nc'
'sdfopen L:/DZQDL/nc/erainterim/temp.nc'
'sdfopen L:/DZQDL/nc/erainterim/z.nc'
'set grid off'
'set grads off'
'set csmooth on'



'set t 21 28'
'set lon 96 108'
'set lat 20 30'


*------------计算比湿  单位:kg/kg--------------*
'define es=(6.11*exp((17.2693882*(t.4-273.16))/(t.4-35.86)))'
'define qs=0.62197*es/(z.5-0.378*es)'
*-------计算水汽通量(单层)  单位:(g*s)/kg=g/(s*hpa*cm)------*
'define qu=u.1*q.3*1e3/9.8'
'define qv=v.2*q.3*1e3/9.8'
*------------计算水汽通量散度------*
'define qdiv=hdivg(qu,qv)*1e5'
'set grads off'
'set parea 1 9.0 2 8'
'set xlopts 1 6 0.2'
'set ylopts 1 6 0.2'
'set ylint 2'
'set xlint 2'
'set t 21 28'
*'set lon 102.04'
*'set lat 25.09'
*'set lon 102.10'
*'set lat 24.40'
'set lon 102.4'
'set lat 24.1'
'set gxout contour'
'd mag(qu,qv)'

'set gxout contour'
'set cthick 6'
'set ccolor 1'
'd qdiv'

'C:\grads\Contents\Cygwin\Versions\2.0.1.oga.1\i686\cbar_matlab 5 1 0 5'
*'run axis.gs -type b -position i -label on -interval 3 -sinterval 2 -lsize 0.2  -voffset -0.2'
*'run axis.gs -type l -position i -label on -interval 100 -lsize 0.2'
'printim l:/DZQDL/2014datafiles/divergence_time/div_t2_sanchahe.wmf white'
'disable gxprint'
*
搜狗截图20220420173454.jpg

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2022-10-15 09:59:50 | 显示全部楼层
我也遇到这个问题了,请问楼主解决了吗
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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