- 积分
- 147
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-4-29
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
本人画出了等熵位涡在设定等压面后的剖面图 现在想画等熵位涡南北向的剖面图 修改后出来的全是一条一条直线 求教啊
*等熵位涡分析
'reinit'
'open f:/zhang/fnl_20130710_18_00.00.ctl'
'set lon 80 130'
'set lat 15 45'
'set lev 1000 100'
"define t=TMPprs"
"define u=ugrdprs"
"define v=vgrdprs"
"define pp=lev"
"define coriol=2*7.29e-5*sin(lat*3.1415/180)"
"define dudy=cdiff(u,y)/(111177*cdiff(lat,y))"
"define dvdx=cdiff(v,x)/(111177*cdiff(lon,x)*cos(lat*3.1415/180))"
"define dt=t(z-1)*pow(1000/PP(z-1),0.286)-t(z+1)*pow(1000/PP(z+1),0.286)"
"define dp=100*(PP(z-1)-PP(z+1))"
"define dtdp=dt/dp"
"define part1="isen(dvdx,t,PP,320)
"define part2="isen(dudy,t,PP,320)
"define part3="isen(dtdp,t,PP,320)
"define pv320=-9.8*(coriol+part1-part2)*part3"
"set z 0"
'set lon 80 130'
'set lat 15 40'
"d pv320*1000000"
'printim f:/zhang/1000ipv1018.jpg x1024 y768 white'
按这个可以出来图 可是修改SET Z 和 SET LON 之后 就变成了一条一条直线
|
-
1000HPa的等熵位涡
-
修改 SET LEVE后
|