用Grads利用nc资料话画温度平流

Eegle · 发表于 2014-12-23 16:40:37

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

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

x

以下的gs文件是我在论坛里找的方法，可是为啥画出的图，跟别人文章里面的图是不一样的呢？（左边是我自己画的，右边是别人文章的）
请大家帮我看看。谢谢了

登录/注册后可看大图

'reinit'
'enable print e:\zhuang\wendu\T_adv_25N_1608.gmf'
'sdfopen E:\zhuang\air.2014.nc'
'sdfopen E:\zhuang\uwnd.2014.nc'
'sdfopen E:\zhuang\vwnd.2014.nc'
'set lev 1000 100'
'set t 665'
'define pi=3.1415726/180.0'
'define vor=hcurl(uwnd.2,vwnd.3)'
'define r=6.356e6'
'define dtx = cdiff(air.1,x)'
'define dty = cdiff(air.1,y)'
'define dx=cdiff(lon,x)*cos(lat*pi/180)*pi*r/180'
'define dy=cdiff(lat,y)*pi*r'
'define advT=-(uwnd.2*dtx/dx+vwnd.3*dty/dy)'

'set grads off'
'set grid off'
'set csmooth on'
'set map 15 1 9'
'set parea 1 10.2 1.5 8'
'set xlopts 1 6 0.14'
'set ylopts 1 6 0.14'
'set font 0'
'set lon 112 122'
'set lat 25'
'set lev 1000 100'
'set zlog on'
'set gxout contour'
'set cthick 9'
'set cint 0.5'
'd advT*1e3'

nzy718d · 发表于 2015-1-12 06:45:37

'define pi=3.1415726/180.0'
cos(lat*pi/180)*

是不是重复了，/180?

笑哈乐呵 · 发表于 2015-4-16 16:36:38

挣不到钱

听雨残荷 · 发表于 2015-4-17 17:21:33

'define dx=cdiff(lon,x)*cos(lat*pi/180)*pi*r/180'这句中cos(lat*pi/180）是不是多除了一个180？？

djcj520 · 发表于 2016-6-22 19:26:31

好，谢谢分享

djcj520 · 发表于 2016-6-22 19:26:57

有fnl资料做的吗

djcj520 · 发表于 2016-6-22 19:27:25

现在比较多用fnl

李李李李李 · 发表于 2017-4-24 15:58:35

楼主，请问这里为什么要定义一个vor呢

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

[秀图] 用Grads利用nc资料话画温度平流

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

相关帖子