elongation.f90

Go to the documentation of this file.

!> \file
!> \brief Compute Waist thickness and height in \f$\varphi = 0, \pi\f$ symmetry planes.
 
!> \brief Compute Waist thickness and height in \f$\varphi = 0, \pi\f$ symmetry planes.
!>
!> @param r1 \f$R\f$
!> @param z1 \f$Z\f$
!> @param waist
!> @param height
SUBROUTINE elongation (r1, z1, waist, height)
  USE vmec_main
  IMPLICIT NONE
 
  REAL(rprec), DIMENSION(ns,nzeta,ntheta3,0:1), INTENT(in) :: r1, z1
  REAL(rprec), INTENT(out) :: waist(2), height(2)
  INTEGER :: nv, n1
 
  n1 = 0
  DO nv = 1, nzeta/2+1
     IF (nv.ne.1 .and. nv.ne.nzeta/2+1) cycle
     n1 = n1+1
     waist(n1) = (r1(ns,nv,1,0) + r1(ns,nv,1,1)) - (r1(ns,nv,ntheta2,0) + r1(ns,nv,ntheta2,1))
     height(n1) = 2*maxval(abs(z1(ns,nv,:,0) + z1(ns,nv,:,1)))
  END DO
 
END SUBROUTINE elongation