<div dir="ltr"><div><div>Dear All GrADs users, <br><br></div>Please help me to solve this problem.<br> 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. <br><br></div>Thank you in advance<br><div><br>'reinit'<br>'open SSt12.precip.lon-lat.allssws.ctl'<br>'open SSt12.lat-lon.precip.allssws.climatology.ctl'<br><br>'set lat -90 90'<br>'set lon <span style="background-color:rgb(0,255,0)">180</span>'<br>'set t 1 21'<br>'set xyrev on'<br>'set gxout contour'<br>'set csmooth on'<br>'set grads off'<br>'set grid off'<br>'x1=1000*24*60*60*(aprl.1+aprc.1)'<br>'x2=1000*24*60*60*(aprl.2+aprc.2)'<br>'xx1=tloop(ave(x1,lon=120,lon=180))'<br>'xx2=tloop(ave(x2,lon=120,lon=180))'<br>'set clab on'<br>'set xlabs 10|9| 8| 7| 6| 5| 4| 3| 2| 1| 0| 1| 2| 3| 4| 5| 6| 7| 8| 9| 10'<br>'d xx1-xx2'<br><br><br>'draw title zonal precipitation (120-180E) for all warmings (SSt12)'<br><br>'set string 1 tr 6'<br>'set strsiz 0.1'<br>'draw string 10 0.7 [mm/day]'<br><br>'set string 1 c 6'<br>'set strsiz 0.1'<br>'draw string 5 0.7 days before and after SSW event'<br><br>'enable print figure1.zonalpre120-180E.mf'<br>'print'<br>'disable print'<br>'!./gxeps -c -i figure1.zonalpre120-180E.mf -o SSt12.preci.120-180E.all.ssws.eps'<br><br><br clear="all"><div><div><div><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr">Dr. Sourabh Bal</div><div dir="ltr">Post-Doc Fellow </div><div dir="ltr">Institute of Meteorology, Freie Universitat, Berlin<div><br></div><div><br><br></div></div></div></div></div></div>
</div></div></div></div></div>