**pspeed
int=**ccd-intensity **vel=**ccd-velocity [parameter=value]

In the following (x,y) are
defined along the major and minor axis of the galaxy resp., and pspeed will
then compute the numerical sums

V(y) = SUM ( vel(x,y) * den(x,y)) / SUM ( den(x,y) ) X(y) = SUM ( x * den(x,y)) / SUM ( den(x,y) )where the summation is along

The output will be a tableV = Omega * SIN(inc) * X

**int=***int-file*- input file, must
be the intensity in
*image(5NEMO)*format [no default]. **vel=***vel-file*- input
file, must be the radial velocity in
*image(5NEMO)*format [no default]. **debugy=***y*- Will produce a table how the numerical integration along only one specific
y value resulted in a pattern speed. For debugging only. [default:
**-1**, i.e. not used].

~/src/pjt/imageauthors original sources

17-Aug-87V1.0: Original createdPJT 1-Nov-87V1.2: more debug output PJT