[ncl-talk] Coherence analysis

Dennis Shea shea at ucar.edu
Thu Mar 21 07:28:21 MDT 2019


*http://www.ncl.ucar.edu/Applications/spec.shtml*
<http://www.ncl.ucar.edu/Applications/spec.shtml>
*Example: CrosSpc_TimeLon_1.ncl*

===
The datasets used were some you have.
The example illustrates the 'mechanics' of doing the freq vs longitude.


On Thu, Mar 21, 2019 at 3:56 AM Soma Roy <somaroy892 at gmail.com> wrote:

> Please kindly inform me..
>
> Soma
>
> On Wed, Mar 20, 2019, 22:45 Soma Roy <somaroy892 at gmail.com wrote:
>
>> I am trying with the attached script as you said. But getting constant
>> values..
>>
>> Summary of the input data along with intermediate steps are as follows;
>>
>> Variable: dslp
>> Type: double
>> Total Size: 93312000 bytes
>>             11664000 values
>> Number of Dimensions: 3
>> Dimensions and sizes:   [time | 45] x [lat | 360] x [lon | 720]
>> Coordinates:
>>             time: [25748..41819]
>>             lat: [-89.75..89.75]
>>             lon: [-179.75..179.75]
>> Number Of Attributes: 5
>>   long_name :   diurnal temperature range
>>   units :       degrees Celsius
>>   _FillValue :  9.969209968386869e+36
>>   missing_value :       9.969209968386869e+36
>>   correlation_decay_distance :   750
>>
>> Variable: dslp2
>> Type: double
>> Total Size: 259200 bytes
>>             32400 values
>> Number of Dimensions: 2
>> Dimensions and sizes:   [time | 45] x [lon | 720]
>> Coordinates:
>>             time: [25748..41819]
>>             lon: [-179.75..179.75]
>> Number Of Attributes: 6
>>   _FillValue :  9.969209968386869e+36
>>   long_name :   diurnal temperature range
>>   units :       degrees Celsius
>>   correlation_decay_distance :   750
>>   missing_value :       9.969209968386869e+36
>>   average_op_ncl :      dim_avg_n over dimension(s): lat
>>
>> Variable: tslp2
>> Type: double
>> Total Size: 259200 bytes
>>             32400 values
>> Number of Dimensions: 2
>> Dimensions and sizes:   [time | 45] x [lon | 720]
>> Coordinates:
>>             time: [25748..41819]
>>             lon: [-179.75..179.75]
>> Number Of Attributes: 6
>>   _FillValue :  9.969209968386869e+36
>>   long_name :   precipitation
>>   units :       mm/month
>>   correlation_decay_distance :   450
>>   missing_value :       9.969209968386869e+36
>>   average_op_ncl :      dim_avg_n over dimension(s): lat
>>
>> Variable: pr
>> Type: double
>> Total Size: 5760 bytes
>>             720 values
>> Number of Dimensions: 1
>> Dimensions and sizes:   [720]
>> Coordinates:
>> Number Of Attributes: 20
>>   yslope :      1e+36
>>   ylag1 :       1e+36
>>   yvaro :       1e+36
>>   yvari :       1e+36
>>   yavei :       1e+36
>>   xslope :      1e+36
>>   xlag1 :       1e+36
>>   xvaro :       1e+36
>>   xvari :       1e+36
>>   xavei :       1e+36
>>   coher_probability :   (    0,    0,    0,    0 )
>>   bw :  1e+36
>>   phase :       <ARRAY of 22 elements>
>>   coher :       <ARRAY of 22 elements>
>>   quspc :       <ARRAY of 22 elements>
>>   cospc :       <ARRAY of 22 elements>
>>   frq : <ARRAY of 22 elements>
>>   spcy :        <ARRAY of 22 elements>
>>   spcx :        <ARRAY of 22 elements>
>>   _FillValue :  9.969209968386869e+36
>>
>> Variable: h1
>> Type: double
>> Total Size: 176 bytes
>>             22 values
>> Number of Dimensions: 1
>> Dimensions and sizes:   [22]
>> Coordinates:
>>
>> Variable: h2
>> Type: double
>> Total Size: 176 bytes
>>             22 values
>> Number of Dimensions: 1
>> Dimensions and sizes:   [22]
>> Coordinates:
>>
>>
>> All values are constant like;
>>  yslope :      1e+36
>>   ylag1 :       1e+36
>>   yvaro :       1e+36
>>   yvari :       1e+36
>> and so on..
>>
>> Please let me know where I am doing wrong/missing something.
>>
>> Thanking you,
>> Soma
>>
>> On Wed, Mar 20, 2019 at 8:07 PM Soma Roy <somaroy892 at gmail.com> wrote:
>>
>>> Can you kindly give me a hint how to use the function over a variable
>>> (say 3D variable as time, lat, lon) as you said?
>>>
>>> Thanking you,
>>> Soma
>>>
>>> On Wed, Mar 20, 2019, 20:00 Dennis Shea <shea at ucar.edu wrote:
>>>
>>>> For sure, these were created with NCL. However, I am not sure how these
>>>> were created.
>>>>
>>>> NCL's *specxy_anal*
>>>> <http://www.ncl.ucar.edu/Document/Functions/Built-in/specxy_anal.shtml>
>>>>
>>>> could have been applied *at each longitude* (loop over longitudes) and
>>>> the results stored in arrays:
>>>>
>>>>    COHER(frq,lon), PHASE(frq,lon), COSPC(frq,lon),  QUSPC(frq,lon)
>>>>
>>>> ===
>>>> Other, possibly useful, links:
>>>>
>>>> *http://www.ncl.ucar.edu/Applications/space_time.shtml*
>>>> <http://www.ncl.ucar.edu/Applications/space_time.shtml>
>>>>
>>>> *http://www.ncl.ucar.edu/Applications/mjoclivar.shtml*
>>>> <http://www.ncl.ucar.edu/Applications/mjoclivar.shtml>
>>>> *Example 11*
>>>> ====
>>>>
>>>> Good luck
>>>>
>>>> On Wed, Mar 20, 2019 at 2:42 AM Soma Roy <somaroy892 at gmail.com> wrote:
>>>>
>>>>> Hello,
>>>>>
>>>>> I am trying to calculate Coherence, Cospectrum, Quadrature and phase
>>>>> between two variables (two signals).
>>>>>
>>>>> Is there any function available in ncl to plot like below attached
>>>>> figure for the above?
>>>>>
>>>>> Please kindly inform me.
>>>>>
>>>>> Thanking you,
>>>>> Soma
>>>>> _______________________________________________
>>>>> ncl-talk mailing list
>>>>> ncl-talk at ucar.edu
>>>>> List instructions, subscriber options, unsubscribe:
>>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>>>
>>>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190321/c164011f/attachment.html>


More information about the ncl-talk mailing list