请问有朋友做相对散度的吗？求指导

ypliu2015

请问有朋友做相对散度的吗？
在单个层次散度输出都是正常的，但是我想做300hpa到850hpa的散度差，出来的就是恒定场
找不到问题出在哪里了？
各位路过的大神就指导！

'reinit'
'open D:\fnl\fnl_060326_h.ctl '
'set lon 90 140'
'set lat 25 65'
'set t 1 16'
'set lev 1000 100'
'set display color white'

'define div=hdivg(ugrdprs,vgrdprs)*1e-05'

'set mpdset cnworld'
'set map 1 1 2'

tt=1
while(tt<17)
'set lon 100 130'
'set lat 35 55'
     'q time'
     res = subwrd(result,3)
     year= substr(res,1,12)
     'q time'
     res = subwrd(result,3)
     hour = substr(res,1,2)
     say hour
     'q time'
     res = subwrd(result,3)
     day=substr(res,4,2)
     say day
     'q time'
     res = subwrd(result,3)
     month=substr(res,6,3)
     say month
'set t 'tt
'set grads off'
'set grid off'
'set csmooth on'

*跳过-0.1到0.1（指0线）不画。
'set black -0.1 0.1'   

*每间隔2条等值线作标注。
'set clskip 2'
'set ccols  1'
'set cthick  4'

*散度
'd div(lev=300)-div(lev=850)'

res=subwrd(result,3)
'draw title div 300-850hPa'hour'Z.'day'.'month'.2006'
'printim d:\fnl\czdiv300-850hPa2006.'month'.'day'.'hour'.png white '
'c'
tt=tt+1
endwhile

来自群组: 南气院校友总会

kaiofchina

{:5_275:}