PetThermoTools 0.2.38__py3-none-any.whl → 0.2.39__py3-none-any.whl

PetThermoTools/Melting.py

@@ -49,14 +49,15 @@ def AdiabaticDecompressionMelting(cores = multiprocessing.cpu_count(),
                                   P_path_bar = None, Frac = False, prop = None, 
                                   fO2_buffer = None, fO2_offset = None, Fe3Fet = None, MELTS_filter = True):
     
-    try:
-        import pyMelt as m
-        Lithologies = {'KLB-1': m.lithologies.matthews.klb1(),
-                    'KG1': m.lithologies.matthews.kg1(),
-                    'G2': m.lithologies.matthews.eclogite(),
-                    'hz': m.lithologies.shorttle.harzburgite()}
-    except ImportError:
-        raise RuntimeError('You havent installed pyMelt or there is an error when importing pyMelt. pyMelt is currently required to estimate the starting point for the melting calculations.')
+    if Tp_Method == "pyMelt":
+        try:
+            import pyMelt as m
+            Lithologies = {'KLB-1': m.lithologies.matthews.klb1(),
+                        'KG1': m.lithologies.matthews.kg1(),
+                        'G2': m.lithologies.matthews.eclogite(),
+                        'hz': m.lithologies.shorttle.harzburgite()}
+        except ImportError:
+            raise RuntimeError('You havent installed pyMelt or there is an error when importing pyMelt. pyMelt is currently required to estimate the starting point for the melting calculations.')
 
     if bulk is not None and comp_lith_1 is None:
         comp_lith_1 = bulk
@@ -155,7 +156,9 @@ def AdiabaticDecompressionMelting(cores = multiprocessing.cpu_count(),
 
     return Results
 
-def AdiabaticMelt(q, index, *, Model = None, comp_1 = None, comp_2 = None, comp_3 = None, Tp_C = None, P_start_bar = None, P_end_bar = None, dp_bar = None, P_path_bar = None, Frac = None, fO2_buffer = None, fO2_offset = None, prop = None):
+def AdiabaticMelt(q, index, *, Model = None, comp_1 = None, comp_2 = None, comp_3 = None, 
+                  Tp_C = None, P_start_bar = None, P_end_bar = None, dp_bar = None, P_path_bar = None, 
+                  Frac = None, fO2_buffer = None, fO2_offset = None, prop = None):
     '''
     Melting calculations to be performed in parallel.
 
@@ -163,21 +166,16 @@ def AdiabaticMelt(q, index, *, Model = None, comp_1 = None, comp_2 = None, comp_
     Results = {}
     if "MELTS" in Model:
         try:
-            Results = AdiabaticDecompressionMelting_MELTS(Model = Model, comp = comp_1, Tp_C = Tp_C, P_path_bar = P_path_bar, P_start_bar = P_start_bar, P_end_bar = P_end_bar, dp_bar = dp_bar, fO2_buffer = fO2_buffer, fO2_offset = fO2_offset)
+            Results = AdiabaticDecompressionMelting_MELTS(Model = Model, comp = comp_1, Tp_C = Tp_C, 
+                                                          P_path_bar = P_path_bar, P_start_bar = P_start_bar, P_end_bar = P_end_bar, dp_bar = dp_bar, 
+                                                          fO2_buffer = fO2_buffer, fO2_offset = fO2_offset)
             q.put([Results, index])
         except:
             q.put([])
 
         return
-
-    if Model == "Holland":
-        # import pyMAGEMINcalc as MM
-        # Results = MM.AdiabaticDecompressionMelting(comp = comp_1, T_p_C = Tp_C, P_start_kbar = P_start_bar/1000, P_end_kbar = P_end_bar/1000, dp_kbar = dp_bar/1000, Frac = 0)
-        print('Note that the ability to use MAGEMin to performed adiabatic decompression melting in PetThermoTools has been temporarily disabled. The underlying issue will be fixed soon and this funciton will once again become available.')
-        q.put([Results, index])
-        return
     
-    if Model == "pyMelt":
+    elif Model == "pyMelt":
         try:
             import pyMelt as m
             Lithologies = {'KLB-1': m.lithologies.matthews.klb1(),
@@ -239,5 +237,44 @@ def AdiabaticMelt(q, index, *, Model = None, comp_1 = None, comp_2 = None, comp_
         
         q.put([Results, index])
         return
+
+    else:
+        # import pyMAGEMINcalc as MM
+        # Results = MM.AdiabaticDecompressionMelting(comp = comp_1, T_p_C = Tp_C, P_start_kbar = P_start_bar/1000, P_end_kbar = P_end_bar/1000, dp_kbar = dp_bar/1000, Frac = 0)
+        # print('Note that the ability to use MAGEMin to performed adiabatic decompression melting in PetThermoTools has been temporarily disabled. The underlying issue will be fixed soon and this funciton will once again become available.')
+        
+        try:
+            import pyMelt as m
+            Lithologies = {'KLB-1': m.lithologies.matthews.klb1(),
+                        'KG1': m.lithologies.matthews.kg1(),
+                        'G2': m.lithologies.matthews.eclogite(),
+                        'hz': m.lithologies.shorttle.harzburgite()}
+        except ImportError:
+            raise RuntimeError('You havent installed pyMelt or there is an error when importing pyMelt. pyMelt is currently required to estimate the starting point for the melting calculations.')
+
+        lz = m.lithologies.matthews.klb1()
+        mantle = m.mantle([lz], [1], ['Lz'])
+        T_start_C = mantle.adiabat(P_start_bar/10000.0, Tp_C)
+        
+        from juliacall import Main as jl, convert as jlconvert
+
+        jl.seval("using MAGEMinCalc")
+
+        comp_1['O'] = comp_1['Fe3Fet_Liq']*(((159.59/2)/71.844)*comp_1['FeOt_Liq'] - comp_1['FeOt_Liq'])
+
+        if type(comp_1) == dict:
+            comp_julia = jl.seval("Dict")(comp_1) 
+        else:
+            comp_new = comp_1.loc[i].to_dict()
+            comp_julia = jl.seval("Dict")(comp_new)
+
+        Output_jl = jl.MAGEMinCalc.AdiabaticDecompressionMelting(comp = comp_julia, P_start_kbar = P_start_bar/1000.0, 
+                                                              P_end_kbar = P_end_bar/1000.0, dp_kbar = dp_bar/1000.0,
+                                                              T_start_C = T_start_C, fo2_buffer = fO2_buffer, fo2_offset = fO2_offset, Model = Model)
+        Results = dict(Output_jl)
+        # Results = stich(Results, Model = Model)
+        
+        q.put([Results, index])
+        return
         
 

PetThermoTools/Plotting.py

@@ -359,7 +359,7 @@ def residualT_plot(Results = None, P_bar = None, phases = None, H2O_Liq = None,
                         a[i][j].set_zlim([0,50])
 
 def phase_plot(Results = None, y_axis = None, x_axis = None, 
-               phases = ['Liq','Ol', 'Opx', 'Cpx', 'Sp', 'Grt'], cmap = "Reds",
+               phases = ['Liq', 'Ol', 'Opx', 'Cpx', 'Sp', 'Grt'], cmap = "Reds",
                title = None, figsize = None):
     
     if type(Results) != list:
@@ -373,10 +373,11 @@ def phase_plot(Results = None, y_axis = None, x_axis = None,
 
         Stop = np.zeros(len(Results['All']['P_bar']))
         for idx, p in enumerate(phases):
-            a.fill_betweenx(Results['All']['P_bar'], Stop, 
-                            x2= Stop + Results['All']['Mass_' + p], alpha = 0.75, color = PhaseList[p], lw = 0)
+            if 'mass_' + p in Results['All'].keys():
+                a.fill_betweenx(Results['All']['P_bar'], Stop, 
+                                x2= Stop + Results['All']['mass_' + p], alpha = 0.75, color = PhaseList[p], lw = 0)
 
-            Stop = Stop + Results['All']['Mass_'+p]
+                Stop = Stop + Results['All']['mass_' + p]
 
         a.set_ylabel('Pressure (bars)')
         a.set_xlabel('Mass (g)')

PetThermoTools/_version.py

@@ -5,4 +5,4 @@
 # 1) we don't load dependencies by storing it in __init__.py
 # 2) we can import it in setup.py for the same reason
 # 3) we can import it into your module
-__version__ = '0.2.38'
+__version__ = '0.2.39'

{PetThermoTools-0.2.38.dist-info → PetThermoTools-0.2.39.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: PetThermoTools
-Version: 0.2.38
+Version: 0.2.39
 Summary: PetThermoTools
 Home-page: https://github.com/gleesonm1/PetThermoTools
 Author: Matthew Gleeson

{PetThermoTools-0.2.38.dist-info → PetThermoTools-0.2.39.dist-info}/RECORD

@@ -5,16 +5,16 @@ PetThermoTools/Holland.py,sha256=udBFeVUyTBpSfLIhx7Hy6o0I8ApNCDvwU_gZa0diY5w,725
 PetThermoTools/Installation.py,sha256=UfVOW1NZFdzMWPyID5u7t0KwvpJA0AqYohzidXIAwYs,6098
 PetThermoTools/Liq.py,sha256=2kdU8r00Y20khwvPvCvLr6P7y2HEOvCH2cwrF_qw5u8,35541
 PetThermoTools/MELTS.py,sha256=DeTM8j2fJ5dAw_2DVCGfgh1R_lQB4Z5rLopKSJ1Q5xQ,73852
-PetThermoTools/Melting.py,sha256=h1-KchPFw5tI2xoOyNJUfVy_TJqTtcYiHtGuInIwzio,11017
+PetThermoTools/Melting.py,sha256=6EXjDJi5ZEZkQdoZbClRvnm3la_9c1tqurwyagcp5ts,12808
 PetThermoTools/Path.py,sha256=ey1nKnmdtriQ_0W6nG5Hbqhhk4EYqY64y1bMwhoXTW0,33050
 PetThermoTools/Path_wrappers.py,sha256=_0pBs_cK2hICxAHkYxKXICUnUEBSiUg07-qhgBeuTdc,26555
 PetThermoTools/PhaseDiagrams.py,sha256=8S_BcqggBzfUbiCPcsJRWFBenGL4tcCevte3-ATQjQI,30884
-PetThermoTools/Plotting.py,sha256=IMYYusZy6REQATG9NI3dGlGs6h3vQ48aJHkZnHmSDDY,28924
+PetThermoTools/Plotting.py,sha256=uvL_j2emMveGumLQ-IeJqyMXGUQT_PyInOpGnsWziAI,28992
 PetThermoTools/Saturation.py,sha256=XXY6fKVouQM3RLgQgXur4xSq7_uGp7bCw_k7NNlWYi8,14095
 PetThermoTools/__init__.py,sha256=PbiwQj_mNNEwuIZOLETmtMMshiXa50wjCA6mfvpOpOs,2393
-PetThermoTools/_version.py,sha256=av7IcIpYHZK1tH6QJlIbcK-mn65OcisD_rhAr-ihfZo,296
-PetThermoTools-0.2.38.dist-info/LICENSE.txt,sha256=-mkx4iEw8Pk1RZUvncBhGLW87Uur5JB7FBQtOmX-VP0,1752
-PetThermoTools-0.2.38.dist-info/METADATA,sha256=X2PlaPAPaopRnSHQS9gfV8kS0nJ4-XI76CldRlFla5A,796
-PetThermoTools-0.2.38.dist-info/WHEEL,sha256=tZoeGjtWxWRfdplE7E3d45VPlLNQnvbKiYnx7gwAy8A,92
-PetThermoTools-0.2.38.dist-info/top_level.txt,sha256=IqK8iYBR3YJozzMOTRZ8x8mU2k6x8ycoMBxZTm-I06U,15
-PetThermoTools-0.2.38.dist-info/RECORD,,
+PetThermoTools/_version.py,sha256=mHkADyHdAemaOsQ81BakYanLTpx-qI1XT0aama1lkCw,296
+PetThermoTools-0.2.39.dist-info/LICENSE.txt,sha256=-mkx4iEw8Pk1RZUvncBhGLW87Uur5JB7FBQtOmX-VP0,1752
+PetThermoTools-0.2.39.dist-info/METADATA,sha256=EaGuwfuigrt061ZQbwUN0JEYA0jg-D7YIQd-pXPwQe0,796
+PetThermoTools-0.2.39.dist-info/WHEEL,sha256=tZoeGjtWxWRfdplE7E3d45VPlLNQnvbKiYnx7gwAy8A,92
+PetThermoTools-0.2.39.dist-info/top_level.txt,sha256=IqK8iYBR3YJozzMOTRZ8x8mU2k6x8ycoMBxZTm-I06U,15
+PetThermoTools-0.2.39.dist-info/RECORD,,