Re: Removing the first pattern

[1] I think it should be:
     ENSOSST=conform(sst,eof(0,:,:),(/1,2/))*conform(sst,eof_ts(0,:),0)

[2] It is more efficient to use
     backgsst = sst - ENSOSST

     rather than:
     backgsst(:,:,:)=sst(:,:,:)-ENSOSST(:,:,:)

Good luck

On 9/19/12 8:52 PM, Li Yang wrote:
> Hi,
>
> I am trying to removing the first EOF pattern in SST field. My procedure is
>
> neof = 2 ; number of EOFs
> optEOF = True
> optEOF@jopt = 0 ; This is the default; most commonly used; no
> need to specify.
> optETS = False
>
> eof = eofunc_Wrap(x, neof, optEOF)
> eof_ts = eofunc_ts_Wrap (x, eof, optETS)
> ENSOSST=conform(sst,eof(0,:,:),(/1,2/))*conform(sst,eof_ts(0,:),1)
> backgsst(:,:,:)=sst(:,:,:)-ENSOSST(:,:,:)
> copy_VarCoords(sst,backgsst)
> printVarSummary(backgsst)
>
> Computed eof and eof_ts is right. And the fatal is :
>
> 内嵌图片 1
>
> The line 126 is
> "ENSOSST=conform(sst,eof(0,:,:),(/1,2/))*conform(sst,eof_ts(0,:),1)" in
> the above.
>
> How to make it correct?
> Thanks for any suggestions or info anyone might have!
>
> Yang
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Sep 20 07:27:28 2012