# SVD QUESTION

From: shinn wong <shinnshinnwong_at_nyahnyahspammersnyahnyah>
Date: Sun Jul 24 2011 - 10:32:29 MDT

Dear all,
I am currently performing the SVD analysis of sst(time,lat,lon) and zonal
wind(time,lat,lon). What can I do if I want to get the time series
associated with the mode? Is the following code correct? Thanks in advance
for any available help.
Shinn

P = onedtond( ndtooned(sst), (/30 ,dimsizes(sst&lat)*dimsizes(sst&lon)/))
S = onedtond( ndtooned(wd) , (/30 ,dimsizes(wd&lat)*dimsizes(wd&lon)/))
P!0 = "time"
P!1 = "space"
S!0 = "time"
S!1 = "space"
PP = P(space|:,time|:)
SS = S(space|:,time|:)

svLeft = new((/nsvd,dimsizes(sst&lat)*dimsizes(sst&lon)/) ,float)
svRight = new((/nsvd,dimsizes(wd&lat)*dimsizes(wd&lon)/),float)

svdec = svdstd_sv(PP,SS,2,svLeft,svRight)
Le = svLeft
Le!0 = "time"
Le!1 = "space"
Left = Le(space|:,time|:)

transsvLeft = transpose(Left)
xL =
onedtond(ndtooned(transsvLeft),(/2,dimsizes(sst&lat),dimsizes(sst&lon)/))
copy_VarCoords(sst(0,:,:),xL(0,:,:))

sst_ts = eofcov_ts(sst(lat|:,lon|:,year|:),xL)

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sun Jul 24 10:32:36 2011

This archive was generated by hypermail 2.1.8 : Tue Aug 02 2011 - 14:31:37 MDT