de/d3f/wrout_8f90_source.html

SUBROUTINE wrout(bsq, gsqrt, bsubu, bsubv, bsubs, bsupv, bsupu, rzl_array, gc_array, ier_flag)

  USE vmec_main

  USE vparams, p5 => cp5, two => c2p0

  USE vmec_input, ONLY: ns_array, ftol_array

  USE vmec_params

  USE vmercier

  USE vmec_persistent

  USE xstuff

  USE vmec_io

  USE realspace, ONLY: phip, gsqrta=>z1, z1=>z1

  USE vforces, ONLY: bsupua=>brmn_e, bsupva=>czmn_o, bsqa=>bzmn_e,  &

                     bsubsa=>armn_e, bsubua=>azmn_e, bsubva=>armn_o

  USE vacmod, ONLY: potvac

  USE ezcdf

  USE read_wout_mod, ONLY: vn_version, vn_extension, vn_mgrid,      &

    vn_magen, vn_therm, vn_gam, vn_maxr, vn_minr, vn_maxz, vn_fp,   &

    vn_radnod, vn_polmod, vn_tormod, vn_maxmod, vn_maxit,           &

    vn_asym, vn_free, vn_error, vn_aspect, vn_beta,                 &

    vn_pbeta, vn_tbeta, vn_abeta, vn_b0, vn_rbt0, vn_maxmod_nyq,    &

    vn_rbt1, vn_sgs, vn_lar, vn_modb, vn_ctor, vn_amin, vn_rmaj,    &

    vn_vol, vn_ac, vn_ai, vn_am,                                    &

    vn_pmass_type, vn_pcurr_type, vn_piota_type,                    &

    vn_am_aux_s, vn_am_aux_f, vn_ac_aux_s, vn_ac_aux_f,             &

    vn_ai_aux_s, vn_ai_aux_f,                                       &

    vn_ftolv, vn_fsqr, vn_fsqz, vn_fsql,                            &

    vn_pmod, vn_tmod, vn_pmod_nyq, vn_tmod_nyq,                     &

    vn_racc, vn_zacs, vn_racs, vn_zacc, vn_iotaf, vn_qfact,         &

    vn_presf, vn_phi, vn_phipf, vn_jcuru, vn_jcurv, vn_iotah,       &

    vn_chi, vn_chipf,                                               &

    vn_mass, vn_presh, vn_betah, vn_buco, vn_bvco, vn_vp, vn_specw, &

    vn_phip, vn_jdotb, vn_overr, vn_bgrv, vn_merc, vn_mshear,       &

    vn_mwell, vn_mcurr, vn_mgeo, vn_equif, vn_fsq, vn_wdot,         &

    vn_extcur, vn_curlab, vn_rmnc, vn_zmns, vn_lmns, vn_gmnc,       &

    vn_bmnc, vn_bsubumnc, vn_bsubvmnc, vn_bsubsmns,                 &

    vn_bsupumnc, vn_bsupvmnc, vn_rmns, vn_zmnc, vn_lmnc, vn_gmns,   &

    vn_bmns, vn_bsubumns, vn_bsubvmns, vn_bsubsmnc, vn_bsupumns,    &

    vn_bsupvmns, vn_rbc, vn_zbs, vn_rbs, vn_zbc, vn_potvac,         &

    ln_version, ln_extension, ln_mgrid,                             &

    ln_magen, ln_therm, ln_gam, ln_maxr, ln_minr, ln_maxz, ln_fp,   &

    ln_radnod, ln_polmod, ln_tormod, ln_maxmod, ln_maxit,           &

    ln_asym, ln_free, ln_error, ln_aspect, ln_beta,                 &

    ln_pbeta, ln_tbeta, ln_abeta, ln_b0, ln_rbt0, ln_maxmod_nyq,    &

    ln_rbt1, ln_sgs, ln_lar, ln_modb, ln_ctor, ln_amin, ln_rmaj,    &

    ln_mse, ln_thom, ln_next,                                       &

    ln_pmod, ln_tmod, ln_pmod_nyq, ln_tmod_nyq, ln_racc, ln_zacs,   &

    ln_racs, ln_zacc, ln_iotaf, ln_qfact, ln_am, ln_ac, ln_ai,      &

    ln_pmass_type, ln_pcurr_type, ln_piota_type,                    &

    ln_am_aux_s, ln_am_aux_f, ln_ac_aux_s, ln_ac_aux_f,             &

    ln_ai_aux_s, ln_ai_aux_f, ln_chi, ln_chipf,                     &

    ln_presf, ln_phi, ln_phipf, ln_jcuru, ln_jcurv, ln_iotah,       &

    ln_mass, ln_presh, ln_betah, ln_buco, ln_bvco, ln_vp, ln_specw, &

    ln_vol, ln_phip, ln_jdotb, ln_bgrv, ln_merc, ln_mshear,         &

    ln_mwell, ln_mcurr, ln_mgeo, ln_equif, ln_fsq, ln_wdot,         &

    ln_extcur, ln_curlab, ln_rmnc, ln_zmns, ln_lmns, ln_gmnc,       &

    ln_bmnc, ln_bsubumnc, ln_bsubvmnc, ln_bsubsmns,                 &

    ln_bsupumnc, ln_bsupvmnc, ln_rmns, ln_zmnc, ln_lmnc, ln_gmns,   &

    ln_bmns, ln_bsubumns, ln_bsubvmns, ln_bsubsmnc, ln_bsupumns,    &

    ln_bsupvmns, ln_rbc, ln_zbs, ln_rbs, ln_zbc, ln_potvac


  USE safe_open_mod

  USE mgrid_mod


  IMPLICIT NONE


  INTEGER, INTENT(in) :: ier_flag


  ! reverse ns, mnmax for backwards compatibility

  REAL(rprec), DIMENSION(mnmax,ns,3*MAX(ntmax/2,1)), INTENT(inout), TARGET :: rzl_array, gc_array ! only used as temporary storage !

  REAL(rprec), DIMENSION(ns,nznt), INTENT(inout) :: bsq, gsqrt, bsubu, bsubv, bsubs, bsupv, bsupu


  REAL(rprec) :: qfact(ns)


  CHARACTER(LEN=*), PARAMETER, DIMENSION(1) :: &

               r1dim = (/'radius'/),           &

               mn1dim = (/'mn_mode'/),         &

               mn2dim = (/'mn_mode_nyq'/),     &

               mnpddim = (/'mnpd'/),           &

               currg = (/'ext_current'/),      &

               currl = (/'current_label'/)

  CHARACTER(LEN=*), PARAMETER, DIMENSION(2) ::          &

               r2dim = (/'mn_mode',    'radius '    /), &

               r3dim = (/'mn_mode_nyq','radius     '/)


  INTEGER :: j, js, jlk, mn, lk,                              &

             m, n, k, iwout0, n1, nwout, istat, indx1(1)

  REAL(rprec) :: dmult, tcosi, tsini, vversion, sgn, tmult,         &

                 ftolx1

  REAL(rprec), POINTER, DIMENSION(:,:) :: rmnc, rmns, zmns, zmnc, lmns, lmnc

  REAL(rprec), ALLOCATABLE, DIMENSION(:,:) :: gmnc, bmnc, gmns, bmns, &

    bsubumnc, bsubvmnc, bsubsmns, bsubumns, bsubvmns, bsubsmnc

  REAL(rprec), DIMENSION(mnmax) :: rmnc1, zmns1, lmns1,             &

     rmns1, zmnc1, lmnc1, bmodmn, bmodmn1

  REAL(rprec), DIMENSION(:), ALLOCATABLE :: gmn, bmn,               &

     bsubumn, bsubvmn, bsubsmn, bsupumn, bsupvmn

  CHARACTER(LEN=120) :: wout_file

  REAL(rprec), DIMENSION(:), ALLOCATABLE :: xfinal


  ! ELIMINATE THESE EVENTUALLY

  REAL(rprec), ALLOCATABLE, DIMENSION(:,:) :: bsupumnc, bsupumns, bsupvmnc, bsupvmns


  ! THIS SUBROUTINE CREATES THE FILE WOUT.

  ! IT CONTAINS THE CYLINDRICAL COORDINATE SPECTRAL COEFFICIENTS

  ! RMN,ZMN (full), LMN (half_mesh - CONVERTED FROM INTERNAL full REPRESENTATION),

  ! AS WELL AS COEFFICIENTS (ON NYQ MESH) FOR COMPUTED QUANTITIES:

  ! BSQ, BSUPU,V, BSUBU,V, GSQRT (HALF); BSUBS (FULL-CONVERTED IN JXBFORCE)


  ! Pointer assignments for storage arrays

  n1 = max(1,ntmax/2)

  rmnc => rzl_array(:,:,1)            ! store COS(mu-nv) components

  zmns => rzl_array(:,:,1+n1)         ! store SIN(mu-nv)

  lmns => rzl_array(:,:,1+2*n1)       ! store SIN(mu-nv)


  IF (lasym) THEN

     rmns => gc_array(:,:,1)          ! store SIN(mu-nv)

     zmnc => gc_array(:,:,1+n1)       ! store COS(mu-nv)

     lmnc => gc_array(:,:,1+2*n1)     ! store COS(mu-nv)

  END IF


  ALLOCATE (gmn(mnmax_nyq), bmn(mnmax_nyq),                       &

     bsubumn(mnmax_nyq), bsubvmn(mnmax_nyq), bsubsmn(mnmax_nyq), &

     bsupumn(mnmax_nyq), bsupvmn(mnmax_nyq), stat=istat)


  ALLOCATE (gmnc(mnmax_nyq,ns), bmnc(mnmax_nyq,ns),               &

            bsubumnc(mnmax_nyq,ns), bsubvmnc(mnmax_nyq,ns),       &

            bsubsmns(mnmax_nyq,ns), bsupumnc(mnmax_nyq,ns),       &

            bsupvmnc(mnmax_nyq,ns), stat=istat)

  IF (lasym) THEN

    ALLOCATE (gmns(mnmax_nyq,ns), bmns(mnmax_nyq,ns),             &

              bsubumns(mnmax_nyq,ns), bsubvmns(mnmax_nyq,ns),     &

              bsubsmnc(mnmax_nyq,ns), bsupumns(mnmax_nyq,ns),     &

              bsupvmns(mnmax_nyq,ns), stat=istat)

  END IF

  IF (istat .ne. 0) stop 'Error allocating arrays in VMEC WROUT'


  ! IF (nextcur .eq. 0) THEN

  !    DO j = SIZE(extcur), 1, -1

  !      IF (extcur(j) .ne. zero) THEN

  !          nextcur = j

  !          EXIT

  !       END IF

  !    END DO

  ! END IF


  ! ftol info evaluated here!

  indx1=maxloc(ns_array)

  ftolx1=ftol_array(indx1(1))


  ! NYQUIST FREQUENCY REQUIRES FACTOR OF 1/2

  IF (mnyq .ne. 0) cosmui(:,mnyq) = p5*cosmui(:,mnyq)

  IF (nnyq .ne. 0) cosnv(:,nnyq) = p5*cosnv(:,nnyq)


  ! use wout_file as temporary storage for parsing version number to REAL

  wout_file = version_

  READ (wout_file, *) vversion


  wout_file = 'wout_' // trim(input_extension) // '.nc'

  CALL cdf_open(nwout, wout_file, 'w', iwout0)

  IF (iwout0 .ne. 0) stop 'Error opening wout.nc file VMEC WROUT'


  !================================

  ! Define Variables

  !================================

  ! Scalars

  CALL cdf_define(nwout, vn_version, vversion)

  CALL cdf_define(nwout, vn_extension, input_extension)

  CALL cdf_define(nwout, vn_mgrid, mgrid_file)

  CALL cdf_define(nwout, vn_pcurr_type, pcurr_type)

  CALL cdf_define(nwout, vn_pmass_type, pmass_type)

  CALL cdf_define(nwout, vn_piota_type, piota_type)

  CALL cdf_define(nwout, vn_magen, wb)

  CALL cdf_define(nwout, vn_therm, wp)

  CALL cdf_define(nwout, vn_gam, gamma)

  CALL cdf_define(nwout, vn_maxr, rmax_surf)

  CALL cdf_define(nwout, vn_minr, rmin_surf)

  CALL cdf_define(nwout, vn_maxz, zmax_surf)

  CALL cdf_define(nwout, vn_fp, nfp)

  CALL cdf_define(nwout, vn_radnod, ns)

  CALL cdf_define(nwout, vn_polmod, mpol)

  CALL cdf_define(nwout, vn_tormod, ntor)

  CALL cdf_define(nwout, vn_maxmod, mnmax)

  CALL cdf_define(nwout, vn_maxmod_nyq, mnmax_nyq)

  CALL cdf_define(nwout, vn_maxit, iter2)

  CALL cdf_define(nwout, vn_asym, lasym)

  CALL cdf_define(nwout, vn_free, lfreeb)

  CALL cdf_define(nwout, vn_error, ier_flag)

  CALL cdf_define(nwout, vn_aspect, aspect)

  CALL cdf_define(nwout, vn_beta, betatot)

  CALL cdf_define(nwout, vn_pbeta, betapol)

  CALL cdf_define(nwout, vn_tbeta, betator)

  CALL cdf_define(nwout, vn_abeta, betaxis)

  CALL cdf_define(nwout, vn_b0, b0)

  CALL cdf_define(nwout, vn_rbt0, rbtor0)

  CALL cdf_define(nwout, vn_rbt1, rbtor)

  CALL cdf_define(nwout, vn_sgs, nint(signgs))

  CALL cdf_define(nwout, vn_lar, ionlarmor)

  CALL cdf_define(nwout, vn_modb, volavgb)

  CALL cdf_define(nwout, vn_ctor, ctor)

  CALL cdf_define(nwout, vn_amin, aminor_p)

  CALL cdf_define(nwout, vn_rmaj, rmajor_p)

  CALL cdf_define(nwout, vn_vol, volume_p)

  CALL cdf_define(nwout, vn_ftolv, ftolx1)

  CALL cdf_define(nwout, vn_fsqr, fsqr)

  CALL cdf_define(nwout, vn_fsqz, fsqz)

  CALL cdf_define(nwout, vn_fsql, fsql)


  CALL cdf_define(nwout, vn_nextcur, nextcur)

  CALL cdf_define(nwout, vn_extcur, extcur(1:nextcur), dimname=currg)

  CALL cdf_define(nwout, vn_mgmode, mgrid_mode)


  ! 1D Arrays

  CALL cdf_define(nwout, vn_pmod, xm, dimname=mn1dim)

  CALL cdf_setatt(nwout, vn_pmod, ln_pmod)

  CALL cdf_define(nwout, vn_tmod, xn, dimname=mn1dim)

  CALL cdf_setatt(nwout, vn_tmod, ln_tmod)

  CALL cdf_define(nwout, vn_pmod_nyq, xm_nyq, dimname=mn2dim)

  CALL cdf_setatt(nwout, vn_pmod_nyq, ln_pmod_nyq)

  CALL cdf_define(nwout, vn_tmod_nyq, xn_nyq, dimname=mn2dim)

  CALL cdf_setatt(nwout, vn_tmod_nyq, ln_tmod_nyq)


  CALL cdf_define(nwout, vn_racc, raxis_cc(0:ntor), dimname=(/'n_tor'/))

  CALL cdf_setatt(nwout, vn_racc, ln_racc)

  CALL cdf_define(nwout, vn_zacs, zaxis_cs(0:ntor), dimname=(/'n_tor'/))

  CALL cdf_setatt(nwout, vn_zacs, ln_zacs)

  IF (lasym) THEN

     CALL cdf_define(nwout, vn_racs, raxis_cs(0:ntor), dimname=(/'n_tor'/))

     CALL cdf_setatt(nwout, vn_racs, ln_racs)

     CALL cdf_define(nwout, vn_zacc, zaxis_cc(0:ntor), dimname=(/'n_tor'/))

     CALL cdf_setatt(nwout, vn_zacc, ln_zacc)

  END IF


  j = SIZE(am)-1

  CALL cdf_define(nwout, vn_am, am(0:j), dimname=(/'preset'/))

  j = SIZE(ac)-1

  CALL cdf_define(nwout, vn_ac, ac(0:j), dimname=(/'preset'/))

  j = SIZE(ai)-1

  CALL cdf_define(nwout, vn_ai, ai(0:j), dimname=(/'preset'/))


  j = SIZE(am_aux_s)

  CALL cdf_define(nwout, vn_am_aux_s, am_aux_s(1:j), dimname=(/'ndfmax'/))

  j = SIZE(am_aux_f)

  CALL cdf_define(nwout, vn_am_aux_f, am_aux_f(1:j), dimname=(/'ndfmax'/))

  j = SIZE(ai_aux_s)

  CALL cdf_define(nwout, vn_ai_aux_s, ai_aux_s(1:j), dimname=(/'ndfmax'/))

  j = SIZE(ai_aux_f)

  CALL cdf_define(nwout, vn_ai_aux_f, ai_aux_f(1:j), dimname=(/'ndfmax'/))

  j = SIZE(ac_aux_s)

  CALL cdf_define(nwout, vn_ac_aux_s, ac_aux_s(1:j), dimname=(/'ndfmax'/))

  j = SIZE(ac_aux_f)

  CALL cdf_define(nwout, vn_ac_aux_f, ac_aux_f(1:j), dimname=(/'ndfmax'/))


  CALL cdf_define(nwout, vn_iotaf, iotaf(1:ns), dimname=r1dim)

  CALL cdf_setatt(nwout, vn_iotaf, ln_iotaf)


  qfact=huge(qfact)

  WHERE (iotaf(1:ns) .NE. zero) qfact=one/iotaf(1:ns)


  CALL cdf_define(nwout, vn_qfact, qfact(1:ns), dimname=r1dim)

  CALL cdf_setatt(nwout, vn_iotaf, ln_qfact)


  CALL cdf_define(nwout, vn_presf, presf,       dimname=r1dim)

  CALL cdf_setatt(nwout, vn_presf, ln_presf, units='Pa')

  CALL cdf_define(nwout, vn_phi,   phi,         dimname=r1dim)

  CALL cdf_setatt(nwout, vn_phi,   ln_phi, units='wb')

  CALL cdf_define(nwout, vn_phipf, phipf,       dimname=r1dim)

  CALL cdf_setatt(nwout, vn_phipf, ln_phipf)

  CALL cdf_define(nwout, vn_chi,   chi,         dimname=r1dim)

  CALL cdf_setatt(nwout, vn_chi,   ln_chi, units='wb')

  !CALL cdf_define(nwout, vn_chipf, phipf,       dimname=r1dim) ! TODO: wrong quantity !!!

  CALL cdf_define(nwout, vn_chipf, chipf,       dimname=r1dim)

  CALL cdf_setatt(nwout, vn_chipf, ln_chipf)

  CALL cdf_define(nwout, vn_jcuru, jcuru,       dimname=r1dim)

  CALL cdf_define(nwout, vn_jcurv, jcurv,       dimname=r1dim)


  CALL cdf_define(nwout, vn_iotah, iotas(1:ns), dimname=r1dim)

  CALL cdf_setatt(nwout, vn_iotah, ln_iotah)

  CALL cdf_define(nwout, vn_mass,  mass,        dimname=r1dim)

  CALL cdf_setatt(nwout, vn_mass,  ln_mass)

  CALL cdf_define(nwout, vn_presh, pres(1:ns),  dimname=r1dim)

  CALL cdf_setatt(nwout, vn_presh, ln_presh, units='Pa')

  CALL cdf_define(nwout, vn_betah, beta_vol,    dimname=r1dim)

  CALL cdf_define(nwout, vn_buco,  buco,        dimname=r1dim)

  CALL cdf_define(nwout, vn_bvco,  bvco,        dimname=r1dim)

  CALL cdf_define(nwout, vn_vp,    vp(1:ns),    dimname=r1dim)

  CALL cdf_define(nwout, vn_specw, specw,       dimname=r1dim) ! NOTE: specw on full grid

  CALL cdf_define(nwout, vn_phip,  phips(1:ns), dimname=r1dim)

  CALL cdf_define(nwout, vn_overr, overr(1:ns), dimname=r1dim)


  CALL cdf_define(nwout, vn_jdotb, jdotb, dimname=r1dim)

  CALL cdf_define(nwout, vn_bgrv, bdotgradv, dimname=r1dim)


  CALL cdf_define(nwout, vn_merc,   dmerc,  dimname=r1dim)

  CALL cdf_define(nwout, vn_mshear, dshear, dimname=r1dim)

  CALL cdf_define(nwout, vn_mwell,  dwell,  dimname=r1dim)

  CALL cdf_define(nwout, vn_mcurr,  dcurr,  dimname=r1dim)

  CALL cdf_define(nwout, vn_mgeo,   dgeod,  dimname=r1dim)

  CALL cdf_define(nwout, vn_equif,  equif,  dimname=r1dim)


  IF (lfreeb .and. nextcur.gt.0 .and. ALLOCATED(curlabel)) THEN

     CALL cdf_define(nwout, vn_curlab, curlabel(1:nextcur), dimname=currl)

     ! SAL for potvac

     CALL cdf_define(nwout, vn_potvac, potvac, dimname=mnpddim)

  ENDIF


  ! 2D Arrays

  CALL cdf_define(nwout, vn_rmnc, rmnc, dimname=r2dim)

  CALL cdf_setatt(nwout, vn_rmnc, ln_rmnc, units='m')

  CALL cdf_define(nwout, vn_zmns, zmns, dimname=r2dim)

  CALL cdf_setatt(nwout, vn_zmns, ln_zmns, units='m')

  CALL cdf_define(nwout, vn_lmns, lmns, dimname=r2dim)

  CALL cdf_setatt(nwout, vn_lmns, ln_lmns)

  CALL cdf_define(nwout, vn_gmnc, gmnc, dimname=r3dim)

  CALL cdf_setatt(nwout, vn_gmnc, ln_gmnc)

  CALL cdf_define(nwout, vn_bmnc, bmnc, dimname=r3dim)

  CALL cdf_setatt(nwout, vn_bmnc, ln_bmnc)

  CALL cdf_define(nwout, vn_bsubumnc, bsubumnc, dimname=r3dim)

  CALL cdf_setatt(nwout, vn_bsubumnc, ln_bsubumnc)

  CALL cdf_define(nwout, vn_bsubvmnc, bsubvmnc, dimname=r3dim)

  CALL cdf_setatt(nwout, vn_bsubvmnc, ln_bsubvmnc)

  CALL cdf_define(nwout, vn_bsubsmns, bsubsmns, dimname=r3dim)

  CALL cdf_setatt(nwout, vn_bsubsmns, ln_bsubsmns)


  ! ELIMINATE THESE EVENTUALLY: DON'T NEED THEM - CAN COMPUTE FROM GSQRT

  CALL cdf_define(nwout, vn_bsupumnc, bsupumnc, dimname=r3dim)

  CALL cdf_define(nwout, vn_bsupvmnc, bsupvmnc, dimname=r3dim)


  IF (lasym) then

     CALL cdf_define(nwout, vn_rmns, rmns, dimname=r2dim)

     CALL cdf_setatt(nwout, vn_rmns, ln_rmns, units='m')

     CALL cdf_define(nwout, vn_zmnc, zmnc, dimname=r2dim)

     CALL cdf_setatt(nwout, vn_zmnc, ln_zmnc, units='m')

     CALL cdf_define(nwout, vn_lmnc, lmnc, dimname=r2dim)

     CALL cdf_setatt(nwout, vn_lmnc, ln_lmnc)

     CALL cdf_define(nwout, vn_gmns, gmns, dimname=r3dim)

     CALL cdf_setatt(nwout, vn_gmns, ln_gmns)

     CALL cdf_define(nwout, vn_bmns, bmns, dimname=r3dim)

     CALL cdf_setatt(nwout, vn_bmns, ln_bmns)

     CALL cdf_define(nwout, vn_bsubumns, bsubumns, dimname=r3dim)

     CALL cdf_setatt(nwout, vn_bsubumns, ln_bsubumns)

     CALL cdf_define(nwout, vn_bsubvmns, bsubvmns, dimname=r3dim)

     CALL cdf_setatt(nwout, vn_bsubvmns, ln_bsubvmns)

     CALL cdf_define(nwout, vn_bsubsmnc, bsubsmnc, dimname=r3dim)

     CALL cdf_setatt(nwout, vn_bsubsmnc, ln_bsubsmnc)


     ! ELIMINATE THESE EVENTUALLY: DON'T NEED THEM

     CALL cdf_define(nwout, vn_bsupumns, bsupumns, dimname=r3dim)

     CALL cdf_define(nwout, vn_bsupvmns, bsupvmns, dimname=r3dim)

  end if


  !================================

  ! Write Variables

  !================================

  ! Scalars

  CALL cdf_write(nwout, vn_version, vversion)

  CALL cdf_write(nwout, vn_extension, input_extension)

  CALL cdf_write(nwout, vn_mgrid, mgrid_file)

  CALL cdf_write(nwout, vn_pcurr_type, pcurr_type)

  CALL cdf_write(nwout, vn_piota_type, piota_type)

  CALL cdf_write(nwout, vn_pmass_type, pmass_type)

  CALL cdf_write(nwout, vn_magen, wb)

  CALL cdf_write(nwout, vn_therm, wp)

  CALL cdf_write(nwout, vn_gam, gamma)

  CALL cdf_write(nwout, vn_maxr, rmax_surf)

  CALL cdf_write(nwout, vn_minr, rmin_surf)

  CALL cdf_write(nwout, vn_maxz, zmax_surf)

  CALL cdf_write(nwout, vn_fp, nfp)

  CALL cdf_write(nwout, vn_radnod, ns)

  CALL cdf_write(nwout, vn_polmod, mpol)

  CALL cdf_write(nwout, vn_tormod, ntor)

  CALL cdf_write(nwout, vn_maxmod, mnmax)

  CALL cdf_write(nwout, vn_maxmod_nyq, mnmax_nyq)

  CALL cdf_write(nwout, vn_maxit, iter2)

  CALL cdf_write(nwout, vn_asym, lasym)

  CALL cdf_write(nwout, vn_free, lfreeb)

  CALL cdf_write(nwout, vn_error, ier_flag)


  CALL cdf_write(nwout, vn_aspect, aspect)

  CALL cdf_write(nwout, vn_beta, betatot)

  CALL cdf_write(nwout, vn_pbeta, betapol)

  CALL cdf_write(nwout, vn_tbeta, betator)

  CALL cdf_write(nwout, vn_abeta, betaxis)

  CALL cdf_write(nwout, vn_b0, b0)

  CALL cdf_write(nwout, vn_rbt0, rbtor0)

  CALL cdf_write(nwout, vn_rbt1, rbtor)

  CALL cdf_write(nwout, vn_sgs, nint(signgs))

  CALL cdf_write(nwout, vn_lar, ionlarmor)

  CALL cdf_write(nwout, vn_modb, volavgb)

  CALL cdf_write(nwout, vn_ctor, ctor/mu0)

  CALL cdf_write(nwout, vn_amin, aminor_p)

  CALL cdf_write(nwout, vn_rmaj, rmajor_p)

  CALL cdf_write(nwout, vn_vol, volume_p)

  CALL cdf_write(nwout, vn_ftolv, ftolx1)

  CALL cdf_write(nwout, vn_fsql, fsql)

  CALL cdf_write(nwout, vn_fsqr, fsqr)

  CALL cdf_write(nwout, vn_fsqz, fsqz)


  CALL cdf_write(nwout, vn_nextcur, nextcur)

  IF (nextcur .gt. 0) THEN

     CALL cdf_write(nwout, vn_extcur, extcur(1:nextcur))

     CALL cdf_write(nwout, vn_mgmode, mgrid_mode)

  ENDIF

  IF (lfreeb) THEN

     ! TODO: write current labels

     CALL cdf_write(nwout, vn_potvac, potvac)

  END IF


  ! 1D Arrays

  CALL cdf_write(nwout, vn_pmod, xm)

  CALL cdf_write(nwout, vn_tmod, xn)

  CALL cdf_write(nwout, vn_pmod_nyq, xm_nyq)

  CALL cdf_write(nwout, vn_tmod_nyq, xn_nyq)


  ! here start the convert-to-regular-Fourier-output part


  ALLOCATE (xfinal(neqs), stat=js) ! re-use js for allocation return code

  IF (js .ne. 0) stop 'Allocation error for xfinal in WROUT!'

  xfinal = xc ! ignore passed rzl_array !!!

  ! TODO: scalxc ???

  ! --> rzl_array is re-used to store temporary  symmetric rmnc, ... !

  ! -->  gc_array is re-used to store temporary asymmetric rmns, ... ! (for lasym)


  ! MUST CONVERT m=1 MODES... FROM INTERNAL TO PHYSICAL FORM

  ! Extrapolation of m=0 Lambda (cs) modes, which are not evolved at j=1, done in CONVERT


  IF (lthreed) CALL convert_sym  (xfinal(1+mns*(rss-1)), xfinal(1+irzloff+mns*(zcs-1)))

  IF (lasym)   CALL convert_asym (xfinal(1+mns*(rsc-1)), xfinal(1+irzloff+mns*(zcc-1)))


  ! CONVERT TO rmnc, zmns, lmns, etc EXTERNAL representation (without internal mscale, nscale)

  ! IF B^v ~ phip + lamu, MUST DIVIDE BY phipf(js) below to maintain old-style format

  radius1: DO js = 1, ns


     CALL convert (rmnc1, zmns1, lmns1, rmns1, zmnc1, lmnc1, xfinal, js)


     rmnc(:,js) = rmnc1(:)

     zmns(:,js) = zmns1(:)

     lmns(:,js) = (lmns1(:)/phipf(js)) * lamscale

     IF (lasym) THEN

        rmns(:,js) = rmns1(:)

        zmnc(:,js) = zmnc1(:)

        lmnc(:,js) = (lmnc1(:)/phipf(js)) * lamscale

     END IF


  END DO radius1


  DEALLOCATE (xfinal)


  ! INTERPOLATE LAMBDA ONTO HALF-MESH FOR BACKWARDS CONSISTENCY WITH EARLIER VERSIONS OF VMEC

  ! AND SMOOTHS POSSIBLE UNPHYSICAL "WIGGLE" ON RADIAL MESH

  WHERE (nint(xm) .le. 1) lmns(:,1) = lmns(:,2) ! constant extrapolation of lambda

  DO js = ns,2,-1

     WHERE (mod(nint(xm),2) .eq. 0)

        lmns(:,js) = p5*(lmns(:,js) + lmns(:,js-1))

     ELSEWHERE

        lmns(:,js) = p5*(sm(js)*lmns(:,js) + sp(js-1)*lmns(:,js-1))

     END WHERE

  END DO


  lmns(:,1) = 0 ! NOTE: needed constant extrapolation above only in averaging step!

  raxis_cc(0:ntor) = rmnc(1:ntor+1,1)

  zaxis_cs(0:ntor) = zmns(1:ntor+1,1)


  IF (lasym) then

     WHERE (nint(xm) .le. 1) lmnc(:,1) = lmnc(:,2)

     DO js = ns,2,-1

        WHERE (mod(nint(xm),2) .eq. 0)

           lmnc(:,js) = p5*(lmnc(:,js) + lmnc(:,js-1))

        ELSEWHERE

           lmnc(:,js) = p5*(sm(js)*lmnc(:,js) + sp(js-1)*lmnc(:,js-1))

        END WHERE

     END DO


     lmnc(:,1) = 0 ! NOTE: needed constant extrapolation above only in averaging step!

     raxis_cs(0:ntor) = rmns(1:ntor+1,1)

     zaxis_cc(0:ntor) = zmnc(1:ntor+1,1)

  end if


  ! here comes the realspace-to-Fourier transform part...


  ! COMPUTE |B| = SQRT(|B|**2) and store in bsq, bsqa

  DO js = 2, ns

     bsq(js,:nznt) = sqrt(2*abs(bsq(js,:nznt)-pres(js)))

  END DO


  tmult = p5/r0scale**2

  !SPH: FIXED THIS 03-05-07 TO CALL symmetrization routine

  IF (lasym) THEN

     ! Changed integration norm in fixaray, SPH012314

     tmult = 2*tmult

     bsubs(1,:) = 0

     CALL symoutput (bsq,   gsqrt,  bsubu,  bsubv,  bsupu,   &

                     bsupv,  bsubs,                          &

                     bsqa,  gsqrta, bsubua, bsubva, bsupua,  &

                     bsupva, bsubsa )

  END IF


  ! Fourier-transform derived quantities for each surface individually

  radius2: DO js = 2, ns

     gmn = 0

     bmn = 0


     bsubumn = 0

     bsubvmn = 0

     bsubsmn = 0


     bsupumn = 0

     bsupvmn = 0


     mn2: DO mn = 1, mnmax_nyq

        n = nint(xn_nyq(mn))/nfp

        m = nint(xm_nyq(mn))

        n1 = abs(n)


        dmult = mscale(m)*nscale(n1)*tmult

        IF (m.eq.0 .or. n.eq.0) dmult = 2*dmult


        sgn = sign(1, n)


        lk = 0

        DO j = 1, ntheta2

           DO k = 1, nzeta

              lk = lk + 1


              ! cos(mu - nv)

              tcosi = dmult*(cosmui(j,m)*cosnv(k,n1) + sgn*sinmui(j,m)*sinnv(k,n1))


              ! sin(mu - nv)

              tsini = dmult*(sinmui(j,m)*cosnv(k,n1) - sgn*cosmui(j,m)*sinnv(k,n1))


              bmn(mn)     = bmn(mn)     + tcosi*bsq(js,lk)

              gmn(mn)     = gmn(mn)     + tcosi*gsqrt(js,lk)

              bsubumn(mn) = bsubumn(mn) + tcosi*bsubu(js,lk)

              bsubvmn(mn) = bsubvmn(mn) + tcosi*bsubv(js,lk)

              bsubsmn(mn) = bsubsmn(mn) + tsini*bsubs(js,lk)

              bsupumn(mn) = bsupumn(mn) + tcosi*bsupu(js,lk)

              bsupvmn(mn) = bsupvmn(mn) + tcosi*bsupv(js,lk)

           END DO

        END DO

     END DO mn2


     ! needed for freeb_data below

     IF (js .eq. ns/2) bmodmn  = bmn(1:mnmax)

     IF (js .eq. ns  ) bmodmn1 = bmn(1:mnmax)


     gmnc(:,js) = gmn(:)

     bmnc(:,js) = bmn(:)

     bsubumnc(:,js) = bsubumn(:)

     bsubvmnc(:,js) = bsubvmn(:)

     bsubsmns(:,js) = bsubsmn(:)

     bsupumnc(:,js) = bsupumn(:)

     bsupvmnc(:,js) = bsupvmn(:)

  END DO radius2


  ! endpoint values at magnetic axis

  gmnc(:,1) = 0

  bmnc(:,1) = 0


  bsubumnc(:,1) = 0

  bsubvmnc(:,1) = 0


  ! NOTE: This assumes that bsubs (from which bsubsmns is computed) is on the full grid.

  ! However, since crmn_o is passed as bsubs in fileout(),

  ! HERE, bsubs is ACTUALLY on the HALF-grid !!! WTF ???

  ! (The full-grid bsubs array is a LOCAL variable in jxbforce().)

  bsubsmns(:,1) = 2*bsubsmns(:,2) - bsubsmns(:,3) ! extrapolation on full grid


  bsupumnc(:,1) = 0

  bsupvmnc(:,1) = 0


  IF (lasym) then

     radius3: DO js = 2, ns

        gmn = 0

        bmn = 0

        bsubumn = 0

        bsubvmn = 0

        bsubsmn = 0

        bsupumn = 0

        bsupvmn = 0

        mn3: DO mn = 1, mnmax_nyq

           n = nint(xn_nyq(mn))/nfp

           m = nint(xm_nyq(mn))

           n1 = abs(n)


           dmult = mscale(m)*nscale(n1)*tmult

           IF (m.eq.0 .or. n.eq.0) dmult = 2*dmult


           sgn = sign(1, n)


           lk = 0

           jlk = js

           DO j = 1, ntheta2

              DO k = 1, nzeta

                 lk = lk + 1


                 tcosi = dmult*(cosmui(j,m)*cosnv(k,n1) + sgn*sinmui(j,m)*sinnv(k,n1))

                 tsini = dmult*(sinmui(j,m)*cosnv(k,n1) - sgn*cosmui(j,m)*sinnv(k,n1))


                 bmn(mn)     = bmn(mn)     + tsini*bsqa(jlk)

                 gmn(mn)     = gmn(mn)     + tsini*gsqrta(jlk,0) ! TODO: what is this zero index?


                 bsubumn(mn) = bsubumn(mn) + tsini*bsubua(jlk)

                 bsubvmn(mn) = bsubvmn(mn) + tsini*bsubva(jlk)

                 bsubsmn(mn) = bsubsmn(mn) + tcosi*bsubsa(jlk)


                 bsupumn(mn) = bsupumn(mn) + tsini*bsupua(jlk)

                 bsupvmn(mn) = bsupvmn(mn) + tsini*bsupva(jlk)

                 jlk = jlk+ns

              END DO

           END DO

        END DO mn3


        gmns(:,js)     = gmn(:)

        bmns(:,js)     = bmn(:)

        bsubumns(:,js) = bsubumn(:)

        bsubvmns(:,js) = bsubvmn(:)

        bsubsmnc(:,js) = bsubsmn(:)

        bsupumns(:,js) = bsupumn(:)

        bsupvmns(:,js) = bsupvmn(:)

     END DO radius3


     gmns(:,1) = 0

     bmns(:,1) = 0

     bsubumns(:,1) = 0

     bsubvmns(:,1) = 0

     bsubsmnc(:,1) = 2*bsubsmnc(:,2) - bsubsmnc(:,3)

     bsupumns(:,1) = 0

     bsupvmns(:,1) = 0

  end if


  ! WRITE OUT ARRAYS

  CALL cdf_write(nwout, vn_racc, raxis_cc(0:ntor))

  CALL cdf_write(nwout, vn_zacs, zaxis_cs(0:ntor))

  CALL cdf_write(nwout, vn_rmnc, rmnc)

  CALL cdf_write(nwout, vn_zmns, zmns)

  CALL cdf_write(nwout, vn_lmns, lmns)

  CALL cdf_write(nwout, vn_gmnc, gmnc)              !Half mesh

  CALL cdf_write(nwout, vn_bmnc, bmnc)              !Half mesh

  CALL cdf_write(nwout, vn_bsubumnc, bsubumnc)      !Half mesh

  CALL cdf_write(nwout, vn_bsubvmnc, bsubvmnc)      !Half mesh

  CALL cdf_write(nwout, vn_bsubsmns, bsubsmns)      !Full mesh


  ! GET RID OF THESE EVENTUALLY: DON'T NEED THEM (can express in terms of lambdas)

  CALL cdf_write(nwout, vn_bsupumnc, bsupumnc)

  CALL cdf_write(nwout, vn_bsupvmnc, bsupvmnc)


  ! FULL-MESH quantities

  ! NOTE: jdotb is in units_of_A (1/mu0 incorporated in jxbforce...)

  ! prior to version 6.00, this was output in internal VMEC units...

  j = SIZE(am)-1

  CALL cdf_write(nwout, vn_am, am(0:j))

  j = SIZE(ac)-1

  CALL cdf_write(nwout, vn_ac, ac(0:j))

  j = SIZE(ai)-1

  CALL cdf_write(nwout, vn_ai, ai(0:j))


  j = SIZE(am_aux_s)

  CALL cdf_write(nwout, vn_am_aux_s, am_aux_s(1:j))

  j = SIZE(am_aux_f)

  CALL cdf_write(nwout, vn_am_aux_f, am_aux_f(1:j))

  j = SIZE(ac_aux_s)

  CALL cdf_write(nwout, vn_ac_aux_s, ac_aux_s(1:j))

  j = SIZE(ac_aux_f)

  CALL cdf_write(nwout, vn_ac_aux_f, ac_aux_f(1:j))

  j = SIZE(ai_aux_s)

  CALL cdf_write(nwout, vn_ai_aux_s, ai_aux_s(1:j))

  j = SIZE(ai_aux_f)

  CALL cdf_write(nwout, vn_ai_aux_f, ai_aux_f(1:j))


  CALL cdf_write(nwout, vn_iotaf, iotaf(1:ns))

  CALL cdf_write(nwout, vn_qfact, qfact(1:ns))

  CALL cdf_write(nwout, vn_presf, presf/mu0) ! NOTE: scaling !!!

  CALL cdf_write(nwout, vn_phi, phi)

  CALL cdf_write(nwout, vn_phipf, twopi*signgs*phipf) ! NOTE: scaling !!!

  CALL cdf_write(nwout, vn_chi, chi)

  CALL cdf_write(nwout, vn_chipf, twopi*signgs*chipf) ! NOTE: scaling !!!

  CALL cdf_write(nwout, vn_jcuru, jcuru/mu0) ! NOTE: scaling !!!

  CALL cdf_write(nwout, vn_jcurv, jcurv/mu0) ! NOTE: scaling !!!

  CALL cdf_write(nwout, vn_jdotb, jdotb)

  CALL cdf_write(nwout, vn_bgrv, bdotgradv)


  ! HALF-MESH quantities

  iotas(1) = 0

  mass(1) = 0

  pres(1) = 0

  phip(1) = 0

  buco(1) = 0

  bvco(1) = 0

  vp(1) = 0

  overr(1) = 0

  specw(1) = 1

  beta_vol(1) = 0


  CALL cdf_write(nwout, vn_iotah, iotas(1:ns))

  CALL cdf_write(nwout, vn_mass, mass/mu0) ! NOTE: scaling !!!

  CALL cdf_write(nwout, vn_presh, pres(1:ns)/mu0) ! NOTE: scaling !!!

  CALL cdf_write(nwout, vn_betah, beta_vol)

  CALL cdf_write(nwout, vn_buco, buco)

  CALL cdf_write(nwout, vn_bvco, bvco)

  CALL cdf_write(nwout, vn_vp, vp(1:ns))

  CALL cdf_write(nwout, vn_specw, specw)

  CALL cdf_write(nwout, vn_phip, phips(1:ns))

  CALL cdf_write(nwout, vn_overr, overr(1:ns))


  ! MERCIER_CRITERION

  CALL cdf_write(nwout, vn_merc, dmerc(1:ns))

  CALL cdf_write(nwout, vn_mshear, dshear(1:ns))

  CALL cdf_write(nwout, vn_mwell, dwell(1:ns))

  CALL cdf_write(nwout, vn_mcurr, dcurr(1:ns))

  CALL cdf_write(nwout, vn_mgeo, dgeod(1:ns))


  CALL cdf_write(nwout, vn_equif, equif(1:ns))


  IF (lasym) THEN

     CALL cdf_write(nwout, vn_racs, raxis_cs(0:ntor))

     CALL cdf_write(nwout, vn_zacc, zaxis_cc(0:ntor))


     CALL cdf_write(nwout, vn_rmns, rmns)

     CALL cdf_write(nwout, vn_zmnc, zmnc)

     CALL cdf_write(nwout, vn_lmnc, lmnc)


     CALL cdf_write(nwout, vn_gmns, gmns)

     CALL cdf_write(nwout, vn_bmns, bmns)


     CALL cdf_write(nwout, vn_bsubumns, bsubumns)

     CALL cdf_write(nwout, vn_bsubvmns, bsubvmns)

     CALL cdf_write(nwout, vn_bsubsmnc, bsubsmnc)


     ! GET RID OF THESE EVENTUALLY: DON'T NEED THEM

     CALL cdf_write(nwout, vn_bsupumns, bsupumns)

     CALL cdf_write(nwout, vn_bsupvmns, bsupvmns)

  END IF


  CALL cdf_close(nwout)


  ! RESTORE nyq ENDPOINT VALUES

  IF (mnyq .ne. 0) cosmui(:,mnyq) = 2*cosmui(:,mnyq)

  IF (nnyq .ne. 0) cosnv(:,nnyq) = 2*cosnv(:,nnyq)


  ! DEALLOCATIONS

  IF (ALLOCATED(gmnc)) DEALLOCATE(gmnc, bmnc, bsubumnc, bsubvmnc, &

                                  bsubsmns, bsupumnc, bsupvmnc )

  IF (ALLOCATED(gmns)) DEALLOCATE(gmns, bmns, bsubumns, bsubvmns, &

                                  bsubsmnc, bsupumns, bsupvmns )

  IF (ALLOCATED(gmn))  DEALLOCATE (gmn, bmn, bsubumn, bsubvmn, &

                                   bsubsmn, bsupumn, bsupvmn, stat=istat)


  ! FREE BOUNDARY DATA

  ! rmnc1, ... contain the LCFS values at this point

  CALL freeb_data(rmnc1, zmns1, rmns1, zmnc1, bmodmn, bmodmn1)


  rzl_array = 0


END SUBROUTINE wrout

convert
subroutine convert(rmnc, zmns, lmns, rmns, zmnc, lmnc, rzl_array, js)
Convert internal mode representation to standard form for output (coefficients of cos(mu-nv),...
Definition convert.f90:17

freeb_data
subroutine freeb_data(rmnc, zmns, rmns, zmnc, bmodmn, bmodmn1)
Write out edge values of fields.
Definition freeb_data.f90:13

mgrid_mod
Definition mgrid_mod.f:1

mgrid_mod::vn_mgmode
character(len= *), parameter vn_mgmode
Definition mgrid_mod.f:16

mgrid_mod::curlabel
character(len=30), dimension(:), allocatable curlabel
Definition mgrid_mod.f:95

mgrid_mod::mgrid_mode
character(len=1) mgrid_mode
Definition mgrid_mod.f:102

mgrid_mod::nextcur
integer nextcur
Definition mgrid_mod.f:78

mgrid_mod::vn_nextcur
character(len= *), parameter vn_nextcur
Definition mgrid_mod.f:16

read_wout_mod
Reading of wout VMEC output file.
Definition read_wout_mod.f:5

read_wout_mod::vn_bmns
character(len= *), parameter vn_bmns
Definition read_wout_mod.f:32

read_wout_mod::vn_piota_type
character(len= *), parameter vn_piota_type
Definition read_wout_mod.f:32

read_wout_mod::vn_aspect
character(len= *), parameter vn_aspect
Definition read_wout_mod.f:32

read_wout_mod::vn_bsupumnc
character(len= *), parameter vn_bsupumnc
Definition read_wout_mod.f:32

read_wout_mod::ln_ctor
character(len= *), parameter ln_ctor
Definition read_wout_mod.f:94

read_wout_mod::ln_tmod
character(len= *), parameter ln_tmod
Definition read_wout_mod.f:94

read_wout_mod::ln_vp
character(len= *), parameter ln_vp
Definition read_wout_mod.f:94

read_wout_mod::vn_rmns
character(len= *), parameter vn_rmns
Definition read_wout_mod.f:32

read_wout_mod::ln_beta
character(len= *), parameter ln_beta
Definition read_wout_mod.f:94

read_wout_mod::vn_radnod
character(len= *), parameter vn_radnod
Definition read_wout_mod.f:32

read_wout_mod::ln_curlab
character(len= *), parameter ln_curlab
Definition read_wout_mod.f:94

read_wout_mod::ln_pbeta
character(len= *), parameter ln_pbeta
Definition read_wout_mod.f:94

read_wout_mod::vn_zacs
character(len= *), parameter vn_zacs
Definition read_wout_mod.f:32

read_wout_mod::ln_am_aux_f
character(len= *), parameter ln_am_aux_f
Definition read_wout_mod.f:94

read_wout_mod::vn_rbt1
character(len= *), parameter vn_rbt1
Definition read_wout_mod.f:32

read_wout_mod::ln_maxmod_nyq
character(len= *), parameter ln_maxmod_nyq
Definition read_wout_mod.f:94

read_wout_mod::vn_bsupumns
character(len= *), parameter vn_bsupumns
Definition read_wout_mod.f:32

read_wout_mod::vn_bsubsmns
character(len= *), parameter vn_bsubsmns
Definition read_wout_mod.f:32

read_wout_mod::ln_extcur
character(len= *), parameter ln_extcur
Definition read_wout_mod.f:94

read_wout_mod::vn_pmass_type
character(len= *), parameter vn_pmass_type
Definition read_wout_mod.f:32

read_wout_mod::ln_zbs
character(len= *), parameter ln_zbs
Definition read_wout_mod.f:94

read_wout_mod::ln_fsq
character(len= *), parameter ln_fsq
Definition read_wout_mod.f:94

read_wout_mod::ln_amin
character(len= *), parameter ln_amin
Definition read_wout_mod.f:94

read_wout_mod::ln_bsubvmns
character(len= *), parameter ln_bsubvmns
Definition read_wout_mod.f:94

read_wout_mod::ln_chipf
character(len= *), parameter ln_chipf
Definition read_wout_mod.f:94

read_wout_mod::ln_maxit
character(len= *), parameter ln_maxit
Definition read_wout_mod.f:94

read_wout_mod::vn_lmnc
character(len= *), parameter vn_lmnc
Definition read_wout_mod.f:32

read_wout_mod::ln_bsubumns
character(len= *), parameter ln_bsubumns
Definition read_wout_mod.f:94

read_wout_mod::ln_jcurv
character(len= *), parameter ln_jcurv
Definition read_wout_mod.f:94

read_wout_mod::vn_tbeta
character(len= *), parameter vn_tbeta
Definition read_wout_mod.f:32

read_wout_mod::ln_bsubvmnc
character(len= *), parameter ln_bsubvmnc
Definition read_wout_mod.f:94

read_wout_mod::ln_presh
character(len= *), parameter ln_presh
Definition read_wout_mod.f:94

read_wout_mod::ln_merc
character(len= *), parameter ln_merc
Definition read_wout_mod.f:94

read_wout_mod::ln_zacc
character(len= *), parameter ln_zacc
Definition read_wout_mod.f:94

read_wout_mod::ln_thom
character(len= *), parameter ln_thom
Definition read_wout_mod.f:94

read_wout_mod::ln_rmns
character(len= *), parameter ln_rmns
Definition read_wout_mod.f:94

read_wout_mod::vn_overr
character(len= *), parameter vn_overr
Definition read_wout_mod.f:32

read_wout_mod::vn_ftolv
character(len= *), parameter vn_ftolv
Definition read_wout_mod.f:32

read_wout_mod::ln_therm
character(len= *), parameter ln_therm
Definition read_wout_mod.f:94

read_wout_mod::vn_curlab
character(len= *), parameter vn_curlab
Definition read_wout_mod.f:32

read_wout_mod::vn_specw
character(len= *), parameter vn_specw
Definition read_wout_mod.f:32

read_wout_mod::ln_lmns
character(len= *), parameter ln_lmns
Definition read_wout_mod.f:94

read_wout_mod::vn_vp
character(len= *), parameter vn_vp
Definition read_wout_mod.f:32

read_wout_mod::vn_vol
character(len= *), parameter vn_vol
Definition read_wout_mod.f:32

read_wout_mod::vn_magen
character(len= *), parameter vn_magen
Definition read_wout_mod.f:32

read_wout_mod::vn_maxz
character(len= *), parameter vn_maxz
Definition read_wout_mod.f:32

read_wout_mod::ln_racc
character(len= *), parameter ln_racc
Definition read_wout_mod.f:94

read_wout_mod::vn_abeta
character(len= *), parameter vn_abeta
Definition read_wout_mod.f:32

read_wout_mod::vn_ac
character(len= *), parameter vn_ac
Definition read_wout_mod.f:32

read_wout_mod::vn_presf
character(len= *), parameter vn_presf
Definition read_wout_mod.f:32

read_wout_mod::vn_chipf
character(len= *), parameter vn_chipf
Definition read_wout_mod.f:32

read_wout_mod::vn_jcuru
character(len= *), parameter vn_jcuru
Definition read_wout_mod.f:32

read_wout_mod::ln_ai_aux_s
character(len= *), parameter ln_ai_aux_s
Definition read_wout_mod.f:94

read_wout_mod::vn_free
character(len= *), parameter vn_free
Definition read_wout_mod.f:32

read_wout_mod::vn_version
character(len= *), parameter vn_version
Definition read_wout_mod.f:32

read_wout_mod::vn_buco
character(len= *), parameter vn_buco
Definition read_wout_mod.f:32

read_wout_mod::vn_phipf
character(len= *), parameter vn_phipf
Definition read_wout_mod.f:32

read_wout_mod::ln_bsubsmns
character(len= *), parameter ln_bsubsmns
Definition read_wout_mod.f:94

read_wout_mod::ln_pmod
character(len= *), parameter ln_pmod
Definition read_wout_mod.f:94

read_wout_mod::vn_ac_aux_f
character(len= *), parameter vn_ac_aux_f
Definition read_wout_mod.f:32

read_wout_mod::ln_radnod
character(len= *), parameter ln_radnod
Definition read_wout_mod.f:94

read_wout_mod::vn_bgrv
character(len= *), parameter vn_bgrv
Definition read_wout_mod.f:32

read_wout_mod::vn_mgrid
character(len= *), parameter vn_mgrid
Definition read_wout_mod.f:32

read_wout_mod::ln_extension
character(len= *), parameter ln_extension
Definition read_wout_mod.f:94

read_wout_mod::vn_zmnc
character(len= *), parameter vn_zmnc
Definition read_wout_mod.f:32

read_wout_mod::ln_chi
character(len= *), parameter ln_chi
Definition read_wout_mod.f:94

read_wout_mod::vn_iotah
character(len= *), parameter vn_iotah
Definition read_wout_mod.f:32

read_wout_mod::vn_lmns
character(len= *), parameter vn_lmns
Definition read_wout_mod.f:32

read_wout_mod::ln_rbs
character(len= *), parameter ln_rbs
Definition read_wout_mod.f:94

read_wout_mod::ln_ai
character(len= *), parameter ln_ai
Definition read_wout_mod.f:94

read_wout_mod::vn_bvco
character(len= *), parameter vn_bvco
Definition read_wout_mod.f:32

read_wout_mod::vn_fsqr
character(len= *), parameter vn_fsqr
Definition read_wout_mod.f:32

read_wout_mod::vn_zacc
character(len= *), parameter vn_zacc
Definition read_wout_mod.f:32

read_wout_mod::ln_mse
character(len= *), parameter ln_mse
Definition read_wout_mod.f:94

read_wout_mod::vn_gam
character(len= *), parameter vn_gam
Definition read_wout_mod.f:32

read_wout_mod::ln_ac_aux_f
character(len= *), parameter ln_ac_aux_f
Definition read_wout_mod.f:94

read_wout_mod::vn_betah
character(len= *), parameter vn_betah
Definition read_wout_mod.f:32

read_wout_mod::ln_maxr
character(len= *), parameter ln_maxr
Definition read_wout_mod.f:94

read_wout_mod::vn_bsupvmnc
character(len= *), parameter vn_bsupvmnc
Definition read_wout_mod.f:32

read_wout_mod::ln_mgeo
character(len= *), parameter ln_mgeo
Definition read_wout_mod.f:94

read_wout_mod::vn_mshear
character(len= *), parameter vn_mshear
Definition read_wout_mod.f:32

read_wout_mod::ln_bgrv
character(len= *), parameter ln_bgrv
Definition read_wout_mod.f:94

read_wout_mod::vn_merc
character(len= *), parameter vn_merc
Definition read_wout_mod.f:32

read_wout_mod::vn_b0
character(len= *), parameter vn_b0
Definition read_wout_mod.f:32

read_wout_mod::ln_zmnc
character(len= *), parameter ln_zmnc
Definition read_wout_mod.f:94

read_wout_mod::mnmax_nyq
integer mnmax_nyq
Definition read_wout_mod.f:201

read_wout_mod::ln_ac_aux_s
character(len= *), parameter ln_ac_aux_s
Definition read_wout_mod.f:94

read_wout_mod::vn_ai_aux_s
character(len= *), parameter vn_ai_aux_s
Definition read_wout_mod.f:32

read_wout_mod::ln_free
character(len= *), parameter ln_free
Definition read_wout_mod.f:94

read_wout_mod::ln_zbc
character(len= *), parameter ln_zbc
Definition read_wout_mod.f:94

read_wout_mod::vn_rbc
character(len= *), parameter vn_rbc
Definition read_wout_mod.f:32

read_wout_mod::vn_zbs
character(len= *), parameter vn_zbs
Definition read_wout_mod.f:32

read_wout_mod::ln_phip
character(len= *), parameter ln_phip
Definition read_wout_mod.f:94

read_wout_mod::ln_bsupvmns
character(len= *), parameter ln_bsupvmns
Definition read_wout_mod.f:94

read_wout_mod::ln_buco
character(len= *), parameter ln_buco
Definition read_wout_mod.f:94

read_wout_mod::vn_gmnc
character(len= *), parameter vn_gmnc
Definition read_wout_mod.f:32

read_wout_mod::ln_potvac
character(len= *), parameter ln_potvac
Definition read_wout_mod.f:94

read_wout_mod::vn_qfact
character(len= *), parameter vn_qfact
Definition read_wout_mod.f:32

read_wout_mod::ln_bmns
character(len= *), parameter ln_bmns
Definition read_wout_mod.f:94

read_wout_mod::vn_maxr
character(len= *), parameter vn_maxr
Definition read_wout_mod.f:32

read_wout_mod::vn_bsubsmnc
character(len= *), parameter vn_bsubsmnc
Definition read_wout_mod.f:32

read_wout_mod::vn_jcurv
character(len= *), parameter vn_jcurv
Definition read_wout_mod.f:32

read_wout_mod::vn_zmns
character(len= *), parameter vn_zmns
Definition read_wout_mod.f:32

read_wout_mod::vn_jdotb
character(len= *), parameter vn_jdotb
Definition read_wout_mod.f:32

read_wout_mod::ln_bsupvmnc
character(len= *), parameter ln_bsupvmnc
Definition read_wout_mod.f:94

read_wout_mod::ln_equif
character(len= *), parameter ln_equif
Definition read_wout_mod.f:94

read_wout_mod::ln_am
character(len= *), parameter ln_am
Definition read_wout_mod.f:94

read_wout_mod::vn_am_aux_s
character(len= *), parameter vn_am_aux_s
Definition read_wout_mod.f:32

read_wout_mod::vn_tmod_nyq
character(len= *), parameter vn_tmod_nyq
Definition read_wout_mod.f:32

read_wout_mod::ln_gmns
character(len= *), parameter ln_gmns
Definition read_wout_mod.f:94

read_wout_mod::ln_version
character(len= *), parameter ln_version
Definition read_wout_mod.f:94

read_wout_mod::vn_mwell
character(len= *), parameter vn_mwell
Definition read_wout_mod.f:32

read_wout_mod::ln_jcuru
character(len= *), parameter ln_jcuru
Definition read_wout_mod.f:94

read_wout_mod::ln_bvco
character(len= *), parameter ln_bvco
Definition read_wout_mod.f:94

read_wout_mod::vn_fsql
character(len= *), parameter vn_fsql
Definition read_wout_mod.f:32

read_wout_mod::ln_ac
character(len= *), parameter ln_ac
Definition read_wout_mod.f:94

read_wout_mod::vn_amin
character(len= *), parameter vn_amin
Definition read_wout_mod.f:32

read_wout_mod::vn_rmnc
character(len= *), parameter vn_rmnc
Definition read_wout_mod.f:32

read_wout_mod::mnmax
integer mnmax
Definition read_wout_mod.f:201

read_wout_mod::vn_mgeo
character(len= *), parameter vn_mgeo
Definition read_wout_mod.f:32

read_wout_mod::vn_bsubumns
character(len= *), parameter vn_bsubumns
Definition read_wout_mod.f:32

read_wout_mod::vn_bsupvmns
character(len= *), parameter vn_bsupvmns
Definition read_wout_mod.f:32

read_wout_mod::vn_mass
character(len= *), parameter vn_mass
Definition read_wout_mod.f:32

read_wout_mod::vn_modb
character(len= *), parameter vn_modb
Definition read_wout_mod.f:32

read_wout_mod::ln_vol
character(len= *), parameter ln_vol
Definition read_wout_mod.f:94

read_wout_mod::ln_racs
character(len= *), parameter ln_racs
Definition read_wout_mod.f:94

read_wout_mod::vn_gmns
character(len= *), parameter vn_gmns
Definition read_wout_mod.f:32

read_wout_mod::vn_am
character(len= *), parameter vn_am
Definition read_wout_mod.f:32

read_wout_mod::vn_fsq
character(len= *), parameter vn_fsq
Definition read_wout_mod.f:32

read_wout_mod::ln_rbt1
character(len= *), parameter ln_rbt1
Definition read_wout_mod.f:94

read_wout_mod::vn_am_aux_f
character(len= *), parameter vn_am_aux_f
Definition read_wout_mod.f:32

read_wout_mod::vn_lar
character(len= *), parameter vn_lar
Definition read_wout_mod.f:32

read_wout_mod::vn_rmaj
character(len= *), parameter vn_rmaj
Definition read_wout_mod.f:32

read_wout_mod::ln_specw
character(len= *), parameter ln_specw
Definition read_wout_mod.f:94

read_wout_mod::vn_bmnc
character(len= *), parameter vn_bmnc
Definition read_wout_mod.f:32

read_wout_mod::vn_ai
character(len= *), parameter vn_ai
Definition read_wout_mod.f:32

read_wout_mod::ln_phipf
character(len= *), parameter ln_phipf
Definition read_wout_mod.f:94

read_wout_mod::ln_polmod
character(len= *), parameter ln_polmod
Definition read_wout_mod.f:94

read_wout_mod::ln_bmnc
character(len= *), parameter ln_bmnc
Definition read_wout_mod.f:94

read_wout_mod::ln_mcurr
character(len= *), parameter ln_mcurr
Definition read_wout_mod.f:94

read_wout_mod::vn_racc
character(len= *), parameter vn_racc
Definition read_wout_mod.f:32

read_wout_mod::ln_bsupumns
character(len= *), parameter ln_bsupumns
Definition read_wout_mod.f:94

read_wout_mod::vn_pmod
character(len= *), parameter vn_pmod
Definition read_wout_mod.f:32

read_wout_mod::ln_iotah
character(len= *), parameter ln_iotah
Definition read_wout_mod.f:94

read_wout_mod::ln_error
character(len= *), parameter ln_error
Definition read_wout_mod.f:94

read_wout_mod::ln_piota_type
character(len= *), parameter ln_piota_type
Definition read_wout_mod.f:94

read_wout_mod::vn_wdot
character(len= *), parameter vn_wdot
Definition read_wout_mod.f:32

read_wout_mod::vn_asym
character(len= *), parameter vn_asym
Definition read_wout_mod.f:32

read_wout_mod::vn_extension
character(len= *), parameter vn_extension
Definition read_wout_mod.f:32

read_wout_mod::ln_maxz
character(len= *), parameter ln_maxz
Definition read_wout_mod.f:94

read_wout_mod::vn_tmod
character(len= *), parameter vn_tmod
Definition read_wout_mod.f:32

read_wout_mod::ln_mwell
character(len= *), parameter ln_mwell
Definition read_wout_mod.f:94

read_wout_mod::ln_mshear
character(len= *), parameter ln_mshear
Definition read_wout_mod.f:94

read_wout_mod::vn_beta
character(len= *), parameter vn_beta
Definition read_wout_mod.f:32

read_wout_mod::vn_pcurr_type
character(len= *), parameter vn_pcurr_type
Definition read_wout_mod.f:32

read_wout_mod::ln_rmaj
character(len= *), parameter ln_rmaj
Definition read_wout_mod.f:94

read_wout_mod::ln_rmnc
character(len= *), parameter ln_rmnc
Definition read_wout_mod.f:94

read_wout_mod::ln_lar
character(len= *), parameter ln_lar
Definition read_wout_mod.f:94

read_wout_mod::vn_bsubvmns
character(len= *), parameter vn_bsubvmns
Definition read_wout_mod.f:32

read_wout_mod::ln_pmass_type
character(len= *), parameter ln_pmass_type
Definition read_wout_mod.f:94

read_wout_mod::vn_tormod
character(len= *), parameter vn_tormod
Definition read_wout_mod.f:32

read_wout_mod::ln_magen
character(len= *), parameter ln_magen
Definition read_wout_mod.f:94

read_wout_mod::ln_minr
character(len= *), parameter ln_minr
Definition read_wout_mod.f:94

read_wout_mod::ln_presf
character(len= *), parameter ln_presf
Definition read_wout_mod.f:94

read_wout_mod::ln_b0
character(len= *), parameter ln_b0
Definition read_wout_mod.f:94

read_wout_mod::vn_equif
character(len= *), parameter vn_equif
Definition read_wout_mod.f:32

read_wout_mod::ln_rbt0
character(len= *), parameter ln_rbt0
Definition read_wout_mod.f:94

read_wout_mod::ln_fp
character(len= *), parameter ln_fp
Definition read_wout_mod.f:94

read_wout_mod::vn_error
character(len= *), parameter vn_error
Definition read_wout_mod.f:32

read_wout_mod::vn_minr
character(len= *), parameter vn_minr
Definition read_wout_mod.f:32

read_wout_mod::vn_extcur
character(len= *), parameter vn_extcur
Definition read_wout_mod.f:32

read_wout_mod::vn_pmod_nyq
character(len= *), parameter vn_pmod_nyq
Definition read_wout_mod.f:32

read_wout_mod::ln_mass
character(len= *), parameter ln_mass
Definition read_wout_mod.f:94

read_wout_mod::vn_bsubumnc
character(len= *), parameter vn_bsubumnc
Definition read_wout_mod.f:32

read_wout_mod::vn_presh
character(len= *), parameter vn_presh
Definition read_wout_mod.f:32

read_wout_mod::vn_sgs
character(len= *), parameter vn_sgs
Definition read_wout_mod.f:32

read_wout_mod::ln_pmod_nyq
character(len= *), parameter ln_pmod_nyq
Definition read_wout_mod.f:94

read_wout_mod::vn_phi
character(len= *), parameter vn_phi
Definition read_wout_mod.f:32

read_wout_mod::ln_rbc
character(len= *), parameter ln_rbc
Definition read_wout_mod.f:94

read_wout_mod::ln_aspect
character(len= *), parameter ln_aspect
Definition read_wout_mod.f:94

read_wout_mod::vn_ctor
character(len= *), parameter vn_ctor
Definition read_wout_mod.f:32

read_wout_mod::ln_ai_aux_f
character(len= *), parameter ln_ai_aux_f
Definition read_wout_mod.f:94

read_wout_mod::vn_racs
character(len= *), parameter vn_racs
Definition read_wout_mod.f:32

read_wout_mod::vn_pbeta
character(len= *), parameter vn_pbeta
Definition read_wout_mod.f:32

read_wout_mod::vn_chi
character(len= *), parameter vn_chi
Definition read_wout_mod.f:32

read_wout_mod::vn_polmod
character(len= *), parameter vn_polmod
Definition read_wout_mod.f:32

read_wout_mod::vn_mcurr
character(len= *), parameter vn_mcurr
Definition read_wout_mod.f:32

read_wout_mod::ln_zmns
character(len= *), parameter ln_zmns
Definition read_wout_mod.f:94

read_wout_mod::ln_bsubumnc
character(len= *), parameter ln_bsubumnc
Definition read_wout_mod.f:94

read_wout_mod::vn_maxmod_nyq
character(len= *), parameter vn_maxmod_nyq
Definition read_wout_mod.f:32

read_wout_mod::vn_bsubvmnc
character(len= *), parameter vn_bsubvmnc
Definition read_wout_mod.f:32

read_wout_mod::ln_bsupumnc
character(len= *), parameter ln_bsupumnc
Definition read_wout_mod.f:94

read_wout_mod::ln_sgs
character(len= *), parameter ln_sgs
Definition read_wout_mod.f:94

read_wout_mod::ln_wdot
character(len= *), parameter ln_wdot
Definition read_wout_mod.f:94

read_wout_mod::ln_pcurr_type
character(len= *), parameter ln_pcurr_type
Definition read_wout_mod.f:94

read_wout_mod::ln_tbeta
character(len= *), parameter ln_tbeta
Definition read_wout_mod.f:94

read_wout_mod::ln_maxmod
character(len= *), parameter ln_maxmod
Definition read_wout_mod.f:94

read_wout_mod::vn_maxit
character(len= *), parameter vn_maxit
Definition read_wout_mod.f:32

read_wout_mod::vn_ac_aux_s
character(len= *), parameter vn_ac_aux_s
Definition read_wout_mod.f:32

read_wout_mod::vn_rbt0
character(len= *), parameter vn_rbt0
Definition read_wout_mod.f:32

read_wout_mod::ln_lmnc
character(len= *), parameter ln_lmnc
Definition read_wout_mod.f:94

read_wout_mod::ln_tmod_nyq
character(len= *), parameter ln_tmod_nyq
Definition read_wout_mod.f:94

read_wout_mod::vn_fp
character(len= *), parameter vn_fp
Definition read_wout_mod.f:32

read_wout_mod::ln_bsubsmnc
character(len= *), parameter ln_bsubsmnc
Definition read_wout_mod.f:94

read_wout_mod::ln_mgrid
character(len= *), parameter ln_mgrid
Definition read_wout_mod.f:94

read_wout_mod::ln_iotaf
character(len= *), parameter ln_iotaf
Definition read_wout_mod.f:94

read_wout_mod::vn_zbc
character(len= *), parameter vn_zbc
Definition read_wout_mod.f:32

read_wout_mod::vn_potvac
character(len= *), parameter vn_potvac
Definition read_wout_mod.f:32

read_wout_mod::ln_am_aux_s
character(len= *), parameter ln_am_aux_s
Definition read_wout_mod.f:94

read_wout_mod::ln_phi
character(len= *), parameter ln_phi
Definition read_wout_mod.f:94

read_wout_mod::ln_asym
character(len= *), parameter ln_asym
Definition read_wout_mod.f:94

read_wout_mod::ln_tormod
character(len= *), parameter ln_tormod
Definition read_wout_mod.f:94

read_wout_mod::ln_abeta
character(len= *), parameter ln_abeta
Definition read_wout_mod.f:94

read_wout_mod::ln_qfact
character(len= *), parameter ln_qfact
Definition read_wout_mod.f:94

read_wout_mod::vn_maxmod
character(len= *), parameter vn_maxmod
Definition read_wout_mod.f:32

read_wout_mod::vn_iotaf
character(len= *), parameter vn_iotaf
Definition read_wout_mod.f:32

read_wout_mod::ln_betah
character(len= *), parameter ln_betah
Definition read_wout_mod.f:94

read_wout_mod::vn_therm
character(len= *), parameter vn_therm
Definition read_wout_mod.f:32

read_wout_mod::vn_rbs
character(len= *), parameter vn_rbs
Definition read_wout_mod.f:32

read_wout_mod::ln_gmnc
character(len= *), parameter ln_gmnc
Definition read_wout_mod.f:94

read_wout_mod::vn_phip
character(len= *), parameter vn_phip
Definition read_wout_mod.f:32

read_wout_mod::vn_fsqz
character(len= *), parameter vn_fsqz
Definition read_wout_mod.f:32

read_wout_mod::ln_modb
character(len= *), parameter ln_modb
Definition read_wout_mod.f:94

read_wout_mod::ln_gam
character(len= *), parameter ln_gam
Definition read_wout_mod.f:94

read_wout_mod::ln_jdotb
character(len= *), parameter ln_jdotb
Definition read_wout_mod.f:94

read_wout_mod::ln_zacs
character(len= *), parameter ln_zacs
Definition read_wout_mod.f:94

read_wout_mod::vn_ai_aux_f
character(len= *), parameter vn_ai_aux_f
Definition read_wout_mod.f:32

realspace
Definition realspace.f90:2

realspace::z1
real(rprec), dimension(:,:), allocatable, target z1
Definition realspace.f90:11

realspace::phip
real(rprec), dimension(:), allocatable phip
radial derivative of phi/(2*pi) on half-grid
Definition realspace.f90:29

safe_open_mod
fault-tolerant file opening routines
Definition safe_open_mod.f:5

vacmod
Definition vacmod.f90:2

vacmod::potvac
real(rprec), dimension(:), allocatable, target potvac
Definition vacmod.f90:29

vacmod::p5
real(rprec), parameter p5
Definition vacmod.f90:13

vforces
Definition vforces.f90:2

vforces::bzmn_e
real(rprec), dimension(:), pointer bzmn_e
Definition vforces.f90:28

vforces::brmn_e
real(rprec), dimension(:), pointer brmn_e
Definition vforces.f90:21

vforces::azmn_e
real(rprec), dimension(:), pointer azmn_e
Definition vforces.f90:26

vforces::czmn_o
real(rprec), dimension(:), pointer czmn_o
Definition vforces.f90:31

vforces::armn_o
real(rprec), dimension(:), pointer armn_o
Definition vforces.f90:20

vforces::armn_e
real(rprec), dimension(:), pointer armn_e
Definition vforces.f90:19

vmec_input
Definition vmec_input.f90:2

vmec_input::ac
real(rprec), dimension(0:20) ac
array of coefficients in phi-series for the quantity d(Icurv)/ds = toroidal current density * Vprime,...
Definition vmec_input.f90:45

vmec_input::lfreeb
logical lfreeb
Definition vmec_input.f90:63

vmec_input::pcurr_type
character(len=20) pcurr_type
Definition vmec_input.f90:48

vmec_input::ac_aux_f
real(rprec), dimension(ndatafmax) ac_aux_f
Definition vmec_input.f90:56

vmec_input::ai_aux_s
real(rprec), dimension(ndatafmax) ai_aux_s
Definition vmec_input.f90:53

vmec_input::pmass_type
character(len=20) pmass_type
Definition vmec_input.f90:50

vmec_input::gamma
real(rprec) gamma
Definition vmec_input.f90:38

vmec_input::am_aux_s
real(rprec), dimension(ndatafmax) am_aux_s
Definition vmec_input.f90:51

vmec_input::piota_type
character(len=20) piota_type
Definition vmec_input.f90:49

vmec_input::ai_aux_f
real(rprec), dimension(ndatafmax) ai_aux_f
Definition vmec_input.f90:54

vmec_input::am
real(rprec), dimension(0:20) am
array of coefficients in phi-series for mass (NWT/m**2)
Definition vmec_input.f90:43

vmec_input::am_aux_f
real(rprec), dimension(ndatafmax) am_aux_f
Definition vmec_input.f90:52

vmec_input::raxis_cs
real(rprec), dimension(0:ntord) raxis_cs
Definition vmec_input.f90:59

vmec_input::mgrid_file
character(len=200) mgrid_file
Definition vmec_input.f90:80

vmec_input::ac_aux_s
real(rprec), dimension(ndatafmax) ac_aux_s
Definition vmec_input.f90:55

vmec_input::input_extension
character(len=100) input_extension
Definition vmec_input.f90:81

vmec_input::raxis_cc
real(rprec), dimension(0:ntord) raxis_cc
Definition vmec_input.f90:58

vmec_input::lasym
logical lasym
Definition vmec_input.f90:64

vmec_input::extcur
real(rprec), dimension(nigroup) extcur
Definition vmec_input.f90:62

vmec_input::ai
real(rprec), dimension(0:20) ai
array of coefficients in phi-series for iota (ncurr=0)
Definition vmec_input.f90:44

vmec_input::nfp
integer nfp
Definition vmec_input.f90:16

vmec_input::ntor
integer ntor
Definition vmec_input.f90:21

vmec_input::nzeta
integer nzeta
Definition vmec_input.f90:23

vmec_input::mpol
integer mpol
Definition vmec_input.f90:20

vmec_input::zaxis_cs
real(rprec), dimension(0:ntord) zaxis_cs
Definition vmec_input.f90:61

vmec_input::ftol_array
real(rprec), dimension(100) ftol_array
Definition vmec_input.f90:29

vmec_input::zaxis_cc
real(rprec), dimension(0:ntord) zaxis_cc
Definition vmec_input.f90:60

vmec_input::ns_array
integer, dimension(100) ns_array
Definition vmec_input.f90:27

vmec_io
Definition vmec_io.f90:2

vmec_io::ionlarmor
real(rprec) ionlarmor
Definition vmec_io.f90:9

vmec_io::rmax_surf
real(rprec) rmax_surf
Definition vmec_io.f90:22

vmec_io::aminor_p
real(rprec) aminor_p
Definition vmec_io.f90:10

vmec_io::betapol
real(rprec) betapol
Definition vmec_io.f90:13

vmec_io::zmax_surf
real(rprec) zmax_surf
Definition vmec_io.f90:24

vmec_io::betator
real(rprec) betator
Definition vmec_io.f90:14

vmec_io::betaxis
real(rprec) betaxis
Definition vmec_io.f90:15

vmec_io::rmin_surf
real(rprec) rmin_surf
Definition vmec_io.f90:23

vmec_io::betatot
real(rprec) betatot
Definition vmec_io.f90:12

vmec_io::b0
real(rprec) b0
Definition vmec_io.f90:16

vmec_io::volavgb
real(rprec) volavgb
Definition vmec_io.f90:8

vmec_io::volume_p
real(rprec) volume_p
Definition vmec_io.f90:17

vmec_io::rmajor_p
real(rprec) rmajor_p
Definition vmec_io.f90:11

vmec_main
Definition vmec_main.f90:2

vmec_main::equif
real(rprec), dimension(:), allocatable equif
radial force balance error: grad(p) - <j x B>
Definition vmec_main.f90:45

vmec_main::irzloff
integer irzloff
offset in xc array between R,Z,L components
Definition vmec_main.f90:157

vmec_main::rbtor
real(rprec) rbtor
poloidal current at LCFS
Definition vmec_main.f90:135

vmec_main::vp
real(rprec), dimension(:), allocatable vp
radial derivative of enclosed volume
Definition vmec_main.f90:56

vmec_main::bvco
real(rprec), dimension(:), allocatable bvco
enclosed poloidal current profile
Definition vmec_main.f90:43

vmec_main::buco
real(rprec), dimension(:), allocatable buco
enclosed toroidal current profile
Definition vmec_main.f90:42

vmec_main::sp
real(rprec), dimension(:), allocatable sp
shalf(i+1)/sfull(i)
Definition vmec_main.f90:26

vmec_main::fsql
real(rprec) fsql
Definition vmec_main.f90:96

vmec_main::jcuru
real(rprec), dimension(:), allocatable jcuru
poloidal current density
Definition vmec_main.f90:39

vmec_main::wp
real(rprec) wp
kinetic/thermal energy (from pressure)
Definition vmec_main.f90:105

vmec_main::fsqr
real(rprec) fsqr
Definition vmec_main.f90:94

vmec_main::fsqz
real(rprec) fsqz
Definition vmec_main.f90:95

vmec_main::ctor
real(rprec) ctor
toroidal current (?)
Definition vmec_main.f90:137

vmec_main::presf
real(rprec), dimension(:), allocatable presf
pressure profile on full-grid, mass/phip**gamma
Definition vmec_main.f90:66

vmec_main::phi
real(rprec), dimension(:), allocatable phi
toroidal magnetic flux
Definition vmec_main.f90:37

vmec_main::neqs
integer neqs
total number of equations to evolve (size of xc)
Definition vmec_main.f90:156

vmec_main::overr
real(rprec), dimension(:), allocatable overr
Definition vmec_main.f90:54

vmec_main::sm
real(rprec), dimension(:), allocatable sm
shalf(i)/sfull(i)
Definition vmec_main.f90:25

vmec_main::jdotb
real(rprec), dimension(:), allocatable jdotb
Definition vmec_main.f90:41

vmec_main::aspect
real(rprec) aspect
Definition vmec_main.f90:86

vmec_main::bdotgradv
real(rprec), dimension(:), allocatable bdotgradv
Definition vmec_main.f90:44

vmec_main::specw
real(rprec), dimension(:), allocatable specw
spectral width (diagnostic)
Definition vmec_main.f90:46

vmec_main::iotaf
real(rprec), dimension(:), allocatable iotaf
rotational transform (full grid)
Definition vmec_main.f90:34

vmec_main::rbtor0
real(rprec) rbtor0
poloidal current at magnetic axis
Definition vmec_main.f90:136

vmec_main::pres
real(rprec), dimension(:), allocatable pres
pressure profile
Definition vmec_main.f90:55

vmec_main::chipf
real(rprec), dimension(:), allocatable chipf
radial derivative of poloidal magnetic flux (full grid)
Definition vmec_main.f90:36

vmec_main::wb
real(rprec) wb
magnetic energy: volume integral over B^2/2
Definition vmec_main.f90:104

vmec_main::lthreed
logical lthreed
Definition vmec_main.f90:142

vmec_main::chi
real(rprec), dimension(:), allocatable chi
poloidal magnetic flux
Definition vmec_main.f90:64

vmec_main::r0scale
real(rprec) r0scale
Definition vmec_main.f90:90

vmec_main::phips
real(rprec), dimension(:), allocatable phips
toroidal flux (same as phip), one-dimensional array
Definition vmec_main.f90:68

vmec_main::mass
real(rprec), dimension(:), allocatable mass
mass profile on half-grid
Definition vmec_main.f90:71

vmec_main::beta_vol
real(rprec), dimension(:), allocatable beta_vol
Definition vmec_main.f90:38

vmec_main::iter2
integer iter2
total number of iterations
Definition vmec_main.f90:161

vmec_main::iotas
real(rprec), dimension(:), allocatable iotas
rotational transform , on half radial mesh
Definition vmec_main.f90:69

vmec_main::phipf
real(rprec), dimension(:), allocatable phipf
radial derivative of toroidal magnetic flux (full grid)
Definition vmec_main.f90:35

vmec_main::jcurv
real(rprec), dimension(:), allocatable jcurv
toroidal current density
Definition vmec_main.f90:40

vmec_params
Definition vmec_params.f90:2

vmec_params::zcc
integer zcc
Definition vmec_params.f90:51

vmec_params::signgs
real(rprec) signgs
sign of Jacobian : must be =1 (right-handed) or =-1 (left-handed)
Definition vmec_params.f90:62

vmec_params::mscale
real(rprec), dimension(:), allocatable mscale
array for norming theta-trig functions (internal use only) so that the discrete SUM[cos(mu)*cos(m'u)]...
Definition vmec_params.f90:58

vmec_params::mnyq
integer mnyq
Definition vmec_params.f90:54

vmec_params::nnyq
integer nnyq
Definition vmec_params.f90:55

vmec_params::rss
integer rss
Definition vmec_params.f90:46

vmec_params::zcs
integer zcs
Definition vmec_params.f90:50

vmec_params::lamscale
real(rprec) lamscale
Definition vmec_params.f90:64

vmec_params::nscale
real(rprec), dimension(:), allocatable nscale
array for norming zeta -trig functions (internal use only)
Definition vmec_params.f90:60

vmec_params::version_
character(len= *), parameter version_
Definition vmec_params.f90:41

vmec_params::ntmax
integer ntmax
number of contributing Fourier basis function (can be 1, 2 or 4); assigned in read_indata()
Definition vmec_params.f90:43

vmec_params::rsc
integer rsc
Definition vmec_params.f90:47

vmec_persistent
Definition vmec_persistent.f90:2

vmec_persistent::sinmui
real(rprec), dimension(:,:), allocatable sinmui
Definition vmec_persistent.f90:23

vmec_persistent::xm_nyq
real(rprec), dimension(:), allocatable, target xm_nyq
Definition vmec_persistent.f90:29

vmec_persistent::xm
real(rprec), dimension(:), allocatable, target xm
Definition vmec_persistent.f90:27

vmec_persistent::cosmui
real(rprec), dimension(:,:), allocatable cosmui
Definition vmec_persistent.f90:22

vmec_persistent::sinnv
real(rprec), dimension(:,:), allocatable sinnv
Definition vmec_persistent.f90:18

vmec_persistent::xn
real(rprec), dimension(:), allocatable, target xn
Definition vmec_persistent.f90:28

vmec_persistent::xn_nyq
real(rprec), dimension(:), allocatable, target xn_nyq
Definition vmec_persistent.f90:30

vmec_persistent::cosnv
real(rprec), dimension(:,:), allocatable cosnv
Definition vmec_persistent.f90:17

vmercier
Definition vmercier.f90:2

vmercier::dgeod
real(rprec), dimension(nsd) dgeod
Definition vmercier.f90:12

vmercier::dmerc
real(rprec), dimension(nsd) dmerc
Definition vmercier.f90:11

vmercier::dwell
real(rprec), dimension(nsd) dwell
Definition vmercier.f90:9

vmercier::dcurr
real(rprec), dimension(nsd) dcurr
Definition vmercier.f90:10

vmercier::dshear
real(rprec), dimension(nsd) dshear
Definition vmercier.f90:8

vparams
Definition vparams.f90:2

vparams::cp5
real(rprec), parameter cp5
Definition vparams.f90:37

vparams::c2p0
real(rprec), parameter c2p0
Definition vparams.f90:40

xstuff
Definition xstuff.f90:2

xstuff::xc
real(rprec), dimension(:), allocatable, target xc
stacked array of scaled R, Z, Lambda Fourier coefficients (see above for stack order)
Definition xstuff.f90:40

symoutput
subroutine symoutput(bsq, gsqrt, bsubu, bsubv, bsupu, bsupv, bsubs, bsqa, gsqrta, bsubua, bsubva, bsupua, bsupva, bsubsa)
Symmetrize some quantities so that they can be output (?)
Definition symoutput.f90:22

convert_sym
subroutine convert_sym(rmnss, zmncs)
Convert from internal representation to "physical" rmnss, zmncs Fourier form.
Definition totzsp.f90:217

convert_asym
subroutine convert_asym(rmnsc, zmncc)
Convert from internal representation to "physical" rmnsc, zmncc Fourier form.
Definition totzsp.f90:380

wrout
subroutine wrout(bsq, gsqrt, bsubu, bsubv, bsubs, bsupv, bsupu, rzl_array, gc_array, ier_flag)
Write the output files of VMEC.
Definition wrout.f90:17