Re: how to make a composite

Your method looks correct to me...

On Aug 26, 2009, at 6:47 AM, ncl-talk-request_at_ucar.edu wrote:

> Hellooo,
>
> I have taken monthly sst from 1983-2005 and calculated
> monthly
> sst anomaly...From the total sst anomaly(i.e from 1983-2005) I
> selected 4
> years(1983,1990,1999,2004) and found the average of JJAS...So I got 4
> variables of dimension (lat=89,lon=180) (SSTA_1983_JJAS,
> SSTA_1990_JJAS,
> SSTA_1999_JJAS, SSTA_2004_JJAS)..I want to find the composite of
> these 4
> variables..
> The given method is valid or not...
> sum= SSTA_1983_JJAS+SSTA_1990_JJAS+SSTA_1999_JJAS+SSTA_2004_JJAS
> average=sum/4
>
> .Any help would be appreciable...Here I am attaching the script..
>
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> begin
> ;**************************************************************
> in2=addfile("mn_sst_1983-05.nc","r")
> sst=in2->sst
> sstflt=short2flt(sst)
> ;printVarSummary(sst)
> mn_clim_sst=clmMonTLL(sstflt)
> ;printVarSummary(mn_clim_sst)
> ;printVarSummary(mn_clim_sst)
> Mon_sst_anom=calcMonAnomTLL(sstflt,mn_clim_sst)
> ;*****************************************************
> SSTA_1983=Mon_sst_anom(0:11,:,:)
> SSTA_1983_JJAS=dim_avg(SSTA_1983(lat|:,lon|:,time|5:8))
> ;**************************************************
> SSTA_1990=Mon_sst_anom(84:95,:,:)
> SSTA_1990_JJAS=dim_avg(SSTA_1990(lat|:,lon|:,time|5:8))
> ;**************************************************
> SSTA_1999=Mon_sst_anom(192:203,:,:)
> SSTA_1999_JJAS=dim_avg(SSTA_1999(lat|:,lon|:,time|5:8))
> ;********************************************************

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Aug 26 2009 - 08:50:04 MDT