1 SUBROUTINE read_boozer(extension)
8 CHARACTER*(*) :: extension
12 REAL(rprec),
PARAMETER :: zero = 0, one = 1
21 CALL read_boozer_file(extension, istat)
22 IF (istat .ne. 0)
THEN
23 print *,
'Error reading boozer file in BOOTSJ, istat = ', istat
37 aipsi(1:irup) = buco_b(2:ns_b)
38 gpsi(1:irup) = bvco_b(2:ns_b)
39 qsafety(1:irup) = one/(iota_b(2:ns_b) +
40 1 sign(1.0e-14_dp,iota_b(2:ns_b)))
41 pres1(1:irup) = pres_b(2:ns_b)
42 betar(1:irup) = beta_b(2:ns_b)
43 idx(1:irup) = idx_b(2:ns_b)
44 flux(2:ns_b) = phi_b(2:ns_b)
45 phip(1:irup) = phip_b(2:ns_b)
49 IF( gpsi(irup)*phip_b(ns_b) <= zero) sign_jacobian = -one
52 END SUBROUTINE read_boozer