[gradsusr] could not plot latitude vs time (at average over longitude 120 to 180)

[Sourabh Bal]

Dear All GrADs users,

Please help me to solve this problem.
 I have two files (attached), file1=precipitation of ssw yrs and
file2=climatology. I want to calculate the anomalies by taking the
difference. So I started with the following script. In order to use tloop
and plot latitude vs time I had to fix the longitude to 180. When I am
changing the value of the set latitude to say 120, 130, 140..I am getting
different anomalies plot. Finally, Please suggest how to calculate the t
test(standard deviation) over the same plot as shaded.

Thank you in advance

'reinit'
'open SSt12.precip.lon-lat.allssws.ctl'
'open SSt12.lat-lon.precip.allssws.climatology.ctl'

'set lat -90 90'
'set lon 180'
'set t 1 21'
'set xyrev on'
'set gxout contour'
'set csmooth on'
'set grads off'
'set grid off'
'x1=1000*24*60*60*(aprl.1+aprc.1)'
'x2=1000*24*60*60*(aprl.2+aprc.2)'
'xx1=tloop(ave(x1,lon=120,lon=180))'
'xx2=tloop(ave(x2,lon=120,lon=180))'
'set clab on'
'set xlabs 10|9| 8| 7| 6| 5| 4| 3| 2| 1| 0| 1| 2| 3| 4| 5| 6| 7| 8| 9| 10'
'd xx1-xx2'


'draw title zonal precipitation (120-180E) for all warmings (SSt12)'

'set string 1 tr 6'
'set strsiz 0.1'
'draw string 10 0.7 [mm/day]'

'set string 1 c 6'
'set strsiz 0.1'
'draw string 5 0.7 days before and after SSW event'

'enable print figure1.zonalpre120-180E.mf'
'print'
'disable print'
'!./gxeps -c -i figure1.zonalpre120-180E.mf -o
SSt12.preci.120-180E.all.ssws.eps'



-- 
Dr. Sourabh Bal
Post-Doc Fellow
Institute of Meteorology, Freie Universitat, Berlin
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20141120/99d5bb71/attachment-0001.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: SSt12.precip.lon-lat.allssws.grb
Type: application/octet-stream
Size: 390768 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20141120/99d5bb71/attachment-0002.obj 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: SSt12.lat-lon.precip.allssws.climatology.grb
Type: application/octet-stream
Size: 390768 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20141120/99d5bb71/attachment-0003.obj