SLy.py

import numpy as np

SLyGs = np.array([
[ 10.01701764164203, 1.2896183301369382  ],
[ 10.057401977414822, 1.285214691743150  ],
[ 10.097786313187612, 1.283672063043432  ],
[ 10.138170648960402, 1.281869342553345  ],
[ 10.178554984733195, 1.280586805643998  ],
[ 10.218939320505985, 1.278523993363540  ],
[ 10.259323656278777, 1.277111410559008  ],
[ 10.299707992051568, 1.275178644173736  ],
[ 10.340092327824358, 1.272985785998093  ],
[ 10.38047666359715, 1.2707929278224515  ],
[ 10.42086099936994, 1.2688601614371784  ],
[ 10.461245335142733, 1.267187486842275  ],
[ 10.501629670915523, 1.264994628666633  ],
[ 10.542014006688314, 1.262281586910250  ],
[ 10.582398342461106, 1.260218774629793  ],
[ 10.622782678233897, 1.257895870558965  ],
[ 10.663167014006689, 1.255442920592953  ],
[ 10.70355134977948, 1.2525998329413857  ],
[ 10.74393568555227, 1.2500168370801887  ],
[ 10.784320021325062, 1.247303795323805  ],
[ 10.824704357097852, 1.244330661777053  ],
[ 10.865088692870644, 1.241487574125485  ],
[ 10.905473028643435, 1.238644486473918  ],
[ 10.945857364416227, 1.235151169346425  ],
[ 10.986241700189018, 1.231397760428563  ],
[ 11.026626035961808, 1.228684718672180  ],
[ 11.0670103717346, 1.2245411720687631   ],
[ 11.10739470750739, 1.2211779008364556  ],
[ 11.147779043280181, 1.217424491918592  ],
[ 11.188163379052973, 1.212760761734435  ],
[ 11.228547714825764, 1.208877306921387  ],
[ 11.268932050598556, 1.204343622632415  ],
[ 11.309316386371346, 1.199419800657887  ],
[ 11.349700722144139, 1.194886116368914  ],
[ 11.390085057916929, 1.188531789547351  ],
[ 11.43046939368972, 1.184388242943934   ],
[ 11.470853729462512, 1.177253640751260  ],
[ 11.511238065235302, 1.172459864671918  ],
[ 11.540608491251877, 1.165294365470434  ],
[ 11.54427979450395, 1.1269584648464845  ],
[ 11.54427979450395, 1.0954873582117117  ],
[ 11.54427979450395, 1.064016251576939   ],
[ 11.546115446129985, 1.024964283891176  ],
[ 11.546115446129985, 0.993922328710513  ],
[ 11.545503562254641, 0.956728931141502  ],
[ 11.546727330005332, 0.904277629842411  ],
[ 11.546727330005332, 0.847534270910018  ],
[ 11.546115446129985, 0.802949931631324  ],
[ 11.546727330005332, 0.754551455852734  ],
[ 11.54886892356904, 0.707822581260932   ],
[ 11.547951097756021, 0.669675378020847  ],
[ 11.550398633257403, 0.625807648896163  ],
[ 11.552846168758784, 0.581463084822469  ],
[ 11.560800659138273, 0.534021537495085  ],
[ 11.579157175398633, 0.491472164678746  ],
[ 11.606691949789173, 0.458290169507205  ],
[ 11.647076285561964, 0.446994098668612  ],
[ 11.687460621334754, 0.456245279269250  ],
[ 11.727844957107546, 0.475509993799135  ],
[ 11.768229292880337, 0.501927232564195  ],
[ 11.808613628653129, 0.533416261241471  ],
[ 11.845326661173848, 0.563791045761127  ],
[ 11.878368390442494, 0.594482090970201  ],
[ 11.909574468085106, 0.624655750346216  ],
[ 11.940780545727717, 0.656657277026777  ],
[ 11.971986623370329, 0.689056166164847  ],
[ 12.001357049386904, 0.720897933224108  ],
[ 12.03072747540348, 0.7532761396010068  ],
[ 12.060097901420054, 0.786190785295544  ],
[ 12.08946832743663, 0.8189266178842018  ],
[ 12.118838753453206, 0.851662450472859  ],
[ 12.14820917946978, 0.8840406568497583  ],
[ 12.177579605486358, 0.916776489438416  ],
[ 12.206950031502933, 0.948260630285918  ],
[ 12.236320457519508, 0.979387144921661  ],
[ 12.26752653516212, 1.010772759793082   ],
[ 12.300568264430765, 1.043847979747214  ],
[ 12.333609993699413, 1.074856914922295  ],
[ 12.368487374594096, 1.105866664735672  ],
[ 12.405200407114815, 1.135923559289321  ],
[ 12.443749091261571, 1.166761543852376  ],
[ 12.484133427034362, 1.195909746416322  ],
[ 12.524517762807154, 1.222066893391011  ],
[ 12.564902098579944, 1.245232984776446  ],
[ 12.605286434352735, 1.266188295943735  ],
[ 12.645670770125527, 1.284412643312140  ],
[ 12.686055105898319, 1.300166118672029  ],
[ 12.72643944167111, 1.313448722023403   ],
[ 12.7668237774439, 1.324390499261447    ],
[ 12.807208113216692, 1.333771725757270  ],
[ 12.847592448989483, 1.340942172034949  ],
[ 12.887976784762273, 1.346422021675221  ],
[ 12.928361120535065, 1.350341320573275  ],
[ 12.968745456307856, 1.352439976938737  ],
[ 13.009129792080648, 1.353628312037905  ],
[ 13.049514127853438, 1.353646234080408  ],
[ 13.08989846362623, 1.351843513590321   ],
[ 13.130282799399021, 1.350430930785788  ],
[ 13.170667135171811, 1.347327751343851  ],
[ 13.211051470944604, 1.343964480111543  ],
[ 13.251435806717394, 1.340471162984051  ],
[ 13.291820142490185, 1.336067524590263  ],
[ 13.332204478262977, 1.332184069777215  ],
[ 13.372588814035767, 1.328040523173797  ],
[ 13.41297314980856, 1.3248072978366756  ],
[ 13.45335748558135, 1.3220942560802929  ],
[ 13.493741821354142, 1.319901397904650  ],
[ 13.534126157126932, 1.320049365842337  ],
[ 13.574510492899723, 1.322538159893355  ],
[ 13.614894828672515, 1.327107688267334  ],
[ 13.655279164445306, 1.334668272230567  ],
[ 13.695663500218096, 1.345870141258981  ],
[ 13.736047835990888, 1.36162361661887   ],
[ 13.77643217176368, 1.3830991113668998  ],
[ 13.81681650753647, 1.4097764419223298  ],
[ 13.855365191683227, 1.441121605476606  ],
[ 13.888406920951873, 1.472571612979523  ],
[ 13.915941695342411, 1.502394531776993  ],
[ 13.941640818106915, 1.532855254171450  ],
[ 13.967339940871419, 1.564950839248233  ],
[ 13.993039063635923, 1.59806821350147   ],
[ 14.018738186400427, 1.631798661260578  ],
[ 14.046272960790965, 1.665044573799167  ],
[ 14.08115034168565, 1.684680726362822   ],
[ 14.107583725100568, 1.649429212215911  ],
[ 14.118322287112878, 2.190429453425919  ],
[ 14.112356419328261, 2.070085585594628  ],
[ 14.112968303203605, 1.955883886852422  ],
[ 14.112356419328261, 1.91702156696187   ],
[ 14.112968303203605, 1.831906800109378  ],
[ 14.112356419328261, 1.793998150116849  ],
[ 14.112968303203605, 1.708883383264357  ],
[ 14.110520767702225, 1.685994219527401  ],
[ 14.112968303203605, 1.614470063360039  ],
[ 14.111744535452914, 1.573938549601847  ],
[ 14.114192070954296, 2.160923458019654  ],
[ 14.114192070954296, 2.129452351384881  ],
[ 14.114192070954296, 2.035039031480563  ],
[ 14.114192070954296, 2.003567924845790  ],
[ 14.114192070954296, 1.877683498306699  ],
[ 14.114192070954296, 1.751799071767608  ],
[ 14.136219890466727, 2.241518340062180  ],
[ 14.147233800222944, 2.280623693710914  ],
[ 14.15824770997916, 2.317344872614589   ],
[ 14.167425968109338, 2.350488974762380  ],
[ 14.176604226239519, 2.381725737114125  ],
[ 14.187618135995734, 2.41558590632373   ],
[ 14.19863204575195, 2.4527839201764166  ],
[ 14.209645955508165, 2.487597759284045  ],
[ 14.220659865264382, 2.520504258595626  ],
[ 14.231673775020598, 2.553410757907207  ],
[ 14.242687684776815, 2.586317257218788  ],
[ 14.255537246159065, 2.621608745913724  ],
[ 14.270222459167353, 2.660596520101796  ],
[ 14.296708289771587, 2.702552719862055  ],
[ 14.31244244656618, 2.7219481520992255  ],
[ 14.341812872582755, 2.746959258513893  ],
[ 14.360169388843115, 2.781871723164506  ],
[ 14.378525905103475, 2.815353682968084  ],
[ 14.396882421363834, 2.848263440832847  ],
[ 14.41707458925023, 2.882032316244128   ],
[ 14.439102408762661, 2.914705285910977  ],
[ 14.46296587990113, 2.9464254003180983  ],
[ 14.49417195754374, 2.9786256082274134  ],
[ 14.532720641690496, 3.007707973205470  ],
[ 14.573104977463288, 3.026452504154615  ],
[ 14.613489313236078, 3.035443592964883  ],
[ 14.653873649008869, 3.036631928064051  ],
[ 14.694257984781661, 3.032358335565448  ],
[ 14.734642320554451, 3.024183366211295  ],
[ 14.775026656327242, 3.012757249477517  ],
[ 14.815410992100034, 3.000030673791888  ],
[ 14.855795327872825, 2.986133685049595  ],
[ 14.896179663645615, 2.972236696307301  ],
[ 14.936563999418407, 2.957949569879453  ],
[ 14.9769483351912, 2.94405258113716     ],
[ 15.01733267096399, 2.9306757759756064  ],
[ 15.05771700673678, 2.917298970814053   ],
[ 15.098101342509572, 2.904442349233239  ],
[ 15.138485678282363, 2.891325635862056  ],
[ 15.178870014055153, 2.878859151966797  ],
[ 15.219254349827946, 2.865352300910059  ],
[ 15.259638685600738, 2.852365633434060  ],
[ 15.300023021373526, 2.838208552901397  ],
[ 15.340407357146319, 2.823661334683178  ],
[ 15.38079169291911, 2.808723978779405   ],
[ 15.421176028691901, 2.792096026238226  ],
[ 15.461560364464692, 2.775207981906677  ],
[ 15.501944700237484, 2.757409616308833  ],
[ 15.542329036010274, 2.738961021235064  ],
[ 15.582713371783065, 2.720382380266111  ],
[ 15.623097707555857, 2.701933785192341  ],
[ 15.661646391702611, 2.684745820663572  ],
])


SLyPs = np.array([
[ 10.06937076805206, 27.971835449218744],
[ 10.104813343274133, 28.0130527565696 ],
[ 10.140258626044474, 28.05874889026988],
[ 10.175703725872365, 28.10414240056818],
[ 10.211148459815355, 28.14893066406249],
[ 10.246593742585699, 28.19462679776278],
[ 10.282038659471139, 28.23971768465908],
[ 10.31748357635658, 28.284808571555395],
[ 10.352928493242022, 28.3298994584517 ],
[ 10.388373593069913, 28.37529296874999],
[ 10.423818144070452, 28.41977860884233],
[ 10.459262878013442, 28.46456687233664],
[ 10.494707794898883, 28.50965775923295],
[ 10.530152528841873, 28.55444602272727],
[ 10.565597079842412, 28.5989316628196 ],
[ 10.60104126495805, 28.64281205610795 ],
[ 10.63648563301614, 28.68699507279829 ],
[ 10.671929818131778, 28.73087546608664],
[ 10.707374735017218, 28.77596635298295],
[ 10.742819103075307, 28.82014936967329],
[ 10.778263471133396, 28.86433238636363],
[ 10.813708205076386, 28.90912064985795],
[ 10.849152573134475, 28.95330366654829],
[ 10.88459584353786, 28.9956709428267  ],
[ 10.920040211595948, 29.03985395951704],
[ 10.955484396711585, 29.08373435280539],
[ 10.990928398884773, 29.12731212269175],
[ 11.026372035173058, 29.17028464577414],
[ 11.061816037346246, 29.21386241566050],
[ 11.097259673634532, 29.25683493874289],
[ 11.132703492865268, 29.30011008522727],
[ 11.168147312096004, 29.34338523171164],
[ 11.20359094838429, 29.38635775479403 ],
[ 11.239034767615028, 29.42963290127840],
[ 11.274478220960862, 29.4723028009588 ],
[ 11.309921308421796, 29.51436745383522],
[ 11.345364029997828, 29.55582685990767],
[ 11.380807666286113, 29.59879938299005],
[ 11.416251119631948, 29.64146928267045],
[ 11.45169365826553, 29.682626065340905],
[ 11.487136379841562, 29.72408547141334],
[ 11.522579101417595, 29.76554487748579],
[ 11.556275911715137, 29.78685654860276],
[ 11.593446616209487, 29.81880659623579],
[ 11.62887397061966, 29.83484563654119 ],
[ 11.664301873857182, 29.85179254705255],
[ 11.699728862382452, 29.86722634055397],
[ 11.735156948562427, 29.88447587446732],
[ 11.77058485179995, 29.90142278497869 ],
[ 11.806013669749728, 29.91988281249999],
[ 11.841442853584406, 29.93894808682527],
[ 11.876873135073788, 29.95982910156249],
[ 11.912303416563173, 29.98071011629971],
[ 11.947734429822358, 30.00280162464488],
[ 11.983167089563599, 30.02761674360795],
[ 12.018599383419941, 30.05182661576704],
[ 12.054032957873437, 30.07815485174005],
[ 12.089467995866539, 30.10690407492897],
[ 12.124903765629442, 30.13686379172584],
[ 12.160339718334798, 30.16712613192471],
[ 12.195776768694857, 30.19920421253550],
[ 12.231214916709622, 30.23309803355823],
[ 12.266653613551739, 30.26789972478693],
[ 12.30209340804856, 30.30451715642755 ],
[ 12.337533568430281, 30.34173983487215],
[ 12.372974826466708, 30.38077825372869],
[ 12.408416267445585, 30.42011929598721],
[ 12.443858623136716, 30.46097345525567],
[ 12.479301344712747, 30.50243286132812],
[ 12.514745163943484, 30.54570800781249],
[ 12.550188434346868, 30.58807528409090],
[ 12.585632436520056, 30.63165305397726],
[ 12.621077353405497, 30.67674394087357],
[ 12.656521721463585, 30.72092695756392],
[ 12.691967004233927, 30.7666230912642 ],
[ 12.727413201716523, 30.81383234197443],
[ 12.76285958214157, 30.861344216086643],
[ 12.798305047854363, 30.90734297318891],
[ 12.83375161122186, 30.95515747070312 ],
[ 12.869197625762006, 31.00206409801136],
[ 12.904643640302151, 31.0489707253196 ],
[ 12.940090935439452, 31.09799571644176],
[ 12.975536949979597, 31.14490234374999],
[ 13.010983330404644, 31.19241421786221],
[ 13.046430076714591, 31.24053133877840],
[ 13.08187682302454, 31.2886484596946  ],
[ 13.117323020507136, 31.33585771040482],
[ 13.152769949759534, 31.38427745472301],
[ 13.18821614724213, 31.431486705433233],
[ 13.223662344724726, 31.47869595614346],
[ 13.259108908092223, 31.52651045365766],
[ 13.294555105574819, 31.57371970436789],
[ 13.330001303057415, 31.62092895507812],
[ 13.365446768770209, 31.66692771218039],
[ 13.400893332137706, 31.7147422096946 ],
[ 13.436338614908049, 31.76043834339488],
[ 13.471784629448194, 31.80734497070312],
[ 13.50723082693079, 31.85455422141335 ],
[ 13.542676475586033, 31.90085560191761],
[ 13.578122124241277, 31.94715698242187],
[ 13.613568504666324, 31.99466885653408],
[ 13.649014336264019, 32.04127286044034],
[ 13.684460716689065, 32.08878473455255],
[ 13.71990691417166, 32.13599398526278 ],
[ 13.75535384342406, 32.18441372958096 ],
[ 13.790800406791558, 32.23222822709516],
[ 13.826247701928857, 32.28125321821732],
[ 13.861696643548214, 32.33300181995738],
[ 13.897145036340218, 32.38384255149147],
[ 13.932594526786929, 32.4364990234375 ],
[ 13.96804474900344, 32.49036598899147 ],
[ 14.003496068874655, 32.54604869495738],
[ 14.038948120515675, 32.60294189453125],
[ 14.0744016356963, 32.66225608132102  ],
[ 14.109852040855262, 32.71642567027698],
[ 14.134215137480306, 32.77567260742187],
[ 14.16626796186943, 32.83163439941406 ],
[ 14.198510708980123, 32.90586791992187],
[ 14.22914519477764, 32.98280151367187 ],
[ 14.25816777464182, 33.056406249999995],
[ 14.287191472487645, 33.13186035156249],
[ 14.314604410331318, 33.20588119506836],
[ 14.342016845083252, 33.27906982421874],
[ 14.369430090371878, 33.35359924316405],
[ 14.39523215641524, 33.42600189208984 ],
[ 14.421033719366859, 33.49757232666015],
[ 14.446836791593697, 33.57163940429687],
[ 14.472639863820536, 33.64570648193359],
[ 14.498443439139114, 33.72060577392577],
[ 14.524246511365954, 33.79467285156249],
[ 14.548438037010701, 33.86600551060267],
[ 14.572629598590572, 33.93739761352538],
[ 14.598433928546758, 34.01354522705078],
[ 14.624238258502945, 34.08969284057617],
[ 14.650041833821525, 34.16459213256836],
[ 14.675845409140102, 34.23949142456054],
[ 14.701648984458682, 34.31439071655273],
[ 14.72745230823139, 34.388873901367184],
[ 14.753255632004098, 34.46335708618164],
[ 14.779057698047458, 34.53575973510742],
[ 14.804860015636688, 34.60857849121093],
[ 14.830661578588309, 34.68014892578124],
[ 14.856463141539928, 34.75171936035156],
[ 14.88226319521633, 34.82079315185546 ],
[ 14.908064003530342, 34.89111526489257],
[ 14.935475487997879, 34.96273193359374],
[ 14.96450008022902, 35.03966552734374 ],
[ 14.993522436496871, 35.11290039062499],
[ 15.02254501636105, 35.18650512695312 ],
[ 15.051566925436246, 35.25900024414062],
[ 15.080587940126126, 35.33001586914062],
[ 15.11122014524109, 35.403176757812496],
[ 15.14346349606187, 35.478408935546874],
[ 15.174093912406203, 35.54861083984375],
[ 15.206336279403136, 35.62221557617187],
[ 15.241798209936496, 35.69545043945312],
[ 15.27565145432365, 35.7706826171875  ],
[ 15.309501552100656, 35.84070967240767],
[ 15.344960372612352, 35.90879993785511],
[ 15.380421754318359, 35.98112693093039],
[ 15.415879843060251, 36.04800670276988],
[ 15.45133646826254, 36.112465487393465],
[ 15.486794008177078, 36.17843738902698],
[ 15.522248438069958, 36.2392646928267 ],
[ 15.557702502077936, 36.29948674982244],
[ 15.593155285488756, 36.35759044300425],
[ 15.626995431196441, 36.41115478515625]
])



#structure and composition of the inner crust
#nb Z A Xn Rp Rn Rcell u
#1.2126e-4, 42.198, 130.076, 0.0000, 5.451,  5.915,  63.503, 0.063,
#1.6241e-4, 42.698, 135.750, 0.0000, 5.518,  6.016,  58.440, 0.084,
#1.9772e-4, 43.019, 139.956, 0.0000, 5.565,  6.089,  55.287, 0.102,
#2.0905e-4, 43.106, 141.564, 0.0000, 5.578,  6.111,  54.470, 0.107,
#2.2059e-4, 43.140, 142.161, 0.0247, 5.585,  6.122,  54.032, 0.110,
#2.3114e-4, 43.163, 142.562, 0.0513, 5.590,  6.128,  53.745, 0.113,
#2.6426e-4, 43.215, 143.530, 0.1299, 5.601,  6.145,  53.020, 0.118,
#3.0533e-4, 43.265, 144.490, 0.2107, 5.612,  6.162,  52.312, 0.123,
#3.5331e-4, 43.313, 145.444, 0.2853, 5.623,  6.179,  51.617, 0.129,
#4.0764e-4, 43.359, 146.398, 0.3512, 5.634,  6.195,  50.937, 0.135,
#4.6800e-4, 43.404, 147.351, 0.4082, 5.645,  6.212,  50.269, 0.142,
#5.3414e-4, 43.447, 148.306, 0.4573, 5.656,  6.228,  49.615, 0.148,
#6.0594e-4, 43.490, 149.263, 0.4994, 5.667,  6.245,  48.974, 0.155,
#7.6608e-4, 43.571, 151.184, 0.5669, 5.690,  6.278,  47.736, 0.169,
#1.0471e-3, 43.685, 154.094, 0.6384, 5.725,  6.328,  45.972, 0.193,
#1.2616e-3, 43.755, 156.055, 0.6727, 5.748,  6.362,  44.847, 0.211,
#1.6246e-3, 43.851, 159.030, 0.7111, 5.784,  6.413,  43.245, 0.239,
#2.0384e-3, 43.935, 162.051, 0.7389, 5.821,  6.465,  41.732, 0.271,
#2.6726e-3, 44.030, 166.150, 0.7652, 5.871,  6.535,  39.835, 0.320,
#3.4064e-3, 44.101, 170.333, 0.7836, 5.923,  6.606,  38.068, 0.377,
#4.4746e-3, 44.155, 175.678, 0.7994, 5.989,  6.698,  36.012, 0.460,
#5.7260e-3, 44.164, 181.144, 0.8099, 6.059,  6.792,  34.122, 0.560,
#7.4963e-3, 44.108, 187.838, 0.8179, 6.146,  6.908,  32.030, 0.706,
#9.9795e-3, 43.939, 195.775, 0.8231, 6.253,  7.048,  29.806, 0.923,
#1.2513e-2, 43.691, 202.614, 0.8250, 6.350,  7.171,  28.060, 1.159,
#1.6547e-2, 43.198, 211.641, 0.8249, 6.488,  7.341,  25.932, 1.566,
#2.1405e-2, 42.506, 220.400, 0.8222, 6.637,  7.516,  24.000, 2.115,
#2.4157e-2, 42.089, 224.660, 0.8200, 6.718,  7.606,  23.106, 2.458,
#2.7894e-2, 41.507, 229.922, 0.8164, 6.825,  7.721,  22.046, 2.967,
#3.1941e-2, 40.876, 235.253, 0.8116, 6.942,  7.840,  21.053, 3.585,
#3.6264e-2, 40.219, 240.924, 0.8055, 7.072,  7.967,  20.128, 4.337,
#3.9888e-2, 39.699, 245.999, 0.7994, 7.187,  8.077,  19.433, 5.058,
#4.4578e-2, 39.094, 253.566, 0.7900, 7.352,  8.231,  18.630, 6.146,
#4.8425e-2, 38.686, 261.185, 0.7806, 7.505,  8.372,  18.038, 7.202,
#5.2327e-2, 38.393, 270.963, 0.7693, 7.685,  8.538,  17.499, 8.470,
#5.6264e-2, 38.281, 283.993, 0.7553, 7.900,  8.737,  17.014, 10.011,
#6.0219e-2, 38.458, 302.074, 0.7381, 8.167,  8.987,  16.598, 11.914,
#6.4183e-2, 39.116, 328.489, 0.7163, 8.513,  9.315,  16.271, 14.323,
#6.7163e-2, 40.154, 357.685, 0.6958, 8.853,  9.642,  16.107, 16.606,
#7.0154e-2, 42.051, 401.652, 0.6699, 9.312,  10.088, 16.058, 19.501,
#7.3174e-2, 45.719, 476.253, 0.6354, 9.990,  10.753, 16.213, 23.393,
#7.5226e-2, 50.492, 566.654, 0.6038, 10.701, 11.456, 16.557, 26.996,
#7.5959e-2, 53.162, 615.840, 0.5898, 11.051, 11.803, 16.772, 28.603,
#
##composition of the liquid core
##nb xp xe xmu
#0.0771, 3.516,  3.516,  0.000,
#0.0800, 3.592,  3.592,  0.000,
#0.0850, 3.717,  3.717,  0.000,
#0.0900, 3.833,  3.833,  0.000,
#0.1000, 4.046,  4.046,  0.000,
#0.1100, 4.233,  4.233,  0.000,
#0.1200, 4.403,  4.398,  0.005,
#0.1300, 4.622,  4.521,  0.101,
#0.1600, 5.270,  4.760,  0.510,
#0.1900, 5.791,  4.896,  0.895,
#0.2200, 6.192,  4.973,  1.219,
#0.2500, 6.499,  5.014,  1.485,
#0.2800, 6.736,  5.031,  1.705,
#0.3100, 6.920,  5.034,  1.887,
#0.3400, 7.066,  5.026,  2.040,
#0.3700, 7.185,  5.014,  2.170,
#0.4000, 7.283,  4.999,  2.283,
#0.4300, 7.368,  4.984,  2.383,
#0.4600, 7.444,  4.971,  2.473,
#0.4900, 7.516,  4.960,  2.556,
#0.5200, 7.587,  4.954,  2.634,
#0.5500, 7.660,  4.952,  2.708,
#0.5800, 7.736,  4.955,  2.781,
#0.6100, 7.818,  4.964,  2.854,
#0.6400, 7.907,  4.979,  2.927,
#0.6700, 8.003,  5.001,  3.002,
#0.7000, 8.109,  5.030,  3.079,
#0.7500, 8.309,  5.094,  3.215,
#0.8000, 8.539,  5.178,  3.361,
#0.8500, 8.803,  5.284,  3.519,
#0.9000, 9.102,  5.410,  3.692,
#0.9500, 9.437,  5.557,  3.880,
#1.0000, 9.808,  5.726,  4.083,
#1.1000, 10.663, 6.124,  4.539,
#1.2000, 11.661, 6.602,  5.060,
#1.3000, 12.794, 7.151,  5.643,
#1.4000, 14.043, 7.762,  6.281,
#1.5000, 15.389, 8.424,  6.965,
#
#
##EoS for the inner crust
##nb rho P Gamma
#2.0905e-4,  3.4951e+11,  6.2150e+29,  1.177,
#2.2059e-4,  3.6883e+11,  6.4304e+29,  0.527,
#2.3114e-4,  3.8650e+11,  6.5813e+29,  0.476,
#2.6426e-4,  4.4199e+11,  6.9945e+29,  0.447,
#3.0533e-4,  5.1080e+11,  7.4685e+29,  0.466,
#3.5331e-4,  5.9119e+11,  8.0149e+29,  0.504,
#4.0764e-4,  6.8224e+11,  8.6443e+29,  0.554,
#4.6800e-4,  7.8339e+11,  9.3667e+29,  0.610,
#5.3414e-4,  8.9426e+11,  1.0191e+30,  0.668,
#6.0594e-4,  1.0146e+12,  1.1128e+30,  0.726,
#7.6608e-4,  1.2831e+12,  1.3370e+30,  0.840,
#1.0471e-3,  1.7543e+12,  1.7792e+30,  0.987,
#1.2616e-3,  2.1141e+12,  2.1547e+30,  1.067,
#1.6246e-3,  2.7232e+12,  2.8565e+30,  1.160,
#2.0384e-3,  3.4178e+12,  3.7461e+30,  1.227,
#2.6726e-3,  4.4827e+12,  5.2679e+30,  1.286,
#3.4064e-3,  5.7153e+12,  7.2304e+30,  1.322,
#4.4746e-3,  7.5106e+12,  1.0405e+31,  1.344,
#5.7260e-3,  9.6148e+12,  1.4513e+31,  1.353,
#7.4963e-3,  1.2593e+13,  2.0894e+31,  1.351,
#9.9795e-3,  1.6774e+13,  3.0720e+31,  1.342,
#1.2513e-2,  2.1042e+13,  4.1574e+31,  1.332,
#1.6547e-2,  2.7844e+13,  6.0234e+31,  1.322,
#2.1405e-2,  3.6043e+13,  8.4613e+31,  1.320,
#2.4157e-2,  4.0688e+13,  9.9286e+31,  1.325,
#2.7894e-2,  4.7001e+13,  1.2023e+32,  1.338,
#3.1941e-2,  5.3843e+13,  1.4430e+32,  1.358,
#3.6264e-2,  6.1153e+13,  1.7175e+32,  1.387,
#3.9888e-2,  6.7284e+13,  1.9626e+32,  1.416,
#4.4578e-2,  7.5224e+13,  2.3024e+32,  1.458,
#4.8425e-2,  8.1738e+13,  2.6018e+32,  1.496,
#5.2327e-2,  8.8350e+13,  2.9261e+32,  1.536,
#5.6264e-2,  9.5022e+13,  3.2756e+32,  1.576,
#6.0219e-2,  1.0173e+14,  3.6505e+32,  1.615,
#6.4183e-2,  1.0845e+14,  4.0509e+32,  1.650,
#6.7163e-2,  1.1351e+14,  4.3681e+32,  1.672,
#7.0154e-2,  1.1859e+14,  4.6998e+32,  1.686,
#7.3174e-2,  1.2372e+14,  5.0462e+32,  1.685,
#7.5226e-2,  1.2720e+14,  5.2856e+32,  1.662,
#7.5959e-2,  1.2845e+14,  5.3739e+32,  1.644,
#
##eos for the liquid crust
##nb rho P Gamma
#0.0771,  1.3038e+14,  5.3739e+32,  2.159,  
#0.0800,  1.3531e+14,  5.8260e+32,  2.217,  
#0.0850,  1.4381e+14,  6.6828e+32,  2.309,  
#0.0900,  1.5232e+14,  7.6443e+32,  2.394,  
#0.1000,  1.6935e+14,  9.9146e+32,  2.539,  
#0.1100,  1.8641e+14,  1.2701e+33,  2.655,  
#0.1200,  2.0350e+14,  1.6063e+33,  2.708,  
#0.1300,  2.2063e+14,  1.9971e+33,  2.746,  
#0.1600,  2.7223e+14,  3.5927e+33,  2.905,  
#0.1900,  3.2424e+14,  5.9667e+33,  2.990,  
#0.2200,  3.7675e+14,  9.2766e+33,  3.025,  
#0.2500,  4.2983e+14,  1.3668e+34,  3.035,  
#0.2800,  4.8358e+14,  1.9277e+34,  3.032,  
#0.3100,  5.3808e+14,  2.6235e+34,  3.023,  
#0.3400,  5.9340e+14,  3.4670e+34,  3.012,  
#0.3700,  6.4963e+14,  4.4702e+34,  2.999,  
#0.4000,  7.0684e+14,  5.6451e+34,  2.987,  
#0.4300,  7.6510e+14,  7.0033e+34,  2.975,  
#0.4600,  8.2450e+14,  8.5561e+34,  2.964,  
#0.4900,  8.8509e+14,  1.0315e+35,  2.953,  
#0.5200,  9.4695e+14,  1.2289e+35,  2.943,  
#0.5500,  1.0102e+15,  1.4491e+35,  2.933,  
#0.5800,  1.0748e+15,  1.6930e+35,  2.924,  
#0.6100,  1.1408e+15,  1.9616e+35,  2.916,  
#0.6400,  1.2085e+15,  2.2559e+35,  2.908,  
#0.6700,  1.2777e+15,  2.5769e+35,  2.900,  
#0.7000,  1.3486e+15,  2.9255e+35,  2.893,  
#0.7500,  1.4706e+15,  3.5702e+35,  2.881,  
#0.8000,  1.5977e+15,  4.2981e+35,  2.869,  
#0.8500,  1.7302e+15,  5.1129e+35,  2.858,  
#0.9000,  1.8683e+15,  6.0183e+35,  2.847,  
#0.9500,  2.0123e+15,  7.0176e+35,  2.836,  
#1.0000,  2.1624e+15,  8.1139e+35,  2.824,  
#1.1000,  2.4820e+15,  1.0609e+36,  2.801,  
#1.2000,  2.8289e+15,  1.3524e+36,  2.778,  
#1.3000,  3.2048e+15,  1.6876e+36,  2.754,  
#1.4000,  3.6113e+15,  2.0679e+36,  2.731,  
#1.5000,  4.0498e+15,  2.4947e+36,  2.708,