-
Notifications
You must be signed in to change notification settings - Fork 103
/
makemaggridpoint.doc
47 lines (39 loc) · 1.63 KB
/
makemaggridpoint.doc
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
Determine the three components of the magnetic field vector at a single point.
Usage
-----
value = MakeMagGridPoint (cilm, a, r, lat, lon, [lmax, dealloc])
Returns
-------
value : float, dimension (3)
Vector components (r, theta, phi) of the magnetic field at (r, lat, lon).
Parameters
----------
cilm : float, dimension (2, lmaxin+1, lmaxin+1)
The real Schmidt semi-normalized spherical harmonic coefficients of the
magnetic potential. The coefficients C0lm and C1lm refer to the cosine (Clm)
and sine (Slm) coefficients, respectively, with Clm=cilm[0,1,m] and
Slm=cilm[1,l,m].
a: float
The reference radius of the spherical harmonic coefficients.
r: float
The radius to evaluate the magnetic field.
lat : float
The latitude of the point in degrees.
lon : float
The longitude of the point in degrees.
lmax : optional, integer, default = lmaxin
The maximum spherical harmonic degree used in evaluating the function.
dealloc : optional, integer, default = 0
0 (default) = Save variables used in the external Legendre function calls.
(1) Deallocate this memory at the end of the funcion call.
Description
-----------
MakeMagGridPoint will compute the three components of the magnetic field vector
at a single point. The input latitude and longitude are in degrees, and the
output components are in spherical coordinates (r, theta, phi). The magnetic
potential is given by
V = a Sum_{l=0}^lmax (a/r)^(l+1) Sum_{m=-l}^l C_{lm} Y_{lm},
and the vector magnetic field is
B = - Grad V.
The coefficients are referenced to a radius a, and the output values have the
same units as the input spherical harmonic coefficients.