求助水汽通量散度

好好好我不说了

做了个画整层水汽通量散度的距平，目标是山东半岛，结果降水多的时候山东半岛散度是正的，少的时候是负的，也不知道哪里出了问题，发出来大家帮忙看看呗。思路是先求月气候态，再求一个偏多月份的，做差得到距平场。
'reinit'
'sdfopen F:\data\rhum.mon.mean.nc'
'sdfopen F:\data\uwnd.mon.mean.nc'
'sdfopen F:\data\vwnd.mon.mean.nc'
'sdfopen F:\data\pres.mon.mean.nc'
'sdfopen F:\data\air.mon.mean.nc'
'sdfopen F:\data\shum.mon.mean.nc'


'set map 1 1 3'
'set grads off'
'set grid off'
'set lon 74.99 145.01'
'set lat -0.01 60.01'
'set z 1 8'
'set xlint 10'
'set ylint 10'

**********************1月平均********************************
ts=828;
i=157;
while(i<=ts)
'set t 'i
'define u00=uwnd.2'
'define v00=vwnd.3'
'define e00=shum.6'
'define p00=lev'
'define q00=0.622*e00/(p00-0.378*e00)'
'define uq00=u00*q00/9.8'
'define vq00=v00*q00/9.8'
'set dfile 4'
'set z 1'
'ps00=pres.4'
'set dfile 2'
'set z 1'
'define uqm00=vint(ps00,uq00,300)'
'set dfile 3'
'set z 1'
'define vqm00=vint(ps00,vq00,300)'
'define qdiv00=hdivg(uqm00,vqm00)'
'uqm0'i'=uqm00'
'vqm0'i'=vqm00'
'qdiv0'i'=qdiv00'
i=i+12;
endwhile
'uqma1=0'
'vqma1=0'
'qdiva1=0'
i=157;
while(i<=ts)
'uqma1=uqma1+uqm0'i
'vqma1=vqma1+vqm0'i
'qdiva1=qdiva1+qdiv0'i
i=i+12;
endwhile
'uqma1=uqma1/56'
'vqma1=vqma1/56'
'qdiva1=qdiva1/56'

************一个偏多的1月的水汽通量和散度****************
'set t 535'
'define u1=uwnd.2'
'define v1=vwnd.3'
'define e1=shum.6'
'define p1=lev'
'define q1=0.622*e1/(p1-0.378*e1)'
'define uq1=u1*q1/9.8'
'define vq1=v1*q1/9.8'
'set dfile 4'
'set z 1'
'ps1=pres.4'
'set dfile 2'
'set z 1'
'define uqm1=vint(ps1,uq1,300)'
'set dfile 3'
'set z 1'
'define vqm1=vint(ps1,vq1,300)'
'define qdiv1=hdivg(uqm1,vqm1)'
************两者做差，得到本月距平************
'define uqmff=uqm1-uqma7'
'define vqmff=vqm1-vqma7'
'define qdivff=qdiv1-qdiva7'

好好好我不说了

除了这个问题，其他的部分画出图是没问题的，等改好了就发出来

高冷lxy

楼主改的怎么样了呢？？我也是在做散度和降水的对应，但是我还区分看FD厄尔尼诺年和SD厄尔尼诺年，可是6月、7月可以稍微对应，8月和夏季就对应的很不好了。