Re: How to get pcvar_varimax when using eofunc_varimax?

Hello

I ran the attached script using 4.3.0. It works fine.
If desired, I will send you ps or png.

Good luck
D
On Tue, 4 Sep 2007, Yi-chao Wu (???) wrote:

> Hi,
>
> I used eofunc_varimax to do a rotated EOF analysis but the variance returned were alwasy the same. I noticed that in your reply to other people who asked the same question in April, you indicated that the new version will solve this problem. So I downloaded the latest version (ncl-4.3.1.Linux_i686_gcc3) and replaced the old NCL library with the new one (all files in the directory lib/ncarg and "libnclapi.a" are updated). However, neither eofunc_varimax nor eof_varimax returnes the correct pcvar_varimax. Is this because I updated the library only but not the entire NCL package? Or is there somthing
> else missing?
>
> Below is my code and the responses. Thanks.
>
> Yi-chao Wu
>
> ------------------------------------------------------------------------------------------
> 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"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
>
> begin
>
>
> f=addfile("strm85.nc","r")
> str85=f->str85
>
> printVarSummary(str85)
>
>
> st85=str85(latitude|:,longitude|:,time|:) ;indian
>
>
> time=ispan(1958,2005,1)
> lat1=st85&latitude
> lon1=st85&longitude
>
> neval =3
> ev_cor = eofunc(st85,neval,1)
> ev_ts = eofcor_ts(st85,ev_cor)
> printVarSummary(ev_cor)
> print(ev_cor_at_pcvar)
> ev_rot = eofunc_varimax(ev_cor,1)
> ev_rot_ts = eofcor_ts(st85,ev_rot)
> ; eof_varimax_reorder ( ev_rot ) ; reorder the information
> ; printVarSummary(ev_rot)
>
>
> ev_cor!1="lat"
> ev_cor!2="lon"
>
> ev_cor&lat=lat1
> ev_cor&lon=lon1
>
> ev_rot!0="Rotated EOF"
> ev_rot!1="lat"
> ev_rot!2="lon"
>
> ev_rot&lat=lat1
> ev_rot&lon=lon1
>
>
> printVarSummary(ev_rot)
> print(ev_rot_at_pcvar_varimax(0))
> print(ev_rot_at_pcvar_varimax(1))
> print(ev_rot_at_pcvar_varimax(2))
>
> ------------------------------------------------------------------------------------------------
>
> Copyright (C) 1995-2006 - All Rights Reserved
> University Corporation for Atmospheric Research
> NCAR Command Language Version 4.2.0.a034
> The use of this software is governed by a License Agreement.
> See http://www.ncl.ucar.edu/ for more details.
>
>
> Variable: str85
> Type: float
> Total Size: 205632 bytes
> 51408 values
> Number of Dimensions: 3
> Dimensions and sizes: [time | 48] x [latitude | 21] x [longitude | 51]
> Coordinates:
> time: [ 0.. 47]
> latitude: [-12.5..37.5]
> longitude: [ 95.. 220]
> Number Of Attributes: 9
> comments : Unknown1 variable comment
> long_name : stream function (850 mb)
> units :
> grid_name : grid-1
> grid_type : linear
> level_description : Earth surface
> time_statistic : instantaneous
> missing_value : -9.99e+33
> _FillValue : -9.99e+33
>
>
> Variable: ev_cor
> Type: float
> Total Size: 12852 bytes
> 3213 values
> Number of Dimensions: 3
> Dimensions and sizes: [3] x [21] x [51]
> Coordinates:
> Number Of Attributes: 5
> _FillValue : -9.99e+33
> method : transpose
> matrix : covariance
> pcvar : <aRRAY>
> eval : <aRRAY>
> (0) 61.65483
> (1) 11.9065
> (2) 8.448109
>
>
> Variable: ev_rot
> Type: float
> Total Size: 12852 bytes
> 3213 values
> Number of Dimensions: 3
> Dimensions and sizes: [Rotated EOF | 3] x [lat | 21] x [lon | 51]
> Coordinates:
> lat: [-12.5..37.5]
> lon: [ 95.. 220]
> Number Of Attributes: 1
> pcvar_varimax : <aRRAY>
> (0) 0.09337068
> (0) 0.09337068
> (0) 0.09337068
>
>

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk