European Space Agency

Royal Belgian Institute for Aeronomy

uc170 - temporary magnetic field line

type  uc170

Synopsis

INCLUDE 'unilib.h'

COMMON /UC170/ nsg, kgp, mlab, mlin, mele

INTEGER*4   nsg, kgp
TYPE(zlbl)  mlab
TYPE(zfln)  mlin
TYPE(zseg)  mele(nx170)

Arguments

nsg

Number of elementary segment already computed

kgp

Indice of the current elementary segment

mlab

Label of the field line segment (zlbl)

mlin

Global description of the field line segment (zfln)

mele

Set of the elementary field line segments (zseg(nx170) )

Description

The common block uc170 provides a way to store a magnetic field line segment. The magnetic field line segment is defined as a set of elementary segments. The elementary segments are stored in the argument mele of the common block. The arguments mlab and mlin are used to store general information on the magnetic field line segments. The description of an elementary segment includes:

• a start point defined by its geographic position, mele()%beg%coord;

• the magnetic field vector at this point, mele()%beg%b;

• the local radius of curvature of the magnetic field line, mele()%beg%rcurv;

• the arc length along the magnetic field line from a fixed point to the start point mele()%arcl;

• the cosine of the local pitch-angle of any particle mirroring at the magnetic field intensity mlab%fbmp, mele()%csalp;

• the intermediate values of the magnetic field intensity obtained during the Runge-Kutta step, mele()%rkstp.

The pointers of the argument mlin correspond to the first and last segments of the field line segment (mlin%ind%jbeg and mlin%ind%jend) and to both mirror points (mlin%ind%jmirpn and mlin%ind%jmirps).

Note that a similar structure is used with the common blocks uc110, uc120, and uc130.

Dependencies

Used by

• ul220 - get information on a magnetic field line segment

  kgp, mlin, mlab, mele()%arcl, mele()%beg%b%dnrm and mele()%beg%coord

• ul230 - evaluate the integral invariant coordinate I

  mlin%ind%jmirpn, mlin%ind%jmirps, mele()%beg%b%dnrm, mele()%rkstp, mele()%csalp, mele()%arcl, mlab%lbmp and mlab%fbmp

• ud317 - trace a magnetic drift shell

  mele()%beg, mlin%ind

• ud319 - transfer a field line segment from common block uc170 to uc130

  nsg, mlab%lbmp, mlab%linv, mlab%fbmp, mlab%finv, mlin and mele

• ud328 - evaluate the second derivative of B

  kgp, nsg and mele()%arcl

• uf411 - search a local magnetic equator

  mlin%equat and mlab

• uf415 - rebuild the labels of the field line

  mele()%beg%b%dnrm

• uf421 - initialize and close a field line segment

  kgp

• uf422 - follow a field line until a boundary condition is reached

  kgp, mele()%beg%coord, mele()%beg%rcurv and mele()%beg%b%dnrm

• uf423 - Runge Kutta step

  kgp, mele()%beg%coord, mele()%beg%b and mele()%arcl

• uf424 - search the lowest B value

  nsg and mele()%beg

• uf426 - interpolating an extremum of B

  kgp and mele

• uf427 - interpolating a value of B

  kgp, mele()%beg%coord and mele()%beg%b

• uf428 - interpolating an altitude

  kgp, mele()%beg%coord and mele()%beg%b

• uf429 - transpose a field line segment

  mele

• ut985 - transfer a field line of the drift shell to IDL

  mlin%ind%jbeg, mlin%ind%jend, mlin%ind%jmirps, mlin%ind%jmirpn and mele

• ut991 - print the magnetic field line

  nsg, kgp, mlab, mlin and mele

• ut998 - evaluate the magnetic field vector and the normal to the field line

  nsg, kgp, mele()%beg%coord and mele()%beg%b

Modified by

• ul220 - get information on a magnetic field line segment

  mlin and mlab

• ul230 - evaluate the integral invariant coordinate I

  mlab%linv , mlab%lkauf , mlab%finv , mlab%fkauf, mlab%lbmp, mlab%fbmp and mele()%csalp

• ud328 - evaluate the second derivative of B

  mele()%beg%coord, mele()%arcl, mele()%beg%b and kgp

• uf415 - rebuild the labels of the field line

  mlab%linv, mlab%lbmp, mlab%lkauf, mlab%llmi, mlab%lalp0, mlab%lphi, mlab%ltim and mlab%fbmp

• uf420 - trace a magnetic field line segment passing through a given position

  kgp, mlab%linv, mlab%llmi, mlab%lkauf, mlab%fbmp, mlab%lbmp, mlin%ind, mlin%keast, mlin%kwest, mlin%equat%coord, mlin%footpn, mlin%footps, mlin%dtdft and mlin%drift

• uf421 - initialize and close a field line segment

  nsg, kgp, mele()%beg%coord, mele()%beg%b, mele()%beg%rcurv, mele()%arcl and mele()%rkstp()

• uf422 - follow a field line until a boundary condition is reached

  kgp and mele()%beg%rcurv

• uf423 - Runge Kutta step

  kgp, mele()%rkstp, mele()%beg%coord, mele()%beg%b, mele()%arcl, mele()%beg%rcurv and mele()%dtbnd

• uf424 - search the lowest B value

  kgp and mele()%beg%rcurv

• uf426 - interpolating an extremum of B

  kgp and mele

• uf427 - interpolating a value of B

  kgp , mele()%beg%coord and mele()%beg%b

• uf428 - interpolating an altitude

  kgp , mele()%beg%coord and mele()%beg%b

• uf429 - transpose a field line segment

  mele, mlin%ind%jmirps and mlin%ind%jmirpn

• ut998 - evaluate the magnetic field vector and the normal to the field line

  kgp, mele()%beg%coord, mele()%beg%b, mele()%beg%rcurv and mele()%arcl

See also

None

Examples

None

---

UNILIB/tags/v3.02