grads计算散度注意事项及其它

'reinit'

'sdfopen e:/data/ncep/data/uwnd.nc'

'sdfopen e:/data/ncep/data/vwnd.nc'

'set gxout fwrite'

'set fwrite e:/data/ncep/data/divmav.grd'

'set lev 1000'

k=1967

while(k<=1998)

'set lon 10 240'

'set lat 20 50'

'define a=ave(uwnd.1,time=01may'k',time=01sep'k')'

'define b=ave(vwnd.2,time=01may'k',time=01sep'k')'

'define hh=hdivg(a,b)'

'set lon 97.5 102.5'

'set lat 37.5 42.5'

'd hh'

k=k+1

endwhile

'disable fwrite'

;

注意：这是一个求散度的grads程序，蓝色为最终要求的区域(的散度)，注意其位置，红色为自行设定的区域，考虑到中央差的特性，这个自行设定的区域必须大于蓝色范围。

 

此外grads 中time表示具体的年月日时间，t只表示时次。

如

k=1950

while（k<=1999)

'define a=ave(uwnd, time=01may'k', time=01sep'k')
'd a'

k=k+1

endwhile

此处的time不能写成t.