Dear Ansley and Izidine,
Thanks a lot for both of your suggestions.
I could almost get the plot except the annual mean field. Hence I need
some help on that.
Please see below the commands I have used.
-----------------------------------------------------------
yes? use sst_mon_81_07.nc
yes? show data
currently SET data sets:
1> /users/cshaji/NOW/sst_mon_81_07.nc (default)
name title I J K L
SST ** monthly from december 1981 ( 1:360 1:180 1:1
1:311
!Get the climatological axis we want to use
!------------------------------------
yes? USE climatological_axes
*** NOTE: regarding /usr/local/ferret/go/climatological_axes.cdf ...
*** NOTE: Climatological axes SEASONAL_REG, MONTH_REG, MONTH_IRREG,
MONTH_GREGORIAN, MONTH_NOLEAP, MONTH_360_DAY defined
yes? CANCEL DATA climatological_axes
!Define monthly SST for the region & period of interest
!------------------------------------
yes? let sst_area_avg=sst[x=30:78@ave,y=5:31@ave]
yes? let sst_of_interest=sst_area_avg[T=1-jan-1982:31-dec-2006]
!Define monthly SST climatology for the region & period of interest
!------------------------------------
yes? let sst_mon_clm = sst_area_avg[GT=month_reg@MOD]
yes? let sst_of_interest_clm=sst_mon_clm[T=1-jan-1982:31-dec-2006]
!Time series plot of sst monthly & monthly climatology
!------------------------------------
yes? plot sst_of_interest
yes? plot/overlay sst_of_interest_clm
-----------------------------------------------------------
I also attach the figure I am getting along with this email.
In the same figure I would like to have the annual mean SST of each year
also plotted (a curve connecting all 25 SST values).
I could not understand the steps how to do this. Can you please give me
some suggestions in this?
Also how can I get the 3-months running mean SST also plotted in the same
time seris plot?
Thanks in advance,
Shaji.
On Tue, May 25, 2010 9:50 pm, Ansley Manke said:
> Hi -
> A good way to find things in the Ferret Users Guide is to load the Index
> into your browser,
> ../../documentation/users-guide/index
> and search for the words you want in there.
>
> To do an average in X and Y, look at the index entry "average,
> transformation @AVE", where there are examples of the kind of thing you
> want to do.
>
> For making an annual mean from data on a monthly time axis, you will
> define a time axis with units of a year, and use a regridding
> transformation to put your area-averaged variable onto that new time
> axis. You could use viewports to show both the monthly and yearly time
> series together (set view upper; plot monthly_variable; set view lower;
> plot yearly_variable). I'll let you work out the details and if you
> have questions once you've made a start, please feel free to write back
> to the list.
>
> Ansley
>
> On 5/24/2010 1:50 PM, C. Shaji wrote:
>> Dear all,
>> I am interested to get a time series plot of monthly and annual mean SST
>> using multi-year SST field. There is only one data file conatining all
>> time levels data.
>>
>> Please see the data info below.
>>
>> yes? show grid
>> Default grid for DEFINE VARIABLE is ABSTRACT
>> Last successful data access was on grid GEU1
>> GRID GEU1
>> name axis # pts start end
>> LON LONGITUDE 360mr 0.5E 0.5W
>> LAT LATITUDE 180 r 89.5S 89.5N
>> LEV HEIGHT (millib 1 r- 1 1
>> TIME TIME 311 i 01-DEC-1981 00:00 01-OCT-2007
>> 00:00
>> yes? show data
>> currently SET data sets:
>> 1> /users/cshaji/NOW/sst_mon_81_07.nc (default)
>> name title I J K
>> L
>> SST ** monthly from december 1981 ( 1:360 1:180 1:1
>> 1:311
>>
>> First I want to get an area (x=30E:78E,y=5N:31N) average of SST. Then
>> plot
>> the time series of the area-averaged monthly as well as annual mean SST
>> for the time period from l=2 to l=301 only (covering January 1982 to
>> December 2006) in the same figure window.
>>
>> I could not find any demonstration on this.
>>
>> I would highly appreciate if somebody please give me some suggestions.
>>
>> Thanks in advance
>>
>> Shaji.
>>
>>
>
Attachment:
fig_sst.pdf
Description: Adobe PDF document