regcoil.f90

! Main program

program regcoil

  use global_variables, only: totalTime, outputFilename, general_option
  use init_plasma_mod

  implicit none

  integer :: tic, toc, countrate

  print *,"This is REGCOIL,"
  print *,"a regularized least-squares method for computing stellarator coils."
  call system_clock(tic,countrate)

  call read_input()
  call validate_input()
  call compute_lambda()

  ! Define the position vector and normal vector at each grid point for the surfaces:
  call init_plasma()
  call init_coil_surface()

  ! Initialize some of the vectors and matrices needed:
  call read_bnorm()
  call build_matrices()

  select case (general_option)
  case (1)
     call solve()
  case (2)
     call compute_diagnostics_for_nescout_potential()
  case (3)
     call svd_scan()
  case default
     print *,"Invalid general_option:",general_option
     stop
  end select

  call system_clock(toc)
  totalTime = real(toc-tic)/countrate

  call write_output()
 
  print *,"REGCOIL complete. Total time=",totalTime,"sec."
  print *,"You can run regcoilPlot ",trim(outputFilename)," to plot results."

end program regcoil