[ncl-talk] PDO EOF_REGRES
Adam Phillips
asphilli at ucar.edu
Fri Mar 29 15:38:22 MDT 2019
Hi Melanie,
Typically what is done for PDO regressions is that you calculate the 1st
EOF, then the associated principal component (PC) timeseries. You then
regress other fields onto the PC timeseries. That way it is irrelevant if
the original SST grid is different from the grid of the variable you are
regressing.
For example:
evecv = eofunc(sst_CW({lat|20:70},{lon|110:260},time|:),2,75)
pcts = eofunc_ts(sst_CW({lat|20:70},{lon|110:260},time|:),evecv,False)
pctsS = dim_standardize(pcts(0,:),0)
delete([/pcts/])
finarr = regCoef(pctsS,sst2) ; returns global SST regressions
finarr_psl = regCoef(pctsS,psl) ; returns global PSL regressions
Adam
On Thu, Mar 28, 2019 at 8:32 AM Melanie O' hanoly <mel.hanoly at gmail.com>
wrote:
> Hi Soma
>
> I can not do this, because I just want to see the influence of the PDO
> (North Pacific) in the rest of the globe or I didnt understand what you mean
>
> Thanks
>
> On Thu, Mar 28, 2019 at 11:23 AM Soma Roy <somaroy892 at gmail.com> wrote:
>
>> Use printVarSummary and print for lat and lon for both the variables and
>> cut out only the common portion from the two..
>>
>> Hope this idea help you..
>>
>> Cheers,
>> Soma
>>
>> On Thu, Mar 28, 2019, 19:48 Melanie O' hanoly <mel.hanoly at gmail.com
>> wrote:
>>
>>> Hi NCL Users
>>>
>>> I'm trying to calculate the relationship between PDO and some variables
>>> of the globe, but I calculated the first EOF of north Pacific sst and then
>>> I'm trying regress analysis between salinity of global word but I got the
>>> error, because they dont have the same size, can anyone help me with this?
>>> The routine is attached
>>>
>>> Many Thanks in Advance
>>>
>>> Variable: x
>>> Type: float
>>> Total Size: 3225600 bytes
>>> 806400 values
>>> Number of Dimensions: 3
>>> Dimensions and sizes: [time | 576] x [latitude | 28] x [longitude | 50]
>>> Coordinates:
>>> time: [ 360.. 935]
>>> latitude: [20.70000076293945..69.30000305175781]
>>> longitude: [113.4000015258789..289.7999877929688]
>>> Number Of Attributes: 7
>>> depth : 17.5
>>> standard_name : sea_water_potential_temperature
>>> long_name : ocean potential temperature
>>> units : C
>>> _FillValue : 9.96921e+36
>>> missing_value : 9.96921e+36
>>> FillValue : 9.96921e+36
>>> (0) ocean potential temperature (C) : min=-1.7993 max=29.3727
>>>
>>> Variable: x2
>>> Type: float
>>> Total Size: 23040000 bytes
>>> 5760000 values
>>> Number of Dimensions: 3
>>> Dimensions and sizes: [time | 576] x [latitude | 100] x [longitude |
>>> 100]
>>> Coordinates:
>>> time: [ 360.. 935]
>>> latitude: [-89.09999847412109..89.09999847412109]
>>> longitude: [1.799999952316284..358.2000122070312]
>>> Number Of Attributes: 7
>>> depth : 17.5
>>> standard_name : sea_water_salinity
>>> long_name : ocean salinity
>>> units : 1e-3
>>> _FillValue : 9.96921e+36
>>> missing_value : 9.96921e+36
>>> FillValue : 9.96921e+36
>>> (0) ocean salinity (1e-3) : min=28.3825 max=37.6995
>>>
>>> Variable: eof_ts
>>> Type: double
>>> Total Size: 4608 bytes
>>> 576 values
>>> Number of Dimensions: 2
>>> Dimensions and sizes: [evn | 1] x [time | 576]
>>> Coordinates:
>>> evn: [1..1]
>>> time: [ 360.. 935]
>>> Number Of Attributes: 4
>>> ts_mean : 2.095727725316148e-07
>>> matrix : covariance
>>> _FillValue : 9.969209968386869e+36
>>> long_name : EOF: Amplitude: ocean potential temperature
>>> fatal:regCoef_n: the dimension sizes indicated by dims_nx and dims_ny
>>> must match
>>>
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
--
Adam Phillips
Associate Scientist, Climate and Global Dynamics Laboratory, NCAR
www.cgd.ucar.edu/staff/asphilli/ 303-497-1726
<http://www.cgd.ucar.edu/staff/asphilli>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190329/fbcd6368/attachment.html>
More information about the ncl-talk
mailing list