找到一个计算温度平流的gs,不知道有没有用，有待验证，发上来先。

realzly · 发表于 2011-9-26 22:42:10

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

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

x

找到一个计算温度平流的gs,不知道有没有用，有待验证，发上来先。
利用cdiff例子中介绍的温度平流的计算方法，计算温度平流和水汽平流，gs文件如下：
open i:/icbc-era40/IC_P2000070100.ctl
set z 1 7
define dtx = cdiff(t,x)
define dty = cdiff(t,y)
define dx = cdiff(lon,x)*3.1416/180
define dy = cdiff(lat,y)*3.1416/180
dit=-1*( (u*dtx)/(cos(lat*3.1416/180)*dx) + v*dty/dy )/6.37e6
d dit
z从1取到7设置在不同的层（1000hpa到100hpa)
来自：http://blog.sina.com.cn/s/blog_4e1bbc8501000aqq.html

传说中的谁 · 发表于 2011-9-26 22:50:36

貌似没错，但这样写dx那一项看起来物理意义不那么明显，如果把计算纬圈半径一起写进dx那一项里会好一些。我原来发过这个的。

晚晴山 · 发表于 2013-10-16 12:31:00

算出来的值有好几千，单位是什么？

zl.87625 · 发表于 2014-1-4 14:31:57

公式里面的t  是指时间变量吗？我是下面那样算的，但是系统显示不认识t ，怎么解决？是要写成循环吗

'set  parea 1 5.5 2.7 5.2'
'set  grads off'
'set  grid    off'
'set  lat 0 90'
'set  lon  0  360'
'set  xlint    30'
'set  ylint    20'
'set xlab    off'

**************************
'define pi=3.1415726'
'define r=6.356e6'
'define dtx=cdiff(t,x)'
'define dty=cdiff(t,y)'
'define dx=cdiff(lon,x)*cos(lat*pi/180)*pi*r/180'
'define dy=cdiff(lat,y)*pi*r'
'define advT=-(grdprs*dtx/dx+v*dty/dy)'

Legendary · 发表于 2014-4-24 10:44:31

zl.87625 发表于 2014-1-4 14:31
公式里面的t 是指时间变量吗？我是下面那样算的，但是系统显示不认识t ，怎么解决？是要写成循环吗

t显然是温度。

风儿 · 发表于 2016-12-1 10:41:36

把t 改成tmprs应该就可以啦

temperature · 发表于 2016-12-2 12:08:23

有用有用，谢谢楼主

大笨熊 · 发表于 2017-7-2 17:11:05

mark 学习学习

		自动登录	找回密码
密码			立即注册

[分享资料] 找到一个计算温度平流的gs,不知道有没有用，有待验证，发上来先。

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