Self-Consistent Field

Variable Definitions

     Basic global variables/parameters (see scf.h for declarations):

          ax,ay,az    : accelerations of bodies.
          clm, dlm,   : radial functions used to evaluate expansions.
          elm, flm
          costh       : cosine of polar angular coordinate of bodies.
          cputime     : cpu time (secs) used during the simulation.
          cputime0    : cumulative cpu time at start of run.
          cputime1    : cumulative cpu time at end of run.
          dplm        : values of derivatives of Legendre functions 
                        for bodies.
          dtime       : the timestep.
          fixacc      : option to force conservation of linear
                        momentum by setting acceleration of c.o.m.=0.
          G           : the gravitational constant.
          headline    : identification string for the run.
          inptcoef    : option to read in expansion coefficients.
          irsort      : vector of sorting indices.
          lmax        : number of angular eigenfunctions.
          mass        : masses of bodies.
          nbodies     : total number of bodies.
          nbodsmax    : maximum number of bodies.
          nmax        : number of radial eigenfunctions.
          noutbod     : frequency of system record outputs.
          noutlog     : frequency of outputs to log file.
          nsteps      : number of time-steps to integrate the system.
          one         : the constant 1.
          onesixth    : the constant 1/6.
          outpcoef    : option to write out expansion coefficients.
          phi         : azimuthal angular coordinate of bodies.
          plm         : values of Legendre functions for bodies.
          pi          : the constant pi.
          pot         : potentials of bodies (self-gravity).
          potext      : potentials of bodies (external field).
          r           : radial coordinate of each body.
          selfgrav    : option to turn off (.FALSE.) system self-
                        gravity.
          temp1,temp2,: temporary arrays used to improve the
          temp3,temp4,  efficiency of the calculation of 
          temp5,temp6   self-gravity.
          tnow        : current system time.
          tpos        : current position time.
          tvel        : current velocity time.
          twoopi      : the constant 2./pi.
          ultrass,    : vectors containing values of ultraspherical
          ultrasp1      polynomials for the bodies.
          vx,vy,vz    : velocities of bodies.
          x,y,z       : positions of bodies.
          xi          : radially transformed coordinate of bodies.
          zeroeven    : option to zero out all even terms in the
                        basis function expansions.
          zeroodd     : option to zero out all odd terms in the
                        basis function expansions.

-----------------------------------------------------------------------

   Definitions specific to input/output.

          uterm, upars, ulog, ubodsin,   : logical i/o unit numbers.
            ubodsout,utermfil,uoutcoef,
            uincoef,ubodsel
          parsfile, logfile, ibodfile,   : character names of files.
            obodfile,termfile,outcfile,
            incfile,elfile


=======================================================================