MeteoInfoLab脚本示例：计算温度平流

river

MeteoInfo 发表于 2015-7-19 16:44
谢谢！找到问题所在了，已在1楼更新。

1、高度层是第一层（1000 hPa），所以 lidx = 0。（之前为：idx ...

我也更新一个和王老师一样的最终版本吧。脚本如下：
'reinit'
'set grid off'
'set grads off'
'sdfopen G:\ncep\2011\air.2011.nc'
'sdfopen G:\ncep\2011\uwnd.2011.nc'
'sdfopen G:\ncep\2011\vwnd.2011.nc'
'set t 174'
'set lev 1000'
'pi=3.1416'
'dtx = cdiff(air, 'x')'
'dty = cdiff(air, 'y')'
'dx = cdiff(lon, x) * pi/ 180'
'dy = cdiff(lat, y) * pi/ 180'
'tadv = -1*((uwnd.2*dtx)/(cos(lat*pi/180)*dx)+vwnd.3*dty/dy)/6.37e6'
'run G:\GrADS-NCL-colormaker-v1.1\output\8colors.gs'
'set gxout shaded'
'd tadv'
'run cbar_interp.gs 1 1 0'
'draw title Temporature advection(2011-06-23)'
'gxprint g:/adtv.png white'
;


最终出图如下：

carlye2002

river 发表于 2015-7-20 21:50
我也更新一个和王老师一样的最终版本吧。脚本如下：
'reinit'
'set grid off'

问一下温度平流的单位是什么？k/s吗

river

carlye2002 发表于 2015-10-7 15:52
问一下温度平流的单位是什么？k/s吗

没仔细看，应该是K/s，所有物理量应该都是国际单位的

3870383

@MeteoInfo 王老师你好，为什么这个脚本我复制在脚本编辑器后，提示
Open data files...
Traceback (most recent call last):
  File "<iostream>", line 2, in <module>
NameError: name 'addfile' is not defined

MeteoInfo

3870383 发表于 2016-1-21 15:09
@MeteoInfo 王老师你好，为什么这个脚本我复制在脚本编辑器后，提示
Open data files...
Traceback (most ...

要在MeteoInfoLab中运行。

3870383

MeteoInfo 发表于 2016-1-21 19:16
要在MeteoInfoLab中运行。

可以了，谢谢~~添麻烦了

好多鱼

很好的帖子，留下来了。

好多鱼

很好的帖子，留下来了。

dlmuqixiang

多谢楼主分享

1638672281

老师，我这边是打算用python做的，但是cdiff这个函数您是怎么实现的