[gradsusr] Computing velocity potential & divergent from wind component

[syairah atiqah atiq]

Hi Grads expert,
Can anyone check my script. Is it correct? I do manage run the script, the problems is when I change the highlight line as shown below I do still get the answer. My question is which one is correct?Below is my full script. Thank you in advance.
   'reinit'*  'open E:\IMPORTANT\Desktop\Latest\+-mh\+.8w.jja.ctl'   'open E:\IMPORTANT\Desktop\Latest\+-mh\-.8w.jja.ctl'** +MH years* =========*   'uav=ave(u,t=1,t=6)'*   'vav=ave(v,t=1,t=6)'** -MH years* ==========  'uav=ave(u,t=1,t=7)'  'vav=ave(v,t=1,t=7)'** stream function* ---------------   'define psi1=fish_psi(uav,vav)'   'define psi=psi1-aave(psi1,global)'** Rotational Wind* ---------------   'define one=1+0*lat'   'define upsi=mvadv(one,psi)'   'define vpsi=-muadv(one,psi)'** Velocity potential* ------------------  'define chi=fish_chi(uav,vav)'     ---->  define chi1=fish_chi(uav,vav)  'psi=psi-aave(psi,global)'          -----> change to chi=chi1-aave(chi1,global)** Divergent Wind* ---------------  'define one=1+0*lat'  'define uchi=-muadv(one,chi)'  'define vchi=-mvadv(one,chi)'** Plotting* --------  'set map 1 1 12'  'set grid off'  'set grads off'** India  'set lon 45 90'  'set lat -10 30'* Peninsular*   'set lon 20 120'*   'set lat -75 20'* big region*  'set lon 20 170'*  'set lat -75 70'** Shaded colour (rainbow color)* ==============  'set rgb  17  130    0  220'  'set rgb  18  110    0  220'   'set rgb  19   30   60  255'  'set rgb  20   0  160  255'  'set rgb  21    0  200  200'*      'set rgb  22   0  220    0'   'set rgb  23 160  230   50'   'set rgb  24 230  175   45'   'set rgb  25 240  130   40'    'set rgb  26 250   60   60'  *   'set gxout shaded'   'set clevs -0.45 -0.35 -0.25 -0.15 -0.05 0.05 0.15 0.25 0.35 0.45'   'set ccols 17 18 19 20 21 0 22 23 24 25 26'    'set csmooth on'   'd chi/1e6'*  'set gxout contour'  'set csmooth on'  'set clab on'  'set cthick 7'  'set ccolor 1'  'set cint 0.03'  'set clopts -1 6 0.12'  'd chi/1e6'*  'set arrscl 0.5 0.5'  'set ccolor 1'  'set arrlab on'  'd skip(uchi,3,3);vchi'*  'cbarn'  'draw map'*  'set strsiz 0.15 0.15'  'draw string 9.7 1.0 1e6m2/s'*  'printim C:\Users\atiqah\Desktop\-mh.dvp.png white' 


Regards,Atiqah
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20170119/2d96a7bd/attachment.html