GrADS计算湿位涡全套程序

champyuan · 发表于 2014-11-9 17:57:53

这个可是好东西

chuanchuan · 发表于 2014-11-11 14:56:56

感谢分享！

萍水相逢 · 发表于 2014-11-28 09:29:39

正想学习这个，谢谢分享，试试

newbluesea · 发表于 2014-12-4 23:05:17

非常感谢楼主

Jaxon · 发表于 2014-12-18 09:28:25

请问楼主，我这样对高度进行循环行不行？个人觉得比单独每层简单（我有27层）
j = 2
while(j <= 26)
'set lon 70 140'
'set lat 10 60'
'set z 'j-1''
'define lev2=lev'
'define u2=u'
'define v2=v'
'define thse2=thse'
'set z 'j+1''
'define lev3=lev'
'define u3=u'
'define v3=v'
'define thse3=thse'
'set z 'j''
'define lev1=lev'
'define u1=u'
'define v1=v'
'define thse1=thse'

'define vo=hcurl(u1,v1)'
'define f=2*7.292*sin(lat*3.14159/180.0)*0.00001'
'define g=9.8'
'define dp=100.0*(lev2-lev3)'
'define dthse=thse2-thse3'
'define du=u2-u3'
'define dv=v2-v3'
'define dx=2.0*6370949.0*cos(lat*3.14159/180.0)*3.14159/180.0'
'define dy=2.0*6370949.0*3.14159/180.0'
'define dtx=cdiff(thse,x)'
'define dty=cdiff(thse,y)'
'define mpv1=-g*(vo+f)*dthse/dp'
'define mpv2=g*((dv/dp)*(dtx/dx)-(du/dp)*(dty/dy))'
'define mpv=mpv1+mpv2'
'd mpv1*1000000'

小怪兽 · 发表于 2014-12-24 17:16:06

谢谢分享，感激不尽。

风浪 · 发表于 2015-1-21 11:41:35

楼主，你太好了，感谢分享！！！

slqx · 发表于 2015-2-3 08:54:01

非常感谢分享着和楼主！

nzy718d · 发表于 2015-2-5 17:36:19

请教老师：按照您的方法做垂直差分，类似du=u（lev=700）-u（lev=850）怎么等于0?也就是u（lev=700）与（lev=850）的值是一样的。我分别显示了u（lev=700）与u（lev=850）的图形，的确是一样的，dp等等也是一样的，均显示的是1000hpa的数值，请问这是怎么回事啊？谢谢！

13240848 · 发表于 2015-2-6 11:30:19

怎么需要金币？

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

[分享资料] GrADS计算湿位涡全套程序