'reinit'
'open ../ctl/r30_cdas.ctl'
'open ../ctl/mask_10lon_5lat.staggered.ctl'
'enable print div_wind'
*
'set grads off'
'set map 1 1 6'
'set looping off' 
'set lat -50 50'
'set lon -20 360'
'define rad = 3.14159/180.'
'define dx = cdiff(lon,x)*rad*6.37e6*cos(lat*rad)'
'define dy = cdiff(lat,y)*rad*6.37e6'
*
'define uchi = cdiff(vpot,x)/dx'
'define vchi = cdiff(vpot,y)/dy'
*
'set vpage 0.0 11.0 0.0 8.5'
'set parea 1.2 10.7 4.8 8.3'
*
* ----- plot velocity potential -----
*
'set gxout contour'
'set grads off'
'set ccolor 1'
'set cthick 2'
'set cint 1'
'set clab off'
'd 1.e-6*vpot'
'set ccolor 1'
'set cthick 2'
'set clevs -10 -8 -6 -4 -2 0 2 4 6 8 10'
'set clab on'
'd 1.e-6*vpot'
'set ccolor 1'
'set cthick 5'
'set clab off'
'set clevs 0'
'd 1.e-6*vpot'
*
* ----- plot divergent wind -----
*
'set gxout vector'
'set cthick 2'
'set arrscl 0.3 5'
'set arrowhead -0.3'
'set arrlab off'
'set ccolor 1'
'define u = maskout(uchi,stag.2(t=1,lev=1))'
'define v = maskout(vchi,stag.2(t=1,lev=1))'
'd u;v'
*'d skip(u,2);v)'
*
*'draw title 200 mb Chi & Div. Wind \'
*'run ../gs/date.gs'
'run ../gs/equat.gs'
*
******** ANOMALIES ********
*
'define uchi = cdiff(vpota,x)/dx'
'define vchi = cdiff(vpota,y)/dy'
*
'set vpage 0.0 11.0 0.0 8.5'
'set parea 1.2 10.7 1.3 4.8'
*
* ----- plot velocity potential anomaly -----
*
'set gxout contour'
'set grads off'
'set ccolor 1'
'set cthick 2'
'set cint 1'
'set clab off'
'd 1.e-6*vpota'
'set ccolor 1'
'set cthick 2'
'set clevs -10 -8 -6 -4 -2 0 2 4 6 8 10'
'set clab on'
'd 1.e-6*vpota'
'set ccolor 1'
'set cthick 5'
'set clab off'
'set clevs 0'
'd 1.e-6*vpota'
*
* ----- plot divergent wind anomaly-----
*
'set gxout vector'
'set cthick 2'
'set arrlab on'
'set arrscl 0.3 5'
'set arrowhead -0.3'
'set ccolor 1'
'define u = maskout(uchi,stag.2(t=1,lev=1))'
'define v = maskout(vchi,stag.2(t=1,lev=1))'
'd u;v'
*'d skip(u,2);v'
*
*'draw title 200 mb Chi & Div. Wind Anomaly\'
*'run ../gs/date.gs'
'run ../gs/equat.gs'
'print'
'quit'