[cig-commits] r4705 - geodyn/3D/MAG/trunk/doc
sue at geodynamics.org
sue at geodynamics.org
Wed Oct 4 17:15:37 PDT 2006
Author: sue
Date: 2006-10-04 17:15:37 -0700 (Wed, 04 Oct 2006)
New Revision: 4705
Modified:
geodyn/3D/MAG/trunk/doc/mag_book.lyx
Log:
fixed formatting
Modified: geodyn/3D/MAG/trunk/doc/mag_book.lyx
===================================================================
--- geodyn/3D/MAG/trunk/doc/mag_book.lyx 2006-10-05 00:15:03 UTC (rev 4704)
+++ geodyn/3D/MAG/trunk/doc/mag_book.lyx 2006-10-05 00:15:37 UTC (rev 4705)
@@ -1,4 +1,4 @@
-#LyX 1.4.2 created this file. For more info see http://www.lyx.org/
+#LyX 1.4.3 created this file. For more info see http://www.lyx.org/
\lyxformat 245
\begin_document
\begin_header
@@ -58,7 +58,7 @@
placement H
wide false
sideways false
-status open
+status collapsed
\begin_layout Standard
\begin_inset Graphics
@@ -102,6 +102,26 @@
Preface
\end_layout
+\begin_layout Standard
+\begin_inset ERT
+status collapsed
+
+\begin_layout Standard
+
+
+\backslash
+raggedbottom
+\end_layout
+
+\begin_layout Standard
+
+\end_layout
+
+\end_inset
+
+
+\end_layout
+
\begin_layout Section*
About This Document
\end_layout
@@ -710,9 +730,8 @@
After unpacking the source, you will find following directories:
\end_layout
-\begin_layout List
-\labelwidthstring 00.00.0000
-~/src: It contains the set of FORTRAN source code files with suffix ".f".
+\begin_layout Description
+~/src It contains the set of FORTRAN source code files with suffix ".f".
This includes sample grid parameter value files with names like "param32s4.f"
for a coarse grid (up to 32 spherical harmonics, 24 radial grid intervals,
and 4-fold symmetry in
@@ -740,16 +759,14 @@
.
\end_layout
-\begin_layout List
-\labelwidthstring 00.00.0000
-~/doc: This is the directory where you will find this manual and other documenta
-tion files.
+\begin_layout Description
+~/doc This is the directory where you will find this manual and other documentat
+ion files.
\end_layout
-\begin_layout List
-\labelwidthstring 00.00.0000
-~/bench-data: Output files "ls.benchX", "l.benchX", "g.benchx", and "d.benchx"
+\begin_layout Description
+~/bench-data Output files "ls.benchX", "l.benchX", "g.benchx", and "d.benchx"
obtained with short runs of benchmark0 and benchmark1 on a Linux workstation.
Explanations of the contents of these files are found in Appendix
\begin_inset LatexCommand \ref{cha:MAG-Output-File}
@@ -761,9 +778,8 @@
your local run of the MAG.
\end_layout
-\begin_layout List
-\labelwidthstring 00.00.0000
-~/idl: This is where the postprosessing IDL routines reside.
+\begin_layout Description
+~/idl This is where the postprosessing IDL routines reside.
\end_layout
\begin_layout Subsection
@@ -774,19 +790,26 @@
Prepare MAG for running
\end_layout
-\begin_layout Itemize
+\begin_layout Enumerate
First you need to create a path for execution of magx (an example; use your
path) :
\end_layout
+\begin_deeper
\begin_layout LyX-Code
printenv PATH $ PATH=$PATH:/your_mag_dir_path $ export PATH
\end_layout
-\begin_layout Itemize
-make, using existing param.f grid and symmetry
+\end_deeper
+\begin_layout Enumerate
+Compile the program using
+\family typewriter
+make
+\family default
+, which by default uses the existing param.f grid and symmetry
\end_layout
+\begin_deeper
\begin_layout LyX-Code
$ make
\end_layout
@@ -797,6 +820,7 @@
and Z=longitudinal symmetry
\end_layout
+\end_deeper
\begin_layout Section
\begin_inset LatexCommand \label{sec:Software-Repository}
@@ -967,7 +991,7 @@
\end_deeper
\begin_layout Enumerate
-mv magx magx32s64 [Renaming, optional]
+mv magx magx32s64 [Renaming is optional]
\end_layout
\begin_layout Enumerate
@@ -1026,9 +1050,10 @@
\end_layout
\begin_layout Standard
-Upon finishing running MAG, you should have a series of output data files.
- To visualize your results, MAG software package provid a set of IDL routines.
- IDL routines designed for MAG can be found in the directory called
+Once you finish running MAG, you should have a series of output data files.
+ To visualize your results, MAG software package provides a set of IDL routines.
+ IDL routines which are designed for MAG can be found in the directory called
+
\family typewriter
PREFIX/idl
\family default
@@ -1056,7 +1081,7 @@
placement H
wide false
sideways false
-status collapsed
+status open
\begin_layout Standard
\align center
@@ -1079,14 +1104,17 @@
\end_layout
\begin_layout Standard
+\noindent
\align center
\begin_inset Float figure
placement H
wide false
sideways false
-status collapsed
+status open
\begin_layout Standard
+\noindent
+\align center
\begin_inset Graphics
filename F6spct.ps
scale 50
@@ -1132,14 +1160,17 @@
\end_layout
\begin_layout Standard
+\noindent
\align center
\begin_inset Float figure
placement H
wide false
sideways false
-status collapsed
+status open
\begin_layout Standard
+\noindent
+\align center
\begin_inset Graphics
filename map.jpeg
lyxscale 30
@@ -1160,6 +1191,7 @@
\end_layout
\begin_layout Standard
+\noindent
\align center
\begin_inset Float table
placement H
@@ -1505,44 +1537,49 @@
\end_layout
\begin_layout Description
-adrke: : axisymmetric toroidal kinetic energy (diagnostic)
+adrke axisymmetric toroidal kinetic energy (diagnostic)
\end_layout
\begin_layout Description
-ai: : imaginary unit = complex(0,1)
+ai imaginary unit = complex(0,1)
\end_layout
\begin_layout Description
-aj:(nlma,nn+1) : poloidal magnetic field potential (spectral form) the second
+aj\InsetSpace ~
+(nlma,nn+1) poloidal magnetic field potential (spectral form) the second
index is either the Chebycheff degree (n) or the radial grid point (kc)
\end_layout
\begin_layout Description
-ajmat(nn,nn,lmax): LU-decomposed matrix from Chebycheff collocation of tor.
+ajmat\InsetSpace ~
+(nn,nn,lmax) LU-decomposed matrix from Chebycheff collocation of tor.
induction equation.
Built in ludc, used in amhd.
\end_layout
\begin_layout Description
-aleg1: (nlma,ni) : Value of associated Legendre function at grid points
+aleg1\InsetSpace ~
+(nlma,ni) Value of associated Legendre function at grid points
\end_layout
\begin_layout Description
-aleg2: (nlma,ni) : Value of associated Legendre function, multiplied with
- Gaussian weight, at grid points
+aleg2\InsetSpace ~
+(nlma,ni) Value of associated Legendre function, multiplied with Gaussian
+ weight, at grid points
\end_layout
\begin_layout Description
-aleg3: (nlma,ni) : Value of derivative of associated Legendre function multiplie
-d with sin(theta) at grid points
+aleg3\InsetSpace ~
+(nlma,ni) Value of derivative of associated Legendre function multiplied
+ with sin(theta) at grid points
\end_layout
\begin_layout Description
-alfilt: INPUT: Filter parameter for B_r in graphics output, see nfilt
+alfilt [INPUT] Filter parameter for B_r in graphics output, see nfilt
\end_layout
\begin_layout Description
-alffac: INPUT: controls the contribution of the (modified) Alfven velocity
+alffac [INPUT] Controls the contribution of the (modified) Alfven velocity
to the Courant time step limit (see under "courfac").
The modified Alfven velocity is given by v_alfven' = (v_a)^2 / {(v_a)^2
+[pi*(eta+nu)/delx]^2} where v_a = B/sqrt(mu*rho) and delx is the Courant
@@ -1550,256 +1587,284 @@
\end_layout
\begin_layout Description
-alpha: INPUT: =0 linear terms in the equations are treated fully explicit,
+alpha [INPUT] =0 linear terms in the equations are treated fully explicit,
=1 linear terms are treated fully implicit, =0.5: Crank-N.
\end_layout
\begin_layout Description
-alumn0: : factor for scaling heat flow in output amcke: : axisymmetric poloidal
+alumn0 Factor for scaling heat flow in output amcke: axisymmetric poloidal
kinetic energy
\end_layout
\begin_layout Description
-amhd: SUBROUTINE: the "workhorse" of the program, advance solution by nstep
+amhd [SUBROUTINE] The "workhorse" of the program: advance solution by nstep
time steps
\end_layout
\begin_layout Description
-amps: INPUT: can be used to re-scale entropy ampj: INPUT: can be used to
+amps [INPUT] can be used to re-scale entropy ampj: INPUT: can be used to
re-scale toroidal magn.
field
\end_layout
\begin_layout Description
-ampb: INPUT: can be used to re-scale poloidal magn.
+ampb [INPUT] can be used to re-scale poloidal magn.
field
\end_layout
\begin_layout Description
-ampw: INPUT: can be used to re-scale poloidal velocity
+ampw [INPUT] can be used to re-scale poloidal velocity
\end_layout
\begin_layout Description
-ampz: INPUT: can be used to re-scale toroidal velocity
+ampz [INPUT] can be used to re-scale toroidal velocity
\end_layout
\begin_layout Description
-anorm: : = sqrt(2/[nn+1])
+anorm = sqrt(2/[nn+1])
\end_layout
\begin_layout Description
-apome: : axisymmetric poloidal magnetic field energy
+apome axisymmetric poloidal magnetic field energy
\end_layout
\begin_layout Description
-atome: : axisymmetric toroidal magnetic field energy
+atome axisymmetric toroidal magnetic field energy
\end_layout
\begin_layout Description
-b: (nlma,nn+1) : poloidal mag.
+b\InsetSpace ~
+(nlma,nn+1) poloidal mag.
field potential (spectral form, see aj)
\end_layout
\begin_layout Description
-bleg1: (lmax) : auxiliary array for calculation of aleg1
+bleg1\InsetSpace ~
+(lmax) auxiliary array for calculation of aleg1
\end_layout
\begin_layout Description
-bleg2: (lmax) : auxiliary array for calculation of aleg2
+bleg2\InsetSpace ~
+(lmax) auxiliary array for calculation of aleg2
\end_layout
\begin_layout Description
-bleg3: (lmax) : auxiliary array for calculation of aleg3
+bleg3\InsetSpace ~
+(lmax) auxiliary array for calculation of aleg3
\end_layout
\begin_layout Description
-bmat(nn,nn,lmax):: LU-decomposed matrix from Chebycheff collocation of pol.
+bmat(nn,nn,lmax) LU-decomposed matrix from Chebycheff collocation of pol.
induction equation.
Built in ludc, used in amhd.
\end_layout
\begin_layout Description
-bnlc1: (nja/2,ni): bnlr1 stored in complex form
+bnlc1\InsetSpace ~
+(nja/2,ni) bnlr1 stored in complex form
\end_layout
\begin_layout Description
-bnlc2: (nja/2,ni): bnlr2 stored in complex form
+bnlc2\InsetSpace ~
+(nja/2,ni) bnlr2 stored in complex form
\end_layout
\begin_layout Description
-bnlc3: (nja/2,ni): bnlr3 stored in complex form
+bnlc3\InsetSpace ~
+(nja/2,ni) bnlr3 stored in complex form
\end_layout
\begin_layout Description
-bnlr1: (nja,ni) : nonlinear products for updating b (on grid points)
+bnlr1\InsetSpace ~
+(nja,ni) nonlinear products for updating b (on grid points)
\end_layout
\begin_layout Description
-bnlr2: (nja,ni) : nonlinear products for updating aj (on grid points)
+bnlr2\InsetSpace ~
+(nja,ni) nonlinear products for updating aj (on grid points)
\end_layout
\begin_layout Description
-bnlr3: (nja,ni) : nonlinear products for updating aj (on grid points)
+bnlr3\InsetSpace ~
+(nja,ni) nonlinear products for updating aj (on grid points)
\end_layout
\begin_layout Description
-bots(0:lmax,0:mmax): INPUT: harmonic coefficients of prescribed temperature
+bots(0:lmax,0:mmax) [INPUT: harmonic coefficients of prescribed temperature
(entropy) on inner boundary
\end_layout
\begin_layout Description
-br: (nja,ni) : = r^2 * B_r on gridpoints
+br\InsetSpace ~
+(nja,ni) = r^2 * B_r on gridpoints
\end_layout
\begin_layout Description
-brc: : br stored as complex array
+brc br stored as complex array
\end_layout
\begin_layout Description
-bscl: : = dt * radtop^2 bt: (nja,ni) : = r * sin(theta) *b_theta
+bscl = dt * radtop^2 bt: (nja,ni) : = r * sin(theta) *b_theta
\end_layout
\begin_layout Description
-btrdt: (ni) : used in movmout to calculate j_phi
+btrdt\InsetSpace ~
+(ni) used in movmout to calculate j_phi
\end_layout
\begin_layout Description
-bts: (ni,3) : used in movmout to calculate j_phi
+bts\InsetSpace ~
+(ni,3) used in movmout to calculate j_phi
\end_layout
\begin_layout Description
-btc: : bt stored as complex array
+btc bt stored as complex array
\end_layout
\begin_layout Description
-bp: (nja,ni) : = r * sin(theta) * b_phi
+bp\InsetSpace ~
+(nja,ni) = r * sin(theta) * b_phi
\end_layout
\begin_layout Description
-bpc: : bp stored as complex array
+bpc bp stored as complex array
\end_layout
\begin_layout Description
-bpeak: INPUT: maximum value of imposed field on boundaries
+bpeak [INPUT] maximum value of imposed field on boundaries
\end_layout
\begin_layout Description
-bpeakbot: : maximum value of imposed field on inner boundary
+bpeakbot maximum value of imposed field on inner boundary
\end_layout
\begin_layout Description
-bpeaktop: : maximum value of imposed field on outer boundary
+bpeaktop maximum value of imposed field on outer boundary
\end_layout
\begin_layout Description
-cbr: (nja,ni) : = r^2 * curl (B) * e_r
+cbr\InsetSpace ~
+(nja,ni) = r^2 * curl (B) * e_r
\end_layout
\begin_layout Description
-cbrc: : cbr stored as complex array
+cbrc cbr stored as complex array
\end_layout
\begin_layout Description
-cbt: (nja,ni) : = r * sin(theta) * curl (B) * e_theta
+cbt\InsetSpace ~
+(nja,ni) = r * sin(theta) * curl (B) * e_theta
\end_layout
\begin_layout Description
-cbtc: : ctr stored as complex array
+cbtc ctr stored as complex array
\end_layout
\begin_layout Description
-cbp: (nja,ni) : = r * sin(theta) * curl (B) * e_phi
+cbp\InsetSpace ~
+(nja,ni) = r * sin(theta) * curl (B) * e_phi
\end_layout
\begin_layout Description
-cbpc: : cpb stored as complex array
+cbpc cpb stored as complex array
\end_layout
\begin_layout Description
-cheb: (nn,nn) : cheb(i,j) = value of Ch.
+cheb\InsetSpace ~
+(nn,nn) cheb(i,j) = value of Ch.
polyn.
i at grid point j
\end_layout
\begin_layout Description
-chebi: SUBROUTINE: initialize subroutine chebtf
+chebi [SUBROUTINE] initialize subroutine chebtf
\end_layout
\begin_layout Description
-chebtf: SUBROUTINE: multiple fast Chebycheff transform
+chebtf [SUBROUTINE] multiple fast Chebycheff transform
\end_layout
\begin_layout Description
-clm: (lmax,mmax) : normalization factors of spherical harmonics
+clm\InsetSpace ~
+(lmax,mmax) normalization factors of spherical harmonics
\end_layout
\begin_layout Description
-cmb: INPUT: integrated conductivity of thin D"-layer
+cmb [INPUT] integrated conductivity of thin D"-layer
\end_layout
\begin_layout Description
-colat: (ni) : vector of colatitudes (Gauss points), local array in subroutine
+colat\InsetSpace ~
+(ni) vector of colatitudes (Gauss points), local array in subroutine
prep
\end_layout
\begin_layout Description
-courfac: INPUT: factor controlling the time step as fraction of courant
+courfac [INPUT] factor controlling the time step as fraction of courant
advection length.
The time step is limited to dt < min( dx/[ courfac * v + alffac * v_alfven'
] )
\end_layout
\begin_layout Description
-cvr: (nja,ni) : =r^2 * curl(v) * e_r
+cvr\InsetSpace ~
+(nja,ni) =r^2 * curl(v) * e_r
\end_layout
\begin_layout Description
-cvrc: : cvr stored as complex array
+cvrc cvr stored as complex array
\end_layout
\begin_layout Description
-db: (nlma,nn+1) : radial deriv.
+db\InsetSpace ~
+(nlma,nn+1) radial deriv.
of pol.
mag.
potential (spectral form, see aj)
\end_layout
\begin_layout Description
-dbdt: (nlma,nn,2) : time derivative of pol.
+dbdt\InsetSpace ~
+(nlma,nn,2) time derivative of pol.
magn.
pot.
b
\end_layout
\begin_layout Description
-dcheb: (nn,nn) : dcheb(i,j) = 1st derivative of Ch.
+dcheb\InsetSpace ~
+(nn,nn) dcheb(i,j) = 1st derivative of Ch.
polyn.
i at grid point j
\end_layout
\begin_layout Description
-d2cheb: (nn,nn) : d2cheb(i,j) = 2nd derivative of Ch.
+d2cheb\InsetSpace ~
+(nn,nn) d2cheb(i,j) = 2nd derivative of Ch.
polyn.
i at grid point j
\end_layout
\begin_layout Description
-d3cheb: (nn,nn) : d2cheb(i,j) = 3rd derivative of Ch.
+d3cheb\InsetSpace ~
+(nn,nn) d2cheb(i,j) = 3rd derivative of Ch.
polyn.
i at grid point j
\end_layout
\begin_layout Description
-ddb:(nlma,nn+1) : 2nd rad.
+ddb\InsetSpace ~
+(nlma,nn+1) 2nd rad.
derivative of pol.
mag.
potential b
\end_layout
\begin_layout Description
-ddj:(nlma,nn+1) : 2nd rad.
+ddj\InsetSpace ~
+(nlma,nn+1) 2nd rad.
deriv.
of tor.
magn.
@@ -1807,334 +1872,367 @@
\end_layout
\begin_layout Description
-ddw:(nlma,nn+1) : 0.25 * 2nd radial derivative of pol.
+ddw\InsetSpace ~
+(nlma,nn+1) 0.25 * 2nd radial derivative of pol.
velocity pot.
w
\end_layout
\begin_layout Description
-ddz:(nlma,nn+1) : 0.25 * 2nd radial derivative of tor.
+ddz\InsetSpace ~
+(nlma,nn+1) 0.25 * 2nd radial derivative of tor.
velocity pot.
z
\end_layout
\begin_layout Description
-djdt:(nlma,nn,2) : time derivative of tor.
+djdt\InsetSpace ~
+(nlma,nn,2) time derivative of tor.
magn.
pot.
z
\end_layout
\begin_layout Description
-dpdt:(nlma,nn,2) : time derivative of pressure
+dpdt\InsetSpace ~
+(nlma,nn,2) time derivative of pressure
\end_layout
\begin_layout Description
-dsdt:(nlma,nn,2) : time derivative of temperature (entropy)
+dsdt\InsetSpace ~
+nlma,nn,2) time derivative of temperature (entropy)
\end_layout
\begin_layout Description
-dw:(nlma,nn+1) : 0.50 * radial derivative of pol.
+dw\InsetSpace ~
+(nlma,nn+1) 0.50 * radial derivative of pol.
velocity pot.
w
\end_layout
\begin_layout Description
-dwdt:(nlma,nn,2) : time derivative of pol.
+dwdt\InsetSpace ~
+(nlma,nn,2) time derivative of pol.
velocity pot.
w
\end_layout
\begin_layout Description
-dz:(nlma,nn+1) : 0.50 * radial derivative of tor.
+dz\InsetSpace ~
+(nlma,nn+1) 0.50 * radial derivative of tor.
velocity pot.
z
\end_layout
\begin_layout Description
-dzdt:(nlma,nn,2) : time derivative of tor.
+dzdt\InsetSpace ~
+(nlma,nn,2) time derivative of tor.
velocity pot.
z
\end_layout
\begin_layout Description
-delxh(nn) : horizontal Courant length squared
+delxh\InsetSpace ~
+(nn) horizontal Courant length squared
\end_layout
\begin_layout Description
-delxr(nn) : radial Courant length
+delxr\InsetSpace ~
+(nn) radial Courant length
\end_layout
\begin_layout Description
-difamp: INPUT: amplitude of hyperdiffusivity D=D*(1 + difamp *[(l+1-ldif)/(lmax+
+difamp [INPUT] amplitude of hyperdiffusivity D=D*(1 + difamp *[(l+1-ldif)/(lmax+
1-ldif)]^ldifexp) when l>ldif
\end_layout
\begin_layout Description
-dipfilt: INPUT: If nfilt>0 multiply axial dipole component of B_r on outer
+dipfilt [INPUT] If nfilt>0 multiply axial dipole component of B_r on outer
surface by dipfilt in graphics file
\end_layout
\begin_layout Description
-dj: (nlma,nn+1) : radial deriv.
+dj (nlma,nn+1) radial deriv.
of tor.
magn.
potential (spectral form, see aj)
\end_layout
\begin_layout Description
-dt: : current time step
+dt current time step
\end_layout
\begin_layout Description
-dtchck: SUBROUTINE: controls time step
+dtchck [SUBROUTINE] controls time step
\end_layout
\begin_layout Description
-dth: : Courant time based on horiz.
+dth Courant time based on horiz.
velocity + Alfven veloc.
\end_layout
\begin_layout Description
-dtmax: INPUT: Upper limit on time step (and initial step)
+dtmax [INPUT] Upper limit on time step (and initial step)
\end_layout
\begin_layout Description
-dtmin: : Lower limit on time step (stop when dt < dtmin)
+dtmin Lower limit on time step (stop when dt < dtmin)
\end_layout
\begin_layout Description
-dtold: : Time step of previous iterative step
+dtold Time step of previous iterative step
\end_layout
\begin_layout Description
-dtr: : Courant time based on radial velocity + Alfven veloc.
+dtr Courant time based on radial velocity + Alfven veloc.
\end_layout
\begin_layout Description
-dtstart: INPUT: Initial time step.
+dtstart [INPUT] Initial time step.
If =0, dtmax, or when beginning from restart file, the old dt is taken
\end_layout
\begin_layout Description
-dvpdr: (nja,ni) : = d [r * sin(theta) * v_phi]/dr on gridpoints
+dvpdr\InsetSpace ~
+(nja,ni) = d [r * sin(theta) * v_phi]/dr on gridpoints
\end_layout
\begin_layout Description
-dvpdrc: : dvpdr stored as complex array
+dvpdrc dvpdr stored as complex array
\end_layout
\begin_layout Description
-dvpdp: (nja,ni) : = d [r * sin(theta) * v_phi]/dphi on gridpoints
+dvpdp\InsetSpace ~
+(nja,ni) = d [r * sin(theta) * v_phi]/dphi on gridpoints
\end_layout
\begin_layout Description
-dvpdpc: : dvpdp stored as complex array
+dvpdpc dvpdp stored as complex array
\end_layout
\begin_layout Description
-dvrdp: (nja,ni) : = d [r^2 * v_r]/dphi on gridpoints
+dvrdp\InsetSpace ~
+(nja,ni) = d [r^2 * v_r]/dphi on gridpoints
\end_layout
\begin_layout Description
-dvrdpc: : dvrdp stored as complex array
+dvrdpc dvrdp stored as complex array
\end_layout
\begin_layout Description
-dvrdr: (nja,ni) : = d [r^2 * v_r]/dr on gridpoints
+dvrdr\InsetSpace ~
+(nja,ni) = d [r^2 * v_r]/dr on gridpoints
\end_layout
\begin_layout Description
-dvrdrc: : dvrdr stored as complex array
+dvrdrc dvrdr stored as complex array
\end_layout
\begin_layout Description
-dvrdt: (nja,ni) : = sin(theta) * d [r^2 * v_r]/dtheta on gridpoints
+dvrdt\InsetSpace ~
+(nja,ni) = sin(theta) * d [r^2 * v_r]/dtheta on gridpoints
\end_layout
\begin_layout Description
-dvrdtc: : dvrdt stored as complex array d
+dvrdtc dvrdt stored as complex array d
\end_layout
\begin_layout Description
-vtdp: (nja,ni) : = d [r *sin(theta) * v_theta]/dphi on gridpoints
+vtdp\InsetSpace ~
+(nja,ni) = d [r *sin(theta) * v_theta]/dphi on gridpoints
\end_layout
\begin_layout Description
-dvtdpc: : dvtdp stored as complex array
+dvtdpc dvtdp stored as complex array
\end_layout
\begin_layout Description
-dvtdr: (nja,ni) : = d [r * sin(theta) * v_theta]/dr on gridpoints
+dvtdr\InsetSpace ~
+(nja,ni) = d [r * sin(theta) * v_theta]/dr on gridpoints
\end_layout
\begin_layout Description
-dvtdrc: : dvtdr stored as complex array
+dvtdrc dvtdr stored as complex array
\end_layout
\begin_layout Description
-dw: (nlma,nn+1) : 0.5*times radial deriv of pol.
+dw\InsetSpace ~
+(nlma,nn+1) 0.5*times radial deriv of pol.
velocity potential w
\end_layout
\begin_layout Description
-dz: (nlma,nn+1) : 0.5*times radial deriv of tor.
+dz\InsetSpace ~
+(nlma,nn+1) 0.5*times radial deriv of tor.
velocity potential z
\end_layout
\begin_layout Description
-escale: : scaling factor for energies in output
+escale scaling factor for energies in output
\end_layout
\begin_layout Description
-ek: INPUT: Ekman number
+ek [INPUT] Ekman number
\end_layout
\begin_layout Description
-enb: OUTPUT: magnetic energy
+enb [OUTPUT] magnetic energy
\end_layout
\begin_layout Description
-ens: OUTPUT: thermal energy
+ens [OUTPUT] thermal energy
\end_layout
\begin_layout Description
-enscale: INPUT: in output listing, energies are multiplied by enscale
+enscale [INPUT] in output listing, energies are multiplied by enscale
\end_layout
\begin_layout Description
-ent: OUTPUT: total energy
+ent [OUTPUT] total energy
\end_layout
\begin_layout Description
-env: OUTPUT: kinetic energy
+env [OUTPUT] kinetic energy
\end_layout
\begin_layout Description
-epsc0: INPUT: internal heating rate
+epsc0 [INPUT] internal heating rate
\end_layout
\begin_layout Description
-flmb1: (nlma+..) : r-component of (v x B) term
+flmb1\InsetSpace ~
+(nlma+..) r-component of (v x B) term
\end_layout
\begin_layout Description
-flmb2: (nlma+..) : theta-component of (v x B) term
+flmb2\InsetSpace ~
+(nlma+..) theta-component of (v x B) term
\end_layout
\begin_layout Description
-flmb3: (nlma+..) : phi-component of (v x B) term
+flmb3\InsetSpace ~
+(nlma+..) phi-component of (v x B) term
\end_layout
\begin_layout Description
-flms1: (nlma+..) : r-component of entropy advection term
+flms1\InsetSpace ~
+(nlma+..) r-component of entropy advection term
\end_layout
\begin_layout Description
-flms2: (nlma+..) : theta-component of entropy advection term
+flms2\InsetSpace ~
+(nlma+..) theta-component of entropy advection term
\end_layout
\begin_layout Description
-flms3: (nlma+..) : phi-component of entropy advection term
+flms3\InsetSpace ~
+(nlma+..) phi-component of entropy advection term
\end_layout
\begin_layout Description
-flmw1: (nlma+..) : r-component of v*grad(v) + Lorentz force term
+flmw1\InsetSpace ~
+(nlma+..) r-component of v*grad(v) + Lorentz force term
\end_layout
\begin_layout Description
-flmw2: (nlma+..) : theta-component of v*grad(v) + Lorentz force term
+flmw2\InsetSpace ~
+(nlma+..) theta-component of v*grad(v) + Lorentz force term
\end_layout
\begin_layout Description
-flmw3: (nlma+..) : phi-component of v*grad(v) + Lorentz force term
+flmw3\InsetSpace ~
+(nlma+..) phi-component of v*grad(v) + Lorentz force term
\end_layout
\begin_layout Description
-gauss: (ni) : vector with Gaussian weighting factors, local array in subroutine
+gauss\InsetSpace ~
+(ni) vector with Gaussian weighting factors, local array in subroutine
prep
\end_layout
\begin_layout Description
-gquad: SUBROUTINE: finds zeros and Gauss.
+gquad [SUBROUTINE] finds zeros and Gauss.
weight of assc.
Legendre fct.
\end_layout
\begin_layout Description
-grafile: CHARACT : file name for data on spatial grid for graphics prefix
+grafile [CHARACT] file name for data on spatial grid for graphics prefix
"g." added to outfile
\end_layout
\begin_layout Description
-grav(nn): : gravity at radial levels
+grav\InsetSpace ~
+(nn) gravity at radial levels
\end_layout
\begin_layout Description
-ib(nn,lmax): : pivot array for LU-decomposition of matrix bmat created in
- sgefa, used in sgesl
+ib\InsetSpace ~
+(nn,lmax) pivot array for LU-decomposition of matrix bmat created in sgefa,
+ used in sgesl
\end_layout
\begin_layout Description
-ic: : stepping variable commonly used for steps in colatid.
+ic :stepping variable commonly used for steps in colatid.
\end_layout
\begin_layout Description
-icour: INPUT: Courant criterion is checked each ICOUR'th time step
+icour [INPUT] Courant criterion is checked each ICOUR'th time step
\end_layout
\begin_layout Description
-idiftype: INPUT: controls radial variation of diffusivity, =0: no var.
+idiftype [INPUT] controls radial variation of diffusivity, =0: no var.
\end_layout
\begin_layout Description
-ifaxc: (13) : auxiliary array (factorization) for Chebycheff transform
+ifaxc [13] auxiliary array (factorization) for Chebycheff transform
\end_layout
\begin_layout Description
-ifaxf: (13) : auxiliary array (factorization) for Fourier transform
+ifaxf [13] auxiliary array (factorization) for Fourier transform
\end_layout
\begin_layout Description
-ifbfrz: INPUT: logical, if .T., do not update magnetic field
+ifbfrz [INPUT] logical, if .T., do not update magnetic field
\end_layout
\begin_layout Description
-ifirst: : =1 before first call of time-step checking routine, =0 thereafter
+ifirst =1 before first call of time-step checking routine, =0 thereafter
\end_layout
\begin_layout Description
-iframes: INPUT: write altogether iframes frames on the movie files (see
+iframes [INPUT] write altogether iframes frames on the movie files (see
description under imovopt)
\end_layout
\begin_layout Description
-ifsfrz: INPUT: logical, if .T., do not update temperature (entropy)
+ifsfrz [INPUT] logical, if .T., do not update temperature (entropy)
\end_layout
\begin_layout Description
-ifvfrz: INPUT: logical, if .T., do not update velocity
+ifvfrz [INPUT] logical, if .T., do not update velocity
\end_layout
\begin_layout Description
-ij(nn,lmax): : pivot array for LU-decomposition of matrix ajmat created
- in sgefa, used in sgesl
+ij\InsetSpace ~
+(nn,lmax) pivot array for LU-decomposition of matrix ajmat created in
+ sgefa, used in sgesl
\end_layout
\begin_layout Description
-imagcon: INPUT: <0 imposed poloidal field (l=1,m=0) at ICB >=0 imposed toroidal
+imagcon [INPUT] <0 imposed poloidal field (l=1,m=0) at ICB >=0 imposed toroidal
field (l=2,m=0) at ICB >=10 additionally imposed field at CMB, field is
of same sign and amplitude if imagcon=10 and of opposite sign if imagcon=11
\end_layout
\begin_layout Description
-imovopt: INPUT: three-digit integer number, controls options for generating
+imovopt [INPUT] three-digit integer number, controls options for generating
movie-files.
Last digit>0 - write B_z, W_z (vortic) and T in the equatorial plane on
file with prefix "me." 2nd last digit>0 - write longitud.
@@ -2146,15 +2244,15 @@
\end_layout
\begin_layout Description
-imovct: : counter variable for movie frames
+imovct counter variable for movie frames
\end_layout
\begin_layout Description
-infile: CHA INPUT: name of input file for initial values (restart)
+infile [CHA INPUT] name of input file for initial values (restart)
\end_layout
\begin_layout Description
-init: INPUT: =0 start from dat-file, =1: random initial cond., =-1: hydro.
+init [INPUT] =0 start from dat-file, =1: random initial cond., =-1: hydro.
condition from dat-file, magnetic random >=100: initial temperature perturbatio
n in a single mode l,m.
Here m is given by the last two digits of init and l by the preceding digits.
@@ -2162,195 +2260,202 @@
\end_layout
\begin_layout Description
-ip0(nn): : pivot array for LU-decomposition of matrix p0mat created in sgefa,
+ip0(nn) pivot array for LU-decomposition of matrix p0mat created in sgefa,
used in sgesl
\end_layout
\begin_layout Description
-iprnt: : counting blocks in time iteration sequence with printed output
- created at completion of block is(nn,lmax): : pivot array for LU-decomposition
+iprnt counting blocks in time iteration sequence with printed output created
+ at completion of block is(nn,lmax): : pivot array for LU-decomposition
of matrix smat created in sgefa, used in sgesl
\end_layout
\begin_layout Description
-is0(nn): : pivot array for LU-decomposition of matrix s0mat created in sgefa,
+is0\InsetSpace ~
+(nn) pivot array for LU-decomposition of matrix s0mat created in sgefa,
used in sgesl
\end_layout
\begin_layout Description
-iscale: INPUT: determines which diffusivity is used for scaling of time,
+iscale [INPUT] determines which diffusivity is used for scaling of time,
velocity, energy.
1=viscous, 2=therm., 3=magn.
\end_layout
\begin_layout Description
-istep: : time step counter (routine amh)
+istep time step counter (routine amh)
\end_layout
\begin_layout Description
-istor: : counting superblocks in time iteration sequence, upon completion
- of superblock disk file with data saved
+istor counting superblocks in time iteration sequence, upon completion of
+ superblock disk file with data saved
\end_layout
\begin_layout Description
-ivfilt: INPUT: Apply filter to v_r at radial level ivfilt and right into
+ivfilt [INPUT] Apply filter to v_r at radial level ivfilt and right into
first radial position in graphics file, see nfilt
\end_layout
\begin_layout Description
-iwp(nn,lmax): : pivot array for LU-decomposition of matrix wpmat created
- in sgefa, used in sgesl
+iwp\InsetSpace ~
+(nn,lmax) pivot array for LU-decomposition of matrix wpmat created in
+ sgefa, used in sgesl
\end_layout
\begin_layout Description
-iz(nn,lmax): : pivot array for LU-decomposition of matrix zmat created in
- sgefa, used in sgesl
+iz\InsetSpace ~
+(nn,lmax) pivot array for LU-decomposition of matrix zmat created in sgefa,
+ used in sgesl
\end_layout
\begin_layout Description
-k2k: (nn1) : auxiliary array for Chebycheff transform
+k2k\InsetSpace ~
+(nn1) auxiliary array for Chebycheff transform
\end_layout
\begin_layout Description
-kc: : stepping variable commonly used for steps in radius
+kc stepping variable commonly used for steps in radius
\end_layout
\begin_layout Description
-kcour: : auxiliary variable for time step checking procedure
+kcour auxiliary variable for time step checking procedure
\end_layout
\begin_layout Description
-kbotb: INPUT: magnetic bot condition =1 insulat., =2 perfect cond.
+kbotb [INPUT] magnetic bot condition =1 insulat., =2 perfect cond.
\end_layout
\begin_layout Description
-kbotv: INPUT: mechan.
+kbotv [INPUT] mechan.
bottom condition =1 free, =2 rigid
\end_layout
\begin_layout Description
-kbots: INPUT: thermal bottom condition =1 fixed entropy, =2 flux
+kbots [INPUT] thermal bottom condition =1 fixed entropy, =2 flux
\end_layout
\begin_layout Description
-kei: SUBROUTINE: calculates kinetic energy
+kei [SUBROUTINE] calculates kinetic energy
\end_layout
\begin_layout Description
-kstep: : global time step counter
+kstep global time step counter
\end_layout
\begin_layout Description
-ktops: INPUT: thermal top condition =1 fixed entropy, =2 flux
+ktops [INPUT] thermal top condition =1 fixed entropy, =2 flux
\end_layout
\begin_layout Description
-ktopb: INPUT: magnetic top condition =1 insulat., =2 perfect cond
+ktopb [INPUT:]magnetic top condition =1 insulat., =2 perfect cond
\end_layout
\begin_layout Description
-ktopv: INPUT: mechan.
+ktopv [INPUT] mechan.
top condition =1 free, =2 rigid
\end_layout
\begin_layout Description
-ldif: INPUT: control parameter for hyperdiffusivity, see difamp
+ldif [INPUT] control parameter for hyperdiffusivity, see difamp
\end_layout
\begin_layout Description
-ldifexp: INPUT: control parameter for hyperdiffusivity, see difamp
+ldifexp [INPUT] control parameter for hyperdiffusivity, see difamp
\end_layout
\begin_layout Description
-lm: : stepping variable used to cover all l and m lm = m*(lmax+1)/minc -
- m*(m-minc)/(2*minc) +l-m+1
+lm stepping variable used to cover all l and m lm = m*(lmax+1)/minc - m*(m-minc)
+/(2*minc) +l-m+1
\end_layout
\begin_layout Description
-lmax: : max.
+lmax :max.
harmonic degree, calculated as (nj-1)/3
\end_layout
\begin_layout Description
-logfile: CHARACT : file name for continuous log of enregies and other data
+logfile [CHARACT] file name for continuous log of enregies and other data
prefix "l." added to outfile
\end_layout
\begin_layout Description
-lot: PARAM: =2*nlma (twice the number of harmonic modes)
+lot [PARAM] =2*nlma (twice the number of harmonic modes)
\end_layout
\begin_layout Description
-lpfile: CHARACT : file name continuous log of specified values pre-fix "lp."
+lpfile [CHARACT] file name continuous log of specified values pre-fix "lp."
added to outfile
\end_layout
\begin_layout Description
-lsfile: CHARACT : file name for power spectra of magnetic and kinetic as
+lsfile: [CHARACT] file name for power spectra of magnetic and kinetic as
function of l and m; pre-fix "ls."
\end_layout
\begin_layout Description
-ludc: SUBROUTINE: Chebychev collocation
+ludc: [SUBROUTINE] Chebychev collocation
\end_layout
\begin_layout Description
-mclm: (nlma) : used to unscramble harmonic order m from variable lm
+mclm\InsetSpace ~
+(nlma) :used to unscramble harmonic order m from variable lm
\end_layout
\begin_layout Description
-mclma: (nlma) : = m/minc+1 for given lm (storage order of m)
+mclma\InsetSpace ~
+(nlma) = m/minc+1 for given lm (storage order of m)
\end_layout
\begin_layout Description
-kei: SUBROUTINE: calculates magnetic energy
+kei [SUBROUTINE] calculates magnetic energy
\end_layout
\begin_layout Description
-minc: PARAM: if >1, minc-fold symmetry in longitude assumed
+minc [PARAM] if >1, minc-fold symmetry in longitude assumed
\end_layout
\begin_layout Description
-mmax: : max.
+mmax max.
harmonic order, is the largest integer <= lmax divisible by minc
\end_layout
\begin_layout Description
-movafile: CHARACT: file name for movie data (longitudinal averages) prefix
+movafile [CHARACT] file name for movie data (longitudinal averages) prefix
"ma."
\end_layout
\begin_layout Description
-movefile: CHARACT: file name for movie data in equatorial plane prefix "me."
+movefile [CHARACT] file name for movie data in equatorial plane prefix "me."
\end_layout
\begin_layout Description
-movmfile: CHARACT: file name for movie data in map views prefix "mm."
+movmfile [CHARACT] file name for movie data in map views prefix "mm."
\end_layout
\begin_layout Description
-n, nc: : stepping variables commonly used for steps over Chebycheff polynomia
+n,\InsetSpace ~
+nc stepping variables commonly used for steps over Chebycheff polynomia
\end_layout
\begin_layout Description
-ncp: PARAM: =nja/2 used for storage of points in phi in complex array
+ncp [PARAM] =nja/2 used for storage of points in phi in complex array
\end_layout
\begin_layout Description
-nfilt: INPUT: Apply filter F(l)=exp(-[l/lfilt]^nfilt) to B_r on outer surface
+nfilt [INPUT] Apply filter F(l)=exp(-[l/lfilt]^nfilt) to B_r on outer surface
in graphics output file (if nfilt>0 and alfilt>0) When nfilt>0, alfilt<0,
apply cos-tapered filtered with cutoff at nfilt and taper width |alfilt|
\end_layout
\begin_layout Description
-ngcolat: INPUT: graphics output on each ngcolat'th point in latitude
+ngcolat [INPUT] graphics output on each ngcolat'th point in latitude
\end_layout
\begin_layout Description
-ngform: INPUT: if .ne.
+ngform [INPUT] if .ne.
0, graphics output is written each time a re- start file is (finally) written.
ngform=1 or -1: formatted graphics file, ngform=2: unformatted for ngform=-1
additional comment lines are inserted (this is to look at the file, not
@@ -2358,485 +2463,527 @@
\end_layout
\begin_layout Description
-nglon: INPUT: graphics output for each nglon'th point in longitude
+nglon [INPUT] graphics output for each nglon'th point in longitude
\end_layout
\begin_layout Description
-ngrad: INPUT: graphics output on each ngrad'th radial level ni: PARAM: #
+ngrad [INPUT] graphics output on each ngrad'th radial level ni: PARAM: #
of grid points in colatide , must be even.
\end_layout
\begin_layout Description
-nip1: PARAM: =ni+1 nj: PARAM: # of grid points in longitude, nj/minc must
+nip1 [PARAM] =ni+1 nj: PARAM: # of grid points in longitude, nj/minc must
be multiple of four.
\end_layout
\begin_layout Description
-nja:: : =nj/minc, # of actually needed grid points in phi
+nja =nj/minc, # of actually needed grid points in phi
\end_layout
\begin_layout Description
-njp1: PARAM: =nj+1 nlaf: PARAM: = lmax+1
+njp1 [PARAM] =nj+1 nlaf: PARAM: = lmax+1
\end_layout
\begin_layout Description
-nlafp1: PARAM: = lmax+2 nlm: PARAM: = (mmax+1)*(mmax+2)/2
+nlafp1 [PARAM] = lmax+2 nlm: PARAM: = (mmax+1)*(mmax+2)/2
\end_layout
\begin_layout Description
-nlma: PARAM: # of angular modes employed nlma = mmax*(lmax+1)/minc - mmax*(mmax-
+nlma [PARAM] # of angular modes employed nlma = mmax*(lmax+1)/minc - mmax*(mmax-
minc)/(2*minc) + lmax-mmax+1.
\end_layout
\begin_layout Description
-nlmpa: PARAM: = nlma + mmax/minc + 1
+nlmpa [PARAM] = nlma + mmax/minc + 1
\end_layout
\begin_layout Description
-nlogstep: INPUT: write data on logfile (prefix l.) after each nlogstep steps.
+nlogstep [INPUT] write data on logfile (prefix l.) after each nlogstep steps.
\end_layout
\begin_layout Description
-nmaf: PARAM: = mmax+1 nmafa: PARAM: = mmax/minc+1
+nmaf [PARAM] = mmax+1 nmafa: PARAM: = mmax/minc+1
\end_layout
\begin_layout Description
-nn: PARAM: # of radial grid points, nn-1 must be multiple of 4, and contain
+nn [PARAM] # of radial grid points, nn-1 must be multiple of 4, and contain
no prime factors larger than 5
\end_layout
\begin_layout Description
-nn1: PARAM: =nn-1 nn2: PARAM: =nn-2
+nn1 [PARAM] =nn-1 nn2: PARAM: =nn-2
\end_layout
\begin_layout Description
-nn3: PARAM: =nn-3
+nn3 [PARAM] =nn-3
\end_layout
\begin_layout Description
-nnp1: PARAM: =nn+1
+nnp1 [PARAM] =nn+1
\end_layout
\begin_layout Description
-nnp2: PARAM: =nn+2
+nnp2 [PARAM] =nn+2
\end_layout
\begin_layout Description
-nnaf: PARAM: # of radial chebychev modes, must be <= nn
+nnaf [PARAM] # of radial chebychev modes, must be <= nn
\end_layout
\begin_layout Description
-nnx2: PARAM: =2*nn
+nnx2 [PARAM] =2*nn
\end_layout
\begin_layout Description
-nplog: INPUT: if >0 write velocity values at specific points of the grid
+nplog [INPUT] if >0 write velocity values at specific points of the grid
on separate logfile (prefix "lp.") after every nplog steps (see for arrays
vrpoint, vppoint, vtpoint in subroutine amhd for details)
\end_layout
\begin_layout Description
-nprnt: INPUT: # of printed output blocks created until next data storage
+nprnt [INPUT] # of printed output blocks created until next data storage
for restart
\end_layout
\begin_layout Description
-nps2: PARAM: =(nn+1)/2
+nps2 [PARAM] =(nn+1)/2
\end_layout
\begin_layout Description
-nrp: PARAM: =nja+2 (# of points in phi +2)
+nrp [PARAM] =nja+2 (# of points in phi +2)
\end_layout
\begin_layout Description
-ns2: PARAM: =(nn-1)/2
+ns2 [PARAM] =(nn-1)/2
\end_layout
\begin_layout Description
-nstep: INPUT: # of time steps done until next printed output (total # of
+nstep [INPUT] # of time steps done until next printed output (total # of
time steps is nstep*nprnt*nstor)
\end_layout
\begin_layout Description
-nstor: INPUT: # of data storages before program termination
+nstor [INPUT] # of data storages before program termination
\end_layout
\begin_layout Description
-ntf: PARAM: =3*nja/2+1, used for Fourier transform array trigsf
+ntf [PARAM] =3*nja/2+1, used for Fourier transform array trigsf
\end_layout
\begin_layout Description
-ocorevol: : volume of spherical shell (outer core)
+ocorevol volume of spherical shell (outer core)
\end_layout
\begin_layout Description
-oek: : = 1.
+oek = 1.
/ Ekman number
\end_layout
\begin_layout Description
-oekpm: : = 1.
+oekpm = 1.
/ (Ekman number * Mag.Prandtl number)
\end_layout
\begin_layout Description
-oodt: : = 1.
+oodt = 1.
/ dt (inverse time step)
\end_layout
\begin_layout Description
-oosscl: : = 1.
+oosscl = 1.
/ dt
\end_layout
\begin_layout Description
-opr: : = 1.
+opr = 1.
/ Prandtl number
\end_layout
\begin_layout Description
-outfile: CH INPUT: Name of output files (pre-fixes d.,l.,ls.,g.,me.,ma.,mm., lp.
+outfile [CHARACT INPUT] Name of output files (pre-fixes d.,l.,ls.,g.,me.,ma.,mm.,
+ lp.
added)
\end_layout
\begin_layout Description
-p0mat(nn,nn): : LU-decomposed matrix from Chebycheff collocation of pol.
+p0mat\InsetSpace ~
+(nn,nn) LU-decomposed matrix from Chebycheff collocation of pol.
equation of motion, l=0-term for pressure.
Constructed in ludc, used in amhd
\end_layout
\begin_layout Description
-pbar: SUBROUTINE: Calculates value of assoc.
+pbar [SUBROUTINE] Calculates value of assoc.
Legendre function
\end_layout
\begin_layout Description
-pscale: : scaling pressure in output
+pscale scaling pressure in output
\end_layout
\begin_layout Description
-pr: INPUT: Prandtl number
+pr [INPUT] Prandtl number
\end_layout
\begin_layout Description
-prmag: INPUT: Magnetic Prandtl number
+prmag [INPUT] Magnetic Prandtl number
\end_layout
\begin_layout Description
-prnt: SUBROUTINE: print diagnostic data
+prnt [SUBROUTINE] print diagnostic data
\end_layout
\begin_layout Description
-pscl: : = radtop^2
+pscl = radtop^2
\end_layout
\begin_layout Description
-qi: (ni,5) : array with various coefficients depending on colatid.
+qi\InsetSpace ~
+(ni,5) array with various coefficients depending on colatid.
(look in subroutine prep, loop "do 32 " for details
\end_layout
\begin_layout Description
-qk: (nn,16) : array with various coefficients depending on radius (look
- in subroutine prep for details)
+qk\InsetSpace ~
+(nn,16) array with various coefficients depending on radius (look in subroutin
+e prep for details)
\end_layout
\begin_layout Description
-ql: (nlma,10) : various expressions depending on l and m (look in subrout.
+ql\InsetSpace ~
+(nlma,10) various expressions depending on l and m (look in subrout.
prep, loop "do 35" for details)
\end_layout
\begin_layout Description
-qn: (nn,6) : Chebycheff integrals
+qn\InsetSpace ~
+(nn,6) Chebycheff integrals
\end_layout
\begin_layout Description
-r: (nn) : vector with radial levels, r(1)=radtop, r(nn)=radbot
+r\InsetSpace ~
+(nn) vector with radial levels, r(1)=radtop, r(nn)=radbot
\end_layout
\begin_layout Description
-ra: INPUT: Rayleigh number
+ra [INPUT] Rayleigh number
\end_layout
\begin_layout Description
-rapr: : = Rayleigh number / Prandtl number radbot: : radius of inner boundary
+rapr = Rayleigh number / Prandtl number radbot: : radius of inner boundary
\end_layout
\begin_layout Description
-radratio: INPUT: ratio of inner radius to outer radius
+radratio [INPUT] ratio of inner radius to outer radius
\end_layout
\begin_layout Description
-radtop: : radius of outer boundary
+radtop radius of outer boundary
\end_layout
\begin_layout Description
-rderiv: SUBROUTINE: ?? radial derivative
+rderiv [SUBROUTINE] ?? radial derivative
\end_layout
\begin_layout Description
-rffti: SUBROUTINE: ??
+rffti [SUBROUTINE] ??
\end_layout
\begin_layout Description
-rstfile: CHARACT : file name for data in spectral form ('restart data')
- prefix "d." or "d0.", "d1." ....
+rstfile [CHARACT] file name for data in spectral form ('restart data') prefix
+ "d." or "d0.", "d1." ....
added to outfile
\end_layout
\begin_layout Description
-runid: CHAR*64 : text identifying the run
+runid [CHAR*64] text identifying the run
\end_layout
\begin_layout Description
-rva(nn): : auxiliary array used in prep
+rva\InsetSpace ~
+(nn) auxiliary array used in prep
\end_layout
\begin_layout Description
-rvap(nn): : auxiliary array used in kei, mei
+rvap\InsetSpace ~
+(nn) auxiliary array used in kei, mei
\end_layout
\begin_layout Description
-rvat(nn): : auxiliary array used in kei, mei
+rvat\InsetSpace ~
+(nn) auxiliary array used in kei, mei
\end_layout
\begin_layout Description
-rvb(nn): : auxiliary array used in prep, kei, mei
+rvb\InsetSpace ~
+(nn) auxiliary array used in prep, kei, mei
\end_layout
\begin_layout Description
-rvc(nn): : auxiliary array used in kei, mei
+rvc\InsetSpace ~
+(nn) auxiliary array used in kei, mei
\end_layout
\begin_layout Description
-p: (nlma,nn+1) : pressure (spectral form)
+p\InsetSpace ~
+(nlma,nn+1) pressure (spectral form)
\end_layout
\begin_layout Description
-p00co: : = 4/sqrt(3)
+p00co = 4/sqrt(3)
\end_layout
\begin_layout Description
-prep: SUBROUTINE: parameter input, set up auxiliary arrays, set initial
+prep [SUBROUTINE] parameter input, set up auxiliary arrays, set initial
conditions, etc
\end_layout
\begin_layout Description
-s: (nlma,nn+1) : entropy perturbation (spectral form)
+s\InsetSpace ~
+(nlma,nn+1) entropy perturbation (spectral form)
\end_layout
\begin_layout Description
-sc: (ncp,ni) : sr stored in complex form
+sc\InsetSpace ~
+(ncp,ni) sr stored in complex form
\end_layout
\begin_layout Description
-snlc1: (ncp,ni) : slnr1 stored in complex form
+snlc1\InsetSpace ~
+(ncp,ni) slnr1 stored in complex form
\end_layout
\begin_layout Description
-snlc2: (ncp,ni) : slnr2 stored in complex form
+snlc2\InsetSpace ~
+(ncp,ni) slnr2 stored in complex form
\end_layout
\begin_layout Description
-snlc3: (ncp,ni) : slnr3 stored in complex form
+snlc3\InsetSpace ~
+(ncp,ni) slnr3 stored in complex form
\end_layout
\begin_layout Description
-snlr1: (nrp,ni) : nonlinear term (radial advection) for updating temperature
+snlr1\InsetSpace ~
+(nrp,ni) nonlinear term (radial advection) for updating temperature
\end_layout
\begin_layout Description
-snlr2: (nrp,ni) : nonlinear term (theta advection) for updating temperature
+snlr2\InsetSpace ~
+(nrp,ni) nonlinear term (theta advection) for updating temperature
\end_layout
\begin_layout Description
-snlr3: (nrp,ni) : nonlinear term (phi advection) for updating temperature
-
+snlr3\InsetSpace ~
+(nrp,ni) nonlinear term (phi advection) for updating temperature
\end_layout
\begin_layout Description
-sr: (nrp,ni) : temperature (entropy) on grid points
+sr\InsetSpace ~
+(nrp,ni) temperature (entropy) on grid points
\end_layout
\begin_layout Description
-s0mat(nn,nn): : LU-decomposed matrix from Chebycheff collocation of temperature
+s0mat\InsetSpace ~
+(nn,nn) LU-decomposed matrix from Chebycheff collocation of temperature
equation, l=0-term.
Constructed in ludc, used in amhd
\end_layout
\begin_layout Description
-samp: INPUT: amplitude of initial entropy perturbation
+samp [INPUT] amplitude of initial entropy perturbation
\end_layout
\begin_layout Description
-smat(nn,nn,lmax):: LU-decomposed matrix from Chebycheff collocation of temperatu
-re equation.
+smat\InsetSpace ~
+(nn,nn,lmax) LU-decomposed matrix from Chebycheff collocation of temperature
+ equation.
Built in ludc, used in amhd.
\end_layout
\begin_layout Description
-sr: (nja,ni) : entropy on gridpoints
+sr\InsetSpace ~
+(nja,ni) entropy on gridpoints
\end_layout
\begin_layout Description
-src: : sr stored as complex array
+src sr stored as complex array
\end_layout
\begin_layout Description
-sscl: : = dt stor: SUBROUTINE: store data in restart file
+sscl = dt stor: [SUBROUTINE] store data in restart file
\end_layout
\begin_layout Description
-tei: SUBROUTINE: calculates thermal energy
+tei [SUBROUTINE] calculates thermal energy
\end_layout
\begin_layout Description
-time: : time timediff: : time
+time time timediff: : time
\end_layout
\begin_layout Description
-tipdipole: INPUT: rotate poloidal dipole term when beginning from restart
+tipdipole [INPUT] rotate poloidal dipole term when beginning from restart
file
\end_layout
\begin_layout Description
-tmovnext: : auxiliary variable (next output time) for movie file generation
+tmovnext auxiliary variable (next output time) for movie file generation
\end_layout
\begin_layout Description
-tmovstart: INPUT: time at which to start writing movie-frames on m.*-file
+tmovstart [INPUT] time at which to start writing movie-frames on m.*-file
\end_layout
\begin_layout Description
-tmovstep: INPUT: time increments for writing movie-frames on m.*-file
+tmovstep [INPUT] time increments for writing movie-frames on m.*-file
\end_layout
\begin_layout Description
-tops(0:lmax,0:mmax): INPUT: harmonic coefficients of prescribed temperature
+tops\InsetSpace ~
+(0:lmax,0:mmax) [INPUT] harmonic coefficients of prescribed temperature
(entropy) on outer boundary
\end_layout
\begin_layout Description
-treset: INPUT: (LOGICAL) if true reset time and step counter to zero when
+treset [INPUT; LOGICAL] if true reset time and step counter to zero when
starting from a stored dataset
\end_layout
\begin_layout Description
-trigsc: (nn) : auxiliary array for Chebycheff transform routine created
- in chebi, used in chebtf
+trigsc\InsetSpace ~
+(nn) auxiliary array for Chebycheff transform routine created in chebi,
+ used in chebtf
\end_layout
\begin_layout Description
-trigsf: (ntf) : auxiliary array for Fourier transform routine created in
- fftrig, used in fourtf
+trigsf\InsetSpace ~
+(ntf) auxiliary array for Fourier transform routine created in fftrig,
+ used in fourtf
\end_layout
\begin_layout Description
-tscale: : scaling of time in output
+tscale scaling of time in output
\end_layout
\begin_layout Description
-up: (nja,3) : phi-component of velocity in equatorial plane for three consecutiv
-e radial levels.
+up\InsetSpace ~
+(nja,3) phi-component of velocity in equatorial plane for three consecutive
+ radial levels.
Used in moveout to calculate vorticity
\end_layout
\begin_layout Description
-urdp: (nja) : derivative dv_r/dphi in equatorial plane, used in moveout
- to calculate vorticity
+urdp\InsetSpace ~
+(nja) :derivative dv_r/dphi in equatorial plane, used in moveout to
+ calculate vorticity
\end_layout
\begin_layout Description
-vr: (nja,ni) : = r^2 * v_r on grid points vrc: : vr stored as complex array
+vr\InsetSpace ~
+(nja,ni) = r^2 * v_r on grid points vrc: : vr stored as complex array
\end_layout
\begin_layout Description
-vp: (nja,ni) : = c * sin(theta) * v_phi on grid points vpc: : vp stored
- as complex array
+vp\InsetSpace ~
+(nja,ni) = c * sin(theta) * v_phi on grid points vpc: : vp stored as complex
+ array
\end_layout
\begin_layout Description
-vscale: : scaling of velocity in output
+vscale scaling of velocity in output
\end_layout
\begin_layout Description
-vt: (nja,ni) : = r * sin(theta) * v_theta on grid points
+vt\InsetSpace ~
+(nja,ni) = r * sin(theta) * v_theta on grid points
\end_layout
\begin_layout Description
-vtc: : vt stored as complex array
+vtc vt stored as complex array
\end_layout
\begin_layout Description
-w: (nlma,nn+1) : poloidal velocity potential (spectral form)
+w\InsetSpace ~
+(nlma,nn+1) poloidal velocity potential (spectral form)
\end_layout
\begin_layout Description
-wpmat(nn,nn,lmax): LU-decomposed matrix from Chebycheff collocation of pol.
+wpmat\InsetSpace ~
+(nn,nn,lmax) LU-decomposed matrix from Chebycheff collocation of pol.
equation of motion.
Built in ludc, used in amhd.
\end_layout
\begin_layout Description
-wnlc1: (nja/2,ni): wnlr1 stored in complex form
+wnlc1\InsetSpace ~
+(nja/2,ni) wnlr1 stored in complex form
\end_layout
\begin_layout Description
-wnlc2: (nja/2,ni): wnlr2 stored in complex form
+wnlc2\InsetSpace ~
+(nja/2,ni) wnlr2 stored in complex form
\end_layout
\begin_layout Description
-wnlc3: (nja/2,ni): wnlr3 stored in complex form
+wnlc3\InsetSpace ~
+(nja/2,ni) wnlr3 stored in complex form
\end_layout
\begin_layout Description
-wnlr1: (nja,ni) : nonlinear products for updating w (on grid points)
+wnlr1\InsetSpace ~
+(nja,ni) nonlinear products for updating w (on grid points)
\end_layout
\begin_layout Description
-wnlr2: (nja,ni) : nonlinear products for updating z (on grid points)
+wnlr2\InsetSpace ~
+(nja,ni) nonlinear products for updating z (on grid points)
\end_layout
\begin_layout Description
-wnlr3: (nja,ni) : nonlinear products for updating z (on grid points)
+wnlr3\InsetSpace ~
+(nja,ni) nonlinear products for updating z (on grid points)
\end_layout
\begin_layout Description
-work: (lot,nnp2) : work array used in Fourier and Chebycheff transforms
-
+work\InsetSpace ~
+(lot,nnp2) work array used in Fourier and Chebycheff transforms
\end_layout
\begin_layout Description
-wsave: (nn) : auxiliary array used for Chebycheff transform
+wsave\InsetSpace ~
+(nn) auxiliary array used for Chebycheff transform
\end_layout
\begin_layout Description
-wscl: : = dt * radtop^2
+wscl = dt * radtop^2
\end_layout
\begin_layout Description
-y00: : = 1/sqrt(4*pi)
+y00 = 1/sqrt(4*pi)
\end_layout
\begin_layout Description
-z: (nlma,nn+1) : toroidal velocity potential (spectral form)
+z (nlma,nn+1) toroidal velocity potential (spectral form)
\end_layout
\begin_layout Description
-zscl: : = dt * radtop^2
+zscl = dt * radtop^2
\end_layout
\begin_layout Description
-zmat(nn,nn,lmax):: LU-decomposed matrix from Chebycheff collocation of tor.
+zmat\InsetSpace ~
+(nn,nn,lmax) LU-decomposed matrix from Chebycheff collocation of tor.
equation of motion.
Built in ludc, used in amhd.
@@ -2883,20 +3030,20 @@
\end_layout
\begin_layout Description
-outfile: Name of output files (pre-fixes d., g., l., ls., me., ma., mm., are added)
+outfile Name of output files (pre-fixes d., g., l., ls., me., ma., mm., are added)
\end_layout
\begin_layout Description
-infile: Complete name of file from which initial values are read (restart-file).
+infile Complete name of file from which initial values are read (restart-file).
\end_layout
\begin_layout Description
-runid: arbitrary text of up to 64 characters to describe the model
+runid arbitrary text of up to 64 characters to describe the model
\end_layout
\begin_layout Description
-init: set 1 to start from scratch (random noise initial condition) set 0
+init set 1 to start from scratch (random noise initial condition) set 0
to start from a previous result obtained on the same grid and has been
written into a file named d[0-9].<name> set to a value >= 100 to start from
an initial temperature perturbation of one given mode l,m.
@@ -2907,19 +3054,18 @@
\end_layout
\begin_layout Description
-samp: amplitude of initial perturbation (whether random or single mode)
-
+samp amplitude of initial perturbation (whether random or single mode)
\end_layout
\begin_layout Description
-nstep: do one block of nstep time step before producing a summary printout
+nstep do one block of nstep time step before producing a summary printout
of some diagnostics standard output.
nstep should be even.
\end_layout
\begin_layout Description
-nprnt: do one 'superblock' consisting of nprnt blocks of nstep time steps
+nprnt do one 'superblock' consisting of nprnt blocks of nstep time steps
each, before saving all data in file 'd[0-9].name'.
If nstor=1 there is no number added after the 'd', if nstor>1 the number
is incremented by one for each new superblock, starting with zero.
@@ -2927,7 +3073,7 @@
\end_layout
\begin_layout Description
-nstor: do nstor 'superblocks' consisting of nstep*nprnt time steps before
+nstor do nstor 'superblocks' consisting of nstep*nprnt time steps before
terminating the process.
The total number of time steps is nstep*nprnt*nstor.
nstor must be <=10.
@@ -2935,8 +3081,8 @@
\end_layout
\begin_layout Description
-ngform: Write data at grid points for graphics processing and other post-process
-ing (programs column.f diagnos.f) into file 'g[0-9].<name>' each time a superblock
+ngform Write data at grid points for graphics processing and other post-processi
+ng (programs column.f diagnos.f) into file 'g[0-9].<name>' each time a superblock
is written.
ngform=2: unformatted file, ngform=1: formatted file ngform=0: nothing
written, ngform=-1: comment lines are included into file for easier reading
@@ -2944,21 +3090,21 @@
\end_layout
\begin_layout Description
-ngrad: Output on graphics file for each ngrad'th radial point.
+ngrad Output on graphics file for each ngrad'th radial point.
\end_layout
\begin_layout Description
-ngcolat: Output on graphics file every ngcolat'th point in colatitude.
+ngcolat Output on graphics file every ngcolat'th point in colatitude.
\end_layout
\begin_layout Description
-nglon: Output on graphics file every nglon'th point in longitude.
+nglon Output on graphics file every nglon'th point in longitude.
\end_layout
\begin_layout Description
-nfilt: If>0 apply filter of type F(l)=exp[-(l/alfilt)^nfil] to the radial
+nfilt If>0 apply filter of type F(l)=exp[-(l/alfilt)^nfil] to the radial
component of the magnetic field on the outer radius (kc=1) before writing
data into graphics file (for alfilt >0).
When alfilt<0 then apply filter F(l)=(1+sin(pi*(l-nfilt)/alfilt) as long
@@ -2967,8 +3113,8 @@
\end_layout
\begin_layout Description
-alfilt: See under nfilt ivfilt: If >0 apply the same filter as above to
- the radial velocity at radial level
+alfilt See under nfilt ivfilt: If >0 apply the same filter as above to the
+ radial velocity at radial level
\end_layout
\begin_layout Description
@@ -2977,68 +3123,68 @@
\end_layout
\begin_layout Description
-dipfilt: If nfilt>0 multiply axial dipole component of B_r on outer surface
+dipfilt If nfilt>0 multiply axial dipole component of B_r on outer surface
by dipfilt in graphics output
\end_layout
\begin_layout Description
-nlogstep: write data on logfile (prefix l.) after each nlogstep steps.
+nlogstep write data on logfile (prefix l.) after each nlogstep steps.
\end_layout
\begin_layout Description
-nplog: if >0 write velocity values at specific points of the grid on separate
+nplog if >0 write velocity values at specific points of the grid on separate
logfile (prefix "lp.") after every nplog steps (see for arrays vrpoint,
vppoint, vtpoint in subroutine amhd for details)
\end_layout
\begin_layout Description
-iscale: determines which diffusivity is used for scaling of time, velocity
+iscale determines which diffusivity is used for scaling of time, velocity
and energy.
1=viscous, 2=therm., 3=magn.
\end_layout
\begin_layout Description
-enscale: in output listings, energies are multiplied by enscale treset:
- (LOGICAL) if true reset time and step counter to zero when starting from
- a stored dataset
+enscale in output listings, energies are multiplied by enscale treset: (LOGICAL)
+ if true reset time and step counter to zero when starting from a stored
+ dataset
\end_layout
\begin_layout Description
-tipdipole: when starting calculation without imposed symmetry (minc=1) from
+tipdipole when starting calculation without imposed symmetry (minc=1) from
a data file with symmetry (minc>1), add an equatorial dipole component
with tipdipole times the magnitude of the polar dipole
\end_layout
\begin_layout Description
-amps: Option for rescaling temperature perturbation (from restart file)
- by factor amps (if not equal 1)
+amps Option for rescaling temperature perturbation (from restart file) by
+ factor amps (if not equal 1)
\end_layout
\begin_layout Description
-ampw: Same for poloidal velocity ampz: Same for toroidal velocity
+ampw Same for poloidal velocity ampz: Same for toroidal velocity
\end_layout
\begin_layout Description
-ampb: Same for poloidal magnetic field
+ampb Same for poloidal magnetic field
\end_layout
\begin_layout Description
-ampj: Same for toroidal magnetic field
+ampj Same for toroidal magnetic field
\end_layout
\begin_layout Description
-ifvfrz: (logical) if true, do not update velocity during iteration
+ifvfrz (logical) if true, do not update velocity during iteration
\end_layout
\begin_layout Description
-ifbfrz: (logical) if true, do not update mag.
+ifbfrz (logical) if true, do not update mag.
field during iteration
\end_layout
\begin_layout Description
-ifsfrz: (logical) if true, do not update temperature during iteration
+ifsfrz (logical) if true, do not update temperature during iteration
\end_layout
\begin_layout Standard
@@ -3046,20 +3192,20 @@
\end_layout
\begin_layout Description
-dtmin: Minimum time step (in sec).
+dtmin Minimum time step (in sec).
If the dynamically determined time step becomes less, the program terminates.
\end_layout
\begin_layout Description
-dtmax: Maximum (and usually initial) time step.
+dtmax Maximum (and usually initial) time step.
This must be less than 0.25*ek.
Between dtmax and dtmin the actual time step is controlled by a Courant
criterion (see below).
\end_layout
\begin_layout Description
-dtstart: Initial time step.
+dtstart Initial time step.
If dtmax=0, dtmax is used for the initial time step when init>0 and the
last time step used in the previous run (stored in the restart file) is
used when init=0.
@@ -3067,17 +3213,17 @@
\end_layout
\begin_layout Description
-courfac: controls the contribution of the fluid velocity to the Courant
- time step limit (a larger value leads to smaller dt)
+courfac controls the contribution of the fluid velocity to the Courant time
+ step limit (a larger value leads to smaller dt)
\end_layout
\begin_layout Description
-alffac: controls the contribution of the (modified) Alfven velocity to the
+alffac controls the contribution of the (modified) Alfven velocity to the
Courant time step limit (a larger value leads to smaller dt)
\end_layout
\begin_layout Description
-icour: check Courant criterion after each icour time steps (even number)
+icour check Courant criterion after each icour time steps (even number)
\end_layout
@@ -3086,34 +3232,34 @@
\end_layout
\begin_layout Description
-ra: Rayleigh number (defined with gravity on outer boundary)
+ra Rayleigh number (defined with gravity on outer boundary)
\end_layout
\begin_layout Description
-ek: Ekman number
+ek Ekman number
\end_layout
\begin_layout Description
-pr: Prandtl number
+pr Prandtl number
\end_layout
\begin_layout Description
-prmag: Magnetic Prandtl number
+prmag Magnetic Prandtl number
\end_layout
\begin_layout Description
-radratio: Ratio of inner to outer radius
+radratio Ratio of inner to outer radius
\end_layout
\begin_layout Description
-bpeak: peak value of magnetic field imposed by bound.
+bpeak peak value of magnetic field imposed by bound.
cond.
at ICB (also when imagcon=0, bpeak controls the initial magnetic field:
toroidal when bpeak>0, poloidal dipole when bpeak<0!)
\end_layout
\begin_layout Description
-epsc0: Volumetric rate of internal heating
+epsc0 Volumetric rate of internal heating
\end_layout
\begin_layout Standard
@@ -3121,42 +3267,42 @@
\end_layout
\begin_layout Description
-ktops: thermal boundary condition at CMB.
+ktops thermal boundary condition at CMB.
1-fixed temp, 2-fixed radial heat flow.
(ktops=2 not tested !).
\end_layout
\begin_layout Description
-kbots: thermal boundary condition at ICB.
+kbots thermal boundary condition at ICB.
As above.
\end_layout
\begin_layout Description
-ktopv: velocity condition at CMB.
+ktopv velocity condition at CMB.
1-free, 2-rigid.
\end_layout
\begin_layout Description
-kbotv: velocity condition at ICB.
+kbotv velocity condition at ICB.
As above.
\end_layout
\begin_layout Description
-kbotb: =1 for insulating inner core =2: ideally conducting inner core
+kbotb =1 for insulating inner core =2: ideally conducting inner core
\end_layout
\begin_layout Description
-ktopb: =1 for insulating mantle =2: not implemented! imagcon: <0 imposed
+ktopb =1 for insulating mantle =2: not implemented! imagcon: <0 imposed
poloidal field (l=1,m=0) at ICB >=0 imposed toroidal field (l=2,m=0) at
ICB >=10 imposed toroidal field (l=2,m=0) at both CMB and ICB (same amplitude
and same sign if =10, opposite sign if =11)
\end_layout
\begin_layout Description
-cmb: If >0, thin conducting layer at bottom of mantle (not tested!)
+cmb If >0, thin conducting layer at bottom of mantle (not tested!)
\end_layout
\begin_layout Standard
@@ -3164,15 +3310,15 @@
\end_layout
\begin_layout Description
-difamp: Amplitude of hyperdiffusivities
+difamp Amplitude of hyperdiffusivities
\end_layout
\begin_layout Description
-ldif: Hyperdiffusivites applied for harmonic degrees l >= ldif
+ldif Hyperdiffusivites applied for harmonic degrees l >= ldif
\end_layout
\begin_layout Description
-ldifexp: Exponent for increase of hyperdiffusities with l (analytical details
+ldifexp Exponent for increase of hyperdiffusities with l (analytical details
see definition of ql(lm,11) in prep.f)
\end_layout
@@ -3181,18 +3327,18 @@
\end_layout
\begin_layout Description
-imovopt: three-digit integer number, options for generating movie files
- Last digit>0 - write B_z, W_z (vorticity) and T in the equatorial plane
- on file with prefix "me." 2nd last digit>0 - write longitudinally averaged
- B_phi, j_phi and v_phi on file with prefix "ma." 3rd last digit>0 - write
- B_r at outer surface and B_r and v_r at mid- depth on file with prefix
- "mm." 4th last digit>0 - write spherical harmonic coeffs for poloidal field
- at outer boundary and for velocity potentials at radial level given by
- this digit on file with prefix "cc."
+imovopt three-digit integer number, options for generating movie files Last
+ digit>0 - write B_z, W_z (vorticity) and T in the equatorial plane on file
+ with prefix "me." 2nd last digit>0 - write longitudinally averaged B_phi,
+ j_phi and v_phi on file with prefix "ma." 3rd last digit>0 - write B_r at
+ outer surface and B_r and v_r at mid- depth on file with prefix "mm." 4th
+ last digit>0 - write spherical harmonic coeffs for poloidal field at outer
+ boundary and for velocity potentials at radial level given by this digit
+ on file with prefix "cc."
\end_layout
\begin_layout Description
-iframes: write altogether iframes frames on the movie files tmovstart: time
+iframes write altogether iframes frames on the movie files tmovstart: time
at which to start writing movie-frames tmovstep: time increments for writing
movie-frames
\end_layout
@@ -3542,14 +3688,11 @@
\end_layout
\begin_layout Description
-l.[outfile]:
+l.[outfile] printed every nlogstep time steps one record is printed that
+ contains:
\end_layout
\begin_layout Standard
-printed every nlogstep time steps one record is printed that contains:
-\end_layout
-
-\begin_layout Standard
1) time 2) mean kinetic energy density 3) mean poloidal kinetic energy density
4) mean magnetic energy density 5) mean poloidal magnetic energy density
6) mean axisymmetric toroidal kinetic energy density 7) mean axisymmetric
@@ -3562,17 +3705,13 @@
\end_layout
\begin_layout Description
-ls.[outfile]:
+ls.[outfile] printed each nprint time steps are four records with time being
+ the first variable followed by the spectral power of kinetic and mag- netic
+ energy, respectively, as a function of harmonic degree l, from l=0 to lmax
+ (first two records in a block) and spectral power as function of harmonic
+ order m in the last two records of a block.
\end_layout
-\begin_layout Standard
-printed each nprint time steps are four records with time being the first
- variable followed by the spectral power of kinetic and mag- netic energy,
- respectively, as a function of harmonic degree l, from l=0 to lmax (first
- two records in a block) and spectral power as function of harmonic order
- m in the last two records of a block.
-\end_layout
-
\begin_layout Chapter*
License
\end_layout
More information about the cig-commits
mailing list