'reinit'
'open pptn.ctl'
'enable print pptn_climo_Hovmuller.gx'
*
*  define climatology (assuming t=1 is Jan, start at t=13 since earlier
*  values are undefined and average to the end of the dataset)...
*
'set z 1'
'set t 1 12'
'define pptnclim = ave(pptn,t+12,t=120,12)'
*
*  do not plot year for date (climatology)...
*
'set lon 270 300'
'set y 1'
'set tlsupp year'
'set gxout shaded'
'd ave(pptnclim,lat=15,lat=17)'
'draw title \\ "pptn" averaged for latitudes 15N to 17N'
'run cbar 0.7'
'print'
'disable print'
