init_plasma Subroutine

public subroutine init_plasma(NA1, ABC, BTOR, RTOR, UPDWN, GP2, AMETR, RHO, SHIF, ELON, TRIA, MU, NE, TE, TI, ZEF, UPL)

variant amy=(btor/q)rho(drho/dr) is a function of "minor radius" r=rh(i). Poloidal magnetic field: B_pol=amy(r)*sqrt(g22/g), where g is determinant of 3D metric tensor and g22 is the (22) element of the tensor, normalized on ABC^4 and ABC^2, correspondingly.

Polinomial approximation of the amy(r):

!!!!!!!!!!!!! spline approximation of plasma profiles !!!!!!!!!!!!!!!!

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: NA1
real(kind=wp), intent(in) :: ABC
real(kind=wp), intent(in) :: BTOR
real(kind=wp), intent(in) :: RTOR
real(kind=wp), intent(in) :: UPDWN
real(kind=wp), intent(in) :: GP2
real(kind=wp), dimension(*) :: AMETR
real(kind=wp), dimension(*) :: RHO
real(kind=wp), dimension(*) :: SHIF
real(kind=wp), dimension(*) :: ELON
real(kind=wp), dimension(*) :: TRIA
real(kind=wp), dimension(*) :: MU
real(kind=wp), dimension(*) :: NE
real(kind=wp), dimension(*) :: TE
real(kind=wp), dimension(*) :: TI
real(kind=wp), dimension(*) :: ZEF
real(kind=wp), dimension(*) :: UPL