diff --git a/pyexamples/d3dHsmNew/README b/pyexamples/d3dHsmNew/README new file mode 100644 index 00000000..e35c896d --- /dev/null +++ b/pyexamples/d3dHsmNew/README @@ -0,0 +1,4 @@ +python -i runcase.py + +Or, in Python do: +exec(open("runcase.py").read()) diff --git a/pyexamples/d3dHsmNew/aeqdsk b/pyexamples/d3dHsmNew/aeqdsk new file mode 100644 index 00000000..86ded5a2 Binary files /dev/null and b/pyexamples/d3dHsmNew/aeqdsk differ diff --git a/pyexamples/d3dHsmNew/d3dHsm.h5 b/pyexamples/d3dHsmNew/d3dHsm.h5 new file mode 100644 index 00000000..653f63ec Binary files /dev/null and b/pyexamples/d3dHsmNew/d3dHsm.h5 differ diff --git a/pyexamples/d3dHsmNew/neqdsk b/pyexamples/d3dHsmNew/neqdsk new file mode 100644 index 00000000..354940d8 --- /dev/null +++ b/pyexamples/d3dHsmNew/neqdsk @@ -0,0 +1,669 @@ + EFITD 09/07/90 # 66832 ,2384ms 3 33 65 + 0.170000005E+01 0.320000005E+01 0.169550002E+01 0.839999974E+00 0.000000000E+00 + 0.172453952E+01 0.910593718E-01-0.380751818E+00 0.164451320E-01-0.206767821E+01 + 0.120704213E+07-0.380751818E+00 0.000000000E+00 0.172453952E+01 0.000000000E+00 + 0.910593718E-01 0.000000000E+00 0.164451320E-01 0.000000000E+00 0.000000000E+00 +-0.376038146E+01-0.373506474E+01-0.371134663E+01-0.368919587E+01-0.366857719E+01 +-0.364945292E+01-0.363178301E+01-0.361552477E+01-0.360063291E+01-0.358705997E+01 +-0.357475591E+01-0.356366873E+01-0.355374527E+01-0.354492903E+01-0.353716350E+01 +-0.353038955E+01-0.352454734E+01-0.351957560E+01-0.351541209E+01-0.351199365E+01 +-0.350925636E+01-0.350713611E+01-0.350556731E+01-0.350448513E+01-0.350382376E+01 +-0.350351667E+01-0.350349855E+01-0.350370312E+01-0.350406408E+01-0.350451517E+01 +-0.350499034E+01-0.350542355E+01-0.350574851E+01 +-0.434976123E+04-0.182651587E+04 0.542709229E+03 0.275791406E+04 0.481909863E+04 + 0.672626270E+04 0.847940625E+04 0.100785293E+05 0.115236318E+05 0.128147148E+05 + 0.139517773E+05 0.149348193E+05 0.157638408E+05 0.164388418E+05 0.169598223E+05 + 0.173267832E+05 0.175397227E+05 0.175986426E+05 0.175035430E+05 0.172544219E+05 + 0.168512813E+05 0.162941191E+05 0.155829375E+05 0.147177354E+05 0.136985127E+05 + 0.125252695E+05 0.111980068E+05 0.971672363E+04 0.808141943E+04 0.629209521E+04 + 0.434875049E+04 0.225138525E+04 0.000000000E+00 +-0.791346884E+01-0.737451935E+01-0.685424232E+01-0.635263681E+01-0.586970234E+01 +-0.540543985E+01-0.495984745E+01-0.453292751E+01-0.412467861E+01-0.373510170E+01 +-0.336419606E+01-0.301196194E+01-0.267839861E+01-0.236350751E+01-0.206728768E+01 +-0.178973937E+01-0.153086221E+01-0.129065645E+01-0.106912267E+01-0.866259813E+00 +-0.682069004E+00-0.516548812E+00-0.369700968E+00-0.241523921E+00-0.132018119E+00 +-0.411849171E-01 0.309779700E-01 0.844686925E-01 0.119288653E+00 0.135435984E+00 + 0.132912546E+00 0.111718789E+00 0.718519539E-01 + 0.209488453E+06 0.197079875E+06 0.184671297E+06 0.172262719E+06 0.159854141E+06 + 0.147445547E+06 0.135036969E+06 0.122628398E+06 0.110219805E+06 0.978112266E+05 + 0.854026406E+05 0.729940625E+05 0.605854766E+05 0.481768984E+05 0.357683203E+05 + 0.233597324E+05 0.109511543E+05-0.145742261E+04-0.138660000E+05-0.262745996E+05 +-0.386831758E+05-0.510917539E+05-0.635003320E+05-0.759089063E+05-0.883175078E+05 +-0.100726086E+06-0.113134664E+06-0.125543242E+06-0.137951813E+06-0.150360391E+06 +-0.162769000E+06-0.175177578E+06-0.187586156E+06 + 0.721479878E-01 0.624937564E-01 0.511088483E-01 0.383830927E-01 0.243493877E-01 + 0.176306255E-01 0.263733789E-02-0.535020605E-02-0.142600462E-01-0.170640498E-01 +-0.114340410E-01-0.876979530E-02-0.116666108E-01-0.195798501E-01-0.345983803E-01 +-0.567315072E-01-0.600479245E-01-0.446733832E-01-0.804638863E-02 0.217667222E-01 + 0.455679297E-01 0.661828071E-01 0.849517584E-01 0.102494366E+00 0.119042829E+00 + 0.134606332E+00 0.149060979E+00 0.162210926E+00 0.173846364E+00 0.183798358E+00 + 0.191983402E+00 0.198422328E+00 0.203235507E+00 0.700871497E-01 0.590926744E-01 + 0.459876247E-01 0.306905545E-01 0.215200409E-01 0.161816850E-01 0.618049502E-02 +-0.171422213E-02-0.809847564E-02-0.110300183E-01-0.975672901E-02-0.923087448E-02 +-0.125951618E-01-0.204607174E-01-0.350873843E-01-0.570527166E-01-0.600595474E-01 +-0.445291996E-01-0.730913877E-02 0.230438411E-01 0.478390008E-01 0.696861222E-01 + 0.898093134E-01 0.108786933E+00 0.126826316E+00 0.143893346E+00 0.159788683E+00 + 0.174213290E+00 0.186846688E+00 0.197431520E+00 0.205845386E+00 0.212131843E+00 + 0.216483817E+00 0.688181370E-01 0.571269281E-01 0.420824103E-01 0.300944224E-01 + 0.238809213E-01 0.189776458E-01 0.123297684E-01 0.579196960E-02 0.604785979E-03 +-0.277537853E-02-0.456120819E-02-0.634963810E-02-0.991243869E-02-0.161561593E-01 +-0.257470980E-01-0.360031873E-01-0.355222821E-01-0.220062584E-01 0.358977914E-02 + 0.294892043E-01 0.531211495E-01 0.751331672E-01 0.960884243E-01 0.116299786E+00 + 0.135839671E+00 0.154563874E+00 0.172138885E+00 0.188093260E+00 0.201917335E+00 + 0.213200182E+00 0.221754372E+00 0.227667257E+00 0.231258318E+00 0.684757158E-01 + 0.571621954E-01 0.434899107E-01 0.346383378E-01 0.290759504E-01 0.235628448E-01 + 0.177345909E-01 0.121091530E-01 0.730229169E-02 0.358518958E-02 0.764682889E-03 +-0.184364617E-02-0.500249863E-02-0.904437155E-02-0.134546012E-01-0.157886669E-01 +-0.121455640E-01-0.650554895E-03 0.175053626E-01 0.386574864E-01 0.603689775E-01 + 0.819932818E-01 0.103474438E+00 0.124831364E+00 0.145975843E+00 0.166629016E+00 + 0.186276823E+00 0.204185113E+00 0.219527468E+00 0.231622711E+00 0.240174800E+00 + 0.245363995E+00 0.247743115E+00 0.687776282E-01 0.586588010E-01 0.482516140E-01 + 0.404306948E-01 0.341601074E-01 0.283303820E-01 0.227091499E-01 0.174457729E-01 + 0.128099024E-01 0.895442814E-02 0.581165403E-02 0.312928110E-02 0.672996044E-03 +-0.148733705E-02-0.268607587E-02-0.146185607E-02 0.395810604E-02 0.145211965E-01 + 0.297434479E-01 0.481281057E-01 0.682940036E-01 0.895029232E-01 0.111456260E+00 + 0.134010330E+00 0.156980976E+00 0.179993287E+00 0.202342734E+00 0.222918943E+00 + 0.240351737E+00 0.253460288E+00 0.261758506E+00 0.265637338E+00 0.266105562E+00 + 0.690192208E-01 0.601205267E-01 0.522296093E-01 0.450656526E-01 0.386598185E-01 + 0.327284522E-01 0.271573998E-01 0.220119022E-01 0.174349211E-01 0.135394558E-01 + 0.103561282E-01 0.786284357E-02 0.609362870E-02 0.528478622E-02 0.599200279E-02 + 0.909225643E-02 0.155083835E-01 0.257546753E-01 0.396791101E-01 0.566278175E-01 + 0.758704692E-01 0.968931988E-01 0.119422033E+00 0.143315077E+00 0.168417528E+00 + 0.194381908E+00 0.220397249E+00 0.244863704E+00 0.265424937E+00 0.279904127E+00 + 0.287460029E+00 0.288978934E+00 0.286423385E+00 0.683522671E-01 0.607914105E-01 + 0.550814942E-01 0.486991629E-01 0.424978770E-01 0.365896709E-01 0.310205370E-01 + 0.258784331E-01 0.212847814E-01 0.173584670E-01 0.141999647E-01 0.119142085E-01 + 0.106694773E-01 0.107598826E-01 0.126118809E-01 0.167776868E-01 0.237944946E-01 + 0.339843705E-01 0.473422110E-01 0.636000410E-01 0.824079588E-01 0.103493713E+00 + 0.126724109E+00 0.152084500E+00 0.179600969E+00 0.209182546E+00 0.240238369E+00 + 0.270769328E+00 0.296457112E+00 0.312964380E+00 0.318753332E+00 0.315872103E+00 + 0.308505028E+00 0.712740123E-01 0.637281016E-01 0.577843972E-01 0.518189669E-01 + 0.458026230E-01 0.399136618E-01 0.342914835E-01 0.290676802E-01 0.243798681E-01 + 0.203693882E-01 0.171847641E-01 0.150029138E-01 0.140682533E-01 0.147226602E-01 + 0.172666162E-01 0.220628679E-01 0.294589177E-01 0.396882147E-01 0.528178215E-01 + 0.687768236E-01 0.874445960E-01 0.108745955E+00 0.132716447E+00 0.159535378E+00 + 0.189533651E+00 0.223166257E+00 0.260802925E+00 0.301590264E+00 0.336497635E+00 + 0.356238335E+00 0.358198136E+00 0.346335530E+00 0.331488043E+00 0.719051659E-01 + 0.651957318E-01 0.604294762E-01 0.546944104E-01 0.487194881E-01 0.427543223E-01 + 0.369887128E-01 0.315881222E-01 0.267113969E-01 0.225211605E-01 0.191947520E-01 + 0.169456378E-01 0.160665736E-01 0.170676559E-01 0.200299099E-01 0.252011791E-01 + 0.328335539E-01 0.431188121E-01 0.561631769E-01 0.720036924E-01 0.906597525E-01 + 0.112198636E+00 0.136801422E+00 0.164835826E+00 0.196953297E+00 0.234274149E+00 + 0.278785735E+00 0.333935350E+00 0.382554889E+00 0.407972217E+00 0.406733125E+00 + 0.378263652E+00 0.353199124E+00 0.757619888E-01 0.688297600E-01 0.633113235E-01 + 0.574831590E-01 0.513507053E-01 0.451630540E-01 0.391329937E-01 0.334433168E-01 + 0.282657892E-01 0.237748027E-01 0.201584883E-01 0.176287964E-01 0.164071321E-01 + 0.177386999E-01 0.209532380E-01 0.263328403E-01 0.340944454E-01 0.444105640E-01 + 0.574081540E-01 0.731832087E-01 0.918348581E-01 0.113516450E+00 0.138498917E+00 + 0.167261302E+00 0.200636253E+00 0.240111977E+00 0.288545549E+00 0.351945877E+00 + 0.413900256E+00 0.445567876E+00 0.444172204E+00 0.404598564E+00 0.369906157E+00 + 0.779870823E-01 0.714537874E-01 0.664044097E-01 0.602587834E-01 0.537447892E-01 + 0.471653640E-01 0.407339633E-01 0.346318781E-01 0.290316045E-01 0.241096653E-01 + 0.200564265E-01 0.170943439E-01 0.155964643E-01 0.168119073E-01 0.200834125E-01 + 0.255374610E-01 0.333374441E-01 0.436370969E-01 0.565705523E-01 0.722599700E-01 + 0.908402354E-01 0.112497166E+00 0.137520805E+00 0.166386262E+00 0.199880227E+00 + 0.239337981E+00 0.287072569E+00 0.347073019E+00 0.417508781E+00 0.453214049E+00 + 0.453254402E+00 0.416859329E+00 0.377476096E+00 0.830354691E-01 0.763174519E-01 + 0.699361339E-01 0.630261749E-01 0.558904521E-01 0.487529151E-01 0.417851880E-01 + 0.351470299E-01 0.290001109E-01 0.235150158E-01 0.188685022E-01 0.151905566E-01 + 0.134715028E-01 0.141325742E-01 0.174078718E-01 0.228512436E-01 0.306144804E-01 + 0.408425406E-01 0.536683947E-01 0.692190528E-01 0.876317695E-01 0.109079979E+00 + 0.133809403E+00 0.162188470E+00 0.194772109E+00 0.232387692E+00 0.276204765E+00 + 0.327719986E+00 0.391630620E+00 0.422646165E+00 0.428285092E+00 0.408697575E+00 + 0.373105764E+00 0.886474848E-01 0.817487314E-01 0.736334622E-01 0.656376556E-01 + 0.577105209E-01 0.498864353E-01 0.422663577E-01 0.349778384E-01 0.281656608E-01 + 0.219970867E-01 0.166804455E-01 0.124590658E-01 0.102801993E-01 0.103594176E-01 + 0.129400119E-01 0.182980746E-01 0.259570405E-01 0.360574126E-01 0.487261117E-01 + 0.640803427E-01 0.822370127E-01 0.103326946E+00 0.127510756E+00 0.154995099E+00 + 0.186030611E+00 0.220866114E+00 0.259552330E+00 0.301282585E+00 0.341756910E+00 + 0.369875699E+00 0.380467087E+00 0.374313891E+00 0.359036595E+00 0.938926265E-01 + 0.864606947E-01 0.767870024E-01 0.678369701E-01 0.590935014E-01 0.505121052E-01 + 0.421514921E-01 0.341128893E-01 0.265198760E-01 0.195198096E-01 0.133587606E-01 + 0.876656920E-02 0.606423616E-02 0.551274046E-02 0.736675784E-02 0.118926540E-01 + 0.193831176E-01 0.293042660E-01 0.417731553E-01 0.568884537E-01 0.747344196E-01 + 0.953841433E-01 0.118900254E+00 0.145324647E+00 0.174635038E+00 0.206640586E+00 + 0.240726352E+00 0.275305241E+00 0.306971252E+00 0.330945134E+00 0.344194055E+00 + 0.347022712E+00 0.342863888E+00 0.956574455E-01 0.885940716E-01 0.791120976E-01 + 0.694731995E-01 0.599569455E-01 0.505859405E-01 0.414200462E-01 0.325504690E-01 + 0.240776502E-01 0.160646737E-01 0.934401527E-02 0.417944416E-02 0.840134919E-03 +-0.414751470E-03 0.670015812E-03 0.435190648E-02 0.108976215E-01 0.205958523E-01 + 0.328371227E-01 0.476977304E-01 0.652244762E-01 0.854319334E-01 0.108292460E+00 + 0.133713782E+00 0.161484331E+00 0.191172019E+00 0.221941575E+00 0.252300054E+00 + 0.279964983E+00 0.302356511E+00 0.317742795E+00 0.326032400E+00 0.328657836E+00 + 0.997702405E-01 0.916645974E-01 0.808631480E-01 0.705471411E-01 0.602759235E-01 + 0.500874482E-01 0.400616340E-01 0.303007066E-01 0.209217407E-01 0.120403655E-01 + 0.457078591E-02-0.132605061E-02-0.542169437E-02-0.746523589E-02-0.720171630E-02 +-0.437980145E-02 0.124763697E-02 0.992981344E-02 0.219345912E-01 0.365553871E-01 + 0.538033545E-01 0.736408532E-01 0.959686860E-01 0.120602056E+00 0.147224039E+00 + 0.175315708E+00 0.204055116E+00 0.232230902E+00 0.258281976E+00 0.280601412E+00 + 0.298033267E+00 0.310228854E+00 0.317624092E+00 0.997893065E-01 0.923687071E-01 + 0.820303261E-01 0.711124912E-01 0.600672998E-01 0.490164198E-01 0.380678438E-01 + 0.273414850E-01 0.169855841E-01 0.725288689E-02-0.106307119E-02-0.781808048E-02 +-0.127868950E-01-0.157126710E-01-0.163323879E-01-0.143931210E-01-0.965905562E-02 +-0.191175193E-02 0.906132907E-02 0.234882981E-01 0.405417942E-01 0.601420440E-01 + 0.821379274E-01 0.106288522E+00 0.132232100E+00 0.159451008E+00 0.187223181E+00 + 0.214601189E+00 0.240467474E+00 0.263694644E+00 0.283363283E+00 0.298915416E+00 + 0.310175806E+00 0.103398286E+00 0.947051421E-01 0.829287097E-01 0.712576136E-01 + 0.593674891E-01 0.473858193E-01 0.354369134E-01 0.236446355E-01 0.121238194E-01 + 0.160373375E-02-0.764080137E-02-0.153765306E-01-0.213423166E-01-0.252558384E-01 +-0.268334392E-01-0.258095358E-01-0.219495837E-01-0.150557868E-01-0.496050715E-02 + 0.849099457E-02 0.254717357E-01 0.450121574E-01 0.669255704E-01 0.909380093E-01 + 0.116675749E+00 0.143663868E+00 0.171302795E+00 0.198864043E+00 0.225522503E+00 + 0.250421137E+00 0.272735417E+00 0.291698247E+00 0.306614608E+00 0.103025205E+00 + 0.950401723E-01 0.835409984E-01 0.710467175E-01 0.582101829E-01 0.452147573E-01 + 0.321843326E-01 0.192545634E-01 0.666677579E-02-0.487890095E-02-0.152111370E-01 +-0.240784045E-01-0.311840251E-01-0.362089463E-01-0.388357416E-01-0.387737341E-01 +-0.357778445E-01-0.296601728E-01-0.202880278E-01-0.756858662E-02 0.858043507E-02 + 0.282715149E-01 0.503811315E-01 0.746127963E-01 0.100581802E+00 0.127876654E+00 + 0.156024739E+00 0.184486389E+00 0.212666586E+00 0.239915669E+00 0.265485555E+00 + 0.288407862E+00 0.307390630E+00 0.106172390E+00 0.974910706E-01 0.841548666E-01 + 0.705044568E-01 0.566085540E-01 0.425186679E-01 0.283145886E-01 0.140657444E-01 + 0.469096005E-03-0.122417379E-01-0.238295272E-01-0.340002626E-01-0.424109362E-01 +-0.486928746E-01-0.524807908E-01-0.534445606E-01-0.513155162E-01-0.459044911E-01 +-0.371042341E-01-0.248753335E-01-0.921531394E-02 0.988516957E-02 0.324927270E-01 + 0.573125258E-01 0.838961527E-01 0.111926541E+00 0.141061634E+00 0.170938715E+00 + 0.201185197E+00 0.231392980E+00 0.261001468E+00 0.288980931E+00 0.313279301E+00 + 0.108466648E+00 0.990995094E-01 0.843003765E-01 0.695417970E-01 0.545542911E-01 + 0.393218845E-01 0.239067450E-01 0.841036066E-02-0.643950701E-02-0.205054525E-01 +-0.335396565E-01-0.452082530E-01-0.551130734E-01-0.628222004E-01-0.679065958E-01 +-0.699803457E-01-0.687367022E-01-0.639737099E-01-0.556014366E-01-0.436284989E-01 +-0.281256549E-01-0.917414948E-02 0.131963454E-01 0.390320942E-01 0.665329024E-01 + 0.956161767E-01 0.126065060E+00 0.157682240E+00 0.190332294E+00 0.223936498E+00 + 0.258377701E+00 0.293115169E+00 0.325654328E+00 0.108109899E+00 0.986868218E-01 + 0.837900564E-01 0.681070387E-01 0.520477034E-01 0.356368758E-01 0.189439915E-01 + 0.217206404E-02-0.140795484E-01-0.296836011E-01-0.443671457E-01-0.577470139E-01 +-0.693580136E-01-0.786890686E-01-0.852295235E-01-0.885194018E-01-0.881976038E-01 +-0.840371996E-01-0.759591609E-01-0.640182048E-01-0.483611822E-01-0.291555244E-01 +-0.648099184E-02 0.198845789E-01 0.483992212E-01 0.787462443E-01 0.110697709E+00 + 0.144182488E+00 0.179300770E+00 0.216368333E+00 0.255974680E+00 0.298988193E+00 + 0.345067650E+00 0.110124052E+00 0.995138884E-01 0.828729048E-01 0.662502125E-01 + 0.491303355E-01 0.314993896E-01 0.133725200E-01-0.465883501E-02-0.224373396E-01 +-0.397600941E-01-0.563040301E-01-0.716234818E-01-0.851719901E-01-0.963421911E-01 +-0.104521342E+00-0.109154873E+00-0.109809093E+00-0.106219709E+00-0.983130857E-01 +-0.861941949E-01-0.700991452E-01-0.503171496E-01-0.270935353E-01-0.537738204E-03 + 0.293026678E-01 0.611412376E-01 0.946811661E-01 0.129969284E+00 0.167304382E+00 + 0.207346037E+00 0.251358569E+00 0.301784098E+00 0.367813706E+00 0.107843198E+00 + 0.975676104E-01 0.814865455E-01 0.640605092E-01 0.458915047E-01 0.270411652E-01 + 0.755168684E-02-0.119622042E-01-0.314403027E-01-0.506719835E-01-0.692926794E-01 +-0.867883116E-01-0.102518238E+00-0.115760878E+00-0.125779390E+00-0.131902009E+00 +-0.133602619E+00-0.130565509E+00-0.122717962E+00-0.110220209E+00-0.934119672E-01 +-0.727236047E-01-0.485703088E-01-0.212485325E-01 0.913102180E-02 0.427533761E-01 + 0.778691843E-01 0.114734814E+00 0.153762266E+00 0.195767418E+00 0.242293373E+00 + 0.296201408E+00 0.362947375E+00 0.110547587E+00 0.982940570E-01 0.800264701E-01 + 0.616705008E-01 0.424271226E-01 0.223324988E-01 0.145780668E-02-0.196594223E-01 +-0.409864187E-01-0.623053722E-01-0.832111612E-01-0.103115492E+00-0.121271670E+00 +-0.136825055E+00-0.148892641E+00-0.156660944E+00-0.159489602E+00-0.156997025E+00 +-0.149106532E+00-0.136038974E+00-0.118251771E+00-0.963338837E-01-0.708771497E-01 +-0.423406847E-01-0.109030697E-01 0.238047652E-01 0.602973625E-01 0.984174386E-01 + 0.138486117E+00 0.181217805E+00 0.227944508E+00 0.281122953E+00 0.343576372E+00 + 0.106664397E+00 0.955226794E-01 0.784057826E-01 0.592022799E-01 0.388338268E-01 + 0.174408052E-01-0.484151952E-02-0.276364498E-01-0.509275012E-01-0.744846612E-01 +-0.978592709E-01-0.120383769E+00-0.141193956E+00-0.159285530E+00-0.173606053E+00 +-0.183176011E+00-0.187218711E+00-0.185270652E+00-0.177245915E+00-0.163432300E+00 +-0.144420102E+00-0.120976977E+00-0.938943475E-01-0.638313591E-01-0.311749950E-01 + 0.406262279E-02 0.420984328E-01 0.812231004E-01 0.121757820E+00 0.164145306E+00 + 0.209107608E+00 0.257787496E+00 0.311996251E+00 0.109487720E+00 0.970690772E-01 + 0.772285536E-01 0.567766465E-01 0.352052227E-01 0.124047082E-01-0.112240259E-01 +-0.357293338E-01-0.610595681E-01-0.869635269E-01-0.112947181E+00-0.138260216E+00 +-0.161912650E+00-0.182735637E+00-0.199487031E+00-0.210998207E+00-0.216333672E+00 +-0.214933515E+00-0.206696823E+00-0.191984087E+00-0.171532929E+00-0.146309987E+00 +-0.117327452E+00-0.854644030E-01-0.513131768E-01-0.150457956E-01 0.237673000E-01 + 0.636264086E-01 0.104258582E+00 0.145700574E+00 0.187971652E+00 0.230632305E+00 + 0.271143347E+00 0.108472645E+00 0.959820673E-01 0.758646429E-01 0.543803833E-01 + 0.316740125E-01 0.760102086E-02-0.174713861E-01-0.437188745E-01-0.711209774E-01 +-0.994223207E-01-0.128090009E+00-0.156291738E+00-0.182907447E+00-0.206593066E+00 +-0.225901440E+00-0.239455447E+00-0.246142432E+00-0.245291844E+00-0.236782327E+00 +-0.221050978E+00-0.198995501E+00-0.171797514E+00-0.140707269E+00-0.106839545E+00 +-0.710103810E-01-0.336268879E-01 0.534822419E-02 0.461801328E-01 0.868630260E-01 + 0.127300501E+00 0.167103320E+00 0.205303162E+00 0.239773110E+00 0.108515389E+00 + 0.954614878E-01 0.744008273E-01 0.520305596E-01 0.283115413E-01 0.306869857E-02 +-0.234130155E-01-0.513666756E-01-0.808048397E-01-0.111475423E+00-0.142813236E+00 +-0.173907802E+00-0.203510046E+00-0.230096906E+00-0.252007663E+00-0.267644614E+00 +-0.275705159E+00-0.275397241E+00-0.266575009E+00-0.249753460E+00-0.225997403E+00 +-0.196711808E+00-0.163393214E+00-0.127395093E+00-0.897499025E-01-0.510552302E-01 +-0.114105865E-01 0.296980441E-01 0.704756677E-01 0.110183023E+00 0.148249403E+00 + 0.183759674E+00 0.215285391E+00 0.108179197E+00 0.946479812E-01 0.728809386E-01 + 0.497963540E-01 0.252081081E-01-0.108141731E-02-0.288661327E-01-0.584158078E-01 +-0.897708982E-01-0.122686312E+00-0.156569883E+00-0.190440744E+00-0.222926378E+00 +-0.252330005E+00-0.276778787E+00-0.294452637E+00-0.303855449E+00-0.304069489E+00 +-0.294920146E+00-0.276998550E+00-0.251534879E+00-0.220156550E+00-0.184602082E+00 +-0.146454975E+00-0.106950618E+00-0.668640137E-01-0.264815744E-01 0.143401735E-01 + 0.558190681E-01 0.952347219E-01 0.132297114E+00 0.166339546E+00 0.196436435E+00 + 0.106608398E+00 0.930255204E-01 0.713842064E-01 0.477811843E-01 0.224769600E-01 +-0.471066032E-02-0.336374566E-01-0.646002814E-01-0.976644531E-01-0.132593021E+00 +-0.168773770E+00-0.205164298E+00-0.240282297E+00-0.272271067E+00-0.299060076E+00 +-0.318618804E+00-0.329265863E+00-0.329961210E+00-0.320499569E+00-0.301541418E+00 +-0.274471670E+00-0.241125047E+00-0.203462705E+00-0.163276836E+00-0.121978439E+00 +-0.804817975E-01-0.391756035E-01 0.204949826E-02 0.437841713E-01 0.830959231E-01 + 0.119678639E+00 0.153022453E+00 0.182497576E+00 0.107437558E+00 0.929754302E-01 + 0.701431930E-01 0.461139940E-01 0.202403273E-01-0.767659815E-02-0.375377163E-01 +-0.696629956E-01-0.104139924E+00-0.140741631E+00-0.178841755E+00-0.217349529E+00 +-0.254691362E+00-0.288875490E+00-0.317661375E+00-0.338836104E+00-0.350554407E+00 +-0.351666987E+00-0.341938704E+00-0.322088629E+00-0.293632448E+00-0.258582622E+00 +-0.219089687E+00-0.177119404E+00-0.134223461E+00-0.914174095E-01-0.491506718E-01 +-0.732490048E-02 0.346972942E-01 0.741077140E-01 0.110596515E+00 0.143738598E+00 + 0.173055679E+00 0.104138426E+00 0.904099643E-01 0.691768229E-01 0.449235886E-01 + 0.186068397E-01-0.985502917E-02-0.404017940E-01-0.733776912E-01-0.108891211E+00 +-0.146726072E+00-0.186248198E+00-0.226333186E+00-0.265340865E+00-0.301178783E+00 +-0.331477791E+00-0.353884429E+00-0.366427362E+00-0.367870301E+00-0.357951909E+00 +-0.337432444E+00-0.307925791E+00-0.271576643E+00-0.230676964E+00-0.187319562E+00 +-0.143155828E+00-0.992680937E-01-0.561413877E-01-0.136849415E-01 0.287349075E-01 + 0.684482381E-01 0.105151512E+00 0.138448402E+00 0.167910203E+00 0.108899273E+00 + 0.934590772E-01 0.691408366E-01 0.442984551E-01 0.176373497E-01-0.111616626E-01 +-0.421073399E-01-0.755712390E-01-0.111680448E+00-0.150228098E+00-0.190578625E+00 +-0.231589690E+00-0.271583557E+00-0.308409154E+00-0.339620888E+00-0.362781525E+00 +-0.375840545E+00-0.377506614E+00-0.367499709E+00-0.346599787E+00-0.316475928E+00 +-0.279348642E+00-0.237590492E+00-0.193365782E+00-0.148379117E+00-0.103745036E+00 +-0.599593557E-01-0.169277657E-01 0.259896405E-01 0.662027448E-01 0.103402480E+00 + 0.137166753E+00 0.167024821E+00 0.106082328E+00 0.919306949E-01 0.691163689E-01 + 0.441152081E-01 0.173228029E-01-0.115621733E-01-0.425863490E-01-0.761404783E-01 +-0.112362713E+00-0.151052535E+00-0.191577435E+00-0.232793614E+00-0.273017675E+00 +-0.310086459E+00-0.341537595E+00-0.364912421E+00-0.378139704E+00-0.379910320E+00 +-0.369932771E+00-0.348984987E+00-0.318741918E+00-0.281435430E+00-0.239450738E+00 +-0.194961742E+00-0.149677008E+00-0.104705624E+00-0.605287291E-01-0.170278214E-01 + 0.264676418E-01 0.673820674E-01 0.105381630E+00 0.139961645E+00 0.170509607E+00 + 0.108138137E+00 0.933433250E-01 0.692931041E-01 0.442848988E-01 0.176206026E-01 +-0.110663623E-01-0.418291278E-01-0.750612095E-01-0.110900365E+00-0.149147764E+00 +-0.189177305E+00-0.229860291E+00-0.269539148E+00-0.306086898E+00-0.337084979E+00 +-0.360119879E+00-0.373158187E+00-0.374911368E+00-0.365086913E+00-0.344437301E+00 +-0.314593583E+00-0.277732998E+00-0.236184001E+00-0.192067310E+00-0.147043630E+00 +-0.102178879E+00-0.579113290E-01-0.140741579E-01 0.300674625E-01 0.719059557E-01 + 0.111082122E+00 0.146956086E+00 0.178638160E+00 0.106306896E+00 0.922275409E-01 + 0.693566874E-01 0.447291993E-01 0.184785984E-01-0.971650705E-02-0.398795642E-01 +-0.723876059E-01-0.107365668E+00-0.144610509E+00-0.183506072E+00-0.222954631E+00 +-0.261354178E+00-0.296659172E+00-0.326551914E+00-0.348725468E+00-0.361238778E+00 +-0.362860888E+00-0.353306592E+00-0.333282769E+00-0.304331034E+00-0.268512011E+00 +-0.228033483E+00-0.184903488E+00-0.140686542E+00-0.963676572E-01-0.523134321E-01 +-0.828482956E-02 0.365518332E-01 0.795566216E-01 0.120403655E+00 0.158303842E+00 + 0.191897869E+00 0.106239870E+00 0.921886265E-01 0.694286823E-01 0.454096459E-01 + 0.198475942E-01-0.757687259E-02-0.368266180E-01-0.682443008E-01-0.101930521E+00 +-0.137673169E+00-0.174870342E+00-0.212469473E+00-0.248950854E+00-0.282389045E+00 +-0.310615510E+00-0.331481069E+00-0.343182236E+00-0.344576001E+00-0.335391998E+00 +-0.316274613E+00-0.288637906E+00-0.254378080E+00-0.215525895E+00-0.173926666E+00 +-0.131006226E+00-0.876326561E-01-0.440712348E-01 0.637210906E-05 0.455151126E-01 + 0.898924693E-01 0.133021131E+00 0.174100116E+00 0.211083844E+00 0.104935400E+00 + 0.912121758E-01 0.693669096E-01 0.463101752E-01 0.216885190E-01-0.472578546E-02 +-0.327950045E-01-0.628125742E-01-0.948474258E-01-0.128676817E+00-0.163717806E+00 +-0.198975623E+00-0.233034581E+00-0.264119387E+00-0.290246308E+00-0.309461504E+00 +-0.320131540E+00-0.321221828E+00-0.312482655E+00-0.294484317E+00-0.268488109E+00 +-0.236188054E+00-0.199400693E+00-0.159770861E+00-0.118550584E+00-0.764582679E-01 +-0.336163379E-01 0.104248095E-01 0.563589595E-01 0.102124706E+00 0.148105279E+00 + 0.194051877E+00 0.237502724E+00 0.101298369E+00 0.886008516E-01 0.693148226E-01 + 0.474952459E-01 0.239776671E-01-0.125063024E-02-0.279332511E-01-0.563121699E-01 +-0.864227563E-01-0.118033417E+00-0.150585666E+00-0.183153600E+00-0.214441061E+00 +-0.242842793E+00-0.266582251E+00-0.283922076E+00-0.293417186E+00-0.294152856E+00 +-0.285903156E+00-0.269158483E+00-0.245012492E+00-0.214938596E+00-0.180517003E+00 +-0.143170193E+00-0.103955798E+00-0.634165630E-01-0.214626752E-01 0.228156522E-01 + 0.682688877E-01 0.115022667E+00 0.163894027E+00 0.216175407E+00 0.273821950E+00 + 0.104736276E+00 0.904621333E-01 0.699748546E-01 0.490563810E-01 0.266818348E-01 + 0.275011361E-02-0.224049557E-01-0.489822924E-01-0.769877285E-01-0.106184460E+00 +-0.136045426E+00-0.165721059E+00-0.194044560E+00-0.219590366E+00-0.240797460E+00 +-0.256152153E+00-0.264401764E+00-0.264754385E+00-0.257008016E+00-0.241572082E+00 +-0.219370648E+00-0.191652402E+00-0.159755021E+00-0.124874152E+00-0.878731683E-01 +-0.491360687E-01-0.845680945E-02 0.351161733E-01 0.799905285E-01 0.127035171E+00 + 0.177729517E+00 0.234691992E+00 0.300893903E+00 0.100594975E+00 0.887555331E-01 + 0.707772151E-01 0.509076715E-01 0.297252629E-01 0.716857053E-02-0.163764730E-01 +-0.410631038E-01-0.668707341E-01-0.935638845E-01-0.120652117E+00-0.147368073E+00 +-0.172677323E+00-0.195334375E+00-0.213990927E+00-0.227350995E+00-0.234347761E+00 +-0.234306440E+00-0.227048039E+00-0.212904528E+00-0.192636862E+00-0.167276770E+00 +-0.137928098E+00-0.105567202E+00-0.708703399E-01-0.340751708E-01 0.509181432E-02 + 0.472190268E-01 0.907884538E-01 0.136912897E+00 0.187317923E+00 0.245030075E+00 + 0.315692067E+00 0.103110529E+00 0.906343907E-01 0.718365312E-01 0.529291704E-01 + 0.330056921E-01 0.119070113E-01-0.100049917E-01-0.327794664E-01-0.563750900E-01 +-0.805661306E-01-0.104904599E+00-0.128706932E+00-0.151069030E+00-0.170918077E+00 +-0.187106431E+00-0.198540539E+00-0.204324633E+00-0.203891098E+00-0.197081000E+00 +-0.184154049E+00-0.165721118E+00-0.142614037E+00-0.115720712E+00-0.858157352E-01 +-0.534013286E-01-0.185437948E-01 0.193651263E-01 0.589775555E-01 0.100314200E+00 + 0.144112334E+00 0.191766635E+00 0.245858535E+00 0.315825790E+00 0.999498740E-01 + 0.890942290E-01 0.726569965E-01 0.549634099E-01 0.364141129E-01 0.169340521E-01 +-0.343187572E-02-0.243300200E-01-0.457613841E-01-0.675247461E-01-0.892180949E-01 +-0.110240169E+00-0.129809871E+00-0.147014320E+00-0.160888076E+00-0.170518965E+00 +-0.175162748E+00-0.174344838E+00-0.167922795E+00-0.156091169E+00-0.139326319E+00 +-0.118281372E+00-0.936524719E-01-0.660396069E-01-0.358166695E-01-0.299849920E-02 + 0.329746492E-01 0.700855479E-01 0.108588278E+00 0.148930013E+00 0.191818595E+00 + 0.238222674E+00 0.287859082E+00 0.100721255E+00 0.898159072E-01 0.733174607E-01 + 0.568919703E-01 0.398034863E-01 0.219017845E-01 0.316727906E-02-0.158836544E-01 +-0.352360532E-01-0.546987094E-01-0.739102289E-01-0.923456997E-01-0.109336182E+00 +-0.124112070E+00-0.135868385E+00-0.143850073E+00-0.147442669E+00-0.146249652E+00 +-0.140140250E+00-0.129251048E+00-0.113939770E+00-0.947024971E-01-0.720653012E-01 +-0.464741364E-01-0.182049293E-01 0.126766395E-01 0.463654585E-01 0.807643831E-01 + 0.116060078E+00 0.152347788E+00 0.189623237E+00 0.227405146E+00 0.263170660E+00 + 0.972504541E-01 0.874627233E-01 0.734389499E-01 0.586630479E-01 0.431428105E-01 + 0.268321354E-01 0.973232463E-02-0.754037127E-02-0.249385126E-01-0.422674604E-01 +-0.591996312E-01-0.752783567E-01-0.899350643E-01-0.102525055E+00-0.112380974E+00 +-0.118880399E+00-0.121515773E+00-0.119954072E+00-0.114071153E+00-0.103951082E+00 +-0.898487717E-01-0.721219704E-01-0.511428565E-01-0.272028875E-01-0.437689945E-03 + 0.290700085E-01 0.598539375E-01 0.913737640E-01 0.123368733E+00 0.155568197E+00 + 0.187510222E+00 0.218197688E+00 0.245525450E+00 0.939041823E-01 0.849240869E-01 + 0.733717680E-01 0.603976846E-01 0.464755893E-01 0.317519121E-01 0.162766837E-01 + 0.665500760E-03-0.149431815E-01-0.303353369E-01-0.452126153E-01-0.591812283E-01 +-0.717609450E-01-0.824141130E-01-0.905884653E-01-0.957707241E-01-0.975392014E-01 +-0.956097990E-01-0.898596421E-01-0.803257003E-01-0.671729222E-01-0.506373234E-01 +-0.309547521E-01-0.825838186E-02 0.176198091E-01 0.449853279E-01 0.733483434E-01 + 0.102211304E+00 0.131120712E+00 0.159584314E+00 0.186962366E+00 0.212308243E+00 + 0.234254614E+00 0.926011205E-01 0.841351226E-01 0.738050714E-01 0.623026378E-01 + 0.498688854E-01 0.366917178E-01 0.229180735E-01 0.875650533E-02-0.527309394E-02 +-0.189427417E-01-0.319967195E-01-0.441022068E-01-0.548570789E-01-0.638139546E-01 +-0.705150291E-01-0.745342970E-01-0.755169615E-01-0.732134357E-01-0.674972832E-01 +-0.583625734E-01-0.458949134E-01-0.302223954E-01-0.114752790E-01 0.102604721E-01 + 0.349461064E-01 0.606157891E-01 0.869987607E-01 0.113595143E+00 0.139827162E+00 + 0.165072918E+00 0.188661546E+00 0.209844992E+00 0.227805942E+00 0.919332653E-01 + 0.840397328E-01 0.747539029E-01 0.644493103E-01 0.533376187E-01 0.415814593E-01 + 0.293333810E-01 0.167650431E-01 0.407250039E-02-0.808288343E-02-0.195351765E-01 +-0.300117601E-01-0.391784720E-01-0.466615483E-01-0.520793237E-01-0.550736561E-01 +-0.553392507E-01-0.526477545E-01-0.468629971E-01-0.379400887E-01-0.258950517E-01 +-0.107400147E-01 0.751061179E-02 0.288750529E-01 0.519508719E-01 0.761689991E-01 + 0.101012290E+00 0.125841498E+00 0.149921343E+00 0.172502339E+00 0.192912653E+00 + 0.210615456E+00 0.225219503E+00 0.921635851E-01 0.846964419E-01 0.760896727E-01 + 0.667977706E-01 0.568593331E-01 0.463921688E-01 0.355256908E-01 0.243915040E-01 + 0.130808773E-01 0.229033642E-02-0.775934383E-02-0.168201216E-01-0.246140156E-01 +-0.308230538E-01-0.351248197E-01-0.372139998E-01-0.368189998E-01-0.337185264E-01 +-0.277574677E-01-0.188579373E-01-0.699315825E-02 0.798624847E-02 0.262024291E-01 + 0.464684516E-01 0.684544370E-01 0.916836113E-01 0.115551338E+00 0.139256135E+00 + 0.161828011E+00 0.182299793E+00 0.199951127E+00 0.214457110E+00 0.225852519E+00 + 0.923170671E-01 0.854138285E-01 0.776394978E-01 0.692670643E-01 0.603910498E-01 + 0.511004031E-01 0.415068530E-01 0.317372233E-01 0.219315346E-01 0.122916438E-01 + 0.344657525E-02-0.438740430E-02-0.110029299E-01-0.161122195E-01-0.194413625E-01 +-0.207271781E-01-0.197190251E-01-0.161853638E-01-0.992959458E-02-0.844377093E-03 + 0.110159619E-01 0.261383615E-01 0.435197428E-01 0.630047694E-01 0.843368024E-01 + 0.107114628E+00 0.130713761E+00 0.154154032E+00 0.176069617E+00 0.195032641E+00 + 0.210137978E+00 0.221321762E+00 0.229225114E+00 0.929623544E-01 0.864743441E-01 + 0.792914182E-01 0.717774704E-01 0.638843626E-01 0.556746572E-01 0.472497083E-01 + 0.387339070E-01 0.302729402E-01 0.220415089E-01 0.142186023E-01 0.746976957E-02 + 0.184745248E-02-0.230868720E-02-0.478107855E-02-0.534883700E-02-0.377413444E-02 + 0.206209719E-03 0.689089950E-02 0.165002421E-01 0.283090305E-01 0.426757559E-01 + 0.593965575E-01 0.783534348E-01 0.993984938E-01 0.122272603E+00 0.146456435E+00 + 0.170860380E+00 0.193427444E+00 0.211647600E+00 0.224156216E+00 0.231288552E+00 + 0.234844133E+00 0.932376012E-01 0.872538686E-01 0.809225291E-01 0.742695108E-01 + 0.673005581E-01 0.600859933E-01 0.527286455E-01 0.453519188E-01 0.380999111E-01 + 0.311425328E-01 0.246825684E-01 0.190066695E-01 0.141379638E-01 0.108320834E-01 + 0.913603883E-02 0.921851024E-02 0.113014737E-01 0.156814698E-01 0.228783265E-01 + 0.321694873E-01 0.436593145E-01 0.574996136E-01 0.737049058E-01 0.922945365E-01 + 0.113296218E+00 0.136709645E+00 0.162377775E+00 0.189567044E+00 0.215297490E+00 + 0.233870775E+00 0.243352041E+00 0.244497404E+00 0.241846934E+00 0.934127048E-01 + 0.879085809E-01 0.825293362E-01 0.767281577E-01 0.706189275E-01 0.643139333E-01 + 0.579251945E-01 0.515749902E-01 0.454015061E-01 0.395632498E-01 0.342438743E-01 + 0.296544917E-01 0.259830467E-01 0.235792007E-01 0.226311330E-01 0.233203135E-01 + 0.258437078E-01 0.304150581E-01 0.372283496E-01 0.461572558E-01 0.572431609E-01 + 0.705918819E-01 0.863100290E-01 0.104546919E+00 0.125533566E+00 0.149620518E+00 + 0.177295357E+00 0.209170133E+00 0.247614697E+00 0.265352756E+00 0.269254446E+00 + 0.259846300E+00 0.248184875E+00 0.938735902E-01 0.888414755E-01 0.842444301E-01 + 0.791730583E-01 0.738311410E-01 0.683448389E-01 0.628281534E-01 0.573980361E-01 + 0.521824248E-01 0.473250411E-01 0.429883748E-01 0.393527821E-01 0.366087258E-01 + 0.349854939E-01 0.346726142E-01 0.358444378E-01 0.386717767E-01 0.433180816E-01 + 0.499127321E-01 0.584922135E-01 0.691294745E-01 0.819419771E-01 0.970919207E-01 + 0.114825450E+00 0.135542631E+00 0.159924403E+00 0.189146370E+00 0.225228056E+00 + 0.261813730E+00 0.281611323E+00 0.283673346E+00 0.265449733E+00 0.250653744E+00 + 0.951518789E-01 0.903707296E-01 0.860929787E-01 0.816035718E-01 0.769256949E-01 + 0.721667632E-01 0.674312562E-01 0.628239736E-01 0.584584512E-01 0.544614457E-01 + 0.509752780E-01 0.481562689E-01 0.461720824E-01 0.452016667E-01 0.454123504E-01 + 0.469605550E-01 0.499940291E-01 0.546491556E-01 0.610405542E-01 0.692533627E-01 + 0.793855861E-01 0.915690809E-01 0.105989881E+00 0.122933246E+00 0.142867371E+00 + 0.166622594E+00 0.195799857E+00 0.233458668E+00 0.265146196E+00 0.281153470E+00 + 0.279981881E+00 0.261761844E+00 0.248058453E+00 0.959481969E-01 0.916335285E-01 + 0.879816189E-01 0.839848518E-01 0.798778459E-01 0.757641420E-01 0.717302039E-01 + 0.678614974E-01 0.642525852E-01 0.610120893E-01 0.582630076E-01 0.561411045E-01 + 0.547906123E-01 0.543604605E-01 0.549943037E-01 0.568278283E-01 0.599883161E-01 + 0.645935312E-01 0.707493573E-01 0.785525665E-01 0.881082043E-01 0.995465815E-01 + 0.113043092E+00 0.128854468E+00 0.147381008E+00 0.169297710E+00 0.195850030E+00 + 0.228817254E+00 0.254490554E+00 0.267166197E+00 0.265053004E+00 0.251658916E+00 + 0.241314426E+00 0.977753550E-01 0.935165137E-01 0.898913965E-01 0.862646773E-01 + 0.826484337E-01 0.791160241E-01 0.757223144E-01 0.725242123E-01 0.695937127E-01 + 0.670209527E-01 0.649120733E-01 0.633853376E-01 0.625659004E-01 0.625802502E-01 + 0.635508597E-01 0.655925572E-01 0.688127205E-01 0.733117759E-01 0.791848749E-01 + 0.865245908E-01 0.954291672E-01 0.106010832E+00 0.118406683E+00 0.132794037E+00 + 0.149406284E+00 0.168543503E+00 0.190507263E+00 0.214868590E+00 0.234376997E+00 + 0.244702861E+00 0.244494453E+00 0.238342777E+00 0.232021347E+00 0.986426398E-01 + 0.948681086E-01 0.916982368E-01 0.883595794E-01 0.851808935E-01 0.821973979E-01 + 0.794086009E-01 0.768324509E-01 0.745169669E-01 0.725361183E-01 0.709835440E-01 + 0.699644685E-01 0.695896447E-01 0.699685067E-01 0.712040067E-01 0.733902678E-01 + 0.766133368E-01 0.809550881E-01 0.864967182E-01 0.933202282E-01 0.101509377E+00 + 0.111149736E+00 0.122327641E+00 0.135125309E+00 0.149596646E+00 0.165695280E+00 + 0.183042526E+00 0.200264364E+00 0.214469865E+00 0.223240018E+00 0.226086140E+00 + 0.224774957E+00 0.221772045E+00 0.101063699E+00 0.969668999E-01 0.933295563E-01 + 0.901342556E-01 0.874001756E-01 0.849867389E-01 0.828014016E-01 0.808180571E-01 + 0.790656880E-01 0.776100010E-01 0.765384957E-01 0.759497136E-01 0.759453550E-01 + 0.766222775E-01 0.780661032E-01 0.803483203E-01 0.835299417E-01 0.876709372E-01 + 0.928368419E-01 0.990984738E-01 0.106525481E+00 0.115176998E+00 0.125092387E+00 + 0.136274114E+00 0.148649365E+00 0.161989242E+00 0.175751299E+00 0.188874647E+00 + 0.199865624E+00 0.207490295E+00 0.211443648E+00 0.212444365E+00 0.211659715E+00 + 0.102823384E+00 0.985840484E-01 0.943550467E-01 0.913522020E-01 0.892316028E-01 + 0.874868408E-01 0.859383792E-01 0.845310390E-01 0.832941979E-01 0.822994187E-01 + 0.816374570E-01 0.814070106E-01 0.817071870E-01 0.826280266E-01 0.842406005E-01 + 0.865896568E-01 0.897015259E-01 0.936069191E-01 0.983567685E-01 0.104018331E+00 + 0.110653989E+00 0.118303977E+00 0.126973853E+00 0.136616021E+00 0.147094697E+00 + 0.158131823E+00 0.169237003E+00 0.179673329E+00 0.188577935E+00 0.195266187E+00 + 0.199521393E+00 0.201640308E+00 0.202211052E+00 0.104448400E+00 0.100064203E+00 + 0.949191377E-01 0.919208899E-01 0.906667933E-01 0.897629112E-01 0.889022723E-01 + 0.880465582E-01 0.872685537E-01 0.866652504E-01 0.863386095E-01 0.863942206E-01 + 0.869364515E-01 0.880570337E-01 0.898199826E-01 0.922362655E-01 0.952724591E-01 + 0.989157781E-01 0.103205882E+00 0.108229548E+00 0.114057198E+00 0.120717928E+00 + 0.128198519E+00 0.136429250E+00 0.145257622E+00 0.154421106E+00 0.163525507E+00 + 0.172060028E+00 0.179485336E+00 0.185383603E+00 0.189591348E+00 0.192224339E+00 + 0.193587840E+00 0.106169961E+00 0.101967186E+00 0.966599509E-01 0.928079486E-01 + 0.918251052E-01 0.920024514E-01 0.918393731E-01 0.914624482E-01 0.910619497E-01 + 0.907695889E-01 0.906957537E-01 0.909585282E-01 0.916767344E-01 0.929544345E-01 + 0.948743895E-01 0.974203050E-01 0.100408182E+00 0.103770822E+00 0.107521877E+00 + 0.111861534E+00 0.116872579E+00 0.122568898E+00 0.128932953E+00 0.135893866E+00 + 0.143307343E+00 0.150946468E+00 0.158503398E+00 0.165614724E+00 0.171919882E+00 + 0.177140638E+00 0.181147471E+00 0.183974266E+00 0.185779780E+00 0.108074628E+00 + 0.104498193E+00 0.100286201E+00 0.955352262E-01 0.937536433E-01 0.946706310E-01 + 0.949395075E-01 0.948705077E-01 0.947446674E-01 0.946772024E-01 0.947553068E-01 + 0.951333866E-01 0.959508494E-01 0.973248780E-01 0.993885323E-01 0.102275424E+00 + 0.105280779E+00 0.108385332E+00 0.111377798E+00 0.115015492E+00 0.119210385E+00 + 0.123966821E+00 0.129286140E+00 0.135111660E+00 0.141315013E+00 0.147703573E+00 + 0.154033780E+00 0.160036266E+00 0.165454432E+00 0.170088589E+00 0.173831776E+00 + 0.176678866E+00 0.178708941E+00 0.110113911E+00 0.107487917E+00 0.104204640E+00 + 0.100623935E+00 0.970070660E-01 0.980412588E-01 0.982342139E-01 0.982439220E-01 + 0.982972234E-01 0.983847529E-01 0.985263139E-01 0.989403874E-01 0.997681990E-01 + 0.101141512E+00 0.103226788E+00 0.106115825E+00 0.109117210E+00 0.112096481E+00 + 0.114776745E+00 0.117790304E+00 0.121160291E+00 0.124978907E+00 0.129312754E+00 + 0.134124309E+00 0.139293358E+00 0.144648418E+00 0.149988338E+00 0.155103385E+00 + 0.159799442E+00 0.163923621E+00 0.167385146E+00 0.170159250E+00 0.172280356E+00 + 0.620565712E+00 0.642535746E+00 0.660853863E+00 0.682174385E+00 0.706194222E+00 + 0.732125759E+00 0.760542989E+00 0.791295230E+00 0.824818015E+00 0.860441208E+00 + 0.900061369E+00 0.943472266E+00 0.991030216E+00 0.104318440E+01 0.110040069E+01 + 0.116347051E+01 0.123315895E+01 0.131042480E+01 0.139618754E+01 0.149200296E+01 + 0.159953809E+01 0.172025740E+01 0.185701382E+01 0.201290226E+01 0.219133496E+01 + 0.239664483E+01 0.263768077E+01 0.292083907E+01 0.325979400E+01 0.367491221E+01 + 0.420350933E+01 0.495448399E+01 0.793196821E+01 + 67 50 + 0.105477178E+01 0.499999523E-01 0.105474961E+01 0.149999976E+00 0.105916715E+01 + 0.250000000E+00 0.106844461E+01 0.350000024E+00 0.108386970E+01 0.449999928E+00 + 0.110562503E+01 0.544746220E+00 0.112120986E+01 0.600000024E+00 0.115830672E+01 + 0.699999928E+00 0.118317974E+01 0.750000119E+00 0.121337843E+01 0.800000072E+00 + 0.126499999E+01 0.868783653E+00 0.131812501E+01 0.921098590E+00 0.137125003E+01 + 0.960451126E+00 0.145275426E+01 0.100000012E+01 0.153062499E+01 0.102186286E+01 + 0.163687503E+01 0.102548373E+01 0.174312508E+01 0.100245833E+01 0.179624999E+01 + 0.979757428E+00 0.184937501E+01 0.949727774E+00 0.192143464E+01 0.899999976E+00 + 0.198017967E+01 0.850000024E+00 0.203086877E+01 0.800000072E+00 0.207510018E+01 + 0.750000119E+00 0.211500001E+01 0.696638465E+00 0.216812491E+01 0.611374080E+00 + 0.219678116E+01 0.550000072E+00 0.222125006E+01 0.489633322E+00 0.225156307E+01 + 0.399999976E+00 0.227437496E+01 0.324134737E+00 0.229351902E+01 0.250000000E+00 + 0.231105804E+01 0.149999976E+00 0.231563663E+01 0.499999523E-01 0.230437565E+01 +-0.499999523E-01 0.227702069E+01-0.149999976E+00 0.225719738E+01-0.200000048E+00 + 0.222125006E+01-0.280016690E+00 0.218536472E+01-0.350000024E+00 0.215680003E+01 +-0.399999976E+00 0.211500001E+01-0.468032271E+00 0.206187510E+01-0.541482925E+00 + 0.201534939E+01-0.600000024E+00 0.197094786E+01-0.650000036E+00 0.192361236E+01 +-0.699999988E+00 0.187413692E+01-0.750000000E+00 0.182329488E+01-0.800000012E+00 + 0.177178347E+01-0.850000024E+00 0.172017288E+01-0.900000036E+00 0.166888404E+01 +-0.949999988E+00 0.161811340E+01-0.100000000E+01 0.156749690E+01-0.105000007E+01 + 0.150382864E+01-0.112791836E+01 0.145070362E+01-0.112266588E+01 0.140682149E+01 +-0.105000007E+01 0.137125003E+01-0.995866239E+00 0.131812501E+01-0.904936910E+00 + 0.129201245E+01-0.850000024E+00 0.126499999E+01-0.794037580E+00 0.122389770E+01 +-0.699999988E+00 0.120315039E+01-0.650000036E+00 0.116678965E+01-0.549999952E+00 + 0.114986026E+01-0.500000000E+00 0.112070811E+01-0.399999976E+00 0.110562503E+01 +-0.339724422E+00 0.108619428E+01-0.250000000E+00 0.107028735E+01-0.149999976E+00 + 0.105979693E+01-0.499999523E-01 0.105477178E+01 0.499999523E-01 + 0.101680005E+01 0.000000000E+00 0.101680005E+01 0.399760008E+00 0.101680005E+01 + 0.414759994E+00 0.101680005E+01 0.121730006E+01 0.115090001E+01 0.135140002E+01 + 0.170959997E+01 0.135140002E+01 0.190550005E+01 0.119299996E+01 0.214059997E+01 + 0.100269997E+01 0.220350003E+01 0.850099981E+00 0.232929993E+01 0.544900000E+00 + 0.239059997E+01 0.392199993E+00 0.237560010E+01 0.392199993E+00 0.237560010E+01 + 0.365500003E+00 0.236169982E+01 0.358399987E+00 0.233959985E+01 0.324900001E+00 + 0.234979987E+01 0.218899995E+00 0.235249996E+01 0.201000005E+00 0.235270000E+01 + 0.200000003E+00 0.236029983E+01 0.150000006E+00 0.236649990E+01 0.109700002E+00 + 0.236689997E+01 0.100000001E+00 0.236879992E+01 0.500000007E-01 0.237079978E+01 + 0.000000000E+00 0.236879992E+01-0.500000007E-01 0.236689997E+01-0.100000001E+00 + 0.236649990E+01-0.109700002E+00 0.236029983E+01-0.150000006E+00 0.235270000E+01 +-0.200000003E+00 0.235249996E+01-0.201000005E+00 0.234979987E+01-0.218899995E+00 + 0.233959985E+01-0.324900001E+00 0.236169982E+01-0.358399987E+00 0.237560010E+01 +-0.365500003E+00 0.237560010E+01-0.392199993E+00 0.239059997E+01-0.392199993E+00 + 0.232890010E+01-0.545799971E+00 0.220239997E+01-0.852900028E+00 0.213910007E+01 +-0.100639999E+01 0.183669996E+01-0.136590004E+01 0.182210004E+01-0.136590004E+01 + 0.170000005E+01-0.136590004E+01 0.156299996E+01-0.136590004E+01 0.139999998E+01 +-0.136590004E+01 0.127499998E+01-0.136590004E+01 0.115390003E+01-0.136590004E+01 + 0.101680005E+01-0.122880006E+01 0.101680005E+01-0.800000012E+00 0.101680005E+01 +-0.414759994E+00 0.101680005E+01-0.399760008E+00 0.101680005E+01 0.000000000E+00 + $OUT1 + ISHOT = 66832, + ITIME = 2384, + BETAP0 = 0.5000000 , + RZERO = 1.695500 , + QENP = 0.0000000E+00, + ENP = 1.000000 , + EMP = 1.000000 , + PLASMA = 1207042. , + EXPMP2 = 0.4708018 , 0.4773179 , 0.3603379 , 0.1731134 , + 8.8836350E-02, 6.6125073E-02, 6.3459434E-02, 7.2555467E-02, 9.4245583E-02, + 0.1759746 , 0.2399994 , 0.2538166 , 0.3425750 , 0.3374638 , + 0.2822209 , 0.4350058 , 0.3136042 , 0.1966812 , 0.1151106 , + 2.1058967E-02, -8.1345916E-02, -7.8260824E-02, -0.1639674 , 0.1066249 , + 0.2524083 , 0.3028223 , 0.2626219 , 0.3335645 , 0.3088680 , + 5.9432302E-02, 7.1395457E-02, 7.2097011E-02, 9.4330192E-02, 0.1296596 , + 0.1851648 , 0.2401672 , 0.2550771 , 0.3443176 , 0.3368228 , + 0.2813148 , 0.3108888 , 0.3329883 , 0.2639162 , 0.3030512 , + 0.2476704 , 0.1063703 , -0.1598371 , -8.0829881E-02, -8.0310188E-02, + 2.2322556E-02, 0.1148047 , 0.1977697 , 0.3164317 , 0.4353944 , + 0.4713073 , 0.4773440 , 0.3592431 , 0.1806585 , 8.5105397E-02, + 6.4306766E-02, + COILS = -9.2823125E-02, -8.9525498E-02, -8.4404118E-02, -9.1545321E-02, + -9.1610327E-02, -0.2315395 , -0.2023126 , -9.4911344E-02, -0.1020187 , + -9.5083944E-02, -9.9068999E-02, -9.0892375E-02, -6.5901682E-02, -3.2564409E-02, + -0.2699479 , -0.3252287 , 6.6496804E-03, 3.9711285E-02, -7.8558229E-02, + -7.7792674E-02, -7.8605026E-02, -8.4278710E-02, -9.1448635E-02, -9.0284504E-02, + -9.7341761E-02, -0.1887342 , -0.1653577 , -0.2570573 , -0.1920957 , + -8.5791036E-02, -8.8832252E-02, -7.4515246E-02, -5.2190609E-02, -2.0975456E-02, + -3.8985893E-02, -0.2692401 , -0.2544583 , -0.3149689 , -0.2491259 , + 9.9443272E-04, 1.5731834E-02, + PREXP = 0.0000000E+00, + BTOR = -2.067678 , + RCENTR = 1.695500 , + BRSP = -110450.9 , -73700.44 , -4188.844 , 13269.04 , + 25485.17 , -198939.5 , -129642.4 , -1059.132 , -3718.815 , + -97422.38 , -75386.45 , -38799.90 , 32434.57 , 61597.44 , + -194932.8 , -259359.5 , 92993.48 , 164767.4 , -556.4534 , + 1054.729 , 16727.29 , -37086.41 , 20207.24 , -0.2777848 , + 0.2569490 , 1.2478198E-02, 0.4994555 , -0.2454668 , -0.3186080 , + -0.1362296 , -0.2273455 , 0.3317215 , 0.5245887 , -0.3367561 , + 0.9300969 , 0.2759959 , 0.4367171 , -0.2228742 , 0.3564245 , + 2.0588748E-02, -0.1974280 , -1.7139224E-02, 7.9159841E-02, 0.3827312 , + 0.2814572 , 0.1769798 , -0.3535235 , -1.2577477E-02, -0.6293046 , + -0.1232804 , -0.3905841 , 7.8122015E-03, 0.1830417 , 9.5276833E-03, + 0.4021844 , -0.3404119 , 0.9754753 , 0.4746936 , -0.2150628 , + -0.1389852 , 0.9889172 , 0.1053845 , 0.1360717 , -9.9650815E-02, + -0.5503032 , -0.4042709 , -3.1082230E-02, 0.7354022 , 0.3024738 , + -0.1116976 , 0.4510566 , -9.1353171E-02, -0.5226294 , 4.8181456E-02, + 9.0212382E-02, 4.5185860E-02, -0.5540689 , -0.3065124 , 6.4070046E-02, + 7.2500899E-02, 0.7516741 , -1.247255 , 0.2785085 , -0.6357061 , + -0.3152615 , 0.4278110 , 0.1217311 , 0.2365918 , -0.2087229 , + -7.3400296E-02, 0.2407893 , 0.3041606 , -3.3600659E-03, -0.2205222 , + -0.7147905 , 0.1511142 , 184*0.0000000E+00, + ICURRT = 2, + RBDRY = 1.054772 , 1.054750 , 1.059167 , 1.068445 , + 1.083870 , 1.105625 , 1.121210 , 1.158307 , 1.183180 , + 1.213378 , 1.265000 , 1.318125 , 1.371250 , 1.452754 , + 1.530625 , 1.636875 , 1.743125 , 1.796250 , 1.849375 , + 1.921435 , 1.980180 , 2.030869 , 2.075100 , 2.115000 , + 2.168125 , 2.196781 , 2.221250 , 2.251563 , 2.274375 , + 2.293519 , 2.311058 , 2.315637 , 2.304376 , 2.277021 , + 2.257197 , 2.221250 , 2.185365 , 2.156800 , 2.115000 , + 2.061875 , 2.015349 , 1.970948 , 1.923612 , 1.874137 , + 1.823295 , 1.771783 , 1.720173 , 1.668884 , 1.618113 , + 1.567497 , 1.503829 , 1.450704 , 1.406821 , 1.371250 , + 1.318125 , 1.292012 , 1.265000 , 1.223898 , 1.203150 , + 1.166790 , 1.149860 , 1.120708 , 1.105625 , 1.086194 , + 1.070287 , 1.059797 , 1.054772 , 43*0.0000000E+00, + ZBDRY = 4.9999952E-02, 0.1500000 , 0.2500000 , 0.3500000 , + 0.4499999 , 0.5447462 , 0.6000000 , 0.6999999 , 0.7500001 , + 0.8000001 , 0.8687837 , 0.9210986 , 0.9604511 , 1.000000 , + 1.021863 , 1.025484 , 1.002458 , 0.9797574 , 0.9497278 , + 0.9000000 , 0.8500000 , 0.8000001 , 0.7500001 , 0.6966385 , + 0.6113741 , 0.5500001 , 0.4896333 , 0.4000000 , 0.3241347 , + 0.2500000 , 0.1500000 , 4.9999952E-02, -4.9999952E-02, -0.1500000 , + -0.2000000 , -0.2800167 , -0.3500000 , -0.4000000 , -0.4680323 , + -0.5414829 , -0.6000000 , -0.6500000 , -0.7000000 , -0.7500000 , + -0.8000000 , -0.8500000 , -0.9000000 , -0.9500000 , -1.000000 , + -1.050000 , -1.127918 , -1.122666 , -1.050000 , -0.9958662 , + -0.9049369 , -0.8500000 , -0.7940376 , -0.7000000 , -0.6500000 , + -0.5500000 , -0.5000000 , -0.4000000 , -0.3397244 , -0.2500000 , + -0.1500000 , -4.9999952E-02, 4.9999952E-02, 43*0.0000000E+00, + NBDRY = 67, + FWTSI = 999.9999 , 556.7620 , 561.7230 , 565.8070 , + 566.2111 , 239.1327 , 0.0000000E+00, 556.3060 , 488.1909 , + 557.0880 , 571.6694 , 529.0601 , 568.2542 , 562.8655 , + 188.7680 , 142.0892 , 337.5253 , 254.2535 , 562.9677 , + 559.1573 , 546.5208 , 569.8792 , 567.6462 , 557.4705 , + 485.5167 , 350.1021 , 375.8311 , 205.7207 , 333.7776 , + 562.3906 , 564.2797 , 558.1378 , 554.5692 , 457.6178 , + 444.1566 , 194.2616 , 207.9348 , 151.7483 , 211.6798 , + 352.6099 , 307.0927 , + FWTCUR = 2.7477012E-05, + MXITER = 25, + NXITER = 1, + LIMITR = 50, + XLIM = 4*1.016800 , 1.150900 , 1.709600 , 1.905500 , + 2.140600 , 2.203500 , 2.329300 , 2.390600 , + 2*2.375600 , 2.361700 , 2.339600 , 2.349800 , 2.352500 , + 2.352700 , 2.360300 , 2.366500 , 2.366900 , 2.368800 , + 2.370800 , 2.368800 , 2.366900 , 2.366500 , 2.360300 , + 2.352700 , 2.352500 , 2.349800 , 2.339600 , 2.361700 , + 2*2.375600 , 2.390600 , 2.328900 , 2.202400 , 2.139100 , + 1.836700 , 1.822100 , 1.700000 , 1.563000 , 1.400000 , + 1.275000 , 1.153900 , 5*1.016800 , 2.392100 , + 2*2.379100 , 2*2.373100 , 2*2.372100 , 2*2.373100 , + 2*2.379100 , 2.392100 , 58*0.0000000E+00, + YLIM = 0.0000000E+00, 0.3997600 , 0.4147600 , 1.217300 , + 2*1.351400 , 1.193000 , 1.002700 , 0.8501000 , 0.5449000 , + 2*0.3922000 , 0.3655000 , 0.3584000 , 0.3249000 , 0.2189000 + , 0.2010000 , 0.2000000 , 0.1500000 , 0.1097000 , + 0.1000000 , 5.0000001E-02, 0.0000000E+00, -5.0000001E-02, -0.1000000 , + -0.1097000 , -0.1500000 , -0.2000000 , -0.2010000 , -0.2189000 , + -0.3249000 , -0.3584000 , -0.3655000 , 2*-0.3922000 , + -0.5458000 , -0.8529000 , -1.006400 , 7*-1.365900 , -1.228800 + , -0.8000000 , -0.4147600 , -0.3997600 , 0.0000000E+00, + 2*0.3922000 , 2*0.0000000E+00, 2*0.2000000 , 2*-0.2000000 , + 2*0.0000000E+00, 2*-0.3922000 , 58*0.0000000E+00, + ERROR = 1.0000000E-03, + ICONVR = 2, + IBUNMN = 3, + PRESSR = 132*0.0000000E+00, + RPRESS = 132*0.0000000E+00 + $END + $CHIOUT + SAISIL = 3.0562438E-05, 5.9179747E-03, 5.0062783E-02, 2.7595712E-02, + 0.1488641 , 2.6582466E-02, 0.0000000E+00, 2.7352499E-02, 2.6084600E-02, + 5.3640924E-02, 2.2481044E-03, 7.6721571E-03, 6.4841054E-02, 0.9075497 , + 1.0222324E-02, 9.7214729E-02, 5.9108265E-02, 0.1335209 , 1.0431458E-04, + 1.0925709E-02, 8.2079172E-02, 2.7492193E-03, 1.0652039E-03, 2.5291166E-03, + 5.5883010E-04, 5.9249125E-02, 3.5640191E-05, 0.2042139 , 4.0082157E-02, + 4.4598214E-06, 1.1861410E-02, 1.1110056E-02, 0.1692510 , 0.2043316 , + 0.1057394 , 0.8772833 , 7.5677007E-02, 0.1414347 , 6.1683334E-02, + 6.2691003E-02, 2.3086492E-05, + SAIMPI = 5.6982312E-02, 1.5273975E-04, 1.9133480E-02, 5.9587371E-02, + 1.3435086E-02, 2*0.0000000E+00, 2.0074593E-03, 0.0000000E+00, 0.1366879 + , 0.0000000E+00, 9.3187485E-03, 0.1921618 , 8.9539858E-03, + 2.4765756E-02, 3.8349528E-03, 8.0762355E-04, 1.4646987E-03, 1.2814005E-02, + 0.0000000E+00, 2.7052544E-02, 0.9126120 , 9.5008463E-02, 7.7443495E-02, + 3.6042441E-02, 0.9349735 , 2.3855878E-02, 7.0651330E-02, 5.1497202E-03, + 9.4857492E-02, 0.1292573 , 6.3288331E-02, 0.3360588 , 5.6843691E-02, + 1.4573755E-02, 0.2399018 , 3.9295730E-04, 0.1100582 , 7.0646964E-02, + 0.1019014 , 4.6731960E-02, 0.1775764 , 6.5976396E-02, 1.7087533E-03, + 4.4726793E-04, 0.4480167 , 1.979766 , 5.8638420E-02, 0.2033015 , + 0.1795525 , 1.2593145E-02, 5.1588966E-03, 7.4997754E-03, 7.8012250E-02, + 1.2977102E-02, 6.0126837E-02, 0.1309072 , 1.6650161E-02, 0.1616113 , + 0.3601822 , + SAIPR = 0.0000000E+00, + SAIIP = 0.0000000E+00 + $END + LAO 10-SEP-1990 11:14:50 066832.02384 diff --git a/pyexamples/d3dHsmNew/rd_d3dHsm_in.py b/pyexamples/d3dHsmNew/rd_d3dHsm_in.py new file mode 100644 index 00000000..e8f43940 --- /dev/null +++ b/pyexamples/d3dHsmNew/rd_d3dHsm_in.py @@ -0,0 +1,92 @@ +# +# +########################################################################### +# DESCRIPTION OF PROBLEM (d3dHsm) from FACETS test suite: +# DIII-D single-null geometry with 5 variables (ni,upi,te,ti,ng) and a +# (16+2)*(8+2)=18x10 [poloidal*radial] mesh yielding 900 variables. +# Solver used is Newton Krylov (svrpkg="nksol") and preconditioner uses a +# direct banded solver for the LU decomposition (premeth="banded"). Iterates +# to steady-state solution from an initial profile file (HF5). +########################################################################### +###import uedge + +from uedge import * +from uedge.hdf5 import * +from uedge.rundt import * + + +# Set the geometry +bbb.mhdgeo = 1 #=1 use MHD equilibrium files +#flx.aeqdskfname = "aeqdskd3d" #name of EFIT 'a' file for flux-surface mesh +#flx.geqdskfname = "neqdskd3d" #name of EFIT 'g' or 'n' file for flux-sur mesh +flx.psi0min1 = 0.98 #normalized flux on core bndry +flx.psi0min2 = 0.98 #normalized flux on pf bndry +flx.psi0sep = 1.00001 #normalized flux at separatrix +flx.psi0max = 1.07 #normalized flux on outer wall bndry +bbb.ngrid = 1 #number of mesh sequenc. (always set to 1) +com.nxleg[0,0] = 4 #pol. mesh pts from inner plate to x-point +com.nxcore[0,0] = 4 #pol. mesh pts from x-point to top on inside +com.nxcore[0,1] = 4 #pol. mesh pts from top to x-point on outside +com.nxleg[0,1] = 4 #pol. mesh pts from x-point to outer plate +com.nysol[0] = 6 #rad. mesh pts in SOL +com.nycore[0] = 2 #rad. mesh pts in core + +# Finite-difference algorithms (upwind, central diff, etc.) +bbb.methn = 33 #ion continuty eqn +bbb.methu = 33 #ion parallel momentum eqn +bbb.methe = 33 #electron energy eqn +bbb.methi = 33 #ion energy eqn +bbb.methg = 33 #neutral gas continuity eqn + +# Boundary conditions +bbb.ncore[0] = 2.5e19 #hydrogen ion density on core +## iflcore = 0 #flag; =0, fixed Te,i; =1, fixed power on core +bbb.tcoree = 100. #core Te +bbb.tcorei = 100. #core Ti +bbb.tedge = 2. #fixed wall,pf Te,i if istewcon=1, etc +bbb.recycp[0] = 0.8 #hydrogen recycling coeff at plates + +# Transport coefficients (m**2/s) +bbb.difni[0] = 1. #D for radial hydrogen diffusion +bbb.kye = 1. #chi_e for radial elec energy diffusion +bbb.kyi = 1. #chi_i for radial ion energy diffusion +bbb.travis[0] = 1. #eta_a for radial ion momentum diffusion + +# Flux limits +bbb.flalfe = 0.21 #electron parallel thermal conduct. coeff +bbb.flalfi = 0.21 #ion parallel thermal conduct. coeff +bbb.flalfv = 1. #ion parallel viscosity coeff +bbb.flalfgx = 1.e20 #neut. gas in poloidal direction +bbb.flalfgy = 1.e20 #neut. gas in radial direction + +# Solver package +bbb.svrpkg = "nksol" #Newton solver using Krylov method +bbb.premeth = "banded" #Solution method for precond. Jacobian matrix + +# Restart from a HDF5 or PDB savefile +bbb.restart = 1 #Begin from savefile, not estimated profiles +bbb.allocate() #allocates storage for arrays + + +# Atomic data switches +com.istabon = 0 #-analytic rates +###com.istabon = 10 #=10 hydrogen data file ehr2.dat + + + +if (1): + hdf5_restore("d3dHsm.h5") + bbb.dtreal = 1e20; bbb.exmain() +else: + #-set up some simple initial state + bbb.ngs=1e14; bbb.ng=1e14 + bbb.nis=1e20; bbb.ni=1e20 + bbb.ups=0.0; bbb.up=0.0 + bbb.tes=bbb.ev; bbb.te=bbb.ev + bbb.tis=bbb.ev; bbb.ti=bbb.ev + + #-short run to initialize everything + bbb.dtreal = 1e-12; bbb.isbcwdt=1; bbb.exmain() + + ##-run to steady state + rundt(dtreal=1e-10) diff --git a/pyexamples/d3dHsmNew/runcase.py b/pyexamples/d3dHsmNew/runcase.py new file mode 100644 index 00000000..57f7e4c7 --- /dev/null +++ b/pyexamples/d3dHsmNew/runcase.py @@ -0,0 +1,14 @@ +import matplotlib.pyplot as plt +import numpy as np +from uedge import * +from uedge.uedgeplots import * + + + +exec(open("rd_d3dHsm_in.py").read()) + +#-show some results +plotmesh(z_min=0., z_max=3., r_min=0.8, r_max=2.5) +plotmeshval(bbb.te/bbb.ev, z_min=0., z_max=3., r_min=0.8, r_max=2.5, title="Te [ev]") +plotmeshval(bbb.ni, z_min=0., z_max=3., r_min=0.8, r_max=2.5, title="Ni [m-3]") +plotmeshval(bbb.ng, z_min=0., z_max=3., r_min=0.8, r_max=2.5, title="Ng [m-3]")