温度平流GS

巫毒公主

新人初来咋到，跟大家分享一段grads中画温度平流的gs，木有进行图形美化啥的，就是基本的计算方法。欢迎大家批评指正！

1. 'reinit'
   
2. 'open F:/MPV/DATA/57.ctl'
   
3. 'set t 1'
   
4. 'set lon 90 125'
   
5. 'set lat 10 35'
   
6. 'set lev 1000 200'
   
7. 'define pi=3.1415726/180.0'
   
8. 'define vor=hcurl(ugrdprs,vgrdprs)'
   
9. 'define r=6.356e6'
   
10. 'define dtx = cdiff(tmpprs,x)'
    
11. 'define dty = cdiff(tmpprs,y)'
    
12. 'define dx=cdiff(lon,x)*cos(lat*pi/180)*pi*r/180'
    
13. 'define dy=cdiff(lat,y)*pi*r'
    
14. 'define advT=-(ugrdprs*dtx/dx+vgrdprs*dty/dy)'
    
15. 'set lev 850'
    
16. 'set xlint 5'
    
17. 'set ylint 5'
    
18. *'set cint 5'
    
19. 'set mpdset cnworld guangdong guangdong_q'
    
20. 'd advT*1e3'
    

复制代码

xsjvn_

楼主，您给的程序有点问题，给出个人意见如下：
define pi = 3.1416
define r = 6.37e6
define dtx = cdiff(t,x)
define dty = cdiff(t,y)
define dx = cdiff(lon,x)*pi/180
define dy = cdiff(lat,y)*pi/180
define advT = -1*( (u*dtx)/(cos(lat*pi/180)*dx) + v*dty/dy )/r
首先Pi不应除以180 ，其次dy 没有除以180
仅供参考！

善人/jw

感谢分享~~谢谢LZ
希望多多分享啊

传说中的谁

哇grads专区来新人了，看起来还狠萌。感谢分享啦，以后有好东西别忘了发到家园来。

尽头的尽头

非常感谢啊，顶起来

言深深

啊，风景如此多娇啊！气象家园欢迎原创，多多益善多多益善

mofangbao

希望家园这样的程序越来越多啦

zsmnantong

越来越多啦！感谢分享，欢迎新朋友！

zhaolanxiai

呵呵，学习了！

喜上眉梢

学用一下，感谢分享！

巫毒公主

@善人/jw@传说中的谁@言深深@mofangbao@artmunich矮油，你们给俺那么多钱，多不好意思呀