Table of Contents

Name

potq - query a NEMO potential in the XY plane

Synopsis

potq [parameter=value]

Description

potq queries a potential and forcefield in the XY plane for a number of radii, and reports some quantities of the potential.

Parameters

The following parameters are recognized in any order if the keyword is also given:
potname=
Name of potential. No default.
potpars=
Parameters for potential (1st one is pattern speed)
potfile=
Any optional data file associated with potential []
r=
Radii to sample to find Qt_max [0:2:0.1]
t=
Time to test potential at, if relevant [0.0]
omega=
Use this instead of any returned pattern speed. Not used.
format=
Format used to print numbers [%g] Not used.
niter=
max #iterations to find qt. [10]
eps=
minimal fractional change in qt to abort iterations [0.001]

Examples

The following example prints the radius, Q_b and an estimate of the uncertainty in Q_b based on the convergence speed (eps=):
% potq pfenniger84 r=1:10:1
 1 0.497721   0           10
 2 0.517117   0.000606655 8
 3 0.426776   0.000922551 8
 4 0.275563   0.000832961 8
 5 0.10667    0.000569431 8
 6 0.0305462  0.000926259 6
 7 0.0177048  0.000573299 6
 8 0.0118711  0.000401105 6
 9 0.00860622 0.000298991 6
10 0.00656493 0.000234145 6
where the last column is the number of iterations.

See Also

potccd(1NEMO) , potrot(1NEMO) , potential(5NEMO)

Files

src/orbit/potential    source code

Author

Peter Teuben

Update History


23-Feb-02    V1.0 Created    PJT
8-feb-05    V1.0c  fixed initialization bug     PJT


Table of Contents