 ***************************************
 *  namelist produced by eqdsk2jorek   *
 ***************************************
 EHTR equilibrium eqdsk# dummy  3
 magentic field :   -2.107 T
 current        :   -1.987 MA
 ***************************************
 
  &in1
  restart = .t.
  rst_format = 0
  regrid  = .f.
  tstep   = 5. 
  nstep   = 2000
  nout    = 10 
  
  time_evol_scheme= 'Gears'

  tgnum(2) = 1.,
  tgnum(5) = 1.,
  tgnum(6) = 1.,
  tgnum(7) = 1.,

  gmres_4     = 1.d3
  iter_precon = 15
  gmres_tol   = 1.d-8
  gmres_m     = 20

  !_____________________________________boundary definition

  mf = 0
  n_boundary =          257
  R_boundary(  1) =  0.39800000E+01, Z_boundary(  1) =  0.10000000E+00, psi_boundary(  1) =  0.11891212E+01,
  R_boundary(  2) =  0.39796746E+01, Z_boundary(  2) =  0.14903337E+00, psi_boundary(  2) =  0.11825941E+01,
  R_boundary(  3) =  0.39786978E+01, Z_boundary(  3) =  0.19803721E+00, psi_boundary(  3) =  0.11770495E+01,
  R_boundary(  4) =  0.39770674E+01, Z_boundary(  4) =  0.24698200E+00, psi_boundary(  4) =  0.11724168E+01,
  R_boundary(  5) =  0.39747795E+01, Z_boundary(  5) =  0.29583825E+00, psi_boundary(  5) =  0.11686233E+01,
  R_boundary(  6) =  0.39718291E+01, Z_boundary(  6) =  0.34457653E+00, psi_boundary(  6) =  0.11655994E+01,
  R_boundary(  7) =  0.39682097E+01, Z_boundary(  7) =  0.39316749E+00, psi_boundary(  7) =  0.11632799E+01,
  R_boundary(  8) =  0.39639134E+01, Z_boundary(  8) =  0.44158185E+00, psi_boundary(  8) =  0.11616051E+01,
  R_boundary(  9) =  0.39589310E+01, Z_boundary(  9) =  0.48979046E+00, psi_boundary(  9) =  0.11605226E+01,
  R_boundary( 10) =  0.39532520E+01, Z_boundary( 10) =  0.53776428E+00, psi_boundary( 10) =  0.11599896E+01,
  R_boundary( 11) =  0.39468650E+01, Z_boundary( 11) =  0.58547440E+00, psi_boundary( 11) =  0.11599746E+01,
  R_boundary( 12) =  0.39397572E+01, Z_boundary( 12) =  0.63289209E+00, psi_boundary( 12) =  0.11604602E+01,
  R_boundary( 13) =  0.39319152E+01, Z_boundary( 13) =  0.67998879E+00, psi_boundary( 13) =  0.11614274E+01,
  R_boundary( 14) =  0.39233245E+01, Z_boundary( 14) =  0.72673612E+00, psi_boundary( 14) =  0.11628078E+01,
  R_boundary( 15) =  0.39139700E+01, Z_boundary( 15) =  0.77310593E+00, psi_boundary( 15) =  0.11645262E+01,
  R_boundary( 16) =  0.39038360E+01, Z_boundary( 16) =  0.81907028E+00, psi_boundary( 16) =  0.11665126E+01,
  R_boundary( 17) =  0.38929065E+01, Z_boundary( 17) =  0.86460150E+00, psi_boundary( 17) =  0.11686993E+01,
  R_boundary( 18) =  0.38811650E+01, Z_boundary( 18) =  0.90967215E+00, psi_boundary( 18) =  0.11710179E+01,
  R_boundary( 19) =  0.38685951E+01, Z_boundary( 19) =  0.95425508E+00, psi_boundary( 19) =  0.11733953E+01,
  R_boundary( 20) =  0.38551803E+01, Z_boundary( 20) =  0.99832344E+00, psi_boundary( 20) =  0.11757511E+01,
  R_boundary( 21) =  0.38409046E+01, Z_boundary( 21) =  0.10418507E+01, psi_boundary( 21) =  0.11779937E+01,
  R_boundary( 22) =  0.38257523E+01, Z_boundary( 22) =  0.10848106E+01, psi_boundary( 22) =  0.11800185E+01,
  R_boundary( 23) =  0.38097085E+01, Z_boundary( 23) =  0.11271773E+01, psi_boundary( 23) =  0.11817070E+01,
  R_boundary( 24) =  0.37927590E+01, Z_boundary( 24) =  0.11689252E+01, psi_boundary( 24) =  0.11829380E+01,
  R_boundary( 25) =  0.37748908E+01, Z_boundary( 25) =  0.12100293E+01, psi_boundary( 25) =  0.11836741E+01,
  R_boundary( 26) =  0.37560924E+01, Z_boundary( 26) =  0.12504648E+01, psi_boundary( 26) =  0.11839323E+01,
  R_boundary( 27) =  0.37363537E+01, Z_boundary( 27) =  0.12902072E+01, psi_boundary( 27) =  0.11837442E+01,
  R_boundary( 28) =  0.37156663E+01, Z_boundary( 28) =  0.13292327E+01, psi_boundary( 28) =  0.11831676E+01,
  R_boundary( 29) =  0.36940240E+01, Z_boundary( 29) =  0.13675178E+01, psi_boundary( 29) =  0.11822719E+01,
  R_boundary( 30) =  0.36714229E+01, Z_boundary( 30) =  0.14050393E+01, psi_boundary( 30) =  0.11811239E+01,
  R_boundary( 31) =  0.36478612E+01, Z_boundary( 31) =  0.14417748E+01, psi_boundary( 31) =  0.11797816E+01,
  R_boundary( 32) =  0.36233401E+01, Z_boundary( 32) =  0.14777020E+01, psi_boundary( 32) =  0.11782884E+01,
  R_boundary( 33) =  0.35978637E+01, Z_boundary( 33) =  0.15127993E+01, psi_boundary( 33) =  0.11766694E+01,
  R_boundary( 34) =  0.35714389E+01, Z_boundary( 34) =  0.15470457E+01, psi_boundary( 34) =  0.11749277E+01,
  R_boundary( 35) =  0.35440761E+01, Z_boundary( 35) =  0.15804203E+01, psi_boundary( 35) =  0.11730427E+01,
  R_boundary( 36) =  0.35157890E+01, Z_boundary( 36) =  0.16129033E+01, psi_boundary( 36) =  0.11709697E+01,
  R_boundary( 37) =  0.34865948E+01, Z_boundary( 37) =  0.16444749E+01, psi_boundary( 37) =  0.11686424E+01,
  R_boundary( 38) =  0.34565145E+01, Z_boundary( 38) =  0.16751162E+01, psi_boundary( 38) =  0.11659870E+01,
  R_boundary( 39) =  0.34255727E+01, Z_boundary( 39) =  0.17048086E+01, psi_boundary( 39) =  0.11629672E+01,
  R_boundary( 40) =  0.33937980E+01, Z_boundary( 40) =  0.17335345E+01, psi_boundary( 40) =  0.11595702E+01,
  R_boundary( 41) =  0.33612227E+01, Z_boundary( 41) =  0.17612763E+01, psi_boundary( 41) =  0.11558030E+01,
  R_boundary( 42) =  0.33278830E+01, Z_boundary( 42) =  0.17880174E+01, psi_boundary( 42) =  0.11516970E+01,
  R_boundary( 43) =  0.32938191E+01, Z_boundary( 43) =  0.18137418E+01, psi_boundary( 43) =  0.11473104E+01,
  R_boundary( 44) =  0.32590748E+01, Z_boundary( 44) =  0.18384338E+01, psi_boundary( 44) =  0.11427269E+01,
  R_boundary( 45) =  0.32236977E+01, Z_boundary( 45) =  0.18620787E+01, psi_boundary( 45) =  0.11380392E+01,
  R_boundary( 46) =  0.31877389E+01, Z_boundary( 46) =  0.18846622E+01, psi_boundary( 46) =  0.11333232E+01,
  R_boundary( 47) =  0.31512528E+01, Z_boundary( 47) =  0.19061706E+01, psi_boundary( 47) =  0.11286342E+01,
  R_boundary( 48) =  0.31142973E+01, Z_boundary( 48) =  0.19265911E+01, psi_boundary( 48) =  0.11240077E+01,
  R_boundary( 49) =  0.30769329E+01, Z_boundary( 49) =  0.19459113E+01, psi_boundary( 49) =  0.11194637E+01,
  R_boundary( 50) =  0.30392232E+01, Z_boundary( 50) =  0.19641196E+01, psi_boundary( 50) =  0.11150187E+01,
  R_boundary( 51) =  0.30012337E+01, Z_boundary( 51) =  0.19812050E+01, psi_boundary( 51) =  0.11106886E+01,
  R_boundary( 52) =  0.29630325E+01, Z_boundary( 52) =  0.19971573E+01, psi_boundary( 52) =  0.11064895E+01,
  R_boundary( 53) =  0.29246891E+01, Z_boundary( 53) =  0.20119668E+01, psi_boundary( 53) =  0.11024379E+01,
  R_boundary( 54) =  0.28862744E+01, Z_boundary( 54) =  0.20256246E+01, psi_boundary( 54) =  0.10985521E+01,
  R_boundary( 55) =  0.28478603E+01, Z_boundary( 55) =  0.20381224E+01, psi_boundary( 55) =  0.10948516E+01,
  R_boundary( 56) =  0.28095195E+01, Z_boundary( 56) =  0.20494529E+01, psi_boundary( 56) =  0.10913538E+01,
  R_boundary( 57) =  0.27713243E+01, Z_boundary( 57) =  0.20596090E+01, psi_boundary( 57) =  0.10880650E+01,
  R_boundary( 58) =  0.27333471E+01, Z_boundary( 58) =  0.20685847E+01, psi_boundary( 58) =  0.10849915E+01,
  R_boundary( 59) =  0.26956596E+01, Z_boundary( 59) =  0.20763747E+01, psi_boundary( 59) =  0.10821412E+01,
  R_boundary( 60) =  0.26583321E+01, Z_boundary( 60) =  0.20829741E+01, psi_boundary( 60) =  0.10795221E+01,
  R_boundary( 61) =  0.26214336E+01, Z_boundary( 61) =  0.20883791E+01, psi_boundary( 61) =  0.10771412E+01,
  R_boundary( 62) =  0.25850308E+01, Z_boundary( 62) =  0.20925863E+01, psi_boundary( 62) =  0.10750041E+01,
  R_boundary( 63) =  0.25491883E+01, Z_boundary( 63) =  0.20955933E+01, psi_boundary( 63) =  0.10731139E+01,
  R_boundary( 64) =  0.25139679E+01, Z_boundary( 64) =  0.20973982E+01, psi_boundary( 64) =  0.10714704E+01,
  R_boundary( 65) =  0.24794282E+01, Z_boundary( 65) =  0.20980000E+01, psi_boundary( 65) =  0.10700663E+01,
  R_boundary( 66) =  0.24456244E+01, Z_boundary( 66) =  0.20973982E+01, psi_boundary( 66) =  0.10688917E+01,
  R_boundary( 67) =  0.24126081E+01, Z_boundary( 67) =  0.20955933E+01, psi_boundary( 67) =  0.10679355E+01,
  R_boundary( 68) =  0.23804268E+01, Z_boundary( 68) =  0.20925863E+01, psi_boundary( 68) =  0.10671857E+01,
  R_boundary( 69) =  0.23491239E+01, Z_boundary( 69) =  0.20883791E+01, psi_boundary( 69) =  0.10666297E+01,
  R_boundary( 70) =  0.23187382E+01, Z_boundary( 70) =  0.20829741E+01, psi_boundary( 70) =  0.10662544E+01,
  R_boundary( 71) =  0.22893043E+01, Z_boundary( 71) =  0.20763747E+01, psi_boundary( 71) =  0.10660464E+01,
  R_boundary( 72) =  0.22608518E+01, Z_boundary( 72) =  0.20685847E+01, psi_boundary( 72) =  0.10659921E+01,
  R_boundary( 73) =  0.22334059E+01, Z_boundary( 73) =  0.20596090E+01, psi_boundary( 73) =  0.10660775E+01,
  R_boundary( 74) =  0.22069869E+01, Z_boundary( 74) =  0.20494529E+01, psi_boundary( 74) =  0.10662851E+01,
  R_boundary( 75) =  0.21816105E+01, Z_boundary( 75) =  0.20381224E+01, psi_boundary( 75) =  0.10665959E+01,
  R_boundary( 76) =  0.21572876E+01, Z_boundary( 76) =  0.20256246E+01, psi_boundary( 76) =  0.10669925E+01,
  R_boundary( 77) =  0.21340247E+01, Z_boundary( 77) =  0.20119668E+01, psi_boundary( 77) =  0.10674596E+01,
  R_boundary( 78) =  0.21118235E+01, Z_boundary( 78) =  0.19971573E+01, psi_boundary( 78) =  0.10679845E+01,
  R_boundary( 79) =  0.20906819E+01, Z_boundary( 79) =  0.19812050E+01, psi_boundary( 79) =  0.10685569E+01,
  R_boundary( 80) =  0.20705933E+01, Z_boundary( 80) =  0.19641196E+01, psi_boundary( 80) =  0.10691692E+01,
  R_boundary( 81) =  0.20515473E+01, Z_boundary( 81) =  0.19459113E+01, psi_boundary( 81) =  0.10698166E+01,
  R_boundary( 82) =  0.20335299E+01, Z_boundary( 82) =  0.19265911E+01, psi_boundary( 82) =  0.10704969E+01,
  R_boundary( 83) =  0.20165234E+01, Z_boundary( 83) =  0.19061706E+01, psi_boundary( 83) =  0.10712099E+01,
  R_boundary( 84) =  0.20005074E+01, Z_boundary( 84) =  0.18846622E+01, psi_boundary( 84) =  0.10719569E+01,
  R_boundary( 85) =  0.19854582E+01, Z_boundary( 85) =  0.18620787E+01, psi_boundary( 85) =  0.10727404E+01,
  R_boundary( 86) =  0.19713497E+01, Z_boundary( 86) =  0.18384338E+01, psi_boundary( 86) =  0.10735629E+01,
  R_boundary( 87) =  0.19581535E+01, Z_boundary( 87) =  0.18137418E+01, psi_boundary( 87) =  0.10744264E+01,
  R_boundary( 88) =  0.19458392E+01, Z_boundary( 88) =  0.17880174E+01, psi_boundary( 88) =  0.10753313E+01,
  R_boundary( 89) =  0.19343745E+01, Z_boundary( 89) =  0.17612763E+01, psi_boundary( 89) =  0.10762757E+01,
  R_boundary( 90) =  0.19237261E+01, Z_boundary( 90) =  0.17335345E+01, psi_boundary( 90) =  0.10772545E+01,
  R_boundary( 91) =  0.19138592E+01, Z_boundary( 91) =  0.17048086E+01, psi_boundary( 91) =  0.10782589E+01,
  R_boundary( 92) =  0.19047383E+01, Z_boundary( 92) =  0.16751162E+01, psi_boundary( 92) =  0.10792759E+01,
  R_boundary( 93) =  0.18963274E+01, Z_boundary( 93) =  0.16444749E+01, psi_boundary( 93) =  0.10802882E+01,
  R_boundary( 94) =  0.18885901E+01, Z_boundary( 94) =  0.16129033E+01, psi_boundary( 94) =  0.10812758E+01,
  R_boundary( 95) =  0.18814898E+01, Z_boundary( 95) =  0.15804203E+01, psi_boundary( 95) =  0.10822241E+01,
  R_boundary( 96) =  0.18749903E+01, Z_boundary( 96) =  0.15470457E+01, psi_boundary( 96) =  0.10831188E+01,
  R_boundary( 97) =  0.18690556E+01, Z_boundary( 97) =  0.15127993E+01, psi_boundary( 97) =  0.10839444E+01,
  R_boundary( 98) =  0.18636502E+01, Z_boundary( 98) =  0.14777020E+01, psi_boundary( 98) =  0.10846839E+01,
  R_boundary( 99) =  0.18587396E+01, Z_boundary( 99) =  0.14417748E+01, psi_boundary( 99) =  0.10853200E+01,
  R_boundary(100) =  0.18542897E+01, Z_boundary(100) =  0.14050393E+01, psi_boundary(100) =  0.10858344E+01,
  R_boundary(101) =  0.18502680E+01, Z_boundary(101) =  0.13675178E+01, psi_boundary(101) =  0.10862088E+01,
  R_boundary(102) =  0.18466425E+01, Z_boundary(102) =  0.13292327E+01, psi_boundary(102) =  0.10864245E+01,
  R_boundary(103) =  0.18433830E+01, Z_boundary(103) =  0.12902072E+01, psi_boundary(103) =  0.10864630E+01,
  R_boundary(104) =  0.18404602E+01, Z_boundary(104) =  0.12504648E+01, psi_boundary(104) =  0.10863061E+01,
  R_boundary(105) =  0.18378464E+01, Z_boundary(105) =  0.12100293E+01, psi_boundary(105) =  0.10859357E+01,
  R_boundary(106) =  0.18355151E+01, Z_boundary(106) =  0.11689252E+01, psi_boundary(106) =  0.10853343E+01,
  R_boundary(107) =  0.18334416E+01, Z_boundary(107) =  0.11271773E+01, psi_boundary(107) =  0.10844849E+01,
  R_boundary(108) =  0.18316024E+01, Z_boundary(108) =  0.10848106E+01, psi_boundary(108) =  0.10833837E+01,
  R_boundary(109) =  0.18299756E+01, Z_boundary(109) =  0.10418507E+01, psi_boundary(109) =  0.10820533E+01,
  R_boundary(110) =  0.18285406E+01, Z_boundary(110) =  0.99832344E+00, psi_boundary(110) =  0.10805230E+01,
  R_boundary(111) =  0.18272784E+01, Z_boundary(111) =  0.95425508E+00, psi_boundary(111) =  0.10788262E+01,
  R_boundary(112) =  0.18261715E+01, Z_boundary(112) =  0.90967215E+00, psi_boundary(112) =  0.10770000E+01,
  R_boundary(113) =  0.18252036E+01, Z_boundary(113) =  0.86460150E+00, psi_boundary(113) =  0.10750848E+01,
  R_boundary(114) =  0.18243599E+01, Z_boundary(114) =  0.81907028E+00, psi_boundary(114) =  0.10731248E+01,
  R_boundary(115) =  0.18236269E+01, Z_boundary(115) =  0.77310593E+00, psi_boundary(115) =  0.10711667E+01,
  R_boundary(116) =  0.18229922E+01, Z_boundary(116) =  0.72673612E+00, psi_boundary(116) =  0.10692604E+01,
  R_boundary(117) =  0.18224446E+01, Z_boundary(117) =  0.67998879E+00, psi_boundary(117) =  0.10674580E+01,
  R_boundary(118) =  0.18219743E+01, Z_boundary(118) =  0.63289209E+00, psi_boundary(118) =  0.10658139E+01,
  R_boundary(119) =  0.18215722E+01, Z_boundary(119) =  0.58547440E+00, psi_boundary(119) =  0.10643798E+01,
  R_boundary(120) =  0.18212306E+01, Z_boundary(120) =  0.53776428E+00, psi_boundary(120) =  0.10631903E+01,
  R_boundary(121) =  0.18209423E+01, Z_boundary(121) =  0.48979046E+00, psi_boundary(121) =  0.10622749E+01,
  R_boundary(122) =  0.18207014E+01, Z_boundary(122) =  0.44158185E+00, psi_boundary(122) =  0.10616633E+01,
  R_boundary(123) =  0.18205027E+01, Z_boundary(123) =  0.39316749E+00, psi_boundary(123) =  0.10613848E+01,
  R_boundary(124) =  0.18203418E+01, Z_boundary(124) =  0.34457653E+00, psi_boundary(124) =  0.10614685E+01,
  R_boundary(125) =  0.18202149E+01, Z_boundary(125) =  0.29583825E+00, psi_boundary(125) =  0.10619428E+01,
  R_boundary(126) =  0.18201192E+01, Z_boundary(126) =  0.24698200E+00, psi_boundary(126) =  0.10628354E+01,
  R_boundary(127) =  0.18200525E+01, Z_boundary(127) =  0.19803721E+00, psi_boundary(127) =  0.10641730E+01,
  R_boundary(128) =  0.18200130E+01, Z_boundary(128) =  0.14903337E+00, psi_boundary(128) =  0.10659811E+01,
  R_boundary(129) =  0.18200000E+01, Z_boundary(129) =  0.10000000E+00, psi_boundary(129) =  0.10682822E+01,
  R_boundary(130) =  0.18200130E+01, Z_boundary(130) =  0.50966625E-01, psi_boundary(130) =  0.10710588E+01,
  R_boundary(131) =  0.18200525E+01, Z_boundary(131) =  0.19627867E-02, psi_boundary(131) =  0.10742545E+01,
  R_boundary(132) =  0.18201192E+01, Z_boundary(132) = -0.46981998E-01, psi_boundary(132) =  0.10778098E+01,
  R_boundary(133) =  0.18202149E+01, Z_boundary(133) = -0.95838246E-01, psi_boundary(133) =  0.10816650E+01,
  R_boundary(134) =  0.18203418E+01, Z_boundary(134) = -0.14457653E+00, psi_boundary(134) =  0.10857594E+01,
  R_boundary(135) =  0.18205027E+01, Z_boundary(135) = -0.19316749E+00, psi_boundary(135) =  0.10900326E+01,
  R_boundary(136) =  0.18207014E+01, Z_boundary(136) = -0.24158185E+00, psi_boundary(136) =  0.10944241E+01,
  R_boundary(137) =  0.18209423E+01, Z_boundary(137) = -0.28979046E+00, psi_boundary(137) =  0.10988744E+01,
  R_boundary(138) =  0.18212306E+01, Z_boundary(138) = -0.33776428E+00, psi_boundary(138) =  0.11033246E+01,
  R_boundary(139) =  0.18215722E+01, Z_boundary(139) = -0.38547440E+00, psi_boundary(139) =  0.11077171E+01,
  R_boundary(140) =  0.18219743E+01, Z_boundary(140) = -0.43289209E+00, psi_boundary(140) =  0.11119961E+01,
  R_boundary(141) =  0.18224446E+01, Z_boundary(141) = -0.47998879E+00, psi_boundary(141) =  0.11161073E+01,
  R_boundary(142) =  0.18229922E+01, Z_boundary(142) = -0.52673612E+00, psi_boundary(142) =  0.11199989E+01,
  R_boundary(143) =  0.18236269E+01, Z_boundary(143) = -0.57310593E+00, psi_boundary(143) =  0.11236214E+01,
  R_boundary(144) =  0.18243599E+01, Z_boundary(144) = -0.61907028E+00, psi_boundary(144) =  0.11269280E+01,
  R_boundary(145) =  0.18252036E+01, Z_boundary(145) = -0.66460150E+00, psi_boundary(145) =  0.11298747E+01,
  R_boundary(146) =  0.18261715E+01, Z_boundary(146) = -0.70967215E+00, psi_boundary(146) =  0.11324208E+01,
  R_boundary(147) =  0.18272784E+01, Z_boundary(147) = -0.75425508E+00, psi_boundary(147) =  0.11345288E+01,
  R_boundary(148) =  0.18285406E+01, Z_boundary(148) = -0.79832344E+00, psi_boundary(148) =  0.11361645E+01,
  R_boundary(149) =  0.18299756E+01, Z_boundary(149) = -0.84185068E+00, psi_boundary(149) =  0.11372974E+01,
  R_boundary(150) =  0.18316024E+01, Z_boundary(150) = -0.88481059E+00, psi_boundary(150) =  0.11379006E+01,
  R_boundary(151) =  0.18334416E+01, Z_boundary(151) = -0.92717728E+00, psi_boundary(151) =  0.11379554E+01,
  R_boundary(152) =  0.18355151E+01, Z_boundary(152) = -0.96892524E+00, psi_boundary(152) =  0.11374681E+01,
  R_boundary(153) =  0.18378464E+01, Z_boundary(153) = -0.10100293E+01, psi_boundary(153) =  0.11364522E+01,
  R_boundary(154) =  0.18404602E+01, Z_boundary(154) = -0.10504648E+01, psi_boundary(154) =  0.11349208E+01,
  R_boundary(155) =  0.18433830E+01, Z_boundary(155) = -0.10902072E+01, psi_boundary(155) =  0.11328868E+01,
  R_boundary(156) =  0.18466425E+01, Z_boundary(156) = -0.11292327E+01, psi_boundary(156) =  0.11303622E+01,
  R_boundary(157) =  0.18502680E+01, Z_boundary(157) = -0.11675178E+01, psi_boundary(157) =  0.11273575E+01,
  R_boundary(158) =  0.18542897E+01, Z_boundary(158) = -0.12050393E+01, psi_boundary(158) =  0.11238820E+01,
  R_boundary(159) =  0.18587396E+01, Z_boundary(159) = -0.12417748E+01, psi_boundary(159) =  0.11199427E+01,
  R_boundary(160) =  0.18636502E+01, Z_boundary(160) = -0.12777020E+01, psi_boundary(160) =  0.11155447E+01,
  R_boundary(161) =  0.18690556E+01, Z_boundary(161) = -0.13127993E+01, psi_boundary(161) =  0.11106900E+01,
  R_boundary(162) =  0.18749903E+01, Z_boundary(162) = -0.13470457E+01, psi_boundary(162) =  0.11053782E+01,
  R_boundary(163) =  0.18814898E+01, Z_boundary(163) = -0.13804203E+01, psi_boundary(163) =  0.10996057E+01,
  R_boundary(164) =  0.18885901E+01, Z_boundary(164) = -0.14129033E+01, psi_boundary(164) =  0.10933706E+01,
  R_boundary(165) =  0.18963274E+01, Z_boundary(165) = -0.14444749E+01, psi_boundary(165) =  0.10867337E+01,
  R_boundary(166) =  0.19047383E+01, Z_boundary(166) = -0.14751162E+01, psi_boundary(166) =  0.10798024E+01,
  R_boundary(167) =  0.19138592E+01, Z_boundary(167) = -0.15048086E+01, psi_boundary(167) =  0.10726801E+01,
  R_boundary(168) =  0.19237261E+01, Z_boundary(168) = -0.15335345E+01, psi_boundary(168) =  0.10654659E+01,
  R_boundary(169) =  0.19343745E+01, Z_boundary(169) = -0.15612763E+01, psi_boundary(169) =  0.10582553E+01,
  R_boundary(170) =  0.19458392E+01, Z_boundary(170) = -0.15880174E+01, psi_boundary(170) =  0.10511391E+01,
  R_boundary(171) =  0.19581535E+01, Z_boundary(171) = -0.16137418E+01, psi_boundary(171) =  0.10442028E+01,
  R_boundary(172) =  0.19713497E+01, Z_boundary(172) = -0.16384338E+01, psi_boundary(172) =  0.10375149E+01,
  R_boundary(173) =  0.19854582E+01, Z_boundary(173) = -0.16620787E+01, psi_boundary(173) =  0.10311175E+01,
  R_boundary(174) =  0.20005074E+01, Z_boundary(174) = -0.16846622E+01, psi_boundary(174) =  0.10249832E+01,
  R_boundary(175) =  0.20165234E+01, Z_boundary(175) = -0.17061706E+01, psi_boundary(175) =  0.10190468E+01,
  R_boundary(176) =  0.20335299E+01, Z_boundary(176) = -0.17265911E+01, psi_boundary(176) =  0.10132297E+01,
  R_boundary(177) =  0.20515473E+01, Z_boundary(177) = -0.17459113E+01, psi_boundary(177) =  0.10074210E+01,
  R_boundary(178) =  0.20705933E+01, Z_boundary(178) = -0.17641196E+01, psi_boundary(178) =  0.10014640E+01,
  R_boundary(179) =  0.20906819E+01, Z_boundary(179) = -0.17812050E+01, psi_boundary(179) =  0.99514940E+00,
  R_boundary(180) =  0.21118235E+01, Z_boundary(180) = -0.17971573E+01, psi_boundary(180) =  0.98822953E+00,
  R_boundary(181) =  0.21340247E+01, Z_boundary(181) = -0.18119668E+01, psi_boundary(181) =  0.98049549E+00,
  R_boundary(182) =  0.21572876E+01, Z_boundary(182) = -0.18256246E+01, psi_boundary(182) =  0.97182339E+00,
  R_boundary(183) =  0.21816105E+01, Z_boundary(183) = -0.18381224E+01, psi_boundary(183) =  0.96217540E+00,
  R_boundary(184) =  0.22069869E+01, Z_boundary(184) = -0.18494529E+01, psi_boundary(184) =  0.95158115E+00,
  R_boundary(185) =  0.22334059E+01, Z_boundary(185) = -0.18596090E+01, psi_boundary(185) =  0.94010572E+00,
  R_boundary(186) =  0.22608518E+01, Z_boundary(186) = -0.18685847E+01, psi_boundary(186) =  0.92781308E+00,
  R_boundary(187) =  0.22893043E+01, Z_boundary(187) = -0.18763747E+01, psi_boundary(187) =  0.91481351E+00,
  R_boundary(188) =  0.23187382E+01, Z_boundary(188) = -0.18829741E+01, psi_boundary(188) =  0.90130820E+00,
  R_boundary(189) =  0.23491239E+01, Z_boundary(189) = -0.18883791E+01, psi_boundary(189) =  0.88750653E+00,
  R_boundary(190) =  0.23804268E+01, Z_boundary(190) = -0.18925863E+01, psi_boundary(190) =  0.87362195E+00,
  R_boundary(191) =  0.24126081E+01, Z_boundary(191) = -0.18955933E+01, psi_boundary(191) =  0.85988475E+00,
  R_boundary(192) =  0.24456244E+01, Z_boundary(192) = -0.18973982E+01, psi_boundary(192) =  0.84656446E+00,
  R_boundary(193) =  0.24794282E+01, Z_boundary(193) = -0.18980000E+01, psi_boundary(193) =  0.83399867E+00,
  R_boundary(194) =  0.25139679E+01, Z_boundary(194) = -0.18973982E+01, psi_boundary(194) =  0.82262381E+00,
  R_boundary(195) =  0.25491883E+01, Z_boundary(195) = -0.18955933E+01, psi_boundary(195) =  0.81300150E+00,
  R_boundary(196) =  0.25850308E+01, Z_boundary(196) = -0.18925863E+01, psi_boundary(196) =  0.80580296E+00,
  R_boundary(197) =  0.26214336E+01, Z_boundary(197) = -0.18883791E+01, psi_boundary(197) =  0.80170478E+00,
  R_boundary(198) =  0.26583321E+01, Z_boundary(198) = -0.18829741E+01, psi_boundary(198) =  0.80146284E+00,
  R_boundary(199) =  0.26956596E+01, Z_boundary(199) = -0.18763747E+01, psi_boundary(199) =  0.80590194E+00,
  R_boundary(200) =  0.27333471E+01, Z_boundary(200) = -0.18685847E+01, psi_boundary(200) =  0.81583688E+00,
  R_boundary(201) =  0.27713243E+01, Z_boundary(201) = -0.18596090E+01, psi_boundary(201) =  0.83196101E+00,
  R_boundary(202) =  0.28095195E+01, Z_boundary(202) = -0.18494529E+01, psi_boundary(202) =  0.85471116E+00,
  R_boundary(203) =  0.28478603E+01, Z_boundary(203) = -0.18381224E+01, psi_boundary(203) =  0.88412538E+00,
  R_boundary(204) =  0.28862744E+01, Z_boundary(204) = -0.18256246E+01, psi_boundary(204) =  0.91958598E+00,
  R_boundary(205) =  0.29246891E+01, Z_boundary(205) = -0.18119668E+01, psi_boundary(205) =  0.95935581E+00,
  R_boundary(206) =  0.29630325E+01, Z_boundary(206) = -0.17971573E+01, psi_boundary(206) =  0.10012883E+01,
  R_boundary(207) =  0.30012337E+01, Z_boundary(207) = -0.17812050E+01, psi_boundary(207) =  0.10433153E+01,
  R_boundary(208) =  0.30392232E+01, Z_boundary(208) = -0.17641196E+01, psi_boundary(208) =  0.10839324E+01,
  R_boundary(209) =  0.30769329E+01, Z_boundary(209) = -0.17459113E+01, psi_boundary(209) =  0.11230450E+01,
  R_boundary(210) =  0.31142973E+01, Z_boundary(210) = -0.17265911E+01, psi_boundary(210) =  0.11618565E+01,
  R_boundary(211) =  0.31512528E+01, Z_boundary(211) = -0.17061706E+01, psi_boundary(211) =  0.12022437E+01,
  R_boundary(212) =  0.31877389E+01, Z_boundary(212) = -0.16846622E+01, psi_boundary(212) =  0.12459184E+01,
  R_boundary(213) =  0.32236977E+01, Z_boundary(213) = -0.16620787E+01, psi_boundary(213) =  0.12923804E+01,
  R_boundary(214) =  0.32590748E+01, Z_boundary(214) = -0.16384338E+01, psi_boundary(214) =  0.13387359E+01,
  R_boundary(215) =  0.32938191E+01, Z_boundary(215) = -0.16137418E+01, psi_boundary(215) =  0.13806562E+01,
  R_boundary(216) =  0.33278830E+01, Z_boundary(216) = -0.15880174E+01, psi_boundary(216) =  0.14146184E+01,
  R_boundary(217) =  0.33612227E+01, Z_boundary(217) = -0.15612763E+01, psi_boundary(217) =  0.14391572E+01,
  R_boundary(218) =  0.33937980E+01, Z_boundary(218) = -0.15335345E+01, psi_boundary(218) =  0.14541033E+01,
  R_boundary(219) =  0.34255727E+01, Z_boundary(219) = -0.15048086E+01, psi_boundary(219) =  0.14601218E+01,
  R_boundary(220) =  0.34565145E+01, Z_boundary(220) = -0.14751162E+01, psi_boundary(220) =  0.14587144E+01,
  R_boundary(221) =  0.34865948E+01, Z_boundary(221) = -0.14444749E+01, psi_boundary(221) =  0.14522229E+01,
  R_boundary(222) =  0.35157890E+01, Z_boundary(222) = -0.14129033E+01, psi_boundary(222) =  0.14435255E+01,
  R_boundary(223) =  0.35440761E+01, Z_boundary(223) = -0.13804203E+01, psi_boundary(223) =  0.14351687E+01,
  R_boundary(224) =  0.35714389E+01, Z_boundary(224) = -0.13470457E+01, psi_boundary(224) =  0.14285154E+01,
  R_boundary(225) =  0.35978637E+01, Z_boundary(225) = -0.13127993E+01, psi_boundary(225) =  0.14237047E+01,
  R_boundary(226) =  0.36233401E+01, Z_boundary(226) = -0.12777020E+01, psi_boundary(226) =  0.14205263E+01,
  R_boundary(227) =  0.36478612E+01, Z_boundary(227) = -0.12417748E+01, psi_boundary(227) =  0.14185916E+01,
  R_boundary(228) =  0.36714229E+01, Z_boundary(228) = -0.12050393E+01, psi_boundary(228) =  0.14174328E+01,
  R_boundary(229) =  0.36940240E+01, Z_boundary(229) = -0.11675178E+01, psi_boundary(229) =  0.14165791E+01,
  R_boundary(230) =  0.37156663E+01, Z_boundary(230) = -0.11292327E+01, psi_boundary(230) =  0.14156097E+01,
  R_boundary(231) =  0.37363537E+01, Z_boundary(231) = -0.10902072E+01, psi_boundary(231) =  0.14141853E+01,
  R_boundary(232) =  0.37560924E+01, Z_boundary(232) = -0.10504648E+01, psi_boundary(232) =  0.14120617E+01,
  R_boundary(233) =  0.37748908E+01, Z_boundary(233) = -0.10100293E+01, psi_boundary(233) =  0.14090935E+01,
  R_boundary(234) =  0.37927590E+01, Z_boundary(234) = -0.96892524E+00, psi_boundary(234) =  0.14052209E+01,
  R_boundary(235) =  0.38097085E+01, Z_boundary(235) = -0.92717728E+00, psi_boundary(235) =  0.14004290E+01,
  R_boundary(236) =  0.38257523E+01, Z_boundary(236) = -0.88481059E+00, psi_boundary(236) =  0.13947219E+01,
  R_boundary(237) =  0.38409046E+01, Z_boundary(237) = -0.84185068E+00, psi_boundary(237) =  0.13881248E+01,
  R_boundary(238) =  0.38551803E+01, Z_boundary(238) = -0.79832344E+00, psi_boundary(238) =  0.13807164E+01,
  R_boundary(239) =  0.38685951E+01, Z_boundary(239) = -0.75425508E+00, psi_boundary(239) =  0.13725783E+01,
  R_boundary(240) =  0.38811650E+01, Z_boundary(240) = -0.70967215E+00, psi_boundary(240) =  0.13637852E+01,
  R_boundary(241) =  0.38929065E+01, Z_boundary(241) = -0.66460150E+00, psi_boundary(241) =  0.13544079E+01,
  R_boundary(242) =  0.39038360E+01, Z_boundary(242) = -0.61907028E+00, psi_boundary(242) =  0.13445157E+01,
  R_boundary(243) =  0.39139700E+01, Z_boundary(243) = -0.57310593E+00, psi_boundary(243) =  0.13341786E+01,
  R_boundary(244) =  0.39233245E+01, Z_boundary(244) = -0.52673612E+00, psi_boundary(244) =  0.13234681E+01,
  R_boundary(245) =  0.39319152E+01, Z_boundary(245) = -0.47998879E+00, psi_boundary(245) =  0.13124588E+01,
  R_boundary(246) =  0.39397572E+01, Z_boundary(246) = -0.43289209E+00, psi_boundary(246) =  0.13012287E+01,
  R_boundary(247) =  0.39468650E+01, Z_boundary(247) = -0.38547440E+00, psi_boundary(247) =  0.12898596E+01,
  R_boundary(248) =  0.39532520E+01, Z_boundary(248) = -0.33776428E+00, psi_boundary(248) =  0.12784373E+01,
  R_boundary(249) =  0.39589310E+01, Z_boundary(249) = -0.28979046E+00, psi_boundary(249) =  0.12670517E+01,
  R_boundary(250) =  0.39639134E+01, Z_boundary(250) = -0.24158185E+00, psi_boundary(250) =  0.12557963E+01,
  R_boundary(251) =  0.39682097E+01, Z_boundary(251) = -0.19316749E+00, psi_boundary(251) =  0.12447679E+01,
  R_boundary(252) =  0.39718291E+01, Z_boundary(252) = -0.14457653E+00, psi_boundary(252) =  0.12340665E+01,
  R_boundary(253) =  0.39747795E+01, Z_boundary(253) = -0.95838246E-01, psi_boundary(253) =  0.12237946E+01,
  R_boundary(254) =  0.39770674E+01, Z_boundary(254) = -0.46981998E-01, psi_boundary(254) =  0.12140571E+01,
  R_boundary(255) =  0.39786978E+01, Z_boundary(255) =  0.19627867E-02, psi_boundary(255) =  0.12049602E+01,
  R_boundary(256) =  0.39796746E+01, Z_boundary(256) =  0.50966625E-01, psi_boundary(256) =  0.11966119E+01,
  R_boundary(257) =  0.39800000E+01, Z_boundary(257) =  0.10000000E+00, psi_boundary(257) =  0.11891212E+01,

  ellip  =    1.85000000000000     
  tria_u =   0.400000000000000     
  tria_l =   0.400000000000000     
  quad_u =  -0.200000000000000     
  quad_l =  -0.200000000000000     
  xampl  = +0.
  xpoint = .t.
  freeboundary = .f.
  resistive_wall = .f.
  
  R_geo =   2.90000000000000     
  Z_geo =   0.100000000000000     
  F0    =   6.10957500000000     
  amin  = 1.d0 ! scale factor for plasma size only
 
  !_____________________________________grid parameters
  n_R      = 0
  n_Z      = 0
  n_radial = 101
  n_pol    = 151
  n_flux   = 101
  n_tht    = 151
  n_open   = 21
  n_leg    = 25
  n_private = 11

  dPSI_open    = 0.06
  dPSI_private = 0.02

  !_____________________________________physics parameters
  gamma_sheath = 8.d0
  bc_natural_open = .t
  eta   = 1.d-7
  visco = 1.d-6
  visco_par = 1.d-5

  eta_num       = 1.d-12
  visco_num     = 1.d-11
  visco_par_num = 1.d-9
  D_perp_num    = 1.d-11
  Zk_perp_num   = 1.d-10

  D_prof_neg  = 1.d-6
  ZK_prof_neg = 1.d-6
  T_min       = 1.e-5

!  rho_file     = "jorek_density"
!  T_file       = "jorek_temperature"
!  ffprime_file = "jorek_ffprime"

  FF_0 =  1.8
  FF_1 =  0.
  FF_coef(1) =  -0.0183437991116
  FF_coef(2) =  -1.15102946172
  FF_coef(3) =  -0.0105832518754
  FF_coef(4) =  0.0260326109348
  FF_coef(5) =  0.952047027421
  FF_coef(6) = -0.1 !-0.0187117979807 ! * delta_psi
  FF_coef(7) =  0.95 !0.966206603386 
  FF_coef(8) =  0.015! 0.00792706422311

  T_0       =  0.0702
  T_1       =  0.0005
  T_coef(1) =  -1.0
  T_coef(2) =  0.12
  T_coef(3) =  0.0
  T_coef(4) =  0.025
  T_coef(5) =  0.95

  rho_0       =  0.946681184152
  rho_1       =  0.0216272686248
  rho_coef(1) =  0.197554107262
  rho_coef(2) =  -0.95 !-1.05601347437
  rho_coef(3) =  0.634391814096
  rho_coef(4) =  0.0216272686248
  rho_coef(5) =  0.959122286073

  D_par     = 0.d0
  D_perp(1) = 1.d-6
  D_perp(2) = 0.9d0
  D_perp(3) = 0.10d0
  D_perp(4) = 0.01d0
  D_perp(5) = 0.92d0
  D_perp(6) = 2.d-6

  ZK_par     = 8000.d0
  ZK_par_max = 1600.d0
  ZK_perp(1) = 1.d-6
  ZK_perp(2) = 0.9d0
  ZK_perp(3) = 0.10d0
  ZK_perp(4) = 0.01d0
  ZK_perp(5) = 0.92d0
  ZK_perp(6) = 5.d-6

  pellet_R      = 3.78 
  pellet_Z      = 0.002
  pellet_phi    = 0.45
  pellet_radius = 0.01
  pellet_sig    = 0.01
  pellet_length = 0.2
  pellet_ellipse = 5.
  pellet_theta  = 0.0

  pellet_velocity_R = -77.5d0      ! [m/s]
  pellet_velocity_Z = 0.d0
  pellet_particles  = 3.2d0        ! [x10^20 particles]
  use_pellet        = .f.
  pellet_density    = 4.95d8       ! [x10^20 m^-3]

  heatsource     = 1.d-7
  particlesource = 5.d-6
  central_density = 0.65
  &end
