Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] absi_v1_0.gz(45 Kbytes)
Manuscript Title: SEURAT: a Monte Carlo algorithm for calculating neutral gas transport in non-circular axisymmetric toroidal plasmas.
Authors: D.B. Heifetz, D.E. Post
Program title: SEURAT
Catalogue identifier: ABSI_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 29(1983)287
Programming language: Fortran.
Computer: CRAY-1.
RAM: 244K words
Word size: 64
Keywords: Plasma physics, Thermonuclear fusion, Tokamak, Two-dimensional, Monte carlo, Neutral gas, Transport.
Classification: 19.11.

Nature of problem:
Calculation of mass and energy transport by neutral particles in toroidally symmetric plasmas with non-circular poloidal cross sections.

Solution method:
Monte Carlo simulation of neutral particle flights through a fixed plasma, with scoring by a pseudo-collisional algorithm.

SEURAT has been applied to a wide variety of tomamak simulations, and contains a variety of checks which terminate the computation if it goes outside the range of validity of the model. A number of common input errors are checked for at the beginning of each computation.

Unusual features:
SEURAT is a subroutine package designed for easy incorporation into general transport codes. It is written in 1966 ANSI standard FORTRAN. The machine word length is assumed to be >= 32 bits throughout, except for a pseudo-random number generator included for test comparisons, which assumes word length >= 36 bits. SUER AT uses a binary file generated by the auxiliary code SIGMAV described in section 2. OLYMPUS variable name spelling conventions are used, but no other OLYMPUS conventions or utilities are utilized. A driver is included for stand- alone operation.

Running time:
A typical case, CPC2 described in section 7, required 4.451 CPU s on the CTSS.