温度平流GS

兰北

是的，后面多了一个除以180,正确的应该为：
'reinit'

'open F:/MPV/DATA/57.ctl'

'set t 1'

'set lon 90 125'

'set lat 10 35'

'set lev 1000 200'

'define pi=3.1415726'

'define vor=hcurl(ugrdprs,vgrdprs)'

'define r=6.356e6'

'define dtx = cdiff(tmpprs,x)'

'define dty = cdiff(tmpprs,y)'

'define dx=cdiff(lon,x)*cos(lat*pi/180)*pi*r/180'

'define dy=cdiff(lat,y)*pi*r/180'

'define advT=-(ugrdprs*dtx/dx+vgrdprs*dty/dy)'

'set lev 850'

'set xlint 5'

'set ylint 5'

*'set cint 5'

'set mpdset cnworld guangdong guangdong_q'

'd advT*1e3'

水每桥

如果pi是这样定义的话，dx的计算就错了！

跃宝妈咪

xiexie

最后一个匈奴

感谢楼主无私分享！！！

逝者如斯

正在试用，看看成果在进一步回复啊

洋咩咩

程序中的vor好像没有用到啊~

OomaulenOo

谢谢分享

极地冷高压

非常感谢啊，学习学习

夏日炎炎

好像多除以了一个180

猫猫monty

学习了！谢谢