% == first section =============

2    % no of numbers in each line of data, in this case 2, P and T
     
6    % system size - in this case KFMASH with 6 components
     
2 1  % which columns to be x,y in phase diagram: so x = T, y = P in this case

% ------------------------------
% == second section ============

% lines (u) and points (i) can be in any order, and the numbers used
% for them, are your choice

u1   g bi cd liq ksp sill q   univariant    % bi + sill = g + cd

i1 i4  red     % the part of the line to plot is between i1 and i4, in red

4.000 735.41   %  21bi + 26sill + 96q = 10g + 12cd + 27liq + 14ksp               
4.200 739.68   %  21bi + 27sill + 97q = 10g + 12cd + 26liq + 14ksp             
4.400 744.03   %  21bi + 28sill + 99q = 10g + 13cd + 26liq + 15ksp             
4.600 748.47   %  22bi + 30sill + 101q = 10g + 14cd + 26liq + 15ksp             
4.800 753.01   %  23bi + 31sill + 103q = 10g + 14cd + 26liq + 16ksp             
5.000 757.64   %  23bi + 33sill + 106q = 10g + 15cd + 26liq + 16ksp             
5.200 762.39   %  24bi + 34sill + 109q = 10g + 16cd + 26liq + 17ksp             
5.400 767.25   %  25bi + 36sill + 112q = 10g + 17cd + 26liq + 18ksp             
5.600 772.23   %  26bi + 38sill + 116q = 10g + 18cd + 26liq + 18ksp             
5.800 777.32   %  26bi + 40sill + 120q = 10g + 20cd + 26liq + 19ksp             
6.000 782.53   %  27bi + 42sill + 124q = 10g + 21cd + 26liq + 20ksp             
6.200 787.85   %  28bi + 44sill + 128q = 10g + 22cd + 26liq + 21ksp             
6.400 793.28   %  29bi + 47sill + 133q = 10g + 23cd + 26liq + 22ksp             
6.600 798.82   %  30bi + 49sill + 138q = 10g + 25cd + 26liq + 23ksp             
6.800 804.46   %  32bi + 52sill + 143q = 10g + 26cd + 26liq + 24ksp             
7.000 810.20   %  33bi + 55sill + 148q = 10g + 28cd + 26liq + 25ksp 

% ------------------------------
i3   g cd liq ksp q  - bi sill    % ie calc with g cd liq ksp bi sill, with
                                  % the modes of bi & sill set to zero
5.301 764.85  % bi = 0, sill = 0

% ------------------------------
i1   bi cd liq ksp q  - g sill
                                  
4.175 739.15  % g = 0, sill = 0

% ------------------------------
i2   bi liq ksp sill q  - g cd

4.752 751.91  % g = 0, cd = 0

% ------------------------------
i4   g liq ksp sill q  - bi cd
 
6.740 802.75  % bi = 0, cd = 0

% ------------------------------
u8  g cd liq ksp q  - bi   % ie calc with g cd liq ksp q bi, with
                           % the mode of bi set to zero
i9  i3                     % plot line between i9 and i3

3.400 790.60  % bi = 0
3.600 788.66  % bi = 0
3.800 786.51  % bi = 0
4.000 784.16  % bi = 0
4.200 781.62  % bi = 0
4.400 778.91  % bi = 0
4.600 776.03  % bi = 0
4.800 773.01  % bi = 0
5.000 769.84  % bi = 0
5.200 766.56  % bi = 0

% ------------------------------
u7   g cd liq ksp q  - sill

i3  i6

5.296 760.00  % sill = 0
5.302 765.00  % sill = 0
5.306 770.00  % sill = 0
5.310 775.00  % sill = 0
5.313 780.00  % sill = 0
5.316 785.00  % sill = 0
5.318 790.00  % sill = 0
5.319 795.00  % sill = 0
5.320 800.00  % sill = 0
5.319 805.00  % sill = 0
5.319 810.00  % sill = 0
5.318 815.00  % sill = 0
5.316 820.00  % sill = 0
5.313 825.00  % sill = 0
5.311 830.00  % sill = 0
5.307 835.00  % sill = 0
5.304 840.00  % sill = 0
5.300 845.00  % sill = 0
5.296 850.00  % sill = 0
5.291 855.00  % sill = 0
5.287 860.00  % sill = 0
5.282 865.00  % sill = 0
5.277 870.00  % sill = 0

% ------------------------------
u14   g cd liq ksp q  - opx

i9  i11

3.606 780.00  % opx = 0
3.608 785.00  % opx = 0
3.610 790.00  % opx = 0
3.612 795.00  % opx = 0
3.614 800.00  % opx = 0
3.615 805.00  % opx = 0
3.617 810.00  % opx = 0
3.618 815.00  % opx = 0
3.620 820.00  % opx = 0
3.621 825.00  % opx = 0
3.623 830.00  % opx = 0
3.625 835.00  % opx = 0

% ------------------------------
u10   g bi opx cd liq ksp q   univariant  % bi + g = opx + cd

i7 i9  red
            
1.000 744.33   %  27g + 250bi + 1833q = 393opx + 10cd + 739liq + 34ksp             
1.200 747.92   %  26g + 215bi + 1533q = 338opx + 10cd + 619liq + 32ksp             
1.400 751.45   %  25g + 188bi + 1312q = 298opx + 10cd + 531liq + 31ksp             
1.600 754.92   %  24g + 167bi + 1140q = 265opx + 10cd + 461liq + 30ksp             
1.800 758.35   %  23g + 149bi + 1001q = 239opx + 10cd + 405liq + 29ksp             
2.000 761.74   %  23g + 135bi + 887q = 217opx + 10cd + 359liq + 28ksp             
2.200 765.11   %  22g + 122bi + 793q = 199opx + 10cd + 320liq + 27ksp             
2.400 768.45   %  22g + 112bi + 713q = 183opx + 10cd + 287liq + 26ksp             
2.600 771.78   %  22g + 103bi + 646q = 169opx + 10cd + 259liq + 25ksp             
2.800 775.10   %  21g + 95bi + 589q = 158opx + 10cd + 236liq + 24ksp             
3.000 778.42   %  21g + 88bi + 540q = 148opx + 10cd + 215liq + 23ksp             
3.200 781.74   %  21g + 82bi + 499q = 139opx + 10cd + 198liq + 22ksp             
3.400 785.06   %  21g + 77bi + 463q = 131opx + 10cd + 183liq + 22ksp             
3.600 788.40   %  21g + 73bi + 432q = 125opx + 10cd + 170liq + 21ksp             
3.800 791.74   %  20g + 69bi + 405q = 119opx + 10cd + 158liq + 20ksp             
4.000 795.09   %  20g + 66bi + 381q = 114opx + 10cd + 148liq + 20ksp             
4.200 798.46   %  20g + 63bi + 361q = 110opx + 10cd + 140liq + 19ksp             
4.400 801.83   %  20g + 60bi + 344q = 106opx + 10cd + 132liq + 19ksp
4.600 805.21   %  20g + 58bi + 329q = 103opx + 10cd + 126liq + 18ksp             
4.800 808.59   %  20g + 56bi + 315q = 100opx + 10cd + 120liq + 18ksp
5.000 811.98   %  20g + 54bi + 304q = 97opx + 10cd + 115liq + 18ksp 

% ------------------------------
i7   bi cd liq ksp q  - g opx

2.503 770.16  % g = 0, opx = 0

% ------------------------------
i8  opx cd liq ksp q  - g bi

2.613 771.99  % g = 0, bi = 0

% ------------------------------
i9   g cd liq ksp q  - bi opx

% 3.610 788.56  % bi = 0, opx = 0   % as calculated (currently commented out)

crossover u8 u14                    % to illustrate the use of this keyword

% ------------------------------
u9   bi cd liq ksp q  - g

i7  i1

1.200 783.91  % g = 0
1.400 782.55  % g = 0
1.600 780.90  % g = 0
1.800 778.97  % g = 0
2.000 776.78  % g = 0
2.200 774.33  % g = 0
2.400 771.64  % g = 0
2.600 768.71  % g = 0
2.800 765.57  % g = 0
3.000 762.23  % g = 0
3.200 758.70  % g = 0
3.400 754.98  % g = 0
3.600 751.11  % g = 0
3.800 747.08  % g = 0
4.000 742.91  % g = 0

% ------------------------------
u11   bi cd liq ksp q  - opx

i35  i7

0.200 734.04  % opx = 0
0.400 737.88  % opx = 0
0.600 741.56  % opx = 0
0.800 745.09  % opx = 0
1.000 748.48  % opx = 0
1.200 751.73  % opx = 0
1.400 754.86  % opx = 0
1.600 757.87  % opx = 0
1.800 760.77  % opx = 0
2.000 763.57  % opx = 0
2.200 766.26  % opx = 0
2.400 768.86  % opx = 0
2.600 771.37  % opx = 0
2.800 773.79  % opx = 0
3.000 776.13  % opx = 0

% ------------------------------
u12   opx cd liq ksp q  - bi

i33  i8

0.400 738.17  % bi = 0
0.600 741.86  % bi = 0
0.800 745.40  % bi = 0
1.000 748.80  % bi = 0
1.200 752.06  % bi = 0
1.400 755.21  % bi = 0
1.600 758.23  % bi = 0
1.800 761.15  % bi = 0
2.000 763.96  % bi = 0
2.200 766.68  % bi = 0
2.400 769.30  % bi = 0
2.600 771.83  % bi = 0
2.800 774.28  % bi = 0
3.000 776.64  % bi = 0

% ------------------------------
u13   opx cd liq ksp q  - g

i8 i10

2.606 770.00  % g = 0
2.624 775.00  % g = 0
2.642 780.00  % g = 0
2.660 785.00  % g = 0
2.679 790.00  % g = 0
2.698 795.00  % g = 0
2.716 800.00  % g = 0
2.736 805.00  % g = 0
2.755 810.00  % g = 0
2.775 815.00  % g = 0
2.795 820.00  % g = 0
2.816 825.00  % g = 0

% remainder of the lines and points:

% ------------------------------
i5   g liq sill q  - cd ksp

6.367 871.08  % cd = 0, ksp = 0

% ------------------------------
i6   g cd liq q  - ksp sill

5.277 870.34  % ksp = 0, sill = 0

% ------------------------------
i10   opx cd liq q  - g ksp

2.820 826.06  % g = 0, ksp = 0

% ------------------------------
i11   g cd liq q  - opx ksp

3.627 839.46  % opx = 0, ksp = 0

% ------------------------------
i12   bi mu liq sill q  - g ksp

6.690 742.08  % g = 0, ksp = 0

% ------------------------------
i13   bi liq ksp sill q  - g mu

6.700 742.28  % g = 0, mu = 0

% ------------------------------
i14   g bi mu liq q  - ksp sill

7.870 764.26  % ksp = 0, sill = 0

% ------------------------------
i15   g mu liq ksp q  - bi sill

9.323 790.96  % bi = 0, sill = 0

% ------------------------------
i16   g liq ksp sill q  - bi mu

9.465 793.54  % bi = 0, mu = 0

% ------------------------------
i17   g mu liq ksp q  - sill ky

9.810 799.22  % sill = 0, ky = 0

% ------------------------------
i18   g liq ksp sill q  - mu ky

9.815 799.48  % mu = 0, ky = 0

% ------------------------------
i19   bi mu sill q H2O  - g liq

6.645 735.44  % g = 0, liq = 0

% ------------------------------
i20   bi mu liq sill q  - g H2O

6.646 735.45  % g = 0, H2O = 0

% ------------------------------
i21   g bi mu q H2O  - liq sill

7.732 748.81  % liq = 0, sill = 0

% ------------------------------
i22   g bi mu liq q  - sill H2O

7.781 749.39  % sill = 0, H2O = 0

% ------------------------------
i23   bi mu liq sill q  - ksp H2O

5.728 723.94  % ksp = 0, H2O = 0

% ------------------------------
i24   bi liq ksp sill q  - mu H2O

5.728 723.94  % mu = 0, H2O = 0

% ------------------------------
i25   bi mu sill q H2O  - liq ksp

5.728 723.94  % liq = 0, ksp = 0

% ------------------------------
i26   bi ksp sill q H2O  - mu liq

5.728 723.94  % mu = 0, liq = 0

% ------------------------------
i27   bi liq ksp sill q  - cd H2O

4.229 726.84  % cd = 0, H2O = 0

% ------------------------------
i28   bi ksp sill q H2O  - cd liq

4.228 726.84  % cd = 0, liq = 0

% ------------------------------
i29   bi cd liq ksp q  - sill H2O

3.930 727.47  % sill = 0, H2O = 0

% ------------------------------
i30   bi cd ksp q H2O  - liq sill

3.922 727.49  % liq = 0, sill = 0

% ------------------------------
i31  bi cd ksp q H2O  - sill and

2.786 668.74  % sill = 0, and = 0

% ------------------------------
i32   bi ksp sill q H2O  - cd and

2.901 658.19  % cd = 0, and = 0

% ------------------------------
i33   opx cd liq ksp q  - bi H2O

1.500 756.73  % bi = 0, H2O = 0

% ------------------------------
i34   opx cd ksp q H2O  - bi liq

1.502 756.63  % bi = 0, liq = 0

% ------------------------------
i35   bi cd liq ksp q  - opx H2O

1.506 756.46  % opx = 0, H2O = 0

% ------------------------------
i36  bi cd ksp q H2O  - opx liq

1.507 756.40  % opx = 0, liq = 0

% ------------------------------
i37   opx cd liq ksp  - q H2O

1.010 775.41  % q = 0, H2O = 0

% ------------------------------
i38   cd liq q  - g opx

4.231 977.55  % g = 0, opx = 0

% ------------------------------
i39   cd liq  - opx q

3.872 969.66  % opx = 0, q = 0

% ------------------------------
i40  opx cd liq  - ksp q

1.052 777.58  % ksp = 0, q = 0

% ------------------------------
i41   opx cd liq q  - ksp H2O

1.005 776.75  % ksp = 0, H2O = 0

% ------------------------------
i42   cd liq q  - opx H2O

0.882 915.57  % opx = 0, H2O = 0

% ------------------------------
i43   g bi mu q H2O  - sill ky

7.493 689.39  % sill = 0, ky = 0

% ------------------------------
i44   g bi mu st sill ky q H2O   invariant

7.470 688.31 

% ------------------------------
i45   g bi mu q H2O  - st ky

7.502 688.16  % st = 0, ky = 0

% ------------------------------
i46   bi st q  - g sill

6.329 664.37  % g = 0, sill = 0

% ------------------------------
i47   bi sill q  - g st

6.155 660.80  % g = 0, st = 0

% ------------------------------
u2   bi cd liq ksp q  - sill

i29  i1

1.000 584.03  % sill = 0
1.200 593.98  % sill = 0
1.400 603.92  % sill = 0
1.600 613.88  % sill = 0
1.800 623.84  % sill = 0
2.000 633.75  % sill = 0
2.200 643.62  % sill = 0
2.400 653.45  % sill = 0
2.600 663.25  % sill = 0
2.800 673.00  % sill = 0
3.000 682.72  % sill = 0
3.200 692.41  % sill = 0
3.400 702.06  % sill = 0
3.600 711.68  % sill = 0
3.800 721.26  % sill = 0
4.000 730.81  % sill = 0

% ------------------------------
u3   bi liq ksp sill q  - cd

i27  i2

1.000 567.87  % cd = 0
1.200 577.84  % cd = 0
1.600 597.70  % cd = 0
1.800 607.61  % cd = 0
2.000 617.52  % cd = 0
2.200 627.45  % cd = 0
2.400 637.41  % cd = 0
2.600 647.33  % cd = 0
2.800 657.21  % cd = 0
3.000 667.06  % cd = 0
3.200 676.87  % cd = 0
3.400 686.65  % cd = 0
3.600 696.40  % cd = 0
3.800 706.12  % cd = 0
4.000 715.80  % cd = 0
4.200 725.46  % cd = 0
4.400 735.07  % cd = 0
4.600 744.65  % cd = 0
4.800 754.20  % cd = 0
5.000 763.70  % cd = 0

% ------------------------------
u4   bi liq ksp sill q  - g

i2  i13

4.000 753.31  % g = 0
4.200 753.07  % g = 0
4.400 752.74  % g = 0
4.600 752.30  % g = 0
4.800 751.77  % g = 0
5.000 751.15  % g = 0
5.200 750.43  % g = 0
5.400 749.62  % g = 0
5.600 748.72  % g = 0
5.800 747.74  % g = 0
6.000 746.67  % g = 0
6.200 745.52  % g = 0
6.400 744.28  % g = 0
6.600 742.97  % g = 0
6.800 741.57  % g = 0
7.000 740.11  % g = 0
7.200 738.67  % g = 0
7.400 737.21  % g = 0
7.600 735.73  % g = 0
7.800 734.20  % g = 0
8.000 732.64  % g = 0
8.200 731.04  % g = 0
8.400 729.40  % g = 0
8.600 727.71  % g = 0
8.800 725.98  % g = 0
9.000 724.21  % g = 0
9.200 722.39  % g = 0
9.400 720.53  % g = 0
9.600 718.63  % g = 0
9.800 716.69  % g = 0
10.000 714.70  % g = 0

% ------------------------------
u5   g liq ksp sill q  - bi

i4  i16

5.000 803.05  % bi = 0
5.200 803.28  % bi = 0
5.400 803.43  % bi = 0
5.600 803.52  % bi = 0
5.800 803.54  % bi = 0
6.000 803.49  % bi = 0
6.200 803.37  % bi = 0
6.400 803.19  % bi = 0
6.600 802.95  % bi = 0
6.800 802.65  % bi = 0
7.000 802.29  % bi = 0
7.200 801.87  % bi = 0
7.400 801.39  % bi = 0
7.600 800.86  % bi = 0
7.800 800.27  % bi = 0
8.000 799.63  % bi = 0
8.200 798.94  % bi = 0
8.400 798.20  % bi = 0
8.600 797.41  % bi = 0
8.800 796.57  % bi = 0
9.000 795.69  % bi = 0
9.200 794.75  % bi = 0
9.400 793.84  % bi = 0
9.600 792.92  % bi = 0
9.800 792.00  % bi = 0
10.000 791.06  % bi = 0

% ------------------------------
u6   g liq ksp sill q  - cd

i4  i5

6.752 800.00  % cd = 0
6.730 805.00  % cd = 0
6.707 810.00  % cd = 0
6.683 815.00  % cd = 0
6.659 820.00  % cd = 0
6.634 825.00  % cd = 0
6.608 830.00  % cd = 0
6.581 835.00  % cd = 0
6.553 840.00  % cd = 0
6.525 845.00  % cd = 0
6.496 850.00  % cd = 0
6.466 855.00  % cd = 0
6.436 860.00  % cd = 0
6.405 865.00  % cd = 0
6.374 870.00  % cd = 0

% ------------------------------
u15   g opx cd liq q  - ksp

i10 i11

3.000 828.99  % ksp = 0
3.200 832.32  % ksp = 0
3.400 835.68  % ksp = 0
3.600 839.02  % ksp = 0

% ------------------------------
u16   g bi mu liq ksp sill q   univariant  

i12 i16  green

            
6.000 728.82   %  10g + 520mu + 1031q = 25bi + 874liq + 287ksp + 456sill             
6.200 732.69   %  10g + 553mu + 1088q = 27bi + 931liq + 304ksp + 484sill             
6.400 736.54   %  10g + 587mu + 1146q = 28bi + 988liq + 322ksp + 512sill             
6.600 740.36   %  10g + 622mu + 1204q = 29bi + 1048liq + 339ksp + 541sill
6.800 744.17   %  10g + 658mu + 1263q = 31bi + 1108liq + 357ksp + 570sill             
7.000 747.96   %  10g + 694mu + 1322q = 32bi + 1169liq + 374ksp + 599sill             
7.200 751.74   %  10g + 731mu + 1380q = 34bi + 1231liq + 391ksp + 628sill             
7.400 755.50   %  10g + 768mu + 1438q = 36bi + 1293liq + 409ksp + 658sill             
7.600 759.24   %  10g + 804mu + 1496q = 37bi + 1356liq + 425ksp + 687sill             
7.800 762.97   %  10g + 841mu + 1553q = 39bi + 1419liq + 442ksp + 716sill             
8.000 766.68   %  10g + 878mu + 1609q = 41bi + 1482liq + 458ksp + 745sill             
8.200 770.38   %  10g + 915mu + 1664q = 43bi + 1544liq + 474ksp + 774sill             
8.400 774.07   %  10g + 951mu + 1718q = 45bi + 1607liq + 489ksp + 802sill             
8.600 777.75   %  10g + 987mu + 1771q = 47bi + 1669liq + 504ksp + 831sill             
8.800 781.42   %  10g + 1023mu + 1823q = 49bi + 1732liq + 518ksp + 858sill             
9.000 785.07   %  10g + 1060mu + 1874q = 51bi + 1794liq + 532ksp + 886sill             
9.200 788.72   %  10g + 1096mu + 1926q = 53bi + 1857liq + 545ksp + 914sill             
9.400 792.35   %  10g + 1132mu + 1977q = 55bi + 1921liq + 559ksp + 941sill
9.600 795.98   %  10g + 1170mu + 2029q = 57bi + 1986liq + 572ksp + 969sill             
9.800 799.60   %  10g + 1207mu + 2081q = 60bi + 2052liq + 585ksp + 998sill             
10.000 803.2   %  10g + 1246mu + 2135q = 62bi + 2121liq + 598ksp + 1027sill
10.200 806.8   %  10g + 1287mu + 2191q = 65bi + 2192liq + 612ksp + 1057sill
10.400 810.3   %  10g + 1329mu + 2250q = 68bi + 2267liq + 626ksp + 1089sill
10.600 813.9   %  10g + 1374mu + 2312q = 71bi + 2347liq + 640ksp + 1122sill
10.800 817.5   %  10g + 1422mu + 2379q = 74bi + 2432liq + 655ksp + 1158sill
11.000 821.0   %  10g + 1473mu + 2452q = 77bi + 2524liq + 672ksp + 1197sill
11.200 824.6   %  10g + 1530mu + 2532q = 81bi + 2624liq + 690ksp + 1239sill
11.400 828.1   %  10g + 1592mu + 2621q = 86bi + 2735liq + 710ksp + 1286sill
11.600 831.7   %  10g + 1662mu + 2721q = 90bi + 2859liq + 732ksp + 1338sill
11.800 835.2   %  10g + 1740mu + 2835q = 96bi + 2999liq + 758ksp + 1397sill
12.000 838.7   %  10g + 1830mu + 2966q = 102bi + 3158liq + 787ksp + 1464sill 

% ------------------------------
u17  bi liq ksp sill q  - mu

i24  i13

3.000 669.86  % mu = 0
3.200 674.00  % mu = 0
3.400 678.10  % mu = 0
3.600 682.17  % mu = 0
3.800 686.22  % mu = 0
4.000 690.23  % mu = 0
4.200 694.22  % mu = 0
4.400 698.19  % mu = 0
4.600 702.13  % mu = 0
4.800 706.04  % mu = 0
5.000 709.94  % mu = 0
5.200 713.81  % mu = 0
5.400 717.67  % mu = 0
5.600 721.50  % mu = 0
5.800 725.32  % mu = 0
6.000 729.12  % mu = 0
6.200 732.90  % mu = 0
6.400 736.66  % mu = 0
6.600 740.40  % mu = 0
6.800 744.13  % mu = 0
7.000 747.85  % mu = 0
7.200 751.54  % mu = 0
7.400 755.22  % mu = 0
7.600 758.89  % mu = 0
7.800 762.54  % mu = 0
8.000 766.17  % mu = 0
8.200 769.79  % mu = 0
8.400 773.39  % mu = 0
8.600 776.99  % mu = 0
8.800 780.56  % mu = 0
9.000 784.13  % mu = 0
9.200 787.67  % mu = 0
9.400 791.20  % mu = 0
9.600 794.71  % mu = 0
9.800 798.21  % mu = 0
10.000 801.69  % mu = 0

% ------------------------------
u18   bi mu liq sill q  - ksp

i23  i12

3.000 669.86  % ksp = 0
3.200 673.99  % ksp = 0
3.400 678.10  % ksp = 0
3.600 682.17  % ksp = 0
3.800 686.22  % ksp = 0
4.000 690.23  % ksp = 0
4.200 694.22  % ksp = 0
4.400 698.18  % ksp = 0
4.600 702.13  % ksp = 0
4.800 706.04  % ksp = 0
5.000 709.94  % ksp = 0
5.200 713.81  % ksp = 0
5.400 717.67  % ksp = 0
5.600 721.50  % ksp = 0
5.800 725.32  % ksp = 0
6.000 729.11  % ksp = 0
6.200 732.89  % ksp = 0
6.400 736.65  % ksp = 0
6.600 740.40  % ksp = 0
6.800 744.13  % ksp = 0
7.000 747.84  % ksp = 0
7.200 751.54  % ksp = 0
7.400 755.22  % ksp = 0
7.600 758.88  % ksp = 0
7.800 762.53  % ksp = 0
8.000 766.16  % ksp = 0
8.200 769.78  % ksp = 0
8.400 773.39  % ksp = 0
8.600 776.98  % ksp = 0
8.800 780.55  % ksp = 0
9.000 784.12  % ksp = 0
9.200 787.66  % ksp = 0
9.400 791.19  % ksp = 0
9.600 794.70  % ksp = 0
9.800 798.20  % ksp = 0
10.000 801.67  % ksp = 0

% ------------------------------
u19   bi mu liq sill q  - g

i20  i12

6.416 700.00  % g = 0
6.449 705.00  % g = 0
6.481 710.00  % g = 0
6.514 715.00  % g = 0
6.546 720.00  % g = 0
6.578 725.00  % g = 0
6.610 730.00  % g = 0
6.643 735.00  % g = 0

% ------------------------------
u20   g bi mu liq q  - sill

i22  i14

7.554 700.00  % sill = 0
7.576 705.00  % sill = 0
7.598 710.00  % sill = 0
7.620 715.00  % sill = 0
7.642 720.00  % sill = 0
7.664 725.00  % sill = 0
7.687 730.00  % sill = 0
7.710 735.00  % sill = 0
7.734 740.00  % sill = 0
7.758 745.00  % sill = 0
7.784 750.00  % sill = 0
7.811 755.00  % sill = 0
7.841 760.00  % sill = 0
7.875 765.00  % sill = 0
7.915 770.00  % sill = 0

% ------------------------------
u21   g bi mu liq q  - ksp

i14  end

7.800 763.94  % ksp = 0
8.000 764.86  % ksp = 0
8.200 765.80  % ksp = 0
8.400 766.75  % ksp = 0
8.600 767.72  % ksp = 0
8.800 768.68  % ksp = 0
9.000 769.66  % ksp = 0
9.200 770.64  % ksp = 0
9.400 771.63  % ksp = 0
9.600 772.62  % ksp = 0
9.800 773.63  % ksp = 0
10.000 774.64  % ksp = 0
10.200 775.67  % ksp = 0
10.400 776.70  % ksp = 0
10.600 777.75  % ksp = 0
10.800 778.80  % ksp = 0
11.000 779.87  % ksp = 0
11.200 780.95  % ksp = 0
11.400 782.05  % ksp = 0
11.600 783.15  % ksp = 0
11.800 784.28  % ksp = 0
12.000 785.41  % ksp = 0
12.200 786.56  % ksp = 0

% ------------------------------
u22   g mu liq ksp q  - bi

i15  end

9.000 789.83  % bi = 0
9.200 790.52  % bi = 0
9.400 791.23  % bi = 0
9.600 791.95  % bi = 0
9.800 792.67  % bi = 0
10.000 793.39  % bi = 0
10.200 794.12  % bi = 0
10.400 794.86  % bi = 0
10.600 795.59  % bi = 0
10.800 796.34  % bi = 0
11.000 797.09  % bi = 0
11.200 797.84  % bi = 0
11.400 798.60  % bi = 0
11.600 799.37  % bi = 0
11.800 800.14  % bi = 0
12.000 800.92  % bi = 0
12.200 801.71  % bi = 0

% ------------------------------
u23   g mu liq ksp q  - sill

i15  i17

9.000 785.43  % sill = 0
9.200 788.86  % sill = 0
9.400 792.27  % sill = 0
9.600 795.67  % sill = 0
9.800 799.06  % sill = 0
10.000 802.43  % sill = 0
10.200 805.80  % sill = 0

% ------------------------------
u24   g liq ksp sill q  - mu

i16 i18

9.000 785.59  % mu = 0
9.200 789.02  % mu = 0
9.400 792.43  % mu = 0
9.600 795.83  % mu = 0
9.800 799.22  % mu = 0
10.000 802.59  % mu = 0
10.200 805.96  % mu = 0
10.400 809.31  % mu = 0

% ------------------------------
u25   g mu liq ksp q  - ky

i17  end

9.800 799.09  % ky = 0
10.000 801.78  % ky = 0
10.200 804.47  % ky = 0
10.400 807.14  % ky = 0
10.600 809.81  % ky = 0
10.800 812.47  % ky = 0
11.000 815.11  % ky = 0
11.200 817.75  % ky = 0
11.400 820.38  % ky = 0
11.600 823.00  % ky = 0
11.800 825.61  % ky = 0
12.000 828.21  % ky = 0
12.200 830.80  % ky = 0

% ------------------------------
u26   g liq ksp ky q  - mu

i18  end

9.800 799.27  % mu = 0
10.000 801.97  % mu = 0
10.200 804.66  % mu = 0
10.400 807.34  % mu = 0
10.600 810.01  % mu = 0
10.800 812.67  % mu = 0
11.000 815.32  % mu = 0
11.200 817.96  % mu = 0
11.400 820.59  % mu = 0
11.600 823.22  % mu = 0
11.800 825.83  % mu = 0
12.000 828.44  % mu = 0
12.200 831.04  % mu = 0

% ------------------------------
u27   g bi mu liq sill q H2O   univariant

i19 i22   blue
           
6.000 727.18   %  20bi + 1839liq + 367sill = 10g + 443mu + 1673q + 629H2O             
6.200 729.78   %  20bi + 1912liq + 375sill = 10g + 455mu + 1714q + 670H2O             
6.400 732.34   %  20bi + 1979liq + 381sill = 10g + 464mu + 1748q + 710H2O             
6.600 734.87   %  20bi + 2037liq + 386sill = 10g + 472mu + 1774q + 747H2O             
6.800 737.37   %  20bi + 2088liq + 389sill = 10g + 477mu + 1792q + 782H2O             
7.000 739.85   %  20bi + 2130liq + 391sill = 10g + 481mu + 1803q + 814H2O             
7.200 742.32   %  20bi + 2163liq + 391sill = 10g + 483mu + 1807q + 842H2O             
7.400 744.77   %  20bi + 2189liq + 390sill = 10g + 483mu + 1804q + 867H2O             
7.600 747.21   %  20bi + 2208liq + 388sill = 10g + 482mu + 1795q + 890H2O             
7.800 749.62   %  20bi + 2220liq + 385sill = 10g + 480mu + 1782q + 909H2O             
8.000 752.03   %  20bi + 2227liq + 381sill = 10g + 476mu + 1765q + 926H2O 

% ------------------------------
u28   bi mu liq sill q  - H2O

i23  i20

4.000 699.90  % H2O = 0
4.200 702.88  % H2O = 0
4.400 705.80  % H2O = 0
4.600 708.66  % H2O = 0
4.800 711.48  % H2O = 0
5.000 714.24  % H2O = 0
5.200 716.96  % H2O = 0
5.400 719.64  % H2O = 0
5.600 722.28  % H2O = 0
5.800 724.87  % H2O = 0
6.000 727.43  % H2O = 0
6.200 729.95  % H2O = 0
6.400 732.44  % H2O = 0
6.600 734.89  % H2O = 0
6.800 737.31  % H2O = 0
7.000 739.71  % H2O = 0

% ------------------------------
u29   bi mu sill q H2O  - liq

i25  i19

4.000 699.90  % liq = 0
4.200 702.88  % liq = 0
4.400 705.80  % liq = 0
4.600 708.66  % liq = 0
4.800 711.48  % liq = 0
5.000 714.24  % liq = 0
5.200 716.96  % liq = 0
5.400 719.64  % liq = 0
5.600 722.28  % liq = 0
5.800 724.87  % liq = 0
6.000 727.43  % liq = 0
6.200 729.95  % liq = 0
6.400 732.44  % liq = 0
6.600 734.89  % liq = 0
6.800 737.31  % liq = 0
7.000 739.71  % liq = 0

% ------------------------------
u30   bi mu sill q H2O  - g

i47  i19

5.941 630.00  % g = 0
5.976 635.00  % g = 0
6.012 640.00  % g = 0
6.047 645.00  % g = 0
6.081 650.00  % g = 0
6.116 655.00  % g = 0
6.150 660.00  % g = 0
6.184 665.00  % g = 0
6.218 670.00  % g = 0
6.251 675.00  % g = 0
6.284 680.00  % g = 0
6.318 685.00  % g = 0
6.350 690.00  % g = 0
6.383 695.00  % g = 0
6.416 700.00  % g = 0
6.448 705.00  % g = 0
6.481 710.00  % g = 0
6.513 715.00  % g = 0
6.545 720.00  % g = 0
6.577 725.00  % g = 0
6.609 730.00  % g = 0
6.642 735.00  % g = 0
6.675 740.00  % g = 0

% ------------------------------
u31   g bi mu q H2O  - sill

i43  i21

7.238 630.00  % sill = 0
7.261 635.00  % sill = 0
7.283 640.00  % sill = 0
7.305 645.00  % sill = 0
7.326 650.00  % sill = 0
7.348 655.00  % sill = 0
7.369 660.00  % sill = 0
7.391 665.00  % sill = 0
7.412 670.00  % sill = 0
7.433 675.00  % sill = 0
7.454 680.00  % sill = 0
7.475 685.00  % sill = 0
7.495 690.00  % sill = 0
7.516 695.00  % sill = 0
7.536 700.00  % sill = 0
7.556 705.00  % sill = 0
7.577 710.00  % sill = 0
7.597 715.00  % sill = 0
7.617 720.00  % sill = 0
7.637 725.00  % sill = 0
7.657 730.00  % sill = 0
7.677 735.00  % sill = 0
7.697 740.00  % sill = 0
7.717 745.00  % sill = 0
7.737 750.00  % sill = 0

% ------------------------------
u32   g bi mu q H2O  - liq

i21  end

7.000 746.47  % liq = 0
7.200 747.11  % liq = 0
7.400 747.74  % liq = 0
7.600 748.39  % liq = 0
7.800 749.02  % liq = 0
8.000 749.65  % liq = 0
8.200 750.27  % liq = 0
8.400 750.88  % liq = 0
8.600 751.49  % liq = 0
8.800 752.09  % liq = 0
9.000 752.68  % liq = 0
9.200 753.27  % liq = 0
9.400 753.86  % liq = 0
9.600 754.44  % liq = 0
9.800 755.03  % liq = 0
10.000 755.61  % liq = 0
10.200 756.19  % liq = 0
10.400 756.77  % liq = 0
10.600 757.35  % liq = 0
10.800 757.93  % liq = 0
11.000 758.51  % liq = 0
11.200 759.10  % liq = 0
11.400 759.69  % liq = 0
11.600 760.29  % liq = 0
11.800 760.89  % liq = 0
12.000 761.49  % liq = 0
12.200 762.10  % liq = 0
12.400 762.72  % liq = 0
12.600 763.34  % liq = 0
12.800 763.97  % liq = 0
13.000 764.61  % liq = 0

% ------------------------------
u33   g bi mu liq q  - H2O

i22  end

7.000 746.90  % H2O = 0
7.200 747.54  % H2O = 0
7.400 748.18  % H2O = 0
7.600 748.82  % H2O = 0
7.800 749.46  % H2O = 0
8.000 750.08  % H2O = 0
8.200 750.70  % H2O = 0
8.400 751.32  % H2O = 0
8.600 751.92  % H2O = 0
8.800 752.52  % H2O = 0
9.000 753.11  % H2O = 0
9.200 753.70  % H2O = 0
9.400 754.29  % H2O = 0
9.600 754.87  % H2O = 0
9.800 755.45  % H2O = 0
10.000 756.03  % H2O = 0
10.200 756.60  % H2O = 0
10.400 757.18  % H2O = 0
10.600 757.76  % H2O = 0
10.800 758.34  % H2O = 0
11.000 758.92  % H2O = 0
11.200 759.51  % H2O = 0
11.400 760.10  % H2O = 0
11.600 760.69  % H2O = 0
11.800 761.29  % H2O = 0
12.000 761.89  % H2O = 0
12.200 762.50  % H2O = 0
12.400 763.12  % H2O = 0
12.600 763.74  % H2O = 0
12.800 764.37  % H2O = 0
13.000 765.00  % H2O = 0

% ------------------------------
u34   bi mu sill q H2O  - ksp

begin  i25

2.000 613.86  % ksp = 0
2.200 620.51  % ksp = 0
2.400 626.99  % ksp = 0
2.600 633.33  % ksp = 0
2.800 639.56  % ksp = 0
3.000 645.69  % ksp = 0
3.200 651.73  % ksp = 0
3.400 657.70  % ksp = 0
3.600 663.61  % ksp = 0
3.800 669.46  % ksp = 0
4.000 675.26  % ksp = 0
4.200 681.01  % ksp = 0
4.400 686.72  % ksp = 0
4.600 692.40  % ksp = 0
4.800 698.05  % ksp = 0
5.000 703.67  % ksp = 0
5.200 709.27  % ksp = 0
5.400 714.85  % ksp = 0
5.600 720.40  % ksp = 0
5.800 725.94  % ksp = 0
6.000 731.46  % ksp = 0

% ------------------------------
u35   bi ksp sill q H2O  - mu

begin  i26

2.000 613.86  % mu = 0
2.200 620.51  % mu = 0
2.400 626.99  % mu = 0
2.600 633.33  % mu = 0
2.800 639.56  % mu = 0
3.000 645.69  % mu = 0
3.200 651.73  % mu = 0
3.400 657.70  % mu = 0
3.600 663.61  % mu = 0
3.800 669.46  % mu = 0
4.000 675.26  % mu = 0
4.200 681.01  % mu = 0
4.400 686.73  % mu = 0
4.600 692.41  % mu = 0
4.800 698.06  % mu = 0
5.000 703.68  % mu = 0
5.200 709.27  % mu = 0
5.400 714.85  % mu = 0
5.600 720.40  % mu = 0
5.800 725.94  % mu = 0
6.000 731.47  % mu = 0

% ------------------------------
u36  bi ksp sill q H2O  - liq

i28  i26

2.000 739.35  % liq = 0
2.200 737.01  % liq = 0
2.400 735.09  % liq = 0
2.600 733.51  % liq = 0
2.800 732.19  % liq = 0
3.000 731.07  % liq = 0
3.200 730.12  % liq = 0
3.400 729.31  % liq = 0
3.600 728.59  % liq = 0
3.800 727.97  % liq = 0
4.000 727.41  % liq = 0
4.200 726.90  % liq = 0
4.400 726.44  % liq = 0
4.600 726.02  % liq = 0
4.800 725.61  % liq = 0
5.000 725.23  % liq = 0
5.200 724.87  % liq = 0
5.400 724.51  % liq = 0
5.600 724.16  % liq = 0
5.800 723.81  % liq = 0
6.000 723.46  % liq = 0

% ------------------------------
u37   bi liq ksp sill q  - H2O

i27  i24

2.000 739.35  % H2O = 0
2.200 737.01  % H2O = 0
2.400 735.09  % H2O = 0
2.600 733.51  % H2O = 0
2.800 732.19  % H2O = 0
3.000 731.08  % H2O = 0
3.200 730.13  % H2O = 0
3.400 729.31  % H2O = 0
3.600 728.60  % H2O = 0
3.800 727.97  % H2O = 0
4.000 727.41  % H2O = 0
4.200 726.90  % H2O = 0
4.400 726.44  % H2O = 0
4.600 726.02  % H2O = 0
4.800 725.62  % H2O = 0
5.000 725.23  % H2O = 0
5.200 724.87  % H2O = 0
5.400 724.51  % H2O = 0
5.600 724.16  % H2O = 0
5.800 723.81  % H2O = 0
6.000 723.46  % H2O = 0

% ------------------------------
u38   bi cd liq ksp sill q H2O   univariant

i30 i27  blue
           
3.900 727.54   %  10cd + 6115liq = 30bi + 1496ksp + 195sill + 5710q + 3028H2O             
3.950 727.42   %  10cd + 5825liq = 28bi + 1423ksp + 185sill + 5406q + 2897H2O             
4.000 727.31   %  10cd + 5576liq = 27bi + 1360ksp + 176sill + 5144q + 2785H2O             
4.050 727.20   %  10cd + 5359liq = 26bi + 1304ksp + 169sill + 4914q + 2688H2O             
4.100 727.09   %  10cd + 5168liq = 25bi + 1256ksp + 163sill + 4710q + 2604H2O             
4.150 726.99   %  10cd + 4998liq = 24bi + 1213ksp + 157sill + 4528q + 2529H2O             
4.200 726.89   %  10cd + 4846liq = 23bi + 1174ksp + 152sill + 4364q + 2462H2O             
4.250 726.80   %  10cd + 4710liq = 22bi + 1139ksp + 147sill + 4215q + 2402H2O             
4.300 726.70   %  10cd + 4586liq = 21bi + 1107ksp + 142sill + 4080q + 2349H2O 

% ------------------------------
u39   bi ksp sill q H2O  - cd

i32  i28

2.000 609.60  % cd = 0
2.200 620.46  % cd = 0
2.400 631.28  % cd = 0
2.600 642.08  % cd = 0
2.800 652.79  % cd = 0
3.000 663.41  % cd = 0
3.200 673.94  % cd = 0
3.400 684.39  % cd = 0
3.600 694.76  % cd = 0
3.800 705.05  % cd = 0
4.000 715.28  % cd = 0
4.200 725.44  % cd = 0

% ------------------------------
u40   bi cd ksp q H2O  - sill

i31  i30

2.000 626.32  % sill = 0
2.200 637.27  % sill = 0
2.400 648.11  % sill = 0
2.600 658.83  % sill = 0
2.800 669.46  % sill = 0
3.000 679.99  % sill = 0
3.200 690.44  % sill = 0
3.400 700.81  % sill = 0
3.600 711.09  % sill = 0
3.800 721.31  % sill = 0
4.000 731.45  % sill = 0
4.200 741.52  % sill = 0

% ------------------------------
u41   bi cd ksp q H2O  - liq

i36  i30

0.400 796.56  % liq = 0
0.600 790.28  % liq = 0
0.800 781.28  % liq = 0
1.000 772.67  % liq = 0
1.200 765.30  % liq = 0
1.400 759.21  % liq = 0
1.600 754.18  % liq = 0
1.800 749.99  % liq = 0
2.000 746.45  % liq = 0
2.200 743.44  % liq = 0
2.400 740.82  % liq = 0
2.600 738.52  % liq = 0
2.800 736.47  % liq = 0
3.000 734.61  % liq = 0
3.200 732.90  % liq = 0
3.400 731.30  % liq = 0
3.600 729.79  % liq = 0
3.800 728.34  % liq = 0
4.000 726.94  % liq = 0
4.200 725.57  % liq = 0

% ------------------------------
u42   bi cd liq ksp q  - H2O

i35  i29

1.000 772.74  % H2O = 0
1.200 765.35  % H2O = 0
1.400 759.24  % H2O = 0
1.600 754.20  % H2O = 0
1.800 750.01  % H2O = 0
2.000 746.48  % H2O = 0
2.200 743.46  % H2O = 0
2.400 740.85  % H2O = 0
2.600 738.55  % H2O = 0
2.800 736.49  % H2O = 0
3.000 734.63  % H2O = 0
3.200 732.93  % H2O = 0
3.400 731.33  % H2O = 0
3.600 729.82  % H2O = 0
3.800 728.38  % H2O = 0
4.000 726.98  % H2O = 0
4.200 725.61  % H2O = 0

% ------------------------------
u43   bi cd ksp q H2O  - and

begin  i31

1.600 612.38  % and = 0
1.800 622.25  % and = 0
2.000 631.94  % and = 0
2.200 641.48  % and = 0
2.400 650.89  % and = 0
2.600 660.18  % and = 0
2.800 669.36  % and = 0
3.000 678.44  % and = 0

% ------------------------------
u44   bi ksp and q H2O  - cd

begin  i32

1.600 596.54  % cd = 0
1.800 606.31  % cd = 0
2.000 615.94  % cd = 0
2.200 625.47  % cd = 0
2.400 634.93  % cd = 0
2.600 644.29  % cd = 0
2.800 653.53  % cd = 0
3.000 662.68  % cd = 0

% ------------------------------
u45   bi cd ksp q H2O  - opx

begin  i36

0.200 614.09  % opx = 0
0.300 643.03  % opx = 0
0.400 663.85  % opx = 0
0.500 679.98  % opx = 0
0.600 693.04  % opx = 0
0.700 703.94  % opx = 0
0.800 713.26  % opx = 0
0.900 721.38  % opx = 0
1.000 728.59  % opx = 0
1.100 735.06  % opx = 0
1.200 740.96  % opx = 0
1.300 746.37  % opx = 0
1.400 751.39  % opx = 0
1.500 756.07  % opx = 0
1.600 760.47  % opx = 0

% ------------------------------
u46   opx cd ksp q H2O  - bi

begin  i34

0.200 614.33  % bi = 0
0.300 643.31  % bi = 0
0.400 664.15  % bi = 0
0.500 680.30  % bi = 0
0.600 693.38  % bi = 0
0.700 704.30  % bi = 0
0.800 713.64  % bi = 0
0.900 721.78  % bi = 0
1.000 728.99  % bi = 0
1.100 735.48  % bi = 0
1.200 741.39  % bi = 0
1.300 746.82  % bi = 0
1.400 751.85  % bi = 0
1.500 756.54  % bi = 0
1.600 760.95  % bi = 0

% ------------------------------
u47   opx cd ksp q H2O  - liq

begin  i34

0.200 855.12  % liq = 0
0.300 838.96  % liq = 0
0.400 825.55  % liq = 0
0.500 814.08  % liq = 0
0.600 804.15  % liq = 0
0.700 795.52  % liq = 0
0.800 788.03  % liq = 0
0.900 781.53  % liq = 0
1.000 775.88  % liq = 0
1.100 770.98  % liq = 0
1.200 766.69  % liq = 0
1.300 762.94  % liq = 0
1.400 759.62  % liq = 0
1.500 756.68  % liq = 0
1.600 754.06  % liq = 0

% ------------------------------
u48   opx cd liq ksp q  - H2O

i37  i33

1.100 771.02  % H2O = 0
1.200 766.74  % H2O = 0
1.300 762.98  % H2O = 0
1.400 759.66  % H2O = 0
1.500 756.72  % H2O = 0
1.600 754.10  % H2O = 0


% ------------------------------
u49   bi opx cd liq ksp q H2O   univariant

i33 i36

1.500 756.72   % 62bi + 1320liq = 106opx + 10cd + 457ksp + 1559q + 514H2O

% ------------------------------
u50   opx cd liq ksp H2O  - q

begin  i37

0.100 219.64  % q = 0
0.400 825.63  % q = 0
0.500 814.16  % q = 0
0.600 804.22  % q = 0
0.700 795.59  % q = 0
0.800 788.09  % q = 0
0.900 781.58  % q = 0
1.000 775.94  % q = 0

% ------------------------------
u51   opx cd liq q  - ksp

i40  i10

1.200 782.09  % ksp = 0
1.400 788.04  % ksp = 0
1.600 793.83  % ksp = 0
1.800 799.46  % ksp = 0
2.000 804.95  % ksp = 0
2.200 810.30  % ksp = 0
2.400 815.52  % ksp = 0
2.600 820.60  % ksp = 0
2.800 825.56  % ksp = 0
3.000 830.40  % ksp = 0
3.200 835.12  % ksp = 0
3.400 839.73  % ksp = 0
3.600 844.23  % ksp = 0
3.800 848.62  % ksp = 0
4.000 852.89  % ksp = 0

% ------------------------------
u52   opx cd liq q  - g

i10  i38

2.779 820.00  % g = 0
2.813 825.00  % g = 0
2.848 830.00  % g = 0
2.883 835.00  % g = 0
2.918 840.00  % g = 0
2.954 845.00  % g = 0
2.990 850.00  % g = 0
3.027 855.00  % g = 0
3.064 860.00  % g = 0
3.102 865.00  % g = 0
3.141 870.00  % g = 0
3.180 875.00  % g = 0
3.220 880.00  % g = 0
3.261 885.00  % g = 0
3.302 890.00  % g = 0
3.344 895.00  % g = 0
3.388 900.00  % g = 0
3.432 905.00  % g = 0
3.477 910.00  % g = 0
3.523 915.00  % g = 0
3.571 920.00  % g = 0
3.619 925.00  % g = 0
3.669 930.00  % g = 0
3.721 935.00  % g = 0
3.774 940.00  % g = 0
3.828 945.00  % g = 0
3.885 950.00  % g = 0
3.943 955.00  % g = 0
4.003 960.00  % g = 0
4.065 965.00  % g = 0
4.129 970.00  % g = 0
4.196 975.00  % g = 0

% ------------------------------
u53   g cd liq q  - opx

i11  i38

3.555 820.00  % opx = 0
3.573 825.00  % opx = 0
3.592 830.00  % opx = 0
3.610 835.00  % opx = 0
3.629 840.00  % opx = 0
3.647 845.00  % opx = 0
3.666 850.00  % opx = 0
3.685 855.00  % opx = 0
3.704 860.00  % opx = 0
3.723 865.00  % opx = 0
3.742 870.00  % opx = 0
3.761 875.00  % opx = 0
3.780 880.00  % opx = 0
3.800 885.00  % opx = 0
3.819 890.00  % opx = 0
3.839 895.00  % opx = 0
3.859 900.00  % opx = 0
3.880 905.00  % opx = 0
3.900 910.00  % opx = 0
3.921 915.00  % opx = 0
3.942 920.00  % opx = 0
3.964 925.00  % opx = 0
3.986 930.00  % opx = 0
4.009 935.00  % opx = 0
4.032 940.00  % opx = 0
4.056 945.00  % opx = 0
4.080 950.00  % opx = 0
4.105 955.00  % opx = 0
4.131 960.00  % opx = 0
4.158 965.00  % opx = 0
4.186 970.00  % opx = 0
4.215 975.00  % opx = 0

% ------------------------------
u54   g cd liq sill q  - ksp

i6 i5

5.400 870.63  % ksp = 0
5.600 871.04  % ksp = 0
5.800 871.34  % ksp = 0
6.000 871.46  % ksp = 0
6.200 871.36  % ksp = 0

% ------------------------------
u55   g liq sill q  - cd

i5  end

6.349 860.00  % cd = 0
6.357 865.00  % cd = 0
6.365 870.00  % cd = 0
6.373 875.00  % cd = 0
6.381 880.00  % cd = 0
6.389 885.00  % cd = 0
6.397 890.00  % cd = 0
6.406 895.00  % cd = 0
6.414 900.00  % cd = 0
6.422 905.00  % cd = 0
6.430 910.00  % cd = 0
6.439 915.00  % cd = 0
6.447 920.00  % cd = 0
6.456 925.00  % cd = 0
6.465 930.00  % cd = 0
6.474 935.00  % cd = 0
6.483 940.00  % cd = 0
6.492 945.00  % cd = 0
6.502 950.00  % cd = 0
6.511 955.00  % cd = 0
6.522 960.00  % cd = 0
6.532 965.00  % cd = 0
6.543 970.00  % cd = 0
6.554 975.00  % cd = 0
6.566 980.00  % cd = 0
6.578 985.00  % cd = 0
6.591 990.00  % cd = 0
6.604 995.00  % cd = 0
6.618 1000.00  % cd = 0
6.633 1005.00  % cd = 0
6.649 1010.00  % cd = 0

% ------------------------------
u56   g cd liq q  - sill

i6  end

5.168 850.00  % sill = 0
5.195 855.00  % sill = 0
5.222 860.00  % sill = 0
5.248 865.00  % sill = 0
5.275 870.00  % sill = 0
5.302 875.00  % sill = 0
5.329 880.00  % sill = 0
5.356 885.00  % sill = 0
5.383 890.00  % sill = 0
5.410 895.00  % sill = 0
5.437 900.00  % sill = 0
5.464 905.00  % sill = 0
5.491 910.00  % sill = 0
5.518 915.00  % sill = 0
5.545 920.00  % sill = 0
5.573 925.00  % sill = 0
5.601 930.00  % sill = 0
5.629 935.00  % sill = 0
5.657 940.00  % sill = 0
5.685 945.00  % sill = 0
5.714 950.00  % sill = 0
5.743 955.00  % sill = 0
5.772 960.00  % sill = 0
5.802 965.00  % sill = 0
5.833 970.00  % sill = 0
5.864 975.00  % sill = 0
5.895 980.00  % sill = 0
5.928 985.00  % sill = 0
5.961 990.00  % sill = 0
5.995 995.00  % sill = 0
6.030 1000.00  % sill = 0
6.066 1005.00  % sill = 0
6.104 1010.00  % sill = 0

% ------------------------------
u57   g liq sill q  - ksp

i5  end

6.000 865.03  % ksp = 0
6.200 868.35  % ksp = 0
6.400 871.61  % ksp = 0
6.600 874.81  % ksp = 0
6.800 877.96  % ksp = 0
7.000 881.06  % ksp = 0
7.200 884.10  % ksp = 0
7.400 887.09  % ksp = 0
7.600 890.03  % ksp = 0
7.800 892.92  % ksp = 0
8.000 895.77  % ksp = 0
8.200 898.57  % ksp = 0
8.400 901.32  % ksp = 0
8.600 904.03  % ksp = 0
8.800 906.69  % ksp = 0
9.000 909.31  % ksp = 0
9.200 911.89  % ksp = 0
9.400 914.43  % ksp = 0
9.600 916.92  % ksp = 0
9.800 919.38  % ksp = 0
10.000 921.80  % ksp = 0
10.200 924.18  % ksp = 0
10.400 926.53  % ksp = 0
10.600 928.83  % ksp = 0
10.800 931.10  % ksp = 0
11.000 933.34  % ksp = 0
11.200 935.54  % ksp = 0
11.400 937.71  % ksp = 0
11.600 939.84  % ksp = 0
11.800 941.94  % ksp = 0
12.000 944.00  % ksp = 0

% ------------------------------
u58   g cd liq q  - ksp

i11 i6

3.000 825.47  % ksp = 0
3.200 830.08  % ksp = 0
3.400 834.55  % ksp = 0
3.600 838.89  % ksp = 0
3.800 843.11  % ksp = 0
4.000 847.20  % ksp = 0
4.200 851.16  % ksp = 0
4.400 855.00  % ksp = 0
4.600 858.71  % ksp = 0
4.800 862.30  % ksp = 0
5.000 865.76  % ksp = 0
5.200 869.09  % ksp = 0
5.400 872.30  % ksp = 0
5.600 875.38  % ksp = 0
5.800 878.34  % ksp = 0
6.000 881.17  % ksp = 0

% ------------------------------
u59   cd liq  - g

i38  end

3.858 960.00  % g = 0
3.964 965.00  % g = 0
4.070 970.00  % g = 0
4.177 975.00  % g = 0
4.283 980.00  % g = 0
4.389 985.00  % g = 0
4.495 990.00  % g = 0
4.601 995.00  % g = 0
4.708 1000.00  % g = 0
4.814 1005.00  % g = 0
4.920 1010.00  % g = 0

% ------------------------------
u60   cd liq  - opx

i39  i38

4.000 972.39  % opx = 0
4.200 976.84  % opx = 0
4.400 981.52  % opx = 0
4.600 986.43  % opx = 0
4.800 991.63  % opx = 0


% ------------------------------
u61   opx cd liq  - q

i40  i39

0.200 681.94  % q = 0
0.400 708.79  % q = 0
0.600 732.22  % q = 0
0.800 753.33  % q = 0
1.000 772.75  % q = 0
1.200 790.86  % q = 0
1.400 807.91  % q = 0
1.600 824.08  % q = 0
1.800 839.50  % q = 0
2.000 854.26  % q = 0
2.200 868.44  % q = 0
2.400 882.10  % q = 0
2.600 895.28  % q = 0
2.800 908.01  % q = 0
3.000 920.33  % q = 0
3.200 932.25  % q = 0
3.400 943.79  % q = 0
3.600 954.98  % q = 0
3.800 965.82  % q = 0

% ------------------------------
u62   cd liq  - q

i39  end

3.684 960.00  % q = 0
3.781 965.00  % q = 0
3.879 970.00  % q = 0
3.978 975.00  % q = 0
4.079 980.00  % q = 0
4.180 985.00  % q = 0
4.283 990.00  % q = 0
4.387 995.00  % q = 0
4.493 1000.00  % q = 0
4.599 1005.00  % q = 0
4.707 1010.00  % q = 0

% ------------------------------
u63   opx cd liq ksp q  - H2O

i37 i41

1.008 776.00  % H2O = 0

% ------------------------------
u64   opx cd liq q H2O  - ksp

begin  i41

0.100 871.53  % ksp = 0
0.200 852.81  % ksp = 0
0.300 838.23  % ksp = 0
0.400 825.92  % ksp = 0
0.500 815.20  % ksp = 0
0.600 805.72  % ksp = 0
0.700 797.29  % ksp = 0
0.800 789.76  % ksp = 0
0.900 783.04  % ksp = 0
1.000 777.00  % ksp = 0

% ------------------------------
u65   opx cd liq q  - H2O

i42 i41

0.989 750.00  % H2O = 0
0.993 755.00  % H2O = 0
0.996 760.00  % H2O = 0
0.999 765.00  % H2O = 0
1.002 770.00  % H2O = 0
1.004 775.00  % H2O = 0
1.006 780.00  % H2O = 0
1.007 785.00  % H2O = 0
1.008 790.00  % H2O = 0
1.009 795.00  % H2O = 0
1.009 800.00  % H2O = 0
1.009 805.00  % H2O = 0
1.009 810.00  % H2O = 0
1.008 815.00  % H2O = 0
1.006 820.00  % H2O = 0
1.005 825.00  % H2O = 0
1.002 830.00  % H2O = 0
1.000 835.00  % H2O = 0
0.996 840.00  % H2O = 0
0.993 845.00  % H2O = 0
0.988 850.00  % H2O = 0
0.983 855.00  % H2O = 0
0.978 860.00  % H2O = 0
0.972 865.00  % H2O = 0
0.966 870.00  % H2O = 0
0.959 875.00  % H2O = 0
0.951 880.00  % H2O = 0
0.943 885.00  % H2O = 0
0.935 890.00  % H2O = 0
0.925 895.00  % H2O = 0
0.916 900.00  % H2O = 0
0.905 905.00  % H2O = 0
0.895 910.00  % H2O = 0
0.884 915.00  % H2O = 0

% ------------------------------
u66   cd liq H2O  - opx

i42  end

0.780 920.00  % opx = 0
0.676 925.00  % opx = 0
0.582 930.00  % opx = 0
0.497 935.00  % opx = 0
0.420 940.00  % opx = 0
0.350 945.00  % opx = 0
0.287 950.00  % opx = 0
0.230 955.00  % opx = 0
0.179 960.00  % opx = 0
0.135 965.00  % opx = 0

% ------------------------------
u67   cd liq  - H2O

i42  end

0.882 915.00  % H2O = 0
0.881 920.00  % H2O = 0
0.880 925.00  % H2O = 0
0.879 930.00  % H2O = 0
0.877 935.00  % H2O = 0
0.875 940.00  % H2O = 0
0.872 945.00  % H2O = 0
0.870 950.00  % H2O = 0
0.867 955.00  % H2O = 0
0.864 960.00  % H2O = 0
0.860 965.00  % H2O = 0
0.856 970.00  % H2O = 0
0.852 975.00  % H2O = 0
0.848 980.00  % H2O = 0
0.843 985.00  % H2O = 0
0.839 990.00  % H2O = 0
0.833 995.00  % H2O = 0
0.828 1000.00  % H2O = 0
0.822 1005.00  % H2O = 0
0.816 1010.00  % H2O = 0

% ------------------------------
u68   cd liq  - opx

i42 i39

0.200 905.21  % opx = 0
0.400 908.20  % opx = 0
0.600 911.23  % opx = 0
0.800 914.30  % opx = 0
1.000 917.40  % opx = 0
1.200 920.56  % opx = 0
1.400 923.76  % opx = 0
1.600 927.01  % opx = 0
1.800 930.33  % opx = 0
2.000 933.70  % opx = 0
2.200 937.15  % opx = 0
2.400 940.67  % opx = 0
2.600 944.27  % opx = 0
2.800 947.95  % opx = 0
3.000 951.74  % opx = 0
3.200 955.63  % opx = 0
3.400 959.64  % opx = 0
3.600 963.78  % opx = 0
3.800 968.07  % opx = 0
4.000 972.51  % opx = 0

% ------------------------------
u69   opx cd liq ksp  - q

i37 i40

1.000 774.88  % q = 0
1.010 775.40  % q = 0
1.020 775.92  % q = 0
1.030 776.44  % q = 0
1.040 776.96  % q = 0
1.050 777.47  % q = 0

% ------------------------------
u70

i41 i40 connect

% ------------------------------
u71   g bi mu q H2O  - ky

i45  i43

8.004 630.00  % ky = 0
7.956 635.00  % ky = 0
7.909 640.00  % ky = 0
7.863 645.00  % ky = 0
7.818 650.00  % ky = 0
7.773 655.00  % ky = 0
7.730 660.00  % ky = 0
7.687 665.00  % ky = 0
7.646 670.00  % ky = 0
7.605 675.00  % ky = 0
7.565 680.00  % ky = 0
7.526 685.00  % ky = 0
7.488 690.00  % ky = 0

% ------------------------------
u72   g bi mu st sill q H2O   univariant

i47  i44  magenta
            
6.000 657.63   %  10mu + 42st + 161q = 48g + 10bi + 339sill + 85H2O             
6.200 661.72   %  10mu + 40st + 150q = 44g + 10bi + 318sill + 79H2O             
6.400 665.83   %  10mu + 37st + 139q = 41g + 10bi + 299sill + 74H2O             
6.600 669.96   %  10mu + 35st + 130q = 38g + 10bi + 282sill + 70H2O             
6.800 674.13   %  10mu + 33st + 121q = 35g + 10bi + 266sill + 66H2O             
7.000 678.32   %  10mu + 31st + 113q = 33g + 10bi + 251sill + 62H2O             
7.200 682.54   %  10mu + 29st + 106q = 30g + 10bi + 238sill + 58H2O
7.400 686.81   %  10mu + 28st + 99q = 28g + 10bi + 225sill + 55H2O             
7.600 691.11   %  10mu + 26st + 92q = 26g + 10bi + 213sill + 52H2O             
7.800 695.45   %  10mu + 25st + 86q = 24g + 10bi + 202sill + 49H2O             
8.000 699.84   %  10mu + 23st + 81q = 22g + 10bi + 192sill + 47H2O 

% ------------------------------
u73   g bi mu st ky q H2O   univariant

i44 i45   magenta
           
7.450 688.39   %  10mu + 27st + 97q = 28g + 10bi + 221ky + 54H2O             
7.460 688.35   %  10mu + 27st + 97q = 27g + 10bi + 221ky + 54H2O             
7.470 688.31   %  10mu + 27st + 96q = 27g + 10bi + 221ky + 54H2O             
7.480 688.26   %  10mu + 27st + 96q = 27g + 10bi + 220ky + 54H2O             
7.490 688.22   %  10mu + 27st + 96q = 27g + 10bi + 220ky + 54H2O             
7.500 688.17   %  10mu + 27st + 96q = 27g + 10bi + 220ky + 54H2O             
7.510 688.13   %  10mu + 27st + 96q = 27g + 10bi + 220ky + 54H2O             
7.520 688.08   %  10mu + 27st + 96q = 27g + 10bi + 219ky + 53H2O             
7.530 688.04   %  10mu + 27st + 95q = 27g + 10bi + 219ky + 53H2O             
7.540 687.99   %  10mu + 27st + 95q = 27g + 10bi + 219ky + 53H2O             
7.550 687.95   %  10mu + 27st + 95q = 27g + 10bi + 218ky + 53H2O 

% ------------------------------
u74   g bi  - st      % mu q H2O

begin  i45

9.617 600.00  % st = 0
9.484 605.00  % st = 0
9.351 610.00  % st = 0
9.221 615.00  % st = 0
9.092 620.00  % st = 0
8.965 625.00  % st = 0
8.839 630.00  % st = 0
8.715 635.00  % st = 0
8.593 640.00  % st = 0
8.472 645.00  % st = 0
8.353 650.00  % st = 0
8.236 655.00  % st = 0
8.121 660.00  % st = 0
8.007 665.00  % st = 0
7.895 670.00  % st = 0
7.784 675.00  % st = 0
7.676 680.00  % st = 0
7.569 685.00  % st = 0
7.464 690.00  % st = 0
7.360 695.00  % st = 0
7.259 700.00  % st = 0

% ------------------------------
u75   bi st  - g     % mu q H2O

begin  i46

7.794 580.00  % g = 0
7.700 585.00  % g = 0
7.607 590.00  % g = 0
7.514 595.00  % g = 0
7.423 600.00  % g = 0
7.332 605.00  % g = 0
7.242 610.00  % g = 0
7.153 615.00  % g = 0
7.065 620.00  % g = 0
6.979 625.00  % g = 0
6.893 630.00  % g = 0
6.808 635.00  % g = 0
6.724 640.00  % g = 0
6.641 645.00  % g = 0
6.559 650.00  % g = 0
6.478 655.00  % g = 0
6.398 660.00  % g = 0
6.320 665.00  % g = 0
6.242 670.00  % g = 0
6.166 675.00  % g = 0
6.091 680.00  % g = 0
6.017 685.00  % g = 0
5.944 690.00  % g = 0
5.873 695.00  % g = 0
5.804 700.00  % g = 0

% ------------------------------
u76   bi st  - sill   % mu q H2O

begin  i46

2.000 526.68  % sill = 0
2.200 533.23  % sill = 0
2.400 539.74  % sill = 0
2.600 546.20  % sill = 0
2.800 552.64  % sill = 0
3.000 559.05  % sill = 0
3.200 565.44  % sill = 0
3.400 571.81  % sill = 0
3.600 578.16  % sill = 0
3.800 584.50  % sill = 0
4.000 590.84  % sill = 0
4.200 597.16  % sill = 0
4.400 603.48  % sill = 0
4.600 609.79  % sill = 0
4.800 616.10  % sill = 0
5.000 622.41  % sill = 0
5.200 628.72  % sill = 0
5.400 635.03  % sill = 0
5.600 641.34  % sill = 0
5.800 647.65  % sill = 0
6.000 653.97  % sill = 0


% ------------------------------
u77   bi sill  - st  % mu q H2O

begin  i47

2.000 528.91  % st = 0
2.200 535.45  % st = 0
2.400 541.95  % st = 0
2.600 548.41  % st = 0
2.800 554.83  % st = 0
3.000 561.23  % st = 0
3.200 567.60  % st = 0
3.400 573.96  % st = 0
3.600 580.30  % st = 0
3.800 586.63  % st = 0
4.000 592.95  % st = 0
4.200 599.26  % st = 0
4.400 605.56  % st = 0
4.600 611.86  % st = 0
4.800 618.15  % st = 0
5.000 624.45  % st = 0
5.200 630.74  % st = 0
5.400 637.03  % st = 0
5.600 643.32  % st = 0
5.800 649.61  % st = 0
6.000 655.91  % st = 0

% aluminosilicate invariant and univariants - copied in from mode 3 output

% ------------------------------
i678  and ky sill

4.30 536.00

% ------------------------------
u888  and sill

i31 i678 

0.2	926
0.4	905
0.6	884
0.8	863
1.0	843
1.2	823
1.4	802
1.6	782
1.8	763
2.0	743
2.2	724
2.4	705
2.6	686
2.8	667
3.0	649
3.2	631
3.4	613
3.6	596
3.8	579
4.0	562

% ------------------------------
u666  ky and

begin i678

0.02	193.5
0.42	225.2
0.82	256.8
1.22	288.6
1.62	320.4
2.02	352.3
2.42	384.2
2.82	416.3
3.22	448.5
3.62	480.8
4.02	513.2

% ------------------------------
u777  sill ky

i17  end  % i678 end

4.0     522
4.4		541
4.8		560
5.2		579
5.6		599
6.0		618
6.4		637
6.8		656
7.2		675
7.6		695
8.0		714
8.4		733
8.8		752
9.2		770
9.6		789
10.0	808
10.4    827
10.8	846
11.2	864
11.6	883
12.0	901
12.4	920
12.8	938
13.2	956
13.6	975
14.0	993
14.4	1011
14.8	1029
15.2	1047
15.6	1065
16.0	1083

*  % ---------------------------- 
% == third section ==============

% areas (or fields) - manual assembly
%  the first number in each line is how light is its colour 0 <-> 1, 1=white
%  then the lines defining the field are given in sequence around the field

% divariants are white, so do not need colouring

% trivariants 

0.8 u7 u8 u14 u58   % this area is focussed on in the text
0.8 u32 u31 u71 u74
0.8 u77 u30 u29 u34
0.8 u35 u36 u39 u44
0.8 u43 u40 u41 u45
0.8 u46 u47
0.8 u28 u19 u18
0.8 u37 u17 u4 u3
0.8 u42 u2 u9 u11
0.8 u33 u20 u21
0.8 u22 u23 u25
0.8 u26 u24 u5 u6 u57
0.8 u13 u12 u48 u69 u51 u13
0.8 u58 u14 u8 u7
0.8 u26 u24 u5 u6 u57
0.8 u55 u54 u56
0.8 u53 u15 u52
0.8 u50 u63 u64
0.8 u76 u75

% quadrivariants 

0.7 u57 u55
0.7 u56 u58 u53 u59
0.7 u69 u63 u70
0.7 u60 u52 u51 u61
0.7 u66 u65 u64

% quinivariants 

0.6 u59 u60 u62
0.6 u67 u66
0.6 u61 u70 u65 u68 

% 6-variant 

0.5 u62 u68 u67

*  % ---------------------------- 
% == fourth section =============

info part of Fig 4a for KFMASH (White et al., 2014)

window  735 820 2.25 7   % plotting window, in this case Tlow, Thigh, Plow, Phigh        

bigticks  50 750 1 3     % main ticks and labelling, in this case 
                         %      every 50 starting at 750C, and
                         %      every 1 kbar starting at 3 kbar
                        
smallticks 5 0.25        % minor ticks, in this case T,P, starting at Tlow and Plow

darkcolour 56 16 101     % red green blue in 0<->255, as in Illustrator RGB

numbering 14             % numbering the lines on the diagram
docircles 7              % add numbers and circles for points on the diagram

*

% ================================
% ==== storage ===================

% swap in these blocks of scripts into the fourth section, replacing the 
% ones there, to see other views of the data

window  650 900 2 10       % plotting window, in this case Tlow, Thigh, Plow, Phigh        

bigticks  50 650 1 2       % main ticks and labelling, in this case 
                        
smallticks 10 0.25         % minor ticks, in this case T,P, starting at Tlow and Plow

darkcolour 56 16 101       % red green blue in 0<->255, as in Illustrator RGB

numbering no               % do not number the lines on the diagram
docircles 4                % add small circles for points on the diagram

% ================================

window  685 693 7.4 7.6       % plotting window, in this case Tlow, Thigh, Plow, Phigh        

bigticks  2 685 0.05 7.4      % main ticks and labelling 
                        
smallticks no 10 0.25         % minor ticks, in this case T,P, starting at Tlow and Plow

darkcolour 56 16 101          % red green blue in 0<->255, as in Illustrator RGB

numbering 14                  % do not number the lines on the diagram
docircles 7                   % add small circles for points on the diagram

% ================================

window  740 744 6.6 6.75       % plotting window, in this case Tlow, Thigh, Plow, Phigh        

bigticks  0.5 740 0.05 6.6      % main ticks and labelling 
                        
smallticks no 10 0.25         % minor ticks, in this case T,P, starting at Tlow and Plow

darkcolour 56 16 101          % red green blue in 0<->255, as in Illustrator RGB

numbering 14                  % do not number the lines on the diagram
docircles 7                   % add small circles for points on the diagram

