         
		               HELP_MGAMS

			F. Hofmann, January 1997
		       updated: 16 January 2003

1) List of Variables
********************************************************************************

Type 1  : time independent variable (valid for the entire scenario) 
Type 2  : time dependent variable (one value for each equilibrium)


Type A  : variable which can be adjusted by the MGAMS user
Type B  : variable used for development of new MGAMS functions, 
	   			 must not be changed by the MGAMS user !
Type C  : unused variable

********************************************************************************

Variable	Type	Function				  Standard Value
********************************************************************************

aipdel		1A	produces offset in Ip reference trace 
			is used to adjust OH1, OH2 at t=0.
			d(OH)/d(aipdel) = -0.53
					
aipgain		1A	plasma current feedback coefficient		0.01

aipipz		1A	produces Ip-proportional offset 
			in vertical position reference (t>0)		0.028

alpha		1A	defines current moment for radial
			position feedback 
			(only for iscramb = 3,5,6,8,9,13,14,
			 15,16,19,25,26)				0

amagic		1A	scales the rampdown time: 
			rampdown_time = rampup_time*(amagic/timefac)	1.0
			if amagic<=0, the rampdown phase is eliminated
			(and flattop is ignored)

boost		1B	modifies the Ip*Z observer:

			if((midplan.ne.8).and.(nfast.ne.0)) then
		        aaa(20,14+38) =  boost
      			aaa(20,26+38) = -boost
       		 	endif 

			typical value: boost = -4000.			0

botlim		2A	if((ilie.le.4).and.(ilia.le.4)), the initial 
			estimate for the plasma current distribution 
			is set up in the rectangular area:
                        eftlim<r<ritlim, botlim<z<toplim   

brmzerb		1A	radial field at t=0, r=rax(1)-deltar, z=-zax(1)

brmzero		1A	radial field at t=0, r=rax(1)-deltar, z=zax(1)

brpzerb		1A	radial field at t=0, r=rax(1)+deltar, z=-zax(1)

brpzero		1A	radial field at t=0, r=rax(1)+deltar, z=zax(1)

bzero		2A	vacuum toroidal magnetic field at R=0.88m	1.43

bzmzerb		1A	vertical field at t=0,r=rax(1)-deltar,z=-zax(1) 

bzmzero		1A	vertical field at t=0,r=rax(1)-deltar,z=zax(1)

bzpzerb		1A	vertical field at t=0,r=rax(1)+deltar,z=-zax(1)

bzpzero		1A	vertical field at t=0,r=rax(1)+deltar,z=zax(1) 

capaj		2A	parameter to vary the peaking factor of the 
			initially assumed plasma current distribution	0.5 

cappa1		2A	plasma elongation, only used for iansha=1 
			(see definition of analytic plasma boundary)

cappa2		C

curfac		1A	scales the plasma current without changing 
			any of the poloidal field coil currents !	1.0

delipz		1A	produces offset in vertical position 
			reference (t>0)

delta1		2A	plasma triangularity, only used for iansha=1 
			(see definition of analytic plasma boundary)

delta2		C

deltar		1B	defines extent of initial quadrupole field: 
			field is specified at 
			r=rze+deltar and r=rze-deltar			0.01

diohdt		1A	time derivative of OH-current (in kA/sec) 
			at times toft(i)  

dissi		2A	parameter determining trade-off between 
			shape accuracy and power dissipation in 
			shaping coils					1.0e-9

dpsfac		1A	scales the radial position reference trace	0.9

dpszero		1A	produces offset in radial position 
			reference (t>0)				

eftlim		2A	if((ilie.le.4).and.(ilia.le.4)), the initial 
			estimate for the plasma current distribution 
			is set up in the rectangular area:
                        eftlim<r<ritlim, botlim<z<toplim 

egain		1A	E-coil gain in M-matrix				3.0

ell1		2A	exponent in expression for p' 
			(see definition of source functions)			

ell2		C

emgain		1B	scaling coefficient for egain,fgain and ohgain	200. 

emm1		2A	exponent in expression for TT' 
			(see definition of source functions) 

emm2		C

f36same		1A	time (in seconds) at which vacuum fields and 
			flux are evaluated for diagnostics purposes	0 

fastm		1A	fast-coil gain in M-matrix

fgain		1B	F-coil gain in M-matrix				2.0

flattop		1A	flattop time of plasma current
			(ignored if amagic=0)

gain		1A	proportional gain for feedback 
			control of density				6.0

gainext		1B	gain in kappa feedback loop,
			(to be used only with iscramb=12,15,16,
			 19,25,26,36)					0.05

gainr		1A	proportional gain for feedback 
			control of radial position			0.25

gainvz		1A	derivative gain for feedback control 
			of vertical position (F-coils)

gainvze		1B	derivative gain for feedback control 
			of vertical position (E-coils)

gainz		1A	proportional gain for feedback control 
			of vertical position (F-coils)

gainze		1B	proportional gain for feedback control 
			of vertical position (E-coils) 

gapin		1A	flux extrapolation distance on inner wall	0.04 

gapout		1A	flux extrapolation distance on outer wall	0.08

ggain		1A	integral gain for feedback control 
			of density					0 

hgain		1A	derivative gain for feedback control 
			of density					0

hlamd1		1A	plasma squareness, only used for iansha=1 
			(see definition of analytic plasma boundary)

hlamd2		C

hpla		1B	plasma half height for inductance 
			calculation at early times			0.1

iansha		2A	if iansha=0, plasma boundary is given by 
			coordinates of boundary points.
			if iansha=1, plasma boundary is specified 
			analytically

icoilon		2A	vector with 16 elements, one for each PF coil 
			sequence of elements is 
			(E1,E2,...E8,F1,F2,...F8). 
			Elements can be either 0 or 1. 
			If an element is 0, the corresponding coil 
			will have zero current. 
			If an element is 1, the corresponding coil 
			will have non-zero current.
			 
ierat		1B	if(ierat.ne.2), only the first 14 elements of	0
			diohdt are used, additional elements are 
			ignored
			if(ierat.eq.2), diohdt can be specified up to
			50 ms.
			if(ierat.eq.3), output file 'fbtchea.dat' is 
			written, which can be used as CHEASE input.
			Note that 'fbtchea.dat' gives the parameters of 
			the final equilibrium of the scenario only.
ievolv		C

if36fb		1A	if B-phi is positive, if36fb must be  		+1
			if B-phi is negative, if36fb must be 		-1 

ifour		1B	if iscramb.ne.12,15,16,19,25,26,36 only used for 
			  mvloop=8
			  if(ifour.ne.8) E2,E3,E6,E7,F1,F4,F5,F8 currents 
			   are computed from spcified Br and Bz fields
			  if(ifour.eq.8) F1,F2,F3,F4,F5,F6,F7,F8 currents
 			   are computed from spcified Br and Bz fields
			if iscramb.eq.12,15,16,19,25,26 or 36,
			  if(mod(ifour,10).eq.4), kappa observer is sum(I_p*|z|)
			  if(mod(ifour,10).eq.5), kappa observer is sum(I_p*z^2)
			  if(mod(ifour,10).eq.6), kappa observer is 
			    psi(top)+psi(bottom)-psi(left)-psi(right)
			  if ifour>10, corrections for the coil currents are
			    applied to the observer and reference

ikriz		1B	modifies vertical position observer:
			if(ikriz.eq.2) aaa(23,j+38)=a3ipz3(j)
			if(ikriz.ne.2) b5(22,20)=-gainvz*0.1*strang
			if(ikriz.eq.2) b5(22,23)=-gainvz*0.1*strang
			if(ikriz.eq.2) efwave(n,23)  = zerefb(n)
ilarg		C

ilia		2A	number of approximate boundary points		20	

ilie		2A	number of exact boundary points

imeas		C	if imeas=1 output file "fbtmeas" is written
                          (input file for the MGAMS code)

inova		1B	if inova.eq.2 or 3, A-mat output 23 is used for slow 	
			  derivative feedback
			if inova.ne.2 and .ne.3, A-mat output 20 is used for
			  slow derivative feedback 
			if inova.eq.3, proportional vertical observer includes 
			  corrections from coil currents
			if inova.ne.3, proportional vertical observer does not 
			  include corrections from coil currents

iohfb		1A	if iohfb=1, Ip will be positive
			if iohfb=-1, Ip will be negative

ipr16		1A	if(ipr16.eq.0) flux loops # 10,11,12,28,29,30  
			are switched off				0  

iprcinc		1A	if iprcinc=0, no vacuum field diagnostic
			if iprcinc=1, vacuum field diagnostic turned on

ipripz		1A	Wobbling frequency: wobfreq = float(ipripz)
			(only for flattop phase)

iprmax		1A	if(iprmax.eq.1), create a second set 
			  of G matrices with vertical and radial 
			  feedback switched off
                        if(iprmax.eq.2), create coil current control
                          matrices for t<ts (ts = typically 10ms) 

isaddl		2A	if isaddl=0, X-points will not be used to 
			  define the plasma boundary
			if isaddl=1, the X-point with the highest flux
			  value defines the last closed flux surface

iscale		1A	determines the set of flux loops and B-probes
                        used for radial position feedback.

			if zmajo1(1) > 0.11 then 
			inboard loop and probe number  = 3+iscale
			outboard loop and probe number = 18-iscale

			if -0.11 < majo1(1) < 0.11 then
			inboard loop and probe number  = 1+iscale
                        outboard loop and probe number = 20-iscale

			if zmajo1(1) < -0.11 then
			inboard loop and probe number  = 37+iscale
                        outboard loop and probe number = 22-iscale
   
			note that iscale cannot exceed 2

iscramb		1A	defines coil combinations for shape feedback
			
			iscramb=1:	
			     vertical:	+mu3-mu4
					+F2-F5          Z=-0.23
					+F3-F6 		Z=0
					+F4-F7          Z=+0.23
			     radial:	+mu3+mu4
			iscramb=2:
			     vertical:	-moh1-mu3+mu4+moh4 
					-F1-F2+F5+F6    Z=-0.23
					-F2-F3+F6+F7	Z=0
					-F3-F4+F7+F8    Z=+0.23
			     radial:	+moh1+mu3+mu4+moh4 
			iscramb=3:	
			     vertical:	-F2-F3+F7+F8	Z>0, Z_opt=0.115
			     radial:	+F2+F4+F6+F8
			     vertical:	-F1-F2+F6+F7	Z<0, Z_opt=-0.115
			     radial:	+F1+F3+F5+F7
			      far and near coils weighted as (1+alpha)/(1-alpha)
			iscramb=4:	
			     vertical:	-F1-F2+F7+F8	Z_opt=0
			     radial:    +F2+F3+F6+F7
			iscramb=5:	
			     vertical:	-F3-F4+F6+F7	Z_opt=0.115
                             radial:	+F2+F4+F6+F8
			      far and near coils weighted as (1+alpha)/(1-alpha)
			iscramb=6:	
			     vertical:	-F3-F4+F5+F6	Z_opt=0
			     radial:	+F1+F3+F6+F8
			      far and near coils weighted as (1+alpha)/(1-alpha)
			iscramb=7:	
			     vertical:	-F3-F4+F6+F7	Z_opt=0.115
			     radial:    +F3+F4+F6+F7
			iscramb=8:	
			     vertical:	-F4-F5+F6+F7	Z_opt=0.23
			     radial:    +F4+F5+F6+F7
			      far and near coils weighted as (1+alpha)/(1-alpha)
			iscramb=9:	
			     vertical:	-F3-F4+F7+F8	Z_opt=0.23
			     radial:    +F4+F5+F6+F7
			      far and near coils weighted as (1+alpha)/(1-alpha)
			iscramb=10:	
			     vertical:	-F5-F6+F7+F8	Z_opt=0.46
			     radial:    +F5+F6+F7+F8
			      far and near coils weighted as (1+alpha)/(1-alpha)
			iscramb=12:
			     vertical: 	-moh1-mu3+mu4+moh4 
					-F1-F2+F5+F6    	Z=-0.23
					-F2-F3+F6+F7		Z=0
					-F3-F4+F7+F8    	Z=+0.23
			     radial: 	+moh1+mu3+moh2+moh3+mu4+moh4
					+F1+F2+F3+F4+F5+F6	Z=-0.23
					+F2+F3+F4+F5+F6+F7	Z=0
					+F3+F4+F5+F6+F7+F8	Z=+0.23
			     kappa:	+moh1-moh2-moh3+moh4 
					+F1-F3-F4+F6		Z=-0.23
					+F2-F4-F5+F7		Z=0
					+F3-F5-F6+F8		Z=+0.23
			iscramb=15:
			     vertical: 	-F3-F4+F6+F7		Z_opt=0.115
			     radial:	+F2+F3+F4+F6+F7+F8
			      farthest and near coils weighted as 
				(1+alpha)/(1-alpha)
			     kappa:	+F2-F4-F6+F8
			      far and near coils weighted as (1-alpha)/(1+alpha)
			iscramb=16:
			     vertical: 	-F3-F4+F5+F6		Z_opt=0
			     radial:	+F2+F3+F4+F5+F6+F7
			      farthest and near coils weighted as 
				(1+alpha)/(1-alpha)
			     kappa:	+F2-F4-F5+F7
			      far and near coils weighted as (1-alpha)/(1+alpha)
			iscramb=19:
			     vertical: 	-F3-F4+F7+F8		Z_opt=0.23
			     radial:	+F3+F4+F5+F6+F7+F8
			      farthest and near coils weighted as 
				(1+alpha)/(1-alpha)
			     kappa:	+F4-F5-F6+F7
			iscramb=25:
			     vertical: 	-F3-F4+F6+F7		Z_opt=0.115
			     radial:	+F2+F3+F4+F6+F7+F8
			      farthest and near coils weighted as 
				(1+alpha)/(1-alpha)
			     kappa:	+F2-F4-F6+F8
			iscramb=26:
			     vertical: 	-F3-F4+F5+F6		Z_opt=0
			     radial:	+F2+F3+F4+F5+F6+F7
			      farthest and near coils weighted as 
				(1+alpha)/(1-alpha)
			     kappa:	+F2-F4-F5+F7
			iscramb=36:
			     vertical: 	-F3-F4+F5+F6		Z_opt=0
			     radial:	+F3+F4+F5+F6
			     kappa:	+F3-F4-F5+F6

			Use iscramb=12,15,16,19,25,26 or 36 for kappa feedback; 
			  in this case set nelz=3, zup=b, zlp=-b (where b is the
			  vertical plasma half-height); the gain is gainext and 
			  the observer is chosen by ifour

istop			produces diagnostic printout			0

itamax		1A	maximum number of iterations in FBT 
			equilibrium calculation				50

iwrida		C

ixdr		1B	if((mvloop.eq.8).and.(ixdr.eq.8)) then		0
                        the A-mat output # 19 is Ip-top (.86,.40)
                        if((mvloop.eq.8).and.(ixdr.eq.9)) then
                        the A-mat output # 19 is Ip-bottom (.86,-.40)

mei		1B	equilibrium parameter				2000

mvloop		1A	if(mvloop.ne.8) normal operation		7
                        if(mvloop.eq.8) doublet   
			if(mvloop.eq.1) 
			    E1 coil is put into hybrid mode.  
			    the OH current is fixed for ton < t < toff.
			    the OH current reference is ioh
			    ton, toff and ioh must be specified 
				"hard-wired" in fbte18.pro.
			    relovo should be zero
			if(mvloop.eq.2) 
			    E1 and E2 coils are put into hybrid mode.  
			    the OH current is fixed for ton < t < toff.
			    the OH current reference is the sum of ioh and ioh1
			    ton, toff, ioh and ioh1 must be specified 
				"hard-wired" in fbte18.pro.
			mvloop=5, 6, 15, 16 are for ECRH feedback control 
			  (expert users only)

ncosei		1B	equilibrium parameter				200

nelz		1A	number of finite elements along z		2

neqtcv		1A	equilibrium number				

nfast		1A	if (nfast.eq.0) fast coil inactive
			if (nfast.ne.0) fast coil active
-----------------------------------------------------------------------------
            To turn FPS on, use nfast=-3, rampt=0.5, ikriz=2
-----------------------------------------------------------------------------
			if (nfast.eq.-2) normal P,D feedback on 
			slow coils, D feedback on fast coil same as 
			for nfast=-1

			if (nfast.eq.-3) same as nfast=-2, but with 
			IpZ observer which is completely decoupled 
			from the fast coil current

			if (nfast.eq.-4) same as nfast=-2, but with 
			IpZ observer which is completely decoupled 
			from the slow coil current moment

npsi		2A	number of flux surfaces in plasma 		10

nr		1B	nr+1 = number of radial mesh points		30

nruns		1A	number of equilibria computed
			(max=40)

nshafa		(FP)	number of poloidal coils
			if (nfast.eq.0) nshafa = 18
			if (nfast.eq.1) nshafa = 19

ngroup		(FP)	nshafa + nvvel

ntmax		1A	maximum # of time steps				1600

numeq		1A	# of equilibria used to create shot scenario

nvvel		1B	number of vessel elements (must be 19 or 38)	19

nz		1B	(2*nz)+1 = number of axial mesh points		32

nzaxel		1A	finite elements are vertically symmetric with
                          respect to the magnetic axis position of the
                          nzaxel-th equilibrium 			

nzaxre		1A	Ip*Z and the kappa observer are measured with respect 
			  to the magnetic axis position of the nzaxre-th 
			  equilibrium 

ohback		1A	scaling coeff. for backoff currents in 
			F-coils to cancel effect of OH currents  	1.1

ohcorr		1A	can be used to cancel slow vertical drift 
			due to OH stray fields: 
			For a plasma at Z=0.23, Ip=166kA, increasing
			ohcorr by 0.02 will decrease the slope by
			10mm/sec
			
ohgain		1A	OH-coil gain in M-matrix			1.0

ohsame		1B	proportional feedback gain to make  
			I(oh1) = I(oh2)					1.5	

ohsami		1B	integral feedback gain to make  
			I(oh1) = I(oh2)					0

omega		2A	parameter to modify li in computed equilibria
			  (see definition of source functions)

placex		2B	equilibrium parameter

placu1		2A	plasma current (kA)
			max |dI/dt| = 15 MA/s

placu2		C

ppal		2A	pressure profile parameter (0 < ppal < 1)

ppfac		2A	scaling coefficient for p' (ppfac > 1.)

psifac		1B	coeff. used to scale radial position 
			reference and measurement			7000.

psirat		2B	if the plasma is diverted, the plasma 
			boundary is defined at psirat*100% of the 
			separatrix flux					1.0  

qzero		2A	q on magnetic axis

rampt		1A	modifies the psi/B mix in the Ip*Z observer 
			(only used for nfast.ne.0  or  mvloop.eq.8)	0.5

rbro		2A	radial coordinates of Br=0 points

rbzo		2A	radial coordinates of Bz=0 points 

relovo		1A	first guess for resistive loop voltage 		1.5

ri		1B	inner boundary of computational grid 		0.58

ritlim		2A	if((ilie.le.4).and.(ilia.le.4)),the initial 
			estimate for the plasma current distribution 
			is set up in the rectangular area:
                        eftlim < r < ritlim, botlim < z < toplim

rlia1		2A	radial coordinates of approximate 
			boundary points

rlia2		C

rlim1		2A	radial coordinates of exact boundary points

rlim2		C

rmajo1		2A	major plasma radius (for iansha=1)

rmajo2		C

rmino1		2A	minor plasma radius (for iansha=1)

rmino2		C

ro		1B	outer boundary of computational grid 		1.18

rshift1		1A	produces Ip-proportional offset in radial 
			position reference. rshift1 is valid 
			at time t=0.

rshift2		1A	produces Ip-proportional offset in radial 
			position reference. rshift2 is valid at time 
			t=toft(nend), the time when the plasma 
			current disappears.  

rshift3		1A	produces (Ip*Ip) proportional offset in 
			radial position reference
			Example: for Ip=317kA, increasing rshift3 by 0.1
			will increase R (and a) by 3.5 mm (this will 
			increase the radial position reference by ~ 70.)
                        rshift3 should be chosen, to 1st approximation, as
			|-------------------------------------------------|
			|  rshift3 = dpszero/(Ip[MA]*Ip[MA])              |
			|-------------------------------------------------|

strki		2B	parameter to discourage large differences 
			between adjacent coil currents			0

testa		2A	convergence criterion for equilibrium 
			calculation					0.0001

timeeq		2A	defines the time of an equilibrium during rampup
			the time of the n-th equilibrium is given by
			t(n) = 0.01 + timefac*(timeeq(n)-timeeq(1))

timefac		1A	scales the rampup time (see under timeeq)	0.8

toft		1B	times (before t=10ms) at which diohdt is 
			specified

toplim		2A	if((ilie.le.4).and.(ilia.le.4)),the initial 
			estimate for the plasma current distribution 
			is set up in the rectangular area:
                        eftlim < r < ritlim, botlim < z < toplim		

ttfac1		2A	scaling coefficient for TT'

ttfac2		C

uc1		C	defines the coefficients of the #14 and #26 probes
			(under the fast coil) in the fast d(IpZ)/dt observer

uc2		1B	weighting coeff. for Ip measurement 
			(sum of B's)					1.0e-8

uc3		1B	weighting coeff. for Ip measurement 
			(sum of elements)				1.0e-8 

veback		1A	scaling coeff. for backoff currents in 
			F-coils to cancel effect of vessel currents	0.99

vscal		1B	weighting coeff. for B-probe measurements	10000.

vsec		2A	prescribed Volt secs of E and F coils for 
			a given equilibrium. If vsec=0, Volt secs 
			are not prescribed.

weitam		2B	parameter to adjust the weights of the 
			approximate plasma boundary points 
			automatically according to their vertical 
			distance from the magnetic axis 		0 

weitex		2B	parameter to adjust the weights of the 
			approximate plasma boundary points 
			automatically according to their vertical 
			distance from the magnetic axis  		1.0	

wpla		1B      plasma half width for inductance calculation
                          at early times     				0.1

wscal		1B	weighting coeff. for flux loop measurements  	10000.

xip		1A	finite element matrix covers the rectangular 
			space:  xip < r < xop,  zlp < z < zup

xop		1A	finite element matrix covers the rectangular 
			space:  xip < r < xop,  zlp < z < zup 

zbro		2A      Z-coordinates of Br=0 points		

zbzo		2A	Z-coordinates of Bz=0 points

zeecorr		1A	produces Ip-proportional Z-offset
			zeecorr should be chosen, to 1st approximation, as
			|-------------------------------------------------|
			|  zeecorr = -delipz/(Ip[MA]*Ip[MA])              |
			|-------------------------------------------------|

zlia1		2A      Z-coordinates of approximate boundary points

zlia2		C

zlim1		2A      Z-coordinates of exact boundary points 

zlim2		C

zlp		1A      finite element matrix cover the rectangular 
			space:  xip < r < xop,  zlp < z < zup  

zmajo1		2A      Z-position of magmetic axis (for iansha=1) 

zmajo2		C

zshift		1A	defines vertical wobble amplitude
			(only for flattop phase)

zu		1B	upper boundary of computational grid = zu
			lower boundary of computational grid = -zu
			  					   zu=0.794839

zup		1A      finite element matrix cover the rectangular 
			space:  xip < r < xop,  zlp < z < zup 



2) Source Functions
********************************************************************************

	source functions are defined as follows:


	ppri=ppal*phell+(1.-ppal)*(phell-phell1),

		where phell=phi**ell, phell1=phell*phi

	ttpr=phemm+ttbe*(phemm1-phemm+omega*(phemm2-phemm1))

		where phemm=phi**emm, phemm1=phemm*phi,phemm2=phemm1*phi



3) Analytic Plasma Boundary
********************************************************************************

	the analytic plasma boundary is defined by:


	r(i) = rmajo+rmino*cos(w(i)+delta*sin(w(i))-hlamd*sin(2.*w(i)))

 	z(i) = zmajo+rmino*cappa*sin(w(i))


4) MGAMS File Structure
********************************************************************************

 TCV_ROOT:[TCV_OPER.MGAMS]
	ASD_AXP.COM	MGAMS command file, called by JMM's matlab interface
        ASD18_AXP.COM	MGAMS command file, current version
	MGAMS18.CLD	defines logicals
	NUB5.EXE	equilibrium code (executable)
	MGA18_AXP.EXE	mgams (executable)
	CONW18.PRO	converts CONMAT18.XDR file to IDL SAVE file
	FBTOUT.TXT	equilibrium code output text file
	MGAOUT.TXT	mgams output text file
	AAMMLL.DAT	mgams output data file
	CONMAT18.XDR	mgams output for shot preparation

	MAKENUB.COM	compiles the equilibrium code
	NUB5.FOR	equilibrium code (source 1)
	NUB5.INC	equilibrium code (source 2)
        LINK18_AXP.COM  compiles mgams 
	GGG18.FOR	mgams (source 1)
	GGG18.INC	mgams (source 2)
	MGAMSREAD.INC	mgams (source 3)

	MGAMSREAD.FOR	mgams (source 4)

	CONMAT18.SAV	mgams output file for shot preparation, 
			to be read by FBTE18.PRO

 USER:[DUTCH.BACKUP.XDR]
	XDRLIB_AXP.OBJ		XDR object file
	TRANSFER_ALPHA.OBJ	XDR object file

 USER:[DUTCH.PUBLIC]
	LOADXDR		converts CONMAT XDR file to IDL SAVE file


 TCV_ROOT:[RECETTES]
	FBTE18.PRO	shot preparation
	FBTE1.PRO	searches for sign changes in computed coil currents 
			and generates waveforms for power supplies. This
			routine is called by FBTE18.PRO


