 &nisland
 !
 ! input_DIIID  Example of namelist input for the test_island code
 !
 ! Two magnetic islands with mode numbers 3/2 and 2/1
 ! in a simple DIII-D 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  = 1.69,    ! Edge major radius
 edge_minor_radius  = 0.60,    ! Edge minor radius (half-width)
 edge_elongation    = 1.85,    ! Edge elongation
 edge_triangularity = 0.327,    ! Edge triangularity
 !
 b_toroidal = 2.0,     ! 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)  = 1
 cislnd(7)  = 0.30
 cislnd(26) = 0.005
 ! 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.73574E+06,1.73574E+06,1.65263E+06,1.56424E+06,1.48561E+06,1.40788E+06,1.32981E+06,1.25295E+06,1.17938E+06,1.10875E+06,1.04089E+06,9.75525E+05,9.12025E+05,8.50295E+05,7.89122E+05,7.23427E+05,6.65595E+05,6.11540E+05,5.59442E+05,5.08485E+05,4.67046E+05,4.10199E+05,3.59168E+05,3.33374E+05,3.14437E+05,2.96441E+05,2.79345E+05,2.63211E+05,2.47990E+05,2.33480E+05,2.20148E+05,2.07957E+05,1.96714E+05,1.86413E+05,1.76944E+05,1.68207E+05,1.60121E+05,1.52585E+05,1.45556E+05,1.38991E+05,1.32834E+05,1.26994E+05,1.21326E+05,1.16067E+05,1.11143E+05,1.06520E+05,1.02328E+05,9.85922E+04,9.50593E+04,9.25129E+04,8.83979E+04,0.00000E+00
 !
 jrbphi = 3.43046E+00,3.43046E+00,3.43011E+00,3.42939E+00,3.42832E+00,3.42689E+00,3.42511E+00,3.42296E+00,3.42046E+00,3.41760E+00,3.41439E+00,3.41083E+00,3.40692E+00,3.40273E+00,3.39844E+00,3.39423E+00,3.39028E+00,3.38661E+00,3.38325E+00,3.38022E+00,3.37752E+00,3.37515E+00,3.37309E+00,3.37132E+00,3.36980E+00,3.36852E+00,3.36744E+00,3.36653E+00,3.36575E+00,3.36509E+00,3.36453E+00,3.36406E+00,3.36367E+00,3.36334E+00,3.36307E+00,3.36285E+00,3.36266E+00,3.36251E+00,3.36239E+00,3.36228E+00,3.36220E+00,3.36213E+00,3.36207E+00,3.36202E+00,3.36198E+00,3.36195E+00,3.36192E+00,3.36190E+00,3.36184E+00,3.36192E+00,3.36152E+00,3.36089E+00
 !
jt17r2 = 3.33566E-01,3.33566E-01,3.33639E-01,3.33782E-01,3.33993E-01,3.34274E-01,3.34622E-01,3.35038E-01,3.35522E-01,3.36071E-01,3.36687E-01,3.37371E-01,3.38131E-01,3.38974E-01,3.39907E-01,3.40930E-01,3.42042E-01,3.43240E-01,3.44528E-01,3.45913E-01,3.47403E-01,3.49002E-01,3.50714E-01,3.52545E-01,3.54496E-01,3.56554E-01,3.58722E-01,3.61025E-01,3.63476E-01,3.66062E-01,3.68743E-01,3.71557E-01,3.74570E-01,3.77758E-01,3.81027E-01,3.84457E-01,3.88167E-01,3.92053E-01,3.95988E-01,4.00250E-01,4.04811E-01,4.09394E-01,4.14287E-01,4.19575E-01,4.24871E-01,4.30516E-01,4.36567E-01,4.42605E-01,4.49140E-01,4.55937E-01,4.62838E-01,4.70066E-01
 !
 ppprf = 6.46426E+04,6.46426E+04,6.38860E+04,6.23293E+04,6.04855E+04,5.83452E+04,5.62269E+04,5.41614E+04,5.22752E+04,5.04922E+04,4.87502E+04,4.70284E+04,4.52739E+04,4.35014E+04,4.16235E+04,3.97546E+04,3.83618E+04,3.68106E+04,3.53631E+04,3.39851E+04,3.26464E+04,3.00527E+04,2.77149E+04,2.64284E+04,2.51732E+04,2.39333E+04,2.27142E+04,2.15228E+04,2.03722E+04,1.92519E+04,1.81703E+04,1.71265E+04,1.61159E+04,1.51443E+04,1.42121E+04,1.33200E+04,1.24687E+04,1.16580E+04,1.08887E+04,1.01636E+04,9.48339E+03,8.84885E+03,8.25642E+03,7.70864E+03,7.20533E+03,6.74006E+03,6.31537E+03,5.92944E+03,5.57883E+03,5.27787E+03,5.01217E+03,4.76889E+03
 !
 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 = 6.19286E-01,6.11700E-01,6.19286E-01,6.42043E-01,6.66572E-01,6.91156E-01,7.16959E-01,7.44604E-01,7.74246E-01,8.05778E-01,8.39232E-01,8.74651E-01,9.12132E-01,9.51972E-01,9.94515E-01,1.04028E+00,1.09048E+00,1.14463E+00,1.20292E+00,1.26588E+00,1.33412E+00,1.40695E+00,1.48786E+00,1.57704E+00,1.67086E+00,1.76864E+00,1.87070E+00,1.97754E+00,2.08971E+00,2.20760E+00,2.33161E+00,2.46203E+00,2.59967E+00,2.74518E+00,2.89867E+00,3.06043E+00,3.23226E+00,3.41503E+00,3.60822E+00,3.81358E+00,4.03407E+00,4.26902E+00,4.51903E+00,4.78922E+00,5.07922E+00,5.38898E+00,5.72383E+00,6.08427E+00,6.47081E+00,6.88699E+00,7.33552E+00,7.82687E+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 = 1.73307E+00,1.72940E+00,1.72641E+00,1.72338E+00,1.72025E+00,1.71701E+00,1.71368E+00,1.71024E+00,1.70671E+00,1.70311E+00,1.69938E+00,1.69558E+00,1.69172E+00,1.68772E+00,1.68366E+00,1.67954E+00,1.67527E+00,1.67093E+00,1.66652E+00,1.66196E+00,1.65733E+00,1.65263E+00,1.64778E+00,1.64286E+00,1.63789E+00,1.63279E+00,1.62763E+00,1.62242E+00,1.61715E+00,1.61183E+00,1.60648E+00
 !
 rcrmom = 0.00000E+00,1.38060E-01,1.93795E-01,2.35595E-01,2.69956E-01,2.99472E-01,3.25481E-01,3.48669E-01,3.69683E-01,3.88902E-01,4.06654E-01,4.23107E-01,4.38427E-01,4.52791E-01,4.66265E-01,4.78936E-01,4.90915E-01,5.02228E-01,5.12925E-01,5.23092E-01,5.32729E-01,5.41869E-01,5.50589E-01,5.58872E-01,5.66744E-01,5.74282E-01,5.81455E-01,5.88287E-01,5.94874E-01,6.01159E-01,6.07161E-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,1.72627E-01,2.46324E-01,3.04065E-01,3.53715E-01,3.98239E-01,4.39166E-01,4.77561E-01,5.13865E-01,5.48481E-01,5.81627E-01,6.13588E-01,6.44537E-01,6.74543E-01,7.03781E-01,7.32347E-01,7.60261E-01,7.87651E-01,8.14581E-01,8.41036E-01,8.67130E-01,8.92907E-01,9.18321E-01,9.43485E-01,9.68427E-01,9.93066E-01,1.01752E+00,1.04181E+00,1.06579E+00,1.08962E+00,1.11331E+00
 /
