Introduction to the PEST1 code

J. Manickam
Princeton Plasma Physics Laboratory
Princeton NJ

Fri Nov. 18, 1994

Overview

The Pest1 code determines the stability of a given equilibrium using the delta-W method. It uses the correct kinetic energy and can determine the complete spectrum. In practice we seek the unstable eigenvalues only.

The procedure requires a mapping to the PEST co-ordinates and then computation of delta-W in the plasma region, and in the external vacuum-region to set up the potential and kinetic energy matrices. This is followed by an eigenvalue determination using an inverse vector iteration techinque.

MAPPING

This code maps the equilibrium information onto a new grid appropriate for stability computation. The equilibrium is assumed to be from either

JSOLVER or FLUXEQ

the two Princeton equilibrium codes. The equilibrium file expected is fort.17. The input file is MAPLST. The output is written into two binary files, FORT.41 and FORT.42, which contain the metric as required by the stability code. Some of the quantities are written as text in OUTPUT. There are no graphics.


      INPUT	     	    OUTPUT 		GRAPHICS 	

 FILE  	   TYPE	  	FILE	  TYPE	   	FILE

 MAPLST	  Namelist	OUTPUT	  TEXT		None
 FORT.17  Binary	FORT.41	  Binary
			FORT.41	  Binary


STABILITY

The input file is MODIN. The output is written into two files, FORT.12 a truncated version, and MODOUT. The graphics are in PEST1.PLT.

      INPUT	     	    OUTPUT 		GRAPHICS 	

 FILE  	   TYPE	  	FILE	  TYPE	   	FILE

 MODIN	  Namelist	MODOUT	  TEXT		PEST1.PLT
 FORT.41  Binary	FORT.12	  TEXT
 FORT.42  Binary        

The executables are stored in cfs dir=/000447/public_pest1, they are:

Mapper		god10zb
Stability	go10zb
Vacuum		govaczb
Post-processor	goplotzb

A sample TFTR equilibrium "eqtest" is also stored there.
A sample script, runpest1 which exercises the mapping and stability codes, using the equilibrium file eqtest, can be viewed by clicking HERE

This script file is also available in CFS /000447/public_pest1

If you would like to have access to the PEST-1 codes following these instructions.

manickam@pppl.pppl.gov
Fri Sep 2 13:00:45 EDT 1994