getfsq.f90

Go to the documentation of this file.

 
 
SUBROUTINE getfsq(gcr, gcz, gnormr, gnormz, gnorm, medge)
  USE vmec_main, ONLY: rprec, ns, ns1, mnsize, dp
  USE vmec_params, ONLY: ntmax
  IMPLICIT NONE
 
  INTEGER, INTENT(in) :: medge
  REAL(rprec), INTENT(out) :: gnormr, gnormz
  REAL(rprec), INTENT(in)  :: gnorm
  REAL(rprec), DIMENSION(ns,mnsize*ntmax), INTENT(in) :: gcr, gcz
 
  INTEGER :: jsmax
 
  ! ns1 is ns-1
  ! so  if medge==0, exclude contribution from LCFS
  ! and if medge==1, also take contribution from LCFS into account
  jsmax = ns1 + medge
  gnormr = gnorm * sum(gcr(:jsmax,:)**2.0_dp)
  gnormz = gnorm * sum(gcz(:jsmax,:)**2.0_dp)
 
END SUBROUTINE getfsq