1 SUBROUTINE booz_rzhalf(r, z, rodd, zodd, r12, z12, ohs,
4 USE booz_params,
ONLY: nv_boz
9 INTEGER,
INTENT(in) :: js, nznt, nrep
11 REAL(rprec),
DIMENSION(nznt),
INTENT(in) :: r, z, rodd, zodd
12 REAL(rprec),
DIMENSION(nznt),
INTENT(out) :: r12, z12
16 REAL(rprec),
PARAMETER :: one = 1, c1p5 = 1.5_dp
17 LOGICAL,
PARAMETER :: lwrite=.false.
22 REAL(rprec) :: shalf, hs
25 shalf = sqrt(hs*abs(js-c1p5))
38 IF (lwrite)
CALL writesurface(js, nvplane, nrep, r12, z12)
40 END SUBROUTINE booz_rzhalf
42 SUBROUTINE writesurface(js, nvplane, nrep, r12, z12)
44 USE booz_params,
ONLY: nv_boz
45 USE booz_persistent,
ONLY: nu3_b
50 INTEGER :: js, nvplane, nrep
51 REAL(rprec),
DIMENSION(nv_boz, nu3_b),
INTENT(in) :: r12, z12
56 CHARACTER*(256) :: file_name
59 file_name =
"RZ_VMEC_VPLANE"
61 file_name =
"RZ_BOOZER_VPLANE"
64 WRITE (file_name,
'(a,i3.3,a,i3.3,a)') trim(file_name), nvplane,
67 OPEN(unit=33,file=file_name,status=
'REPLACE')
71 WRITE (33,100) i, r12(nvplane,i), z12(nvplane,i)
75 100
FORMAT(i5,1p,2e12.4)
79 END SUBROUTINE writesurface