-
Notifications
You must be signed in to change notification settings - Fork 106
/
shsjkpg.1
96 lines (96 loc) · 3.38 KB
/
shsjkpg.1
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
.\" Automatically generated by Pandoc 2.5
.\"
.TH "shsjkpg" "1" "2017\-11\-28" "Fortran 95" "SHTOOLS 4.4"
.hy
.SH SHSjkPG
.PP
Calculate the expectation of the product of two functions, each
multiplied by a different data taper, for a given spherical harmonic
degree and two different angular orders.
.SH Usage
.PP
\f[C]value\f[R] = SHSjkPG (\f[C]incspectra\f[R], \f[C]l\f[R],
\f[C]m\f[R], \f[C]mprime\f[R], \f[C]hj_real\f[R], \f[C]hk_real\f[R],
\f[C]mj\f[R], \f[C]mk\f[R], \f[C]lwin\f[R], \f[C]hkcc\f[R])
.SH Parameters
.TP
.B \f[C]value\f[R] : output, complex*16
The expectation of the product of two functions, each multiplied by a
different data taper, for a given spherical harmonic degree and two
different angular orders.
.TP
.B \f[C]incspectra\f[R] : input, real*8, dimension (\f[C]l\f[R]+\f[C]lwin\f[R]+1)
The global cross\-power spectrum of \f[C]f\f[R] and \f[C]g\f[R].
.TP
.B \f[C]l\f[R] : input, integer
The spherical harmonic degree for which to calculate the expectation.
.TP
.B \f[C]m\f[R] : input, integer
The angular order of the first localized function, \f[C]Phi\f[R].
.TP
.B \f[C]mprime\f[R] : input, integer
The angular order of the second localized function, \f[C]Gamma\f[R].
.TP
.B \f[C]hj_real\f[R] : input, real*8, dimension (\f[C]lwin\f[R]+1)
The real spherical harmonic coefficients of angular order \f[C]mj\f[R]
used to localize the first function \f[C]f\f[R].
These are obtained by a call to \f[C]SHReturnTapers\f[R].
.TP
.B \f[C]hk_real\f[R] : input, real*8, dimension (\f[C]lwin\f[R]+1)
The real spherical harmonic coefficients of angular order \f[C]mk\f[R]
used to localize the second function \f[C]g\f[R].
These are obtained by a call to \f[C]SHReturnTapers\f[R].
.TP
.B \f[C]mj\f[R] : input, integer
The angular order of the window coefficients \f[C]hj_real\f[R].
.TP
.B \f[C]mk\f[R] : input, integer
The angular order of the window coefficients \f[C]hk_real\f[R].
.TP
.B \f[C]lwin\f[R] : input, integer
the spherical harmonic bandwidth of the localizing windows
\f[C]hj_real\f[R] and \f[C]hk_real\f[R].
.TP
.B \f[C]hkcc\f[R] : input, integer
If 1, the function described in the \f[C]description\f[R] will be
calculated as is.
If 2, the second localized function \f[C]Gamma\f[R] will not have its
complex conjugate taken.
.SH Description
.PP
\f[C]SHSjkPG\f[R] will calculate the expectation of two functions
(\f[C]f\f[R] and \f[C]g\f[R]), each localized by a different data taper
that is a solution of the spherical cap concentration problem, for a
given spherical harmonic degree and two different angular orders.
As described in Wieczorek and Simons (2007), this is the function
.IP
.nf
\f[C]
/ m(j) mprime(k)* \[rs]
| Phi Gamma |
\[rs] l l /
\f[R]
.fi
.PP
The global cross\-power spectrum of \f[C]f\f[R] and \f[C]g\f[R] is input
as \f[C]incspectra\f[R], and the real coefficients of the two data
tapers of angular order \f[C]mj\f[R] and \f[C]mk\f[R] (obtained by a
call to \f[C]SHReturnTapers\f[R]) are specified by \f[C]hj_real\f[R] and
\f[C]hk_real\f[R].
If \f[C]hkcc\f[R] is set to 1, then the above function is calculated as
is.
However, if this is set to 2, then the complex conjugate of the second
localized function is not taken.
.SH References
.PP
Wieczorek, M.
A.
and F.
J.
Simons, Minimum\-variance multitaper spectral estimation on the sphere,
J.
Fourier Anal.
Appl., 13, doi:10.1007/s00041\-006\-6904\-1, 665\-692, 2007.
.SH See also
.PP
shreturntapers, shmtvaropt