vmec_main.f90

Go to the documentation of this file.

 
MODULE vmec_main
 
  USE vmec_dim
  USE vmec_input
  USE vmec_persistent
  USE vmec_params, ONLY: ndamp
  USE vparams
 
  IMPLICIT NONE
 
  ! 1d-preconditioner-related
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: blam
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: clam
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: dlam
  REAL(rprec), DIMENSION(:,:), ALLOCATABLE :: ard
  REAL(rprec), DIMENSION(:,:), ALLOCATABLE :: arm
  REAL(rprec), DIMENSION(:,:), ALLOCATABLE :: brd
  REAL(rprec), DIMENSION(:,:), ALLOCATABLE :: brm
  REAL(rprec), DIMENSION(:,:), ALLOCATABLE :: azd
  REAL(rprec), DIMENSION(:,:), ALLOCATABLE :: azm
  REAL(rprec), DIMENSION(:,:), ALLOCATABLE :: bzd
  REAL(rprec), DIMENSION(:,:), ALLOCATABLE :: bzm
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: crd
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: sm
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: sp
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: bdamp
  REAL(rprec), DIMENSION(:,:,:,:), ALLOCATABLE :: faclam
  REAL(rprec), DIMENSION(:,:,:,:), ALLOCATABLE :: faclam0
 
  REAL(rprec), DIMENSION(:,:), ALLOCATABLE :: bmin
  REAL(rprec), DIMENSION(:,:), ALLOCATABLE :: bmax
 
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: iotaf
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: phipf
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: chipf
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: phi
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: beta_vol
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: jcuru
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: jcurv
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: jdotb
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: buco
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: bvco
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: bdotgradv
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: equif
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: specw
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: tcon
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: psi
  REAL(rprec), DIMENSION(:,:), ALLOCATABLE :: yellip
  REAL(rprec), DIMENSION(:,:), ALLOCATABLE :: yinden
  REAL(rprec), DIMENSION(:,:), ALLOCATABLE :: ytrian
  REAL(rprec), DIMENSION(:,:), ALLOCATABLE :: yshift
  REAL(rprec), DIMENSION(:,:), ALLOCATABLE :: ygeo
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: overr
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: pres
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: vp
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: jpar2
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: jperp2
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: bdotb
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: vpphi
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: presgrad
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: bucof
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: bvcof
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: chi
 
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: presf
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: chips
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: phips
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: iotas
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: icurv
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: mass
 
  REAL(rprec), DIMENSION(:,:), ALLOCATABLE :: bsqsav
 
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: bredge
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: bpedge
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: bzedge
 
  REAL(rprec), ALLOCATABLE :: xcl0(:)
 
  REAL(rprec), DIMENSION(0:mpol1d,3) :: xmpq
  REAL(rprec), DIMENSION(0:mpol1d) :: faccon
 
  REAL(rprec) :: hs
  REAL(rprec) :: currv
  REAL(rprec) :: aspect
  REAL(rprec) :: ohs
  REAL(rprec) :: voli
  REAL(rprec) :: r00
  REAL(rprec) :: r0scale
  REAL(rprec) :: z00
  REAL(rprec) :: fsqsum0
  REAL(rprec) :: fnorm
  REAL(rprec) :: fsqr = 1.0_dp
  REAL(rprec) :: fsqz = 1.0_dp
  REAL(rprec) :: fsql = 1.0_dp
  REAL(rprec) :: fnorm1
  REAL(rprec) :: fnorml
  REAL(rprec) :: fsqr1
  REAL(rprec) :: fsqz1
  REAL(rprec) :: fsql1
  REAL(rprec) :: fsq
  REAL(rprec) :: fedge
  REAL(rprec) :: wb
  REAL(rprec) :: wp
 
  ! force residuals due to constraint force alone
  REAL(rprec) :: fsqr_con = 1.0_dp
  REAL(rprec) :: fsqz_con = 1.0_dp
  REAL(rprec) :: fsqr1_con
  REAL(rprec) :: fsqz1_con
 
  ! force residuals due to MHD force alone
  REAL(rprec) :: fsqr_mhd = 1.0_dp
  REAL(rprec) :: fsqz_mhd = 1.0_dp
  REAL(rprec) :: fsqr1_mhd
  REAL(rprec) :: fsqz1_mhd
 
  REAL(rprec) :: router
  REAL(rprec) :: rinner
 
  REAL(rprec) :: ftolv
 
  REAL(rprec) :: otav
  REAL(rprec), DIMENSION(ndamp) :: otau
 
  REAL(rprec), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: rmn_bdy
  REAL(rprec), DIMENSION(:,:,:), ALLOCATABLE, TARGET :: zmn_bdy
 
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: bsubu0
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: dbsq
  REAL(rprec), DIMENSION(:), ALLOCATABLE :: rbsq
 
  REAL(rprec) :: rbtor
  REAL(rprec) :: rbtor0
  REAL(rprec) :: ctor
  REAL(rprec) :: delbsq
  REAL(rprec) :: res0
  REAL(rprec) :: delt0r
 
  LOGICAL :: lthreed
  LOGICAL :: lconm1
 
  LOGICAL :: lflip
 
  INTEGER, DIMENSION(:), ALLOCATABLE :: ireflect
  INTEGER :: multi_ns_grid
 
  integer :: niterv
 
  integer :: first
 
  integer :: neqs
  integer :: irzloff
  integer :: iequi
  integer :: ijacob
  integer :: iter1
  integer :: iter2
  integer :: ivac
 
  integer :: vacuum_calls
!  integer :: profil3d_calls
!  integer :: funct3d_calls
  integer :: num_eqsolve_retries
 
END MODULE vmec_main