 &nisland
 !
 ! input_JET  Example of namelist input for the test_island code
 !
 ! Two magnetic islands with mode numbers 3/2 and 2/1
 ! in a simple JET geometry
 !
 max_iterates = 50,   ! maximum number of iterations
 !
 num_profile_radial_grid = 52,  ! Number of radial profile grid points
 !
 leqtyp = 11,          ! Toroidal equilibrium with harmonic mixing
 !
 nxiequ = 31,         ! Number of equilibrium surfaces
 !
 nmodes = 2,          ! Number of magnetic islands
 modem  = 3,2,        ! Poloidal modes
 moden  = 2,1,      ! Toroidal modes
 !
 width_of_islands(1) = 0.06,0.08, ! Initial guesses for the island halfwidths
 !
 xoshift = 10*0.001,  ! Initial guesses for
                      ! relative shift between x- and o-points
 !
 nmomeq = 1,          ! Number of moment harmonics
 nhamad = 4,          ! Number of Hamada harmonics
 !
 ! Shape of plasma edge:
 edge_major_radius  = 2.93,    ! Edge major radius
 edge_minor_radius  = 0.838,    ! Edge minor radius (half-width)
 edge_elongation    = 1.695,    ! Edge elongation
 edge_triangularity = 0.265,    ! Edge triangularity
 !
 b_toroidal = 2.936,     ! Vacuum toroidal field edge major radius (tesla)
 q_axis  = 0.85,        ! q-value at the magnetic axis
 !
 pressure_at_axis = 6.5e4,    ! Plasma pressure at magnetic axis (Newtons)
 pressure_at_edge = 4.8e3,      ! Plasma pressure at edge
 pressure_exponential = 3.5,  ! exponential for pressure profile
 !
 expjt  = 2.0,         ! exponential for current density profile
 !
 wall_to_plasma_radius = 1.0,  ! r_wall / r_plasma
 !
 current_peaking_factor(1) = 10*0.0, ! current peaking factor inside islands
 !
 kjbs = 2,             ! Bootstrap current switch
 coeff_bootstrap_current = 1.0, ! coefficient of bootstrap current density
 !
 half_width_minimum = 1.0e-4,      ! Minimum allowed halfwidth
 !
 lislnd(8)  = 0,       ! pass bootstrap current density array to island
 lislnd(9)  = 5,       ! Runmode option (see documentation)
 lislnd(23) = 0,       ! Debug level for Hamada metric elements
 lislnd(27) = 1,       ! Number of times initial conditions computed per iteration
 ! lislnd(31) = 1,     ! Use Shampine and Gordon ODE solver
 lislnd(31) = 0,       ! Use DLSODE ODE solver
 lislnd(28) = 2        ! Choice of Delta Prime
 lislnd(32) = 12       ! Integration method for LSODE
 !
 cislnd(2)  = 0.08
 cislnd(7)  = 0.30
 cislnd(26) = 0.03
 ! cislnd(27) = 0.0
 cislnd(29) = 1.0e-12
 cislnd(30) = 1.0e-12, ! Absolute error in numerical integrator
 cislnd(31) = 1.0e-7, ! Relative error in numerical integrator
 cislnd(32) = 5.0e-5,  ! Relative error in zero finder
 !
 lprint = 3,           ! amount of diagnostic printout
 !
 pwidth_min  = 0.007999999, 0.0799999999 
 pwidth_max  = 0.020000000, 0.120000
 pwidth_step = 0.001, 0.002 
 !
 kjprf = 3,
 kpprf = 1,
 kqprf = 1,
 keq   = 1,
 pxjt =   -1.00E-02,1.00E-02,3.00E-02,5.00E-02,7.00E-02,9.00E-02,1.10E-01,1.30E-01,1.50E-01,1.70E-01,1.90E-01,2.10E-01,2.30E-01,2.50E-01,2.70E-01,2.90E-01,3.10E-01,3.30E-01,3.50E-01,3.70E-01,3.90E-01,4.10E-01,4.30E-01,4.50E-01,4.70E-01,4.90E-01,5.10E-01,5.30E-01,5.50E-01,5.70E-01,5.90E-01,6.10E-01,6.30E-01,6.50E-01,6.70E-01,6.90E-01,7.10E-01,7.30E-01,7.50E-01,7.70E-01,7.90E-01,8.10E-01,8.30E-01,8.50E-01,8.70E-01,8.90E-01,9.10E-01,9.30E-01,9.50E-01,9.70E-01,9.90E-01,1.01E+00
 !
 jtprf = 1.14384E+06,1.14384E+06,1.07453E+06,1.00198E+06,9.42312E+05,8.91313E+05,8.44906E+05,8.02352E+05,7.63368E+05,7.27914E+05,6.95749E+05,6.66933E+05,6.41099E+05,6.17735E+05,5.96317E+05,5.76484E+05,5.57922E+05,5.40077E+05,5.22376E+05,5.04782E+05,4.86472E+05,4.68073E+05,4.51188E+05,4.25488E+05,4.05852E+05,3.84938E+05,3.62947E+05,3.41112E+05,3.18897E+05,2.96458E+05,2.73554E+05,2.51508E+05,2.31961E+05,2.14544E+05,1.98604E+05,1.83885E+05,1.70287E+05,1.57818E+05,1.46534E+05,1.36473E+05,1.27562E+05,1.19687E+05,1.12575E+05,1.06064E+05,9.99526E+04,9.39491E+04,8.81338E+04,8.26506E+04,7.69791E+04,7.24146E+04,6.38119E+04,0.00000E+00
 !
 jrbphi = 9.35802E+00,9.35802E+00,9.35659E+00,9.35355E+00,9.34888E+00,9.34256E+00,9.33457E+00,9.32488E+00,9.31355E+00,9.30086E+00,9.28715E+00,9.27277E+00,9.25805E+00,9.24318E+00,9.22812E+00,9.21281E+00,9.19721E+00,9.18148E+00,9.16577E+00,9.15024E+00,9.13497E+00,9.12009E+00,9.10588E+00,9.09238E+00,9.07899E+00,9.06548E+00,9.05266E+00,9.04085E+00,9.02982E+00,9.01977E+00,9.01082E+00,9.00282E+00,8.99567E+00,8.98939E+00,8.98406E+00,8.97954E+00,8.97571E+00,8.97252E+00,8.96988E+00,8.96769E+00,8.96589E+00,8.96445E+00,8.96327E+00,8.96233E+00,8.96162E+00,8.96108E+00,8.96062E+00,8.96045E+00,8.95994E+00,8.96054E+00,8.95779E+00,8.95332E+00
 !
jt17r2 = 1.10684E-01,1.10684E-01,1.10700E-01,1.10734E-01,1.10784E-01,1.10852E-01,1.10936E-01,1.11038E-01,1.11157E-01,1.11293E-01,1.11447E-01,1.11617E-01,1.11804E-01,1.12007E-01,1.12223E-01,1.12455E-01,1.12702E-01,1.12967E-01,1.13250E-01,1.13550E-01,1.13871E-01,1.14217E-01,1.14593E-01,1.14991E-01,1.15406E-01,1.15833E-01,1.16275E-01,1.16749E-01,1.17257E-01,1.17797E-01,1.18368E-01,1.18974E-01,1.19621E-01,1.20308E-01,1.21031E-01,1.21799E-01,1.22630E-01,1.23508E-01,1.24415E-01,1.25395E-01,1.26442E-01,1.27512E-01,1.28657E-01,1.29891E-01,1.31140E-01,1.32479E-01,1.33917E-01,1.35355E-01,1.36928E-01,1.38559E-01,1.40187E-01,1.41878E-01
 !
 ppprf = 1.64157E+05,1.64157E+05,1.63236E+05,1.60375E+05,1.57111E+05,1.53792E+05,1.50274E+05,1.46765E+05,1.43279E+05,1.39855E+05,1.36499E+05,1.33253E+05,1.30087E+05,1.26989E+05,1.23956E+05,1.20982E+05,1.18068E+05,1.15163E+05,1.12252E+05,1.09364E+05,1.06319E+05,1.03179E+05,1.00393E+05,9.57250E+04,9.32492E+04,9.04539E+04,8.76796E+04,8.49444E+04,8.22207E+04,7.94750E+04,7.66149E+04,7.38009E+04,7.11319E+04,6.84654E+04,6.57658E+04,6.30536E+04,6.03359E+04,5.76326E+04,5.49720E+04,5.23775E+04,4.98241E+04,4.73382E+04,4.49127E+04,4.25322E+04,4.02328E+04,3.80113E+04,3.58812E+04,3.38586E+04,3.19428E+04,3.01955E+04,2.86695E+04,2.72139E+04
 !
 qxjt = -2.00000E-02,0.00000E+00,2.00000E-02,4.00000E-02,6.00000E-02,8.00000E-02,1.00000E-01,1.20000E-01,1.40000E-01,1.60000E-01,1.80000E-01,2.00000E-01,2.20000E-01,2.40000E-01,2.60000E-01,2.80000E-01,3.00000E-01,3.20000E-01,3.40000E-01,3.60000E-01,3.80000E-01,4.00000E-01,4.20000E-01,4.40000E-01,4.60000E-01,4.80000E-01,5.00000E-01,5.20000E-01,5.40000E-01,5.60000E-01,5.80000E-01,6.00000E-01,6.20000E-01,6.40000E-01,6.60000E-01,6.80000E-01,7.00000E-01,7.20000E-01,7.40000E-01,7.60000E-01,7.80000E-01,8.00000E-01,8.20000E-01,8.40000E-01,8.60000E-01,8.80000E-01,9.00000E-01,9.20000E-01,9.40000E-01,9.60000E-01,9.80000E-01,1.00000E+00
 !
 pqprf = 4.95562E-01,4.87708E-01,4.95562E-01,5.19122E-01,5.44238E-01,5.68402E-01,5.91780E-01,6.14968E-01,6.38224E-01,6.61645E-01,6.85223E-01,7.08931E-01,7.32659E-01,7.56297E-01,7.79778E-01,8.03101E-01,8.26342E-01,8.49601E-01,8.73016E-01,8.96775E-01,9.21045E-01,9.46114E-01,9.72203E-01,9.99245E-01,1.02863E+00,1.05957E+00,1.09245E+00,1.12769E+00,1.16558E+00,1.20646E+00,1.25064E+00,1.29855E+00,1.35056E+00,1.40671E+00,1.46686E+00,1.53105E+00,1.60001E+00,1.67403E+00,1.75265E+00,1.83653E+00,1.92692E+00,2.02294E+00,2.12451E+00,2.23456E+00,2.35175E+00,2.47568E+00,2.61128E+00,2.75558E+00,2.90993E+00,3.07962E+00,3.25820E+00,3.45635E+00
 !
 pxieq = 0.00000E+00,1.82574E-01,2.58199E-01,3.16228E-01,3.65148E-01,4.08248E-01,4.47214E-01,4.83046E-01,5.16398E-01,5.47723E-01,5.77350E-01,6.05530E-01,6.32456E-01,6.58281E-01,6.83130E-01,7.07107E-01,7.30297E-01,7.52773E-01,7.74597E-01,7.95822E-01,8.16497E-01,8.36660E-01,8.56349E-01,8.75595E-01,8.94427E-01,9.12871E-01,9.30949E-01,9.48683E-01,9.66092E-01,9.83192E-01,1.00000E+00
 !
 pr0 = 3.00652E+00,3.00242E+00,2.99838E+00,2.99442E+00,2.99047E+00,2.98652E+00,2.98255E+00,2.97844E+00,2.97432E+00,2.97017E+00,2.96591E+00,2.96158E+00,2.95717E+00,2.95259E+00,2.94792E+00,2.94316E+00,2.93816E+00,2.93307E+00,2.92788E+00,2.92242E+00,2.91687E+00,2.91123E+00,2.90531E+00,2.89929E+00,2.89318E+00,2.88680E+00,2.88033E+00,2.87377E+00,2.86702E+00,2.86020E+00,2.85329E+00
 !
 rcrmom = 0.00000E+00,1.76863E-01,2.48188E-01,3.02175E-01,3.47240E-01,3.86612E-01,4.21882E-01,4.53944E-01,4.83442E-01,5.10809E-01,5.36389E-01,5.60371E-01,5.82937E-01,6.04224E-01,6.24358E-01,6.43444E-01,6.61540E-01,6.78752E-01,6.95148E-01,7.10748E-01,7.25641E-01,7.39874E-01,7.53449E-01,7.66441E-01,7.78883E-01,7.90773E-01,8.02168E-01,8.13093E-01,8.23594E-01,8.33671E-01,8.43340E-01
 !
 py0 = 0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00,0.00000E+00
 !
 rsymom = 0.00000E+00,2.32934E-01,3.32228E-01,4.09498E-01,4.75264E-01,5.33616E-01,5.86736E-01,6.35957E-01,6.82146E-01,7.25875E-01,7.67521E-01,8.07468E-01,8.45971E-01,8.83245E-01,9.19447E-01,9.54704E-01,9.89170E-01,1.02290E+00,1.05596E+00,1.08849E+00,1.12050E+00,1.15202E+00,1.18318E+00,1.21395E+00,1.24438E+00,1.27454E+00,1.30444E+00,1.33408E+00,1.36348E+00,1.39267E+00,1.42170E+00
 /
