Re: How to get pcvar_varimax when using eofunc_varimax?

Also, I'd like to point out that it looks like you are using version
4.2.0.a034 and not 4.3.0/1. When you run "ncl", it echoes the version
number, and your case, it says:

        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.

If you have two versions of NCL installed on your system, you need to
make sure that your search path is set up so that it points to the
newer version first. Also, you may need to adjust your NCARG_ROOT
environment variable. Let me know if you need help with this.

--Mary

On Thu, 6 Sep 2007, Dennis Shea wrote:

> 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
Received on Thu Sep 06 2007 - 15:52:59 MDT