VMEC 8.52
3D Equilibrium solver with nested flux surfaces.
Loading...
Searching...
No Matches
wrout.f90
Go to the documentation of this file.
1!> \file
2!> \brief Write the output files of VMEC
3
4!> \brief Write the output files of VMEC
5!>
6!> @param bsq
7!> @param gsqrt
8!> @param bsubu
9!> @param bsubv
10!> @param bsubs
11!> @param bsupv
12!> @param bsupu
13!> @param rzl_array
14!> @param gc_array
15!> @param ier_flag error flag
16SUBROUTINE wrout(bsq, gsqrt, bsubu, bsubv, bsubs, bsupv, bsupu, rzl_array, gc_array, ier_flag)
17 USE vmec_main
18 USE vparams, p5 => cp5, two => c2p0
19 USE vmec_input, ONLY: ns_array, ftol_array
20 USE vmec_params
21 USE vmercier
22 USE vmec_persistent
23 USE xstuff
24 USE vmec_io
25 USE realspace, ONLY: phip, gsqrta=>z1, z1=>z1
26 USE vforces, ONLY: bsupua=>brmn_e, bsupva=>czmn_o, bsqa=>bzmn_e, &
27 bsubsa=>armn_e, bsubua=>azmn_e, bsubva=>armn_o
28 USE vacmod, ONLY: potvac
29 USE ezcdf
30 USE read_wout_mod, ONLY: vn_version, vn_extension, vn_mgrid, &
31 vn_magen, vn_therm, vn_gam, vn_maxr, vn_minr, vn_maxz, vn_fp, &
32 vn_radnod, vn_polmod, vn_tormod, vn_maxmod, vn_maxit, &
33 vn_asym, vn_free, vn_error, vn_aspect, vn_beta, &
34 vn_pbeta, vn_tbeta, vn_abeta, vn_b0, vn_rbt0, vn_maxmod_nyq, &
35 vn_rbt1, vn_sgs, vn_lar, vn_modb, vn_ctor, vn_amin, vn_rmaj, &
36 vn_vol, vn_ac, vn_ai, vn_am, &
37 vn_pmass_type, vn_pcurr_type, vn_piota_type, &
38 vn_am_aux_s, vn_am_aux_f, vn_ac_aux_s, vn_ac_aux_f, &
39 vn_ai_aux_s, vn_ai_aux_f, &
40 vn_ftolv, vn_fsqr, vn_fsqz, vn_fsql, &
41 vn_pmod, vn_tmod, vn_pmod_nyq, vn_tmod_nyq, &
42 vn_racc, vn_zacs, vn_racs, vn_zacc, vn_iotaf, vn_qfact, &
43 vn_presf, vn_phi, vn_phipf, vn_jcuru, vn_jcurv, vn_iotah, &
44 vn_chi, vn_chipf, &
45 vn_mass, vn_presh, vn_betah, vn_buco, vn_bvco, vn_vp, vn_specw, &
46 vn_phip, vn_jdotb, vn_bdotb, vn_overr, vn_bgrv, vn_merc, &
47 vn_mshear, vn_mwell, vn_mcurr, vn_mgeo, vn_equif, vn_fsq, &
48 vn_wdot, vn_extcur, vn_curlab, vn_rmnc, vn_zmns, vn_lmns, &
49 vn_gmnc, vn_bmnc, vn_bsubumnc, vn_bsubvmnc, vn_bsubsmns, &
50 vn_bsupumnc, vn_bsupvmnc, vn_rmns, vn_zmnc, vn_lmnc, vn_gmns, &
51 vn_bmns, vn_bsubumns, vn_bsubvmns, vn_bsubsmnc, vn_bsupumns, &
52 vn_bsupvmns, vn_rbc, vn_zbs, vn_rbs, vn_zbc, vn_potvac, &
53 vn_currumnc, vn_currvmnc, vn_currumns, vn_currvmns, &
54 compute_currents, &
55 ln_version, ln_extension, ln_mgrid, &
56 ln_magen, ln_therm, ln_gam, ln_maxr, ln_minr, ln_maxz, ln_fp, &
57 ln_radnod, ln_polmod, ln_tormod, ln_maxmod, ln_maxit, &
58 ln_asym, ln_free, ln_error, ln_aspect, ln_beta, &
59 ln_pbeta, ln_tbeta, ln_abeta, ln_b0, ln_rbt0, ln_maxmod_nyq, &
60 ln_rbt1, ln_sgs, ln_lar, ln_modb, ln_ctor, ln_amin, ln_rmaj, &
61 ln_mse, ln_thom, ln_next, &
62 ln_pmod, ln_tmod, ln_pmod_nyq, ln_tmod_nyq, ln_racc, ln_zacs, &
63 ln_racs, ln_zacc, ln_iotaf, ln_qfact, ln_am, ln_ac, ln_ai, &
64 ln_pmass_type, ln_pcurr_type, ln_piota_type, &
65 ln_am_aux_s, ln_am_aux_f, ln_ac_aux_s, ln_ac_aux_f, &
66 ln_ai_aux_s, ln_ai_aux_f, ln_chi, ln_chipf, &
67 ln_presf, ln_phi, ln_phipf, ln_jcuru, ln_jcurv, ln_iotah, &
68 ln_mass, ln_presh, ln_betah, ln_buco, ln_bvco, ln_vp, ln_specw, &
69 ln_vol, ln_phip, ln_jdotb, ln_bdotb, ln_bgrv, ln_merc, &
70 ln_mshear, ln_mwell, ln_mcurr, ln_mgeo, ln_equif, ln_fsq, &
71 ln_wdot, ln_extcur, ln_curlab, ln_rmnc, ln_zmns, ln_lmns, &
72 ln_gmnc, ln_bmnc, ln_bsubumnc, ln_bsubvmnc, ln_bsubsmns, &
73 ln_bsupumnc, ln_bsupvmnc, ln_rmns, ln_zmnc, ln_lmnc, ln_gmns, &
74 ln_bmns, ln_bsubumns, ln_bsubvmns, ln_bsubsmnc, ln_bsupumns, &
75 ln_bsupvmns, ln_rbc, ln_zbs, ln_rbs, ln_zbc, ln_potvac, &
76 ln_currumnc, ln_currvmnc, ln_currumns, ln_currvmns
77
78 USE safe_open_mod
79 USE mgrid_mod
80
81 IMPLICIT NONE
82
83 INTEGER, INTENT(in) :: ier_flag
84
85 ! reverse ns, mnmax for backwards compatibility
86 REAL(rprec), DIMENSION(mnmax,ns,3*MAX(ntmax/2,1)), INTENT(inout), TARGET :: rzl_array, gc_array ! only used as temporary storage !
87 REAL(rprec), DIMENSION(ns,nznt), INTENT(inout) :: bsq, gsqrt, bsubu, bsubv, bsubs, bsupv, bsupu
88
89 REAL(rprec) :: qfact(ns)
90
91 CHARACTER(LEN=*), PARAMETER, DIMENSION(1) :: &
92 r1dim = (/'radius'/), &
93 mn1dim = (/'mn_mode'/), &
94 mn2dim = (/'mn_mode_nyq'/), &
95 mnpddim = (/'mnpd'/), &
96 currg = (/'ext_current'/), &
97 currl = (/'current_label'/)
98 CHARACTER(LEN=*), PARAMETER, DIMENSION(2) :: &
99 r2dim = (/'mn_mode', 'radius ' /), &
100 r3dim = (/'mn_mode_nyq','radius '/)
101
102 INTEGER :: j, js, jlk, mn, lk, &
103 m, n, k, iwout0, n1, nwout, istat, indx1(1)
104 REAL(rprec) :: dmult, tcosi, tsini, vversion, sgn, tmult, &
105 ftolx1
106 REAL(rprec), POINTER, DIMENSION(:,:) :: rmnc, rmns, zmns, zmnc, lmns, lmnc
107 REAL(rprec), ALLOCATABLE, DIMENSION(:,:) :: gmnc, bmnc, gmns, bmns, &
108 bsubumnc, bsubvmnc, bsubsmns, bsubumns, bsubvmns, bsubsmnc, &
109 currumnc, currvmnc, currumns, currvmns
110 REAL(rprec), DIMENSION(mnmax) :: rmnc1, zmns1, lmns1, &
111 rmns1, zmnc1, lmnc1, bmodmn, bmodmn1
112 REAL(rprec), DIMENSION(:), ALLOCATABLE :: gmn, bmn, &
113 bsubumn, bsubvmn, bsubsmn, bsupumn, bsupvmn
114 CHARACTER(LEN=120) :: wout_file
115 REAL(rprec), DIMENSION(:), ALLOCATABLE :: xfinal
116
117 ! ELIMINATE THESE EVENTUALLY
118 REAL(rprec), ALLOCATABLE, DIMENSION(:,:) :: bsupumnc, bsupumns, bsupvmnc, bsupvmns
119
120 INTEGER :: ierr_curr
121
122
123 ! THIS SUBROUTINE CREATES THE FILE WOUT.
124 ! IT CONTAINS THE CYLINDRICAL COORDINATE SPECTRAL COEFFICIENTS
125 ! RMN,ZMN (full), LMN (half_mesh - CONVERTED FROM INTERNAL full REPRESENTATION),
126 ! AS WELL AS COEFFICIENTS (ON NYQ MESH) FOR COMPUTED QUANTITIES:
127 ! BSQ, BSUPU,V, BSUBU,V, GSQRT (HALF); BSUBS (FULL-CONVERTED IN JXBFORCE)
128
129
130 ! Pointer assignments for storage arrays
131 n1 = max(1,ntmax/2)
132 rmnc => rzl_array(:,:,1) ! store COS(mu-nv) components
133 zmns => rzl_array(:,:,1+n1) ! store SIN(mu-nv)
134 lmns => rzl_array(:,:,1+2*n1) ! store SIN(mu-nv)
135
136 IF (lasym) THEN
137 rmns => gc_array(:,:,1) ! store SIN(mu-nv)
138 zmnc => gc_array(:,:,1+n1) ! store COS(mu-nv)
139 lmnc => gc_array(:,:,1+2*n1) ! store COS(mu-nv)
140 END IF
141
142 ALLOCATE (gmn(mnmax_nyq), bmn(mnmax_nyq), &
143 bsubumn(mnmax_nyq), bsubvmn(mnmax_nyq), bsubsmn(mnmax_nyq), &
144 bsupumn(mnmax_nyq), bsupvmn(mnmax_nyq), stat=istat)
145
146 ALLOCATE (gmnc(mnmax_nyq,ns), bmnc(mnmax_nyq,ns), &
147 bsubumnc(mnmax_nyq,ns), bsubvmnc(mnmax_nyq,ns), &
148 bsubsmns(mnmax_nyq,ns), bsupumnc(mnmax_nyq,ns), &
149 bsupvmnc(mnmax_nyq,ns), &
150 currumnc(mnmax_nyq,ns), currvmnc(mnmax_nyq,ns), &
151 stat=istat)
152 IF (lasym) THEN
153 ALLOCATE (gmns(mnmax_nyq,ns), bmns(mnmax_nyq,ns), &
154 bsubumns(mnmax_nyq,ns), bsubvmns(mnmax_nyq,ns), &
155 bsubsmnc(mnmax_nyq,ns), bsupumns(mnmax_nyq,ns), &
156 bsupvmns(mnmax_nyq,ns), &
157 currumns(mnmax_nyq,ns), currvmns(mnmax_nyq,ns), &
158 stat=istat)
159 END IF
160 IF (istat .ne. 0) stop 'Error allocating arrays in VMEC WROUT'
161
162 ! IF (nextcur .eq. 0) THEN
163 ! DO j = SIZE(extcur), 1, -1
164 ! IF (extcur(j) .ne. zero) THEN
165 ! nextcur = j
166 ! EXIT
167 ! END IF
168 ! END DO
169 ! END IF
170
171 ! ftol info evaluated here!
172 indx1=maxloc(ns_array)
173 ftolx1=ftol_array(indx1(1))
174
175 ! NYQUIST FREQUENCY REQUIRES FACTOR OF 1/2
176 IF (mnyq .ne. 0) cosmui(:,mnyq) = p5*cosmui(:,mnyq)
177 IF (nnyq .ne. 0) cosnv(:,nnyq) = p5*cosnv(:,nnyq)
178
179 ! use wout_file as temporary storage for parsing version number to REAL
180 wout_file = version_
181 READ (wout_file, *) vversion
182
183 wout_file = 'wout_' // trim(input_extension) // '.nc'
184 CALL cdf_open(nwout, wout_file, 'w', iwout0)
185 IF (iwout0 .ne. 0) stop 'Error opening wout.nc file VMEC WROUT'
186
187 !================================
188 ! Define Variables
189 !================================
190 ! Scalars
191 CALL cdf_define(nwout, vn_version, vversion)
192 CALL cdf_define(nwout, vn_extension, input_extension)
193 CALL cdf_define(nwout, vn_mgrid, mgrid_file)
194 CALL cdf_define(nwout, vn_pcurr_type, pcurr_type)
195 CALL cdf_define(nwout, vn_pmass_type, pmass_type)
196 CALL cdf_define(nwout, vn_piota_type, piota_type)
197 CALL cdf_define(nwout, vn_magen, wb)
198 CALL cdf_define(nwout, vn_therm, wp)
199 CALL cdf_define(nwout, vn_gam, gamma)
200 CALL cdf_define(nwout, vn_maxr, rmax_surf)
201 CALL cdf_define(nwout, vn_minr, rmin_surf)
202 CALL cdf_define(nwout, vn_maxz, zmax_surf)
203 CALL cdf_define(nwout, vn_fp, nfp)
204 CALL cdf_define(nwout, vn_radnod, ns)
205 CALL cdf_define(nwout, vn_polmod, mpol)
206 CALL cdf_define(nwout, vn_tormod, ntor)
207 CALL cdf_define(nwout, vn_maxmod, mnmax)
208 CALL cdf_define(nwout, vn_maxmod_nyq, mnmax_nyq)
209 CALL cdf_define(nwout, vn_maxit, iter2)
210 CALL cdf_define(nwout, vn_asym, lasym)
211 CALL cdf_define(nwout, vn_free, lfreeb)
212 CALL cdf_define(nwout, vn_error, ier_flag)
213 CALL cdf_define(nwout, vn_aspect, aspect)
214 CALL cdf_define(nwout, vn_beta, betatot)
215 CALL cdf_define(nwout, vn_pbeta, betapol)
216 CALL cdf_define(nwout, vn_tbeta, betator)
217 CALL cdf_define(nwout, vn_abeta, betaxis)
218 CALL cdf_define(nwout, vn_b0, b0)
219 CALL cdf_define(nwout, vn_rbt0, rbtor0)
220 CALL cdf_define(nwout, vn_rbt1, rbtor)
221 CALL cdf_define(nwout, vn_sgs, nint(signgs))
222 CALL cdf_define(nwout, vn_lar, ionlarmor)
223 CALL cdf_define(nwout, vn_modb, volavgb)
224 CALL cdf_define(nwout, vn_ctor, ctor)
225 CALL cdf_define(nwout, vn_amin, aminor_p)
226 CALL cdf_define(nwout, vn_rmaj, rmajor_p)
227 CALL cdf_define(nwout, vn_vol, volume_p)
228 CALL cdf_define(nwout, vn_ftolv, ftolx1)
229 CALL cdf_define(nwout, vn_fsqr, fsqr)
230 CALL cdf_define(nwout, vn_fsqz, fsqz)
231 CALL cdf_define(nwout, vn_fsql, fsql)
232
233 CALL cdf_define(nwout, vn_nextcur, nextcur)
234 CALL cdf_define(nwout, vn_extcur, extcur(1:nextcur), dimname=currg)
235 CALL cdf_define(nwout, vn_mgmode, mgrid_mode)
236
237 ! 1D Arrays
238 CALL cdf_define(nwout, vn_pmod, xm, dimname=mn1dim)
239 CALL cdf_setatt(nwout, vn_pmod, ln_pmod)
240 CALL cdf_define(nwout, vn_tmod, xn, dimname=mn1dim)
241 CALL cdf_setatt(nwout, vn_tmod, ln_tmod)
242 CALL cdf_define(nwout, vn_pmod_nyq, xm_nyq, dimname=mn2dim)
243 CALL cdf_setatt(nwout, vn_pmod_nyq, ln_pmod_nyq)
244 CALL cdf_define(nwout, vn_tmod_nyq, xn_nyq, dimname=mn2dim)
245 CALL cdf_setatt(nwout, vn_tmod_nyq, ln_tmod_nyq)
246
247 CALL cdf_define(nwout, vn_racc, raxis_cc(0:ntor), dimname=(/'n_tor'/))
248 CALL cdf_setatt(nwout, vn_racc, ln_racc)
249 CALL cdf_define(nwout, vn_zacs, zaxis_cs(0:ntor), dimname=(/'n_tor'/))
250 CALL cdf_setatt(nwout, vn_zacs, ln_zacs)
251 IF (lasym) THEN
252 CALL cdf_define(nwout, vn_racs, raxis_cs(0:ntor), dimname=(/'n_tor'/))
253 CALL cdf_setatt(nwout, vn_racs, ln_racs)
254 CALL cdf_define(nwout, vn_zacc, zaxis_cc(0:ntor), dimname=(/'n_tor'/))
255 CALL cdf_setatt(nwout, vn_zacc, ln_zacc)
256 END IF
257
258 j = SIZE(am)-1
259 CALL cdf_define(nwout, vn_am, am(0:j), dimname=(/'preset'/))
260 j = SIZE(ac)-1
261 CALL cdf_define(nwout, vn_ac, ac(0:j), dimname=(/'preset'/))
262 j = SIZE(ai)-1
263 CALL cdf_define(nwout, vn_ai, ai(0:j), dimname=(/'preset'/))
264
265 j = SIZE(am_aux_s)
266 CALL cdf_define(nwout, vn_am_aux_s, am_aux_s(1:j), dimname=(/'ndfmax'/))
267 j = SIZE(am_aux_f)
268 CALL cdf_define(nwout, vn_am_aux_f, am_aux_f(1:j), dimname=(/'ndfmax'/))
269 j = SIZE(ai_aux_s)
270 CALL cdf_define(nwout, vn_ai_aux_s, ai_aux_s(1:j), dimname=(/'ndfmax'/))
271 j = SIZE(ai_aux_f)
272 CALL cdf_define(nwout, vn_ai_aux_f, ai_aux_f(1:j), dimname=(/'ndfmax'/))
273 j = SIZE(ac_aux_s)
274 CALL cdf_define(nwout, vn_ac_aux_s, ac_aux_s(1:j), dimname=(/'ndfmax'/))
275 j = SIZE(ac_aux_f)
276 CALL cdf_define(nwout, vn_ac_aux_f, ac_aux_f(1:j), dimname=(/'ndfmax'/))
277
278
279 CALL cdf_define(nwout, vn_iotaf, iotaf(1:ns), dimname=r1dim)
280 CALL cdf_setatt(nwout, vn_iotaf, ln_iotaf)
281
282 qfact=huge(qfact)
283 WHERE (iotaf(1:ns) .NE. zero) qfact=one/iotaf(1:ns)
284
285 CALL cdf_define(nwout, vn_qfact, qfact(1:ns), dimname=r1dim)
286 CALL cdf_setatt(nwout, vn_iotaf, ln_qfact)
287
288 CALL cdf_define(nwout, vn_presf, presf, dimname=r1dim)
289 CALL cdf_setatt(nwout, vn_presf, ln_presf, units='Pa')
290 CALL cdf_define(nwout, vn_phi, phi, dimname=r1dim)
291 CALL cdf_setatt(nwout, vn_phi, ln_phi, units='wb')
292 CALL cdf_define(nwout, vn_phipf, phipf, dimname=r1dim)
293 CALL cdf_setatt(nwout, vn_phipf, ln_phipf)
294 CALL cdf_define(nwout, vn_chi, chi, dimname=r1dim)
295 CALL cdf_setatt(nwout, vn_chi, ln_chi, units='wb')
296 !CALL cdf_define(nwout, vn_chipf, phipf, dimname=r1dim) ! TODO: wrong quantity !!!
297 CALL cdf_define(nwout, vn_chipf, chipf, dimname=r1dim)
298 CALL cdf_setatt(nwout, vn_chipf, ln_chipf)
299 CALL cdf_define(nwout, vn_jcuru, jcuru, dimname=r1dim)
300 CALL cdf_define(nwout, vn_jcurv, jcurv, dimname=r1dim)
301
302 CALL cdf_define(nwout, vn_iotah, iotas(1:ns), dimname=r1dim)
303 CALL cdf_setatt(nwout, vn_iotah, ln_iotah)
304 CALL cdf_define(nwout, vn_mass, mass, dimname=r1dim)
305 CALL cdf_setatt(nwout, vn_mass, ln_mass)
306 CALL cdf_define(nwout, vn_presh, pres(1:ns), dimname=r1dim)
307 CALL cdf_setatt(nwout, vn_presh, ln_presh, units='Pa')
308 CALL cdf_define(nwout, vn_betah, beta_vol, dimname=r1dim)
309 CALL cdf_define(nwout, vn_buco, buco, dimname=r1dim)
310 CALL cdf_define(nwout, vn_bvco, bvco, dimname=r1dim)
311 CALL cdf_define(nwout, vn_vp, vp(1:ns), dimname=r1dim)
312 CALL cdf_define(nwout, vn_specw, specw, dimname=r1dim) ! NOTE: specw on full grid
313 CALL cdf_define(nwout, vn_phip, phips(1:ns), dimname=r1dim)
314 CALL cdf_define(nwout, vn_overr, overr(1:ns), dimname=r1dim)
315
316 CALL cdf_define(nwout, vn_jdotb, jdotb, dimname=r1dim)
317 CALL cdf_define(nwout, vn_bdotb, bdotb, dimname=r1dim)
318 CALL cdf_define(nwout, vn_bgrv, bdotgradv, dimname=r1dim)
319
320 CALL cdf_define(nwout, vn_merc, dmerc, dimname=r1dim)
321 CALL cdf_define(nwout, vn_mshear, dshear, dimname=r1dim)
322 CALL cdf_define(nwout, vn_mwell, dwell, dimname=r1dim)
323 CALL cdf_define(nwout, vn_mcurr, dcurr, dimname=r1dim)
324 CALL cdf_define(nwout, vn_mgeo, dgeod, dimname=r1dim)
325 CALL cdf_define(nwout, vn_equif, equif, dimname=r1dim)
326
327 IF (lfreeb .and. nextcur.gt.0 .and. ALLOCATED(curlabel)) THEN
328 CALL cdf_define(nwout, vn_curlab, curlabel(1:nextcur), dimname=currl)
329 ! SAL for potvac
330 CALL cdf_define(nwout, vn_potvac, potvac, dimname=mnpddim)
331 ENDIF
332
333 ! 2D Arrays
334 CALL cdf_define(nwout, vn_rmnc, rmnc, dimname=r2dim)
335 CALL cdf_setatt(nwout, vn_rmnc, ln_rmnc, units='m')
336 CALL cdf_define(nwout, vn_zmns, zmns, dimname=r2dim)
337 CALL cdf_setatt(nwout, vn_zmns, ln_zmns, units='m')
338 CALL cdf_define(nwout, vn_lmns, lmns, dimname=r2dim)
339 CALL cdf_setatt(nwout, vn_lmns, ln_lmns)
340 CALL cdf_define(nwout, vn_gmnc, gmnc, dimname=r3dim)
341 CALL cdf_setatt(nwout, vn_gmnc, ln_gmnc)
342 CALL cdf_define(nwout, vn_bmnc, bmnc, dimname=r3dim)
343 CALL cdf_setatt(nwout, vn_bmnc, ln_bmnc)
344 CALL cdf_define(nwout, vn_bsubumnc, bsubumnc, dimname=r3dim)
345 CALL cdf_setatt(nwout, vn_bsubumnc, ln_bsubumnc)
346 CALL cdf_define(nwout, vn_bsubvmnc, bsubvmnc, dimname=r3dim)
347 CALL cdf_setatt(nwout, vn_bsubvmnc, ln_bsubvmnc)
348 CALL cdf_define(nwout, vn_bsubsmns, bsubsmns, dimname=r3dim)
349 CALL cdf_setatt(nwout, vn_bsubsmns, ln_bsubsmns)
350
351 ! ELIMINATE THESE EVENTUALLY: DON'T NEED THEM - CAN COMPUTE FROM GSQRT
352 CALL cdf_define(nwout, vn_bsupumnc, bsupumnc, dimname=r3dim)
353 CALL cdf_define(nwout, vn_bsupvmnc, bsupvmnc, dimname=r3dim)
354
355 CALL cdf_define(nwout, vn_currumnc, currumnc, dimname=r3dim)
356 CALL cdf_setatt(nwout, vn_currumnc, ln_currumnc)
357 CALL cdf_define(nwout, vn_currvmnc, currvmnc, dimname=r3dim)
358 CALL cdf_setatt(nwout, vn_currvmnc, ln_currvmnc)
359
360 IF (lasym) then
361 CALL cdf_define(nwout, vn_rmns, rmns, dimname=r2dim)
362 CALL cdf_setatt(nwout, vn_rmns, ln_rmns, units='m')
363 CALL cdf_define(nwout, vn_zmnc, zmnc, dimname=r2dim)
364 CALL cdf_setatt(nwout, vn_zmnc, ln_zmnc, units='m')
365 CALL cdf_define(nwout, vn_lmnc, lmnc, dimname=r2dim)
366 CALL cdf_setatt(nwout, vn_lmnc, ln_lmnc)
367 CALL cdf_define(nwout, vn_gmns, gmns, dimname=r3dim)
368 CALL cdf_setatt(nwout, vn_gmns, ln_gmns)
369 CALL cdf_define(nwout, vn_bmns, bmns, dimname=r3dim)
370 CALL cdf_setatt(nwout, vn_bmns, ln_bmns)
371 CALL cdf_define(nwout, vn_bsubumns, bsubumns, dimname=r3dim)
372 CALL cdf_setatt(nwout, vn_bsubumns, ln_bsubumns)
373 CALL cdf_define(nwout, vn_bsubvmns, bsubvmns, dimname=r3dim)
374 CALL cdf_setatt(nwout, vn_bsubvmns, ln_bsubvmns)
375 CALL cdf_define(nwout, vn_bsubsmnc, bsubsmnc, dimname=r3dim)
376 CALL cdf_setatt(nwout, vn_bsubsmnc, ln_bsubsmnc)
377
378 ! ELIMINATE THESE EVENTUALLY: DON'T NEED THEM
379 CALL cdf_define(nwout, vn_bsupumns, bsupumns, dimname=r3dim)
380 CALL cdf_define(nwout, vn_bsupvmns, bsupvmns, dimname=r3dim)
381
382 CALL cdf_define(nwout, vn_currumns, currumns, dimname=r3dim)
383 CALL cdf_setatt(nwout, vn_currumns, ln_currumns)
384 CALL cdf_define(nwout, vn_currvmns, currvmns, dimname=r3dim)
385 CALL cdf_setatt(nwout, vn_currvmns, ln_currvmns)
386 end if
387
388 !================================
389 ! Write Variables
390 !================================
391 ! Scalars
392 CALL cdf_write(nwout, vn_version, vversion)
393 CALL cdf_write(nwout, vn_extension, input_extension)
394 CALL cdf_write(nwout, vn_mgrid, mgrid_file)
395 CALL cdf_write(nwout, vn_pcurr_type, pcurr_type)
396 CALL cdf_write(nwout, vn_piota_type, piota_type)
397 CALL cdf_write(nwout, vn_pmass_type, pmass_type)
398 CALL cdf_write(nwout, vn_magen, wb)
399 CALL cdf_write(nwout, vn_therm, wp)
400 CALL cdf_write(nwout, vn_gam, gamma)
401 CALL cdf_write(nwout, vn_maxr, rmax_surf)
402 CALL cdf_write(nwout, vn_minr, rmin_surf)
403 CALL cdf_write(nwout, vn_maxz, zmax_surf)
404 CALL cdf_write(nwout, vn_fp, nfp)
405 CALL cdf_write(nwout, vn_radnod, ns)
406 CALL cdf_write(nwout, vn_polmod, mpol)
407 CALL cdf_write(nwout, vn_tormod, ntor)
408 CALL cdf_write(nwout, vn_maxmod, mnmax)
409 CALL cdf_write(nwout, vn_maxmod_nyq, mnmax_nyq)
410 CALL cdf_write(nwout, vn_maxit, iter2)
411 CALL cdf_write(nwout, vn_asym, lasym)
412 CALL cdf_write(nwout, vn_free, lfreeb)
413 CALL cdf_write(nwout, vn_error, ier_flag)
414
415 CALL cdf_write(nwout, vn_aspect, aspect)
416 CALL cdf_write(nwout, vn_beta, betatot)
417 CALL cdf_write(nwout, vn_pbeta, betapol)
418 CALL cdf_write(nwout, vn_tbeta, betator)
419 CALL cdf_write(nwout, vn_abeta, betaxis)
420 CALL cdf_write(nwout, vn_b0, b0)
421 CALL cdf_write(nwout, vn_rbt0, rbtor0)
422 CALL cdf_write(nwout, vn_rbt1, rbtor)
423 CALL cdf_write(nwout, vn_sgs, nint(signgs))
424 CALL cdf_write(nwout, vn_lar, ionlarmor)
425 CALL cdf_write(nwout, vn_modb, volavgb)
426 CALL cdf_write(nwout, vn_ctor, ctor/mu0)
427 CALL cdf_write(nwout, vn_amin, aminor_p)
428 CALL cdf_write(nwout, vn_rmaj, rmajor_p)
429 CALL cdf_write(nwout, vn_vol, volume_p)
430 CALL cdf_write(nwout, vn_ftolv, ftolx1)
431 CALL cdf_write(nwout, vn_fsql, fsql)
432 CALL cdf_write(nwout, vn_fsqr, fsqr)
433 CALL cdf_write(nwout, vn_fsqz, fsqz)
434
435 CALL cdf_write(nwout, vn_nextcur, nextcur)
436 IF (nextcur .gt. 0) THEN
437 CALL cdf_write(nwout, vn_extcur, extcur(1:nextcur))
438 CALL cdf_write(nwout, vn_mgmode, mgrid_mode)
439 ENDIF
440 IF (lfreeb) THEN
441 ! TODO: write current labels
442 CALL cdf_write(nwout, vn_potvac, potvac)
443 END IF
444
445 ! 1D Arrays
446 CALL cdf_write(nwout, vn_pmod, xm)
447 CALL cdf_write(nwout, vn_tmod, xn)
448 CALL cdf_write(nwout, vn_pmod_nyq, xm_nyq)
449 CALL cdf_write(nwout, vn_tmod_nyq, xn_nyq)
450
451
452 ! here start the convert-to-regular-Fourier-output part
453
454
455 ALLOCATE (xfinal(neqs), stat=js) ! re-use js for allocation return code
456 IF (js .ne. 0) stop 'Allocation error for xfinal in WROUT!'
457 xfinal = xc ! ignore passed rzl_array !!!
458 ! TODO: scalxc ???
459 ! --> rzl_array is re-used to store temporary symmetric rmnc, ... !
460 ! --> gc_array is re-used to store temporary asymmetric rmns, ... ! (for lasym)
461
462 ! MUST CONVERT m=1 MODES... FROM INTERNAL TO PHYSICAL FORM
463 ! Extrapolation of m=0 Lambda (cs) modes, which are not evolved at j=1, done in CONVERT
464
465 IF (lthreed) CALL convert_sym (xfinal(1+mns*(rss-1)), xfinal(1+irzloff+mns*(zcs-1)))
466 IF (lasym) CALL convert_asym (xfinal(1+mns*(rsc-1)), xfinal(1+irzloff+mns*(zcc-1)))
467
468 ! CONVERT TO rmnc, zmns, lmns, etc EXTERNAL representation (without internal mscale, nscale)
469 ! IF B^v ~ phip + lamu, MUST DIVIDE BY phipf(js) below to maintain old-style format
470 radius1: DO js = 1, ns
471
472 CALL convert (rmnc1, zmns1, lmns1, rmns1, zmnc1, lmnc1, xfinal, js)
473
474 rmnc(:,js) = rmnc1(:)
475 zmns(:,js) = zmns1(:)
476 lmns(:,js) = (lmns1(:)/phipf(js)) * lamscale
477 IF (lasym) THEN
478 rmns(:,js) = rmns1(:)
479 zmnc(:,js) = zmnc1(:)
480 lmnc(:,js) = (lmnc1(:)/phipf(js)) * lamscale
481 END IF
482
483 END DO radius1
484
485 DEALLOCATE (xfinal)
486
487 ! INTERPOLATE LAMBDA ONTO HALF-MESH FOR BACKWARDS CONSISTENCY WITH EARLIER VERSIONS OF VMEC
488 ! AND SMOOTHS POSSIBLE UNPHYSICAL "WIGGLE" ON RADIAL MESH
489 WHERE (nint(xm) .le. 1) lmns(:,1) = lmns(:,2) ! constant extrapolation of lambda
490 DO js = ns,2,-1
491 WHERE (mod(nint(xm),2) .eq. 0)
492 lmns(:,js) = p5*(lmns(:,js) + lmns(:,js-1))
493 ELSEWHERE
494 lmns(:,js) = p5*(sm(js)*lmns(:,js) + sp(js-1)*lmns(:,js-1))
495 END WHERE
496 END DO
497
498 lmns(:,1) = 0 ! NOTE: needed constant extrapolation above only in averaging step!
499 raxis_cc(0:ntor) = rmnc(1:ntor+1,1)
500 zaxis_cs(0:ntor) = zmns(1:ntor+1,1)
501
502 IF (lasym) then
503 WHERE (nint(xm) .le. 1) lmnc(:,1) = lmnc(:,2)
504 DO js = ns,2,-1
505 WHERE (mod(nint(xm),2) .eq. 0)
506 lmnc(:,js) = p5*(lmnc(:,js) + lmnc(:,js-1))
507 ELSEWHERE
508 lmnc(:,js) = p5*(sm(js)*lmnc(:,js) + sp(js-1)*lmnc(:,js-1))
509 END WHERE
510 END DO
511
512 lmnc(:,1) = 0 ! NOTE: needed constant extrapolation above only in averaging step!
513 raxis_cs(0:ntor) = rmns(1:ntor+1,1)
514 zaxis_cc(0:ntor) = zmnc(1:ntor+1,1)
515 end if
516
517
518
519
520
521 ! here comes the realspace-to-Fourier transform part...
522
523
524 ! COMPUTE |B| = SQRT(|B|**2) and store in bsq, bsqa
525 DO js = 2, ns
526 bsq(js,:nznt) = sqrt(2*abs(bsq(js,:nznt)-pres(js)))
527 END DO
528
529 tmult = p5/r0scale**2
530 !SPH: FIXED THIS 03-05-07 TO CALL symmetrization routine
531 IF (lasym) THEN
532 ! Changed integration norm in fixaray, SPH012314
533 tmult = 2*tmult
534 bsubs(1,:) = 0
535 CALL symoutput (bsq, gsqrt, bsubu, bsubv, bsupu, &
536 bsupv, bsubs, &
537 bsqa, gsqrta, bsubua, bsubva, bsupua, &
538 bsupva, bsubsa )
539 END IF
540
541 ! Fourier-transform derived quantities for each surface individually
542 radius2: DO js = 2, ns
543 gmn = 0
544 bmn = 0
545
546 bsubumn = 0
547 bsubvmn = 0
548 bsubsmn = 0
549
550 bsupumn = 0
551 bsupvmn = 0
552
553 mn2: DO mn = 1, mnmax_nyq
554 n = nint(xn_nyq(mn))/nfp
555 m = nint(xm_nyq(mn))
556 n1 = abs(n)
557
558 dmult = mscale(m)*nscale(n1)*tmult
559 IF (m.eq.0 .or. n.eq.0) dmult = 2*dmult
560
561 sgn = sign(1, n)
562
563 lk = 0
564 DO j = 1, ntheta2
565 DO k = 1, nzeta
566 lk = lk + 1
567
568 ! cos(mu - nv)
569 tcosi = dmult*(cosmui(j,m)*cosnv(k,n1) + sgn*sinmui(j,m)*sinnv(k,n1))
570
571 ! sin(mu - nv)
572 tsini = dmult*(sinmui(j,m)*cosnv(k,n1) - sgn*cosmui(j,m)*sinnv(k,n1))
573
574 bmn(mn) = bmn(mn) + tcosi*bsq(js,lk)
575 gmn(mn) = gmn(mn) + tcosi*gsqrt(js,lk)
576 bsubumn(mn) = bsubumn(mn) + tcosi*bsubu(js,lk)
577 bsubvmn(mn) = bsubvmn(mn) + tcosi*bsubv(js,lk)
578 bsubsmn(mn) = bsubsmn(mn) + tsini*bsubs(js,lk)
579 bsupumn(mn) = bsupumn(mn) + tcosi*bsupu(js,lk)
580 bsupvmn(mn) = bsupvmn(mn) + tcosi*bsupv(js,lk)
581 END DO
582 END DO
583 END DO mn2
584
585 ! needed for freeb_data below
586 IF (js .eq. ns/2) bmodmn = bmn(1:mnmax)
587 IF (js .eq. ns ) bmodmn1 = bmn(1:mnmax)
588
589 gmnc(:,js) = gmn(:)
590 bmnc(:,js) = bmn(:)
591 bsubumnc(:,js) = bsubumn(:)
592 bsubvmnc(:,js) = bsubvmn(:)
593 bsubsmns(:,js) = bsubsmn(:)
594 bsupumnc(:,js) = bsupumn(:)
595 bsupvmnc(:,js) = bsupvmn(:)
596 END DO radius2
597
598 ! endpoint values at magnetic axis
599 gmnc(:,1) = 0
600 bmnc(:,1) = 0
601
602 bsubumnc(:,1) = 0
603 bsubvmnc(:,1) = 0
604
605 ! NOTE: This assumes that bsubs (from which bsubsmns is computed) is on the full grid.
606 ! However, since crmn_o is passed as bsubs in fileout(),
607 ! HERE, bsubs is ACTUALLY on the HALF-grid !!! WTF ???
608 ! (The full-grid bsubs array is a LOCAL variable in jxbforce().)
609 bsubsmns(:,1) = 2*bsubsmns(:,2) - bsubsmns(:,3) ! extrapolation on full grid
610
611 bsupumnc(:,1) = 0
612 bsupvmnc(:,1) = 0
613
614 IF (lasym) then
615 radius3: DO js = 2, ns
616 gmn = 0
617 bmn = 0
618 bsubumn = 0
619 bsubvmn = 0
620 bsubsmn = 0
621 bsupumn = 0
622 bsupvmn = 0
623 mn3: DO mn = 1, mnmax_nyq
624 n = nint(xn_nyq(mn))/nfp
625 m = nint(xm_nyq(mn))
626 n1 = abs(n)
627
628 dmult = mscale(m)*nscale(n1)*tmult
629 IF (m.eq.0 .or. n.eq.0) dmult = 2*dmult
630
631 sgn = sign(1, n)
632
633 lk = 0
634 jlk = js
635 DO j = 1, ntheta2
636 DO k = 1, nzeta
637 lk = lk + 1
638
639 tcosi = dmult*(cosmui(j,m)*cosnv(k,n1) + sgn*sinmui(j,m)*sinnv(k,n1))
640 tsini = dmult*(sinmui(j,m)*cosnv(k,n1) - sgn*cosmui(j,m)*sinnv(k,n1))
641
642 bmn(mn) = bmn(mn) + tsini*bsqa(jlk)
643 gmn(mn) = gmn(mn) + tsini*gsqrta(jlk,0) ! TODO: what is this zero index?
644
645 bsubumn(mn) = bsubumn(mn) + tsini*bsubua(jlk)
646 bsubvmn(mn) = bsubvmn(mn) + tsini*bsubva(jlk)
647 bsubsmn(mn) = bsubsmn(mn) + tcosi*bsubsa(jlk)
648
649 bsupumn(mn) = bsupumn(mn) + tsini*bsupua(jlk)
650 bsupvmn(mn) = bsupvmn(mn) + tsini*bsupva(jlk)
651 jlk = jlk+ns
652 END DO
653 END DO
654 END DO mn3
655
656 gmns(:,js) = gmn(:)
657 bmns(:,js) = bmn(:)
658 bsubumns(:,js) = bsubumn(:)
659 bsubvmns(:,js) = bsubvmn(:)
660 bsubsmnc(:,js) = bsubsmn(:)
661 bsupumns(:,js) = bsupumn(:)
662 bsupvmns(:,js) = bsupvmn(:)
663 END DO radius3
664
665 gmns(:,1) = 0
666 bmns(:,1) = 0
667 bsubumns(:,1) = 0
668 bsubvmns(:,1) = 0
669 bsubsmnc(:,1) = 2*bsubsmnc(:,2) - bsubsmnc(:,3)
670 bsupumns(:,1) = 0
671 bsupvmns(:,1) = 0
672 end if
673
674 ! Compute current density Fourier coefficients from covariant B components.
675 ! Reuses Compute_Currents from read_wout_mod, following PARVMEC wrout.f.
676 ! Allocate dummy asymmetric arrays when lasym=.false. since
677 ! Compute_Currents requires all arguments to be present.
678 IF (.NOT. lasym) THEN
679 ALLOCATE (bsubsmnc(mnmax_nyq,ns), bsubumns(mnmax_nyq,ns), &
680 bsubvmns(mnmax_nyq,ns), &
681 currumns(mnmax_nyq,ns), currvmns(mnmax_nyq,ns), &
682 stat=istat)
683 bsubsmnc = 0; bsubumns = 0; bsubvmns = 0
684 END IF
685 CALL compute_currents(bsubsmnc, bsubsmns, bsubumnc, bsubumns, &
686 bsubvmnc, bsubvmns, &
687 xm_nyq, xn_nyq, mnmax_nyq, lasym, ns, &
688 currumnc, currvmnc, currumns, currvmns)
689
690 ! WRITE OUT ARRAYS
691 CALL cdf_write(nwout, vn_racc, raxis_cc(0:ntor))
692 CALL cdf_write(nwout, vn_zacs, zaxis_cs(0:ntor))
693 CALL cdf_write(nwout, vn_rmnc, rmnc)
694 CALL cdf_write(nwout, vn_zmns, zmns)
695 CALL cdf_write(nwout, vn_lmns, lmns)
696 CALL cdf_write(nwout, vn_gmnc, gmnc) !Half mesh
697 CALL cdf_write(nwout, vn_bmnc, bmnc) !Half mesh
698 CALL cdf_write(nwout, vn_bsubumnc, bsubumnc) !Half mesh
699 CALL cdf_write(nwout, vn_bsubvmnc, bsubvmnc) !Half mesh
700 CALL cdf_write(nwout, vn_bsubsmns, bsubsmns) !Full mesh
701
702 ! GET RID OF THESE EVENTUALLY: DON'T NEED THEM (can express in terms of lambdas)
703 CALL cdf_write(nwout, vn_bsupumnc, bsupumnc)
704 CALL cdf_write(nwout, vn_bsupvmnc, bsupvmnc)
705
706 CALL cdf_write(nwout, vn_currumnc, currumnc) !Full mesh
707 CALL cdf_write(nwout, vn_currvmnc, currvmnc) !Full mesh
708
709 ! FULL-MESH quantities
710 ! NOTE: jdotb is in units_of_A (1/mu0 incorporated in jxbforce...)
711 ! prior to version 6.00, this was output in internal VMEC units...
712 j = SIZE(am)-1
713 CALL cdf_write(nwout, vn_am, am(0:j))
714 j = SIZE(ac)-1
715 CALL cdf_write(nwout, vn_ac, ac(0:j))
716 j = SIZE(ai)-1
717 CALL cdf_write(nwout, vn_ai, ai(0:j))
718
719 j = SIZE(am_aux_s)
720 CALL cdf_write(nwout, vn_am_aux_s, am_aux_s(1:j))
721 j = SIZE(am_aux_f)
722 CALL cdf_write(nwout, vn_am_aux_f, am_aux_f(1:j))
723 j = SIZE(ac_aux_s)
724 CALL cdf_write(nwout, vn_ac_aux_s, ac_aux_s(1:j))
725 j = SIZE(ac_aux_f)
726 CALL cdf_write(nwout, vn_ac_aux_f, ac_aux_f(1:j))
727 j = SIZE(ai_aux_s)
728 CALL cdf_write(nwout, vn_ai_aux_s, ai_aux_s(1:j))
729 j = SIZE(ai_aux_f)
730 CALL cdf_write(nwout, vn_ai_aux_f, ai_aux_f(1:j))
731
732 CALL cdf_write(nwout, vn_iotaf, iotaf(1:ns))
733 CALL cdf_write(nwout, vn_qfact, qfact(1:ns))
734 CALL cdf_write(nwout, vn_presf, presf/mu0) ! NOTE: scaling !!!
735 CALL cdf_write(nwout, vn_phi, phi)
736 CALL cdf_write(nwout, vn_phipf, twopi*signgs*phipf) ! NOTE: scaling !!!
737 CALL cdf_write(nwout, vn_chi, chi)
738 CALL cdf_write(nwout, vn_chipf, twopi*signgs*chipf) ! NOTE: scaling !!!
739 CALL cdf_write(nwout, vn_jcuru, jcuru/mu0) ! NOTE: scaling !!!
740 CALL cdf_write(nwout, vn_jcurv, jcurv/mu0) ! NOTE: scaling !!!
741 CALL cdf_write(nwout, vn_jdotb, jdotb)
742 CALL cdf_write(nwout, vn_bdotb, bdotb)
743 CALL cdf_write(nwout, vn_bgrv, bdotgradv)
744
745 ! HALF-MESH quantities
746 iotas(1) = 0
747 mass(1) = 0
748 pres(1) = 0
749 phip(1) = 0
750 buco(1) = 0
751 bvco(1) = 0
752 vp(1) = 0
753 overr(1) = 0
754 specw(1) = 1
755 beta_vol(1) = 0
756
757 CALL cdf_write(nwout, vn_iotah, iotas(1:ns))
758 CALL cdf_write(nwout, vn_mass, mass/mu0) ! NOTE: scaling !!!
759 CALL cdf_write(nwout, vn_presh, pres(1:ns)/mu0) ! NOTE: scaling !!!
760 CALL cdf_write(nwout, vn_betah, beta_vol)
761 CALL cdf_write(nwout, vn_buco, buco)
762 CALL cdf_write(nwout, vn_bvco, bvco)
763 CALL cdf_write(nwout, vn_vp, vp(1:ns))
764 CALL cdf_write(nwout, vn_specw, specw)
765 CALL cdf_write(nwout, vn_phip, phips(1:ns))
766 CALL cdf_write(nwout, vn_overr, overr(1:ns))
767
768 ! MERCIER_CRITERION
769 CALL cdf_write(nwout, vn_merc, dmerc(1:ns))
770 CALL cdf_write(nwout, vn_mshear, dshear(1:ns))
771 CALL cdf_write(nwout, vn_mwell, dwell(1:ns))
772 CALL cdf_write(nwout, vn_mcurr, dcurr(1:ns))
773 CALL cdf_write(nwout, vn_mgeo, dgeod(1:ns))
774
775 CALL cdf_write(nwout, vn_equif, equif(1:ns))
776
777 IF (lasym) THEN
778 CALL cdf_write(nwout, vn_racs, raxis_cs(0:ntor))
779 CALL cdf_write(nwout, vn_zacc, zaxis_cc(0:ntor))
780
781 CALL cdf_write(nwout, vn_rmns, rmns)
782 CALL cdf_write(nwout, vn_zmnc, zmnc)
783 CALL cdf_write(nwout, vn_lmnc, lmnc)
784
785 CALL cdf_write(nwout, vn_gmns, gmns)
786 CALL cdf_write(nwout, vn_bmns, bmns)
787
788 CALL cdf_write(nwout, vn_bsubumns, bsubumns)
789 CALL cdf_write(nwout, vn_bsubvmns, bsubvmns)
790 CALL cdf_write(nwout, vn_bsubsmnc, bsubsmnc)
791
792 ! GET RID OF THESE EVENTUALLY: DON'T NEED THEM
793 CALL cdf_write(nwout, vn_bsupumns, bsupumns)
794 CALL cdf_write(nwout, vn_bsupvmns, bsupvmns)
795
796 CALL cdf_write(nwout, vn_currumns, currumns) !Full mesh
797 CALL cdf_write(nwout, vn_currvmns, currvmns) !Full mesh
798 END IF
799
800 CALL cdf_close(nwout)
801
802 ! RESTORE nyq ENDPOINT VALUES
803 IF (mnyq .ne. 0) cosmui(:,mnyq) = 2*cosmui(:,mnyq)
804 IF (nnyq .ne. 0) cosnv(:,nnyq) = 2*cosnv(:,nnyq)
805
806 ! DEALLOCATIONS
807 IF (ALLOCATED(gmnc)) DEALLOCATE(gmnc, bmnc, bsubumnc, bsubvmnc, &
808 bsubsmns, bsupumnc, bsupvmnc )
809 IF (ALLOCATED(gmns)) DEALLOCATE(gmns, bmns, bsupumns, bsupvmns)
810 IF (ALLOCATED(bsubumns)) DEALLOCATE(bsubumns, bsubvmns, bsubsmnc)
811 IF (ALLOCATED(currumnc)) DEALLOCATE(currumnc, currvmnc)
812 IF (ALLOCATED(currumns)) DEALLOCATE(currumns, currvmns)
813 IF (ALLOCATED(gmn)) DEALLOCATE (gmn, bmn, bsubumn, bsubvmn, &
814 bsubsmn, bsupumn, bsupvmn, stat=istat)
815
816 ! FREE BOUNDARY DATA
817 ! rmnc1, ... contain the LCFS values at this point
818 CALL freeb_data(rmnc1, zmns1, rmns1, zmnc1, bmodmn, bmodmn1)
819
820 rzl_array = 0
821
822END 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:96
read_wout_mod::ln_tmod
character(len= *), parameter ln_tmod
Definition read_wout_mod.f:96
read_wout_mod::ln_vp
character(len= *), parameter ln_vp
Definition read_wout_mod.f:96
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:96
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:96
read_wout_mod::ln_pbeta
character(len= *), parameter ln_pbeta
Definition read_wout_mod.f:96
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:96
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:96
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:96
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:96
read_wout_mod::ln_fsq
character(len= *), parameter ln_fsq
Definition read_wout_mod.f:96
read_wout_mod::ln_amin
character(len= *), parameter ln_amin
Definition read_wout_mod.f:96
read_wout_mod::ln_bsubvmns
character(len= *), parameter ln_bsubvmns
Definition read_wout_mod.f:96
read_wout_mod::ln_chipf
character(len= *), parameter ln_chipf
Definition read_wout_mod.f:96
read_wout_mod::ln_maxit
character(len= *), parameter ln_maxit
Definition read_wout_mod.f:96
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:96
read_wout_mod::ln_jcurv
character(len= *), parameter ln_jcurv
Definition read_wout_mod.f:96
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:96
read_wout_mod::ln_presh
character(len= *), parameter ln_presh
Definition read_wout_mod.f:96
read_wout_mod::ln_merc
character(len= *), parameter ln_merc
Definition read_wout_mod.f:96
read_wout_mod::ln_zacc
character(len= *), parameter ln_zacc
Definition read_wout_mod.f:96
read_wout_mod::ln_thom
character(len= *), parameter ln_thom
Definition read_wout_mod.f:96
read_wout_mod::ln_rmns
character(len= *), parameter ln_rmns
Definition read_wout_mod.f:96
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:96
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:96
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:96
read_wout_mod::vn_currvmns
character(len= *), parameter vn_currvmns
Definition read_wout_mod.f:32
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:96
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_currumns
character(len= *), parameter vn_currumns
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:96
read_wout_mod::ln_pmod
character(len= *), parameter ln_pmod
Definition read_wout_mod.f:96
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:96
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:96
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:96
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:96
read_wout_mod::ln_ai
character(len= *), parameter ln_ai
Definition read_wout_mod.f:96
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:96
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:96
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:96
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:96
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:96
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:96
read_wout_mod::mnmax_nyq
integer mnmax_nyq
Definition read_wout_mod.f:208
read_wout_mod::ln_ac_aux_s
character(len= *), parameter ln_ac_aux_s
Definition read_wout_mod.f:96
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:96
read_wout_mod::ln_zbc
character(len= *), parameter ln_zbc
Definition read_wout_mod.f:96
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:96
read_wout_mod::ln_bsupvmns
character(len= *), parameter ln_bsupvmns
Definition read_wout_mod.f:96
read_wout_mod::ln_buco
character(len= *), parameter ln_buco
Definition read_wout_mod.f:96
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:96
read_wout_mod::ln_currvmnc
character(len= *), parameter ln_currvmnc
Definition read_wout_mod.f:96
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:96
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_currvmnc
character(len= *), parameter vn_currvmnc
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:96
read_wout_mod::ln_equif
character(len= *), parameter ln_equif
Definition read_wout_mod.f:96
read_wout_mod::ln_am
character(len= *), parameter ln_am
Definition read_wout_mod.f:96
read_wout_mod::ln_currumnc
character(len= *), parameter ln_currumnc
Definition read_wout_mod.f:96
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:96
read_wout_mod::ln_version
character(len= *), parameter ln_version
Definition read_wout_mod.f:96
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:96
read_wout_mod::ln_bvco
character(len= *), parameter ln_bvco
Definition read_wout_mod.f:96
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:96
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:208
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:96
read_wout_mod::ln_racs
character(len= *), parameter ln_racs
Definition read_wout_mod.f:96
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:96
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:96
read_wout_mod::vn_bmnc
character(len= *), parameter vn_bmnc
Definition read_wout_mod.f:32
read_wout_mod::ln_bdotb
character(len= *), parameter ln_bdotb
Definition read_wout_mod.f:96
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:96
read_wout_mod::ln_polmod
character(len= *), parameter ln_polmod
Definition read_wout_mod.f:96
read_wout_mod::ln_bmnc
character(len= *), parameter ln_bmnc
Definition read_wout_mod.f:96
read_wout_mod::ln_mcurr
character(len= *), parameter ln_mcurr
Definition read_wout_mod.f:96
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:96
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:96
read_wout_mod::ln_error
character(len= *), parameter ln_error
Definition read_wout_mod.f:96
read_wout_mod::ln_piota_type
character(len= *), parameter ln_piota_type
Definition read_wout_mod.f:96
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:96
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:96
read_wout_mod::ln_mshear
character(len= *), parameter ln_mshear
Definition read_wout_mod.f:96
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:96
read_wout_mod::vn_currumnc
character(len= *), parameter vn_currumnc
Definition read_wout_mod.f:32
read_wout_mod::ln_rmnc
character(len= *), parameter ln_rmnc
Definition read_wout_mod.f:96
read_wout_mod::ln_lar
character(len= *), parameter ln_lar
Definition read_wout_mod.f:96
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:96
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:96
read_wout_mod::ln_minr
character(len= *), parameter ln_minr
Definition read_wout_mod.f:96
read_wout_mod::ln_presf
character(len= *), parameter ln_presf
Definition read_wout_mod.f:96
read_wout_mod::ln_b0
character(len= *), parameter ln_b0
Definition read_wout_mod.f:96
read_wout_mod::vn_equif
character(len= *), parameter vn_equif
Definition read_wout_mod.f:32
read_wout_mod::ln_currvmns
character(len= *), parameter ln_currvmns
Definition read_wout_mod.f:96
read_wout_mod::ln_rbt0
character(len= *), parameter ln_rbt0
Definition read_wout_mod.f:96
read_wout_mod::ln_fp
character(len= *), parameter ln_fp
Definition read_wout_mod.f:96
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:96
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:96
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:96
read_wout_mod::ln_aspect
character(len= *), parameter ln_aspect
Definition read_wout_mod.f:96
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:96
read_wout_mod::vn_bdotb
character(len= *), parameter vn_bdotb
Definition read_wout_mod.f:32
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:96
read_wout_mod::ln_bsubumnc
character(len= *), parameter ln_bsubumnc
Definition read_wout_mod.f:96
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:96
read_wout_mod::ln_sgs
character(len= *), parameter ln_sgs
Definition read_wout_mod.f:96
read_wout_mod::ln_wdot
character(len= *), parameter ln_wdot
Definition read_wout_mod.f:96
read_wout_mod::ln_pcurr_type
character(len= *), parameter ln_pcurr_type
Definition read_wout_mod.f:96
read_wout_mod::ln_tbeta
character(len= *), parameter ln_tbeta
Definition read_wout_mod.f:96
read_wout_mod::ln_maxmod
character(len= *), parameter ln_maxmod
Definition read_wout_mod.f:96
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:96
read_wout_mod::ln_tmod_nyq
character(len= *), parameter ln_tmod_nyq
Definition read_wout_mod.f:96
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:96
read_wout_mod::compute_currents
subroutine compute_currents(bsubsmnc_i, bsubsmns_i, bsubumnc_i, bsubumns_i, bsubvmnc_i, bsubvmns_i, xm_nyq_i, xn_nyq_i, mnmax_nyq_i, lasym_i, ns_i, currumnc_o, currvmnc_o, currumns_o, currvmns_o)
Definition read_wout_mod.f:703
read_wout_mod::ln_mgrid
character(len= *), parameter ln_mgrid
Definition read_wout_mod.f:96
read_wout_mod::ln_currumns
character(len= *), parameter ln_currumns
Definition read_wout_mod.f:96
read_wout_mod::ln_iotaf
character(len= *), parameter ln_iotaf
Definition read_wout_mod.f:96
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:96
read_wout_mod::ln_phi
character(len= *), parameter ln_phi
Definition read_wout_mod.f:96
read_wout_mod::ln_asym
character(len= *), parameter ln_asym
Definition read_wout_mod.f:96
read_wout_mod::ln_tormod
character(len= *), parameter ln_tormod
Definition read_wout_mod.f:96
read_wout_mod::ln_abeta
character(len= *), parameter ln_abeta
Definition read_wout_mod.f:96
read_wout_mod::ln_qfact
character(len= *), parameter ln_qfact
Definition read_wout_mod.f:96
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:96
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:96
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:96
read_wout_mod::ln_gam
character(len= *), parameter ln_gam
Definition read_wout_mod.f:96
read_wout_mod::ln_jdotb
character(len= *), parameter ln_jdotb
Definition read_wout_mod.f:96
read_wout_mod::ln_zacs
character(len= *), parameter ln_zacs
Definition read_wout_mod.f:96
read_wout_mod::vn_ai_aux_f
character(len= *), parameter vn_ai_aux_f
Definition read_wout_mod.f:32
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::potvac
real(rprec), dimension(:), allocatable, target potvac
Definition vacmod.f90:29
vacmod::p5
real(rprec), parameter p5
Definition vacmod.f90:13
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::bdotb
real(rprec), dimension(:), allocatable bdotb
Definition vmec_main.f90:59
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::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::cp5
real(rprec), parameter cp5
Definition vparams.f90:37
vparams::c2p0
real(rprec), parameter c2p0
Definition vparams.f90:40
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