Melting using PTpath and issues with pressure integral

[nbaxter]

I'm having some issues with Adiabat_1ph version 3. I'm trying to melt depleted mantle using a pressure/temperature path text file and I'm getting strange results. Why does the pressure integral produce melt compositions that are the same at every step (the composition is the same as the final integrated melt composition)? Also, what are the numbers highlighted in bold in the aggregate integral file? I also tried running the model without log fO2 set, but it still produced the same errors. I've probably missed something really simple or I've installed the program incorrectly. Hopefully, you can help me. Here are the input files:

Environmental File:
ADIABAT_VERSION pMELTS
ADIABAT_MODE PTpath
ADIABAT_PTPATH_FILE /Users/nichellebaxter/Desktop/Adbat3/NH_ptpath.txt
ADIABAT_DELTAP -323.7
ADIABAT_DELTAT -.439
ADIABAT_MAXP +12301.7
ADIABAT_MINP +4208.49
ADIABAT_MAXT +1300.587
ADIABAT_MINT +1248.81
ADIABAT_CONTINUOUS_MELTING true
ADIABAT_MINF 0.005
ADIABAT_MASSIN  0.001
ADIABAT_INTEGRATE_FILE NH_ptpath_bu_integrate

PT path file:
12301.7   1300.587
11978      1298.746
11654.3   1296.898
11330.6   1295.042
11006.8   1293.179
10683.1   1291.309
10359.4   1289.432
10035.6   1287.55
9711.9      1285.661
9388.17   1283.765
9064.44   1281.864
8740.71   1279.956
8416.98   1278.042
8093.25   1276.121
7769.52   1274.194
7445.79   1272.258
7122.06   1270.313
6798.33   1268.357
6474.6      1266.385
6150.87   1264.389
5827.14   1262.356
5503.41   1260.261
5179.68   1258.055
4855.95   1255.646
4532.22   1252.85
4208.49   1248.81

Source
Title: Depleted MORB Mantle of Workmann & Hart 2005
Initial Composition: SiO2 44.71
Initial Composition: Al2O3 3.98
Initial Composition: Fe2O3 0.191
Initial Composition: FeO 8.008
Initial Composition: MnO 0.13
Initial Composition: MgO 38.73
Initial Composition: CaO 3.17
Initial Composition: Na2O 0.28
Initial Composition: Cr2O3 0.57
Initial Composition: TiO2 0.13
Initial Composition: NiO 0.24
Initial Composition: K2O 0.006
Initial Composition: P2O5 0.019
Initial Temperature: 1300.587
Initial Pressure: 12301.7
Log fO2 Path: -1FMQ

The model ran to completion and then I selected option 16 to create integrated melt compositions. Here are selected output files:

Pressure integral
P      F   SiO2_2D   TiO2_2D
12625.4   1573.74   0   53.3747   0.562774
12301.7   1573.74   0.00344351   53.3747   0.562774
11978   1571.9   0.0044529   53.3747   0.562774
11654.3   1570.05   0.00558695   53.3747   0.562774
11330.6   1568.19   0.00685825   53.3747   0.562774
11006.8   1566.33   0.00827971   53.3747   0.562774
10683.1   1564.46   0.00986086   53.3747   0.562774
10359.4   1562.58   0.0116116   53.3747   0.562774
10035.6   1560.7   0.0135417   53.3747   0.562774
9711.9   1558.81   0.015657   53.3747   0.562774
9388.17   1556.91   0.0179675   53.3747   0.562774
9064.44   1555.01   0.0204849   53.3747   0.562774
8740.71   1553.11   0.0232258   53.3747   0.562774
8416.98   1551.19   0.0262164   53.3747   0.562774
8093.25   1549.27   0.0294931   53.3747   0.562774
7769.52   1547.34   0.0331083   53.3747   0.562774
7445.79   1545.41   0.0371327   53.3747   0.562774
7122.06   1543.46   0.041668   53.3747   0.562774
6798.33   1541.51   0.0468581   53.3747   0.562774
6474.6   1539.53   0.0498903   53.3747   0.562774
6150.87   1537.54   0.0548304   53.3747   0.562774
5827.14   1535.51   0.0637795   53.3747   0.562774
5503.41   1533.41   0.0750565   53.3747   0.562774
5179.68   1531.2   0.0823631   53.3747   0.562774
4855.95   1528.8   0.0963591   53.3747   0.562774
4532.22   1526   0.113153   53.3747   0.562774

Aggregate Integral
P   T   F   SiO2   TiO2
12301.7   1573.74   0.00344351   4.94066e-323   4.94066e-323
12301.7   1573.74   6.89E-03   26.6873   0.281387
11978   1571.9   0.00789641   30.0988   0.317356
11654.3   1570.05   0.00903047   33.0218   0.348176
11330.6   1568.19   0.0103018   35.5334   0.374659
11006.8   1566.33   0.0117232   37.6967   0.397468
10683.1   1564.46   0.0133044   39.56   0.417114
10359.4   1562.58   0.0150551   41.1664   0.434052
10035.6   1560.7   0.0169852   42.5537   0.448679
9711.9   1558.81   0.0191005   43.7521   0.461315
9388.17   1556.91   0.0214111   44.7905   0.472264
9064.44   1555.01   0.0239284   45.6936   0.481786
8740.71   1553.11   0.0266693   46.483   0.490109
8416.98   1551.19   0.02966   47.1779   0.497436
8093.25   1549.27   0.0329366   47.7944   0.503936
7769.52   1547.34   0.0365518   48.3463   0.509755
7445.79   1545.41   0.0405762   48.845   0.515014
7122.06   1543.46   0.0451115   49.3004   0.519815
6798.33   1541.51   0.0503016   49.7208   0.524248
6474.6   1539.53   0.0533338   49.9285   0.526438
6150.87   1537.54   0.0582739   50.2207   0.529518
5827.14   1535.51   0.067223   50.6406   0.533946
5503.41   1533.41   0.0785   51.0333   0.538087
5179.68   1531.2   0.0858066   51.2327   0.540189
4855.95   1528.8   0.0998026   51.5331   0.543356
4532.22   1526   0.116596   51.7983   0.546153

Again, there must be something I'm missing about the melt integration step. Hopefully you can help. I look forward to your reply! Thanks!


-Nichelle

[Paula]

Hi Nichelle,

It sounds like you are doing everything right but have run into a couple of bugs, sorry.  We discovered them last week but I just managed to track down the causes yesterday.  I am going to put up release 3.1 very soon but I'm just fixing up / adding a couple of other things first.  In the meantime, if you let me know your operating system I can send you a link to the corrected executable.

Thanks for letting us know about the problem(s).

Paula

[nbaxter]

Thanks for your quick reply! I have a Mac OS X 10.6.8 .

-Nichelle