Next: Pressure profile Up: Profiles Previous: Profiles

profile

The functional form of the profile is chosen by setting

IQSAVE=0

where

and

and are input quantities.

Local perturbation of the profile
It is possible to locally flatten the and profiles at specified -surfaces.This is done by adding a local polynomial to the function.

The number of rational surfaces is set equal to and the corresponding -values and layer widths are defined through input arrays and . The degree of flattening is determined by , where has no effect on and will make at the rational surface.

IQSAVE=1
This option is used when the profile is specified to be a sum of Chebyshev polynomials. where is the Chebyshev polynomial of order evaluated at . The array are constants obtained from a disk file . This option is generally used when recreating a TRANSP generated profile.

The coefficients are obtained from a DATA file:


      call link('unit11=(qpoly,open)//') 
      read(11,2000)mpolyq
2000  format(i2)
      read(11,2100)(cq(i),i=1,mpolyq)
2100  format(1e12.5)

IQSAVE=2
This option specifies the -profile by providing an array of and values, in a disk file, . The data is expected in the form:


      call link('unit21=(qdata,open)//') 
      read(21,2000)nsurf
2000  format(i2)
      read(21,2100)(psidk(i),qdk(i),i=1,nsurf)
2100  format(2e13.5)

IQSAVE=1,2
will request terminal input to verify the input .


manickam@pppl.pppl.gov
Tue Nov 22 17:17:59 EST 1994