-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathJUN92d.bs
1657 lines (1651 loc) · 91.1 KB
/
JUN92d.bs
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
25
O AL BA C CA CL CO
CU F FE H K MN MG
NA NI P S SI SR TI
ZN ZR B E
15.99940 26.98154 137.32700 12.01100 40.07800 35.45270 58.93320
63.54600 18.99840 55.84700 1.00794 39.09830 54.93085 24.30500
22.98977 58.69000 30.97362 32.06600 28.08550 87.62000 47.88000
65.39000 91.22400 10.81000 1.00000
0.0 1.5 1.0 2.0 1.0 0.0 1.0
1.0 0.0 1.0 0.5 0.5 1.0 1.0
0.5 1.0 1.0 0.0 2.0 1.0 2.0
1.0 2.0 0.0 0.0
****MINERAL DATA***** from RG Berman (J.Petrology, v.29, 1988) with
! changes by Capi: comment out: ALBITE (HIGH)
! ALBITE (LOW)
! K-FELDSPAR (HIGH)
! K-FELDSPAR (LOW)
! rename: K-FELDSPAR -> POTASSIUM FELDSPAR
! re-define: ORTHOENSTATITE -> 2 X
! FERROSILITE -> 2 X
!
! preliminary/provisional properties for andradite, baddelyite, Fe cordierite,
! geikelite, hedenbergite, hercynite, kalsilite, leucite, nepheline,
! staurolite, zircon,
! Feb,89: modified Almandine, Ilmenite, Phlogopite; added annite consistent with
! proposed Fe-Mg garnet mixing properties (Berman, American Mineralogist,
! v. 75, p. 328-344, 1990)
! Nov,89: modified staurolite, Fe-cordierite, hercynite for consistency with
! almandine
! Feb,90: corrected Hydrogen gas Cp coefficient
! Aug,90: changed annite and biotite solution properties (McMullin et al., 1991
! Canadian Mineralogist (1991, v. 29, 833-855)
! Aug, 91: added pumpellyite, glaucophane from Evans (1990)
! orthoenstatite formula doubled for ease of solid solution modelling
! calcite, magnesite modified by Mader and Berman (1991,Amer Mineral)
! abbreviation of CORDIERITE changed to 'Cd' in order to match entry into
! hydration model routine
! preliminary properties for ferrotremolite, tschermakite, ferrotschermakite,
! pargasite, ferropargasite: fit 210, IGCP-304 conference, May 92, Urs Mader;
! ref: Mader & Berman, 1992, Current Research, Part E, Geological Survey of Canada,
! Paper 92-1E; (revised properties)
!
!========================================================================
! main block of MINERAL DATA
!========================================================================
!
****MINERAL DATA***** from RG Berman (J.Petrology, v.29, 1988) with
AKERMANITE MG(1)CA(2)SI(2)O(7) Ak AKER
ST -3667112.00 -3860441.19 212.0000 9.252 1.00000
C1 387.06396 -2938.769 0.000 -40790480. 0.00000
V1 2.50112887 0.00067224 -0.07851406 0.00000000 0.00000000
T1 358.00 298.15 0.00000000 0.00000000 452.00
ANTIGORITE MG(48)SI(34)O(99)H(62)O(48) Atg ANTI
ST -66151354.00 -71364156.29 3602.9958 174.246 1.00000
C1 7394.51172 0.000 -548363008.000 87284121600. 0.00000
V1 2.49647051 0.00039444 -0.19777757 0.00049442 0.00000000
BADDELYITE ZR(1)O(2) Bd BADD
! All properties are provisional estimates not included in JPet 1988 paper
ST 0.00 -1100560.00 50.3800 2.115 3.00000
C1 89.49913 -374.135 -1602304.000 161287724. 0.00000
T1 1478.00 298.15 0.00000000 0.00000000 330.00
BRUCITE MG(1)O(2)H(2) Br BRUC
ST -834549.00 -925936.60 63.0642 2.468 1.00000
C1 136.83899 -537.130 -4361949.000 552694528. 0.00000
V1 3.28538493 0.00010940 -0.20228525 0.00067261 0.00000000
CHRYSOTILE MG(3)SI(2)O(9)H(4) Chr CRYS
ST -4039429.00 -4363356.18 220.1345 10.720 1.00000
C1 610.02197 -5581.152 -1857300.000 195472000. 0.00000
V1 2.71512687 0.00067351 -0.18102799 0.00045243 0.00000000
DIASPORE AL(1)O(2)H(1) Dsp DIAS
ST -920021.00 -999377.89 35.3082 1.776 1.00000
C1 143.23953 -1540.409 -323098.000 64632256. 0.00000
V1 2.97180367 0.00000000 -0.05992916 0.00000000 0.00000000
GEHLENITE AL(2)CA(2)SI(1)O(7) Ge GEHL
ST -3775355.00 -3988157.80 198.6000 9.033 1.00000
C1 373.08740 -2276.755 -4778466.000 477911808. 0.00000
V1 2.49260884 0.00005664 -0.09964904 0.00024878 0.00000000
D1 -221.74000 0.00000 17291280.00000 0.00000 0.00000
D2 0.36950 -0.00015 698.15000 1600.00000 0.00000
HEMATITE FE(2)O(3) Hm HEMA
ST -743715.00 -825627.39 87.4371 3.027 1.00000
C1 146.85840 0.000 -5576782.000 525627904. 0.00000
V1 3.83095170 0.00001650 -0.04789400 0.00003040 0.00000000
T1 955.65 298.15 -0.07403236 0.00027921 1287.92
KALSILITE K(1)AL(1)SI(1)O(4) Kls KALI
! All properties are provisional estimates not included in JPet 1988 paper
ST 0.00 -2111813.55 133.9653 5.989 3.00000
C1 186.00261 0.000 -13106745.000 2138928128. 0.00000
V1 2.34867257 0.00140257 -0.15670396 0.00000000 0.00000000
T1 800.15 298.15 -0.07096454 0.00021682 1154.00
KAOLINITE AL(2)SI(2)O(9)H(4) Kln KAOL
ST -3797448.00 -4120327.22 203.7000 9.952 2.00000
C1 523.23291 -4426.668 -2244259.000 92308768. 0.00000
V1 3.20000000 0.00000000 -0.12000000 0.00000000 0.00000000
LAWSONITE CA(1)AL(2)SI(2)O(10)H(4) Lw LAWS
ST -4508229.00 -4865665.95 229.1759 10.144 1.00000
C1 728.67310 -8248.055 0.000 850556416. 0.00000
V1 2.62833695 0.00000000 -0.07688091 0.00019223 0.00000000
ALEUCITE K(1)AL(1)SI(2)O(6) aLc ALUC
! All properties are provisional estimates not included in JPet 1988 paper
ST 0.00 -3021141.60 199.2000 8.839 3.00000
C1 271.13965 -944.133 -7857238.000 959202560. 0.00000
V1 1.87804050 0.01047630 -0.17999774 0.00045028 0.00000000
T1 955.15 298.15 -0.09730853 0.00033730 0.00
BLEUCITE K(1)AL(1)SI(2)O(6) bLc BLUC
! All properties are provisional estimates not included in JPet 1988 paper
ST 0.00 -3012026.03 210.7042 9.263 3.00000
C1 271.13965 -944.133 -7857238.000 959202560. 0.00000
V1 1.25167877 0.00000000 -0.15568254 0.00012739 0.00000000
LIME CA(1)O(1) Lm LIME
ST -603418.00 -635090.00 37.7500 1.676 2.00000
C1 58.79117 -133.904 -1147146.000 102978784. 0.00000
V1 3.46103555 0.00067406 -0.10217132 0.00025650 0.00000000
MAGNETITE FE(3)O(4) Mt MAGN
ST -1014343.00 -1117402.53 146.1142 4.452 1.00000
C1 207.93083 0.000 -7243337.000 664356352. 0.00000
V1 3.02910360 0.00138470 -0.05817030 0.00017510 0.00000000
T1 848.15 298.15 -0.19502452 0.00061037 1565.57
MEIONITE CA(4)AL(6)SI(6)O(27)C(1) Me MEIO
ST -13102086.00 -13849722.94 730.0000 34.036 2.00000
C1 1511.34497 -13243.328 0.000 -751613952. 0.00000
V1 0.93400000 0.00000000 -0.11102715 0.00000000 0.00000000
MERWINITE CA(3)MG(1)SI(2)O(8) Mw MERW
ST -4311116.00 -4537497.44 251.7767 9.847 1.00000
C1 453.61963 -3250.015 0.000 -344227072. 0.00000
V1 2.93755865 0.00087235 -0.05505433 0.00013811 0.00000000
MONTICELLITE CA(1)MG(1)SI(1)O(4) Mo MONT
ST -2135224.00 -2250027.48 108.3000 5.148 1.00000
C1 226.34225 -1542.741 -1179709.000 -23285888. 0.00000
V1 2.78631922 0.00076339 -0.09036547 0.00020000 0.00000000
NEPHELINE NA(1)AL(1)SI(1)O(4) Ne NEPH
! All properties are provisional estimates not included in JPet 1988 paper
ST 0.00 -2087975.81 124.2000 5.422 3.00000
C1 205.24190 -759.898 -10838267.000 2081817344. 0.00000
V1 3.15522000 0.00200000 -0.22205720 0.00060000 0.00000000
T1 467.15 298.15 -1.02784061 0.00339448 241.84
T3 1180.15 298.15 0.00000000 0.00000000 2393.00
PERICLASE MG(1)O(1) Pe PERI
ST -569164.00 -601500.00 26.9514 1.125 1.00000
C1 61.10965 -296.199 -621154.000 5844612. 0.00000
V1 3.74774182 0.00003556 -0.06223329 0.00015114 0.00000000
PREHNITE CA(2)AL(2)SI(3)O(12)H(2) Prh PREH COMPROMISE BETW. CP,S
ST -5817852.00 -6198606.33 288.6335 14.016 3.00000
C1 716.05176 -6404.609 -2182456.000 268502784. 0.00000
V1 0.14684075 0.01128353 -0.14269406 0.00000000 0.00000000
PYROPHYLLITE AL(2)SI(4)O(12)H(2) Prl PYPH
ST -5264450.00 -5640780.66 239.4000 12.760 1.00000
C1 665.93066 -5897.422 -4979913.000 661808896. 0.00000
V1 1.26374216 0.00381661 -0.13544906 0.00000000 0.00000000
RUTILE TI(1)O(2) Rt RUTI
ST -891336.00 -944750.00 50.4600 1.882 1.00000
C1 77.83762 0.000 -3367841.000 402940672. 0.00000
V1 2.57156217 0.00015409 -0.04542508 0.00005845 0.00000000
SPHENE CA(1)TI(1)SI(1)O(5) Sph SPHE
ST -2454203.00 -2596652.15 129.2900 5.565 1.00000
C1 234.62025 -1040.294 -5118324.000 591460608. 0.00000
V1 2.52000000 0.00000000 -0.05900000 0.00000000 0.00000000
TALC MG(3)SI(4)O(12)H(2) Tc TALC
ST -5520258.00 -5897386.72 261.2405 13.610 1.00000
C1 664.10522 -5187.172 -2147218.000 -327371776. 0.00000
V1 2.94468112 0.00000000 -0.16988979 0.00056650 0.00000000
WOLLASTONITE CA(1)SI(1)O(3) Wo WOLL
ST -1545999.00 -1631500.00 81.8100 3.983 1.00000
C1 149.07266 -690.295 -3659348.000 484349440. 0.00000
V1 2.81800151 0.00000000 -0.12449159 0.00031132 0.00000000
PS-WOLLASTONITE CA(1)SI(1)O(3) pWo PSWO
ST -1542618.00 -1627427.11 85.2788 4.016 1.00000
C1 141.15611 -417.232 -5857595.000 940734976. 0.00000
V1 2.81800151 0.00000000 -0.12449159 0.00031132 0.00000000
ZIRCON ZR(1)O(4)SI(1) Zc ZIRC
! All properties are provisional estimates not included in JPet 1988 paper
ST 0.00 -2031900.00 84.0300 3.926 3.00000
C1 166.70241 -434.320 -7348579.000 1040502569. 0.00000
!
!========================================================================
! some additions
!========================================================================
!
!PEROVSKITE CA(1)TI(1)O(3) Pvk PVKT
! ST -1575256. -1660630.000 93.640 3.363
! C3 12.49600 0.045156 -6301800.000
! V1 2.72476341 0.00029968 -0.05839748 0.00012303 0.0
!! ST, C3 from Robie et al. (1978) added by Capi
!! V1 is ilmenite
!
! Mg-carpholite dH fitted to Chopin & Schreyer (1983)
! other parameters as in Holland and Powell (1998)
Mg-carpholite SI(2)AL(2)MG(1)O(10)H(4) mcar 1nh
! ST 0.0 -4781240.000 210.0000 10.5900
ST 0.0 -4782240.000 210.0000 10.5900
C3 667.80000 -0.0125590 -1167100.0 -6440.000 0.00
VHP 0.000050000 525.0 0.0000 0.00 0.00000
VH2 10.0000 4.0000 -0.000078750
!
****COMMENTS******
THE NAMES OF THE GASES CAN NOT BE CHANGED AS THEY ARE USED
TO RECOGNIZE WHICH GAS ROUTINE SHOULD BE USED.
ST is the code for standard STATE PROPERTIES - G, H, S, V
UNITS: J, J, J/DEG, J/bar
NOTE: G = H - T*S where S is third law entropy, not
entropy of formation; i.e. its an apparent G of formation
CP calculated with eqn:
Cp = K0 + K1/SQRT(T) + K2/T/T + K3/T/T/T + K4/T + K5*T + K6*T*T
C1 IS THE CODE FOR THE FIRST LINE OF CP TERMS - K0, K1, K2, K3
C2 IS THE CODE FOR THE 2ND LINE OF CP TERMS - K4, K5, K6
C3 IS THE CODE FOR MAIER-KELLY CP EQUATION - K0, K5, K2
UNITS: J AND DEGREES
CP(disorder) calculated by:
Cp(dis) = D0 + D1/SQRT(T) + D2/T/T + D3/T + D4*T + D5*T*T
To get H, S (disorder), this eqn integrated by TMIN and TMAX
D1 IS THE CODE FOR THE FIRST LINE OF DISORDER TERMS - D0, D1, D2, D3
UNITS: J AND DEGREES
D2 IS THE CODE FOR THE SECOND LINE OF DISORDER TERMS - D4, D5, TMIN, TMAX
UNITS: J AND DEGREES
LAMBDA HEAT CAPACITY EQUATION IS:
CP(LAMBDA) = T * (L1 + L2 * T)**2
T1 (T3,T5 FOR MORE TRANSITIONS) IS THE CODE FOR THE FIRST LINE OF TRANSITION
TERMS - T TRANS(K), REFERENCE T FOR INTEGRATION, LAMBDA
CP TERM L1, LAMBDA CP TERM L2, HEAT OF TRANSITION
UNITS: J AND DEGREES
(SEE BERMAN AND BROWN, 1985, CMP, 89, 168-183)
T2 (T4,T6 FOR MORE TRANSITIONS) IS THE CODE FOR THE SECOND LINE OF TRANS
TERMS - DT/DP SLOPE OF TRANSITION, PARAMETER FOR HIGHER
ORDER FIT OF TRANSITION T AS FUNCTION OF P, DV/DT FOR
LOW T POLYMORPH, DV/DP FOR LOW T POLYMORPH
UNITS: JOULES, BARS, DEGREES
VOLUME EQUATION:
V(P,T)/V(1,298) = 1 + V1(T-298) + V2(T-298)**2 + V3(P-1) + V4(P-1)**2
V1 IS THE CODE FOR VOLUME PARAMETERS - V1, V2,V3,V4
NOTE: V1,V2,V3 terms need to be divided by 10**5, V4 by 10**8
!
!========================================================================
! zeolites and very low grade minerals
!========================================================================
!
****MINERAL DATA***** additions by Capi (Very Low Grade)
HEULANDITE CA(1)AL(2)SI(7)H(12)O(24) Heu Nh HEU
! data added by CAPI.
! C1: predicted (B&B 1985), V1: 3 LW + 5 AQZ - 2 AN
! V0: tables of mineralogy
ST 0.00 -10565396.80 734.372937 31.80250
C1 1355.95400 -9155.550 -22391220.00 3124286000.
V1 2.675699696 -0.000266016 -0.079773389 0.000344341
!HEUkreis CA(1)AL(2)SI(7)H(12)O(24) Hek N
!! data added by CAPI. Produces a circle for reaction: Heu=Lmt+3Qz+2w
!! C1: predicted (B&B 1985), V1: 3 LW + 5 AQZ - 2 AN
! ST 0.00 -10558895.34 745.000000 30.84045
! C1 1355.95400 -9155.550 -22391220.00 3124286000.
! V1 2.675699696 -0.000266016 -0.079773389 0.000344341
LAUMONTITE CA(1)AL(2)SI(4)O(16)H(8) Lmt LMT
! data added by CAPI. C1: predicted (B&B 1985), V1: 2 LW + 2 AQZ - AN
! V0: Liou (1970)
!ST -6709847. -7277149.000 434.750 20.3700 0.0
!ST -6709847.00 -7258743.14 470.204082 20.37000
ST -6709847.00 -7261180.50 463.333333 20.37000
C1 917.37700 -6131.500 -14805540.00 2022583000.
V1 2.610028122 -0.000203757 -0.067858953 0.000260199
STILBITE CA(1)AL(2)SI(7)H(14)O(25) Stb STB
! data added by CAPI.
! C1: predicted (B&B 1985), V1: 3.5 LW + 5 AQZ - 2.5 AN
! V0: Liou (1970)
!ST 0.00 0.00 0.000000 32.86740
ST 0.00 -10889880.20 718.666667 32.86740
C1 1443.57100 -9913.690 -22391220.00 3129577000.
V1 2.827265455 -0.000321747 -0.069550624 0.000314168
WAIRAKITE CA(1)AL(2)SI(4)H(4)O(14) Wa WKA
! data added by CAPI. C1: predicted (B&B 1985), V1: LW + 2 AQZ
! V0: Liou (1970)
!ST -6218243.00 -6666476.00 367.600000 19.04640
ST -6218243.00 -6661077.15 377.777042 19.04640
C1 742.14300 -4615.220 -14805540.00 2012001000.
V1 1.969149330 -0.000000000 -0.098971900 0.000343917
!PUMPELLYITE CA(4)AL(5)MG(1)SI(6)O(28)H(7) Pmp PUMP
!! from Evans (1990, Lithos)
! ST 0.00 -14402300.00 584.0000 29.550 3.00000
! C1 1576.00000 -10603.800 -28304200.000 3838351000. 0.00000
! V1 3.46700000 0.00000000 -0.05151600 0.00012880 0.00000000
PUMPELLYITE2 CA(4)AL(5)MG(1)SI(6)O(28)H(7) Pmp2 PUMP
!! Data entered by El-Shazly, Nov, 89, all from Holland.
!ST -13444859. -14384490. 629.000 29.550
!ST -13444859.00 -14371160.00 629.000000 29.55000
ST -13444859.00 -14380034.80 614.757965 29.55000
C1 1720.8 -14620.300 -5998700.0 0.0 0.0
C2 0.0 -0.0249280 0.0 0.0 0.0
V1 2.978003384 0.0 -0.060913706 0.0
!MG-PUMPELLYITE CA(4)AL(4)MG(2)SI(6)O(28)H(8) Mg-Pmp PUMP
!!ST 0.00 0.00 -11.526067 0.00000
! ST 0.00 0.00 20.473933 0.00000
! C1 0.0 0.000 0.0 0.0 0.0
!COM PUMPELLYITE[1]BRUCITE[1]DIASPORE[-1]
!AL-PUMPELLYITE CA(4)AL(6)SI(6)O(28)H(6) Al-Pmp PUMP
!!ST 0.00 0.00 -11.526067 0.00000
! ST 0.00 0.00 -43.526067 0.00000
! C1 0.0 0.000 0.0 0.0 0.0
!COM PUMPELLYITE[1]BRUCITE[-1]DIASPORE[1]
!
!========================================================================
! System Si-Al-O
!========================================================================
!
****MINERAL DATA***** from RG Berman (J.Petrology, v.29, 1988) Si - Al - O
A-QUARTZ SI(1)O(2) aQz zAQTZ
ST -856227.00 -910699.95 41.4600 2.269 1.00000
C1 80.01199 -240.276 -3546684.000 491568384. 0.00000
V1 2.38945698 0.00000000 -0.24339298 0.00101375 0.00000000
T1 848.00 373.00 -0.09186959 0.00024607 0.00
T2 0.023743 0.000000 0.0000000000 0.0000000000 0.00000
B-QUARTZ SI(1)O(2) bQz BQTZ
ST -856227.00 -908626.77 44.2068 2.370 1.00000
C1 80.01199 -240.276 -3546684.000 491568384. 0.00000
V1 0.00000000 0.00000000 -0.12382672 0.00070871 0.00000000
COESITE SI(1)O(2) Cs COES
ST 0.00 -907604.10 39.4244 2.064 1.00000
C1 94.90714 -767.338 -527867.000 26270730. 0.00000
V1 0.73963178 0.00043605 -0.10368217 0.00030000 0.00000000
CORUNDUM AL(2)O(3) Co CORU
ST -1580958.00 -1675699.95 50.8200 2.558 1.00000
C1 155.01888 -828.387 -3861363.000 409083648. 0.00000
V1 2.13415030 0.00047180 -0.03852190 0.00003750 0.00000000
A-CRISTOBALITE SI(1)O(2) aCr ACRS
ST -853448.00 -907753.35 43.3943 2.587 1.00000
C1 83.51360 -374.693 -2455360.000 280072192. 0.00000
V1 2.08240619 0.00000000 -0.25145068 0.00000000 0.00000000
T1 535.15 298.15 -0.14216187 0.00044142 0.00
T2 0.048000 0.000000 0.0000000000 0.0000000000 0.00000
B-CRISTOBALITE SI(1)O(2) bCr BCRS
ST -853448.00 -906377.23 46.0288 2.730 1.00000
C1 83.51360 -374.693 -2455360.000 280072192. 0.00000
V1 0.31892080 0.00000000 -0.10997269 0.00055353 0.00000000
LOW TRIDYMITE SI(1)O(2) lTr LTRD
ST -854407.00 -907749.56 43.7702 2.675 1.00000
C1 75.37267 0.000 -5958095.000 958246144. 0.00000
V1 1.93394983 0.00000000 -0.25084238 0.00000000 0.00000000
T1 390.15 298.15 0.42670490 -0.00144575 130.54
HIGH TRIDYMITE SI(1)O(2) hTr HTRD
ST -854407.00 -907045.12 45.5237 2.737 1.00000
C1 75.37267 0.000 -5958095.000 958246144. 0.00000
V1 0.48286524 0.00000000 -0.07396833 0.00037354 0.00000000
ANDALUSITE AL(2)SI(1)O(5) And ANDA PHASE EQ. FIT
ST -2440136.00 -2589972.17 91.4337 5.147 1.00000
C1 236.47818 -1102.941 -7526810.000 936442368. 0.00000
V1 2.34426462 0.00007189 -0.07700408 0.00019235 0.00000000
KYANITE AL(2)SI(1)O(5) Ky KYAN
ST -2441634.00 -2594220.46 82.4300 4.412 1.00000
C1 262.68478 -2001.407 -1999740.000 -63181880. 0.00000
V1 2.39725295 0.00000000 -0.06459655 0.00000000 0.00000000
SILLIMANITE AL(2)SI(1)O(5) Si zSILL
ST -2437847.00 -2586090.58 95.9300 4.983 1.00000
C1 256.73158 -1887.225 -2977353.000 250959080. 0.00000
V1 1.34307044 0.00000000 -0.07525587 0.00000000 0.00000000
!
!========================================================================
! Olivine
!========================================================================
!
****MINERAL DATA***** from RG Berman (J.Petrology, v.29, 1988) olivine
FAYALITE FE(2)SI(1)O(4) Fa FAYA
ST -1380699.00 -1479360.11 150.9300 4.630 1.00000
C1 248.92812 -1923.855 0.000 -139104096. 0.00000
V1 2.65457019 0.00079482 -0.07300000 0.00000000 0.00000000
FORSTERITE MG(2)SI(1)O(4) Fo FORS
ST -2057238.00 -2174419.92 94.0100 4.366 1.00000
C1 238.64137 -2001.261 0.000 -116243280. 0.00000
V1 2.94636649 0.00088633 -0.07913336 0.00013513 0.00000000
!FeMgOLIVINE SI(1)FE(1)MG(1)O(4) FeMgOl mN
! ST 0.00 0.00 0.000000000 0.0000
! C1 0.00000 0.0000 0.000 0.
! COM FAYALITE[0.5] FORSTERITE[0.5]
!MgFeOLIVINE SI(1)FE(1)MG(1)O(4) MgFeOl mN
! ST 0.00 0.00 0.000000000 0.0000
! C1 0.00000 0.0000 0.000 0.
! COM FAYALITE[0.5] FORSTERITE[0.5]
!FeMgOLI SI(1)FE(1)MG(1)O(4) MgFeOli mN
! ST 0.00 0.00 5.763033600 0.0000
! C1 0.00000 0.0000 0.000 0.
! COM FAYALITE[0.5] FORSTERITE[0.5]
TEPHROITE MN(2)SI(1)O(4) Teph CN
ST 0.000 -1538156.320 400.250070 4.8610000
C1 222.466500 -1280.89000 -2765238.00 334119400.0
!
***** SOLUTION DATA *****
OLIVINEi (IDEAL)2 M(2):Mg,Fe
FAYALITE Fe,Fe
FORSTERITE Mg,Mg
!
!========================================================================
! Spinel
!========================================================================
!
****MINERAL DATA***** from RG Berman (J.Petrology, v.29, 1988) spinel
!HERCYNITE FE(1)AL(2)O(4) Hc HERC
!! Provisional properties not included in JPet 1988 paper; derived Nov., 1989
! ST 0.00 -1960924.20 107.0130 4.075 1.00000
!!C1 288.04547 -2888.494 -0.0 65342480. 0.0
! C1 204.26948 -258.256 -9807509.000 1208308510. 0.00000
! V1 1.58188633 0.00096276 -0.05922832 0.00014699 0.00000000
!!
! recalibrated hercynite by Philipp
!
HERCYNITE FE(1)AL(2)O(4) Hc HERC
ST 0.00 -1949310.55 119.901744 4.08000
C1 251.76740 -2044.441 -1348330.000 131500400. 0.00000
V1 1.58188633 0.00096276 -0.05922832 0.00014699 0.00000000
SPINEL MG(1)AL(2)O(4) Sp SPIN
ST -2178401.00 -2300312.74 84.5346 3.977 2.00000
C1 235.89993 -1766.578 -1710415.000 40616928. 0.00000
V1 2.16913776 0.00050528 -0.04894168 0.00000000 0.00000000
!
***** SOLUTION DATA ***** Solution model by Philipp (provisional...)
SPIN (IDEAL,MARGULES)
SPINEL
HERCYNITE
***** MARGULES PARAMETER *****
SPINEL - HERCYNITE
12 -4661.1144 0.0000 0.0000 0.0000 0.0000
!
!========================================================================
! Ilmenite-Geikelite
!========================================================================
!
****MINERAL DATA***** from RG Berman (J.Petrology, v.29, 1988) ilmenite
ILMENITE FE(1)TI(1)O(3) Ilm ILME
ST 0.00 -1232448.00 108.6275 3.170 1.00000
C1 149.99950 -441.622 -3323694.000 348151040. 0.00000
V1 2.72476341 0.00029968 -0.05839748 0.00012303 0.00000000
GEIKELITE MG(1)TI(1)O(3) Gk GEIK
! All properties are provisional estimates not included in JPet 1988 paper
ST 0.00 -1570700.00 74.5600 3.086 3.00000
C1 131.15910 -18.045 -4728271.000 455266700. 0.00000
V1 2.72476341 0.00029968 -0.05839748 0.00012303 0.00000000
!
***** SOLUTION DATA *****
ILM (IDEAL) M(1):Fe,Mg
ILMENITE Fe
GEIKELITE Mg
!
!========================================================================
! Garnet
!========================================================================
!
****MINERAL DATA***** from RG Berman (J.Petrology, v.29, 1988) garnet
GROSSULAR CA(3)AL(2)SI(3)O(12) Gr gGROS
ST -6270465.00 -6632859.38 255.1500 12.538 1.00000
C1 573.43042 -2039.405 -18887168.000 2319311872. 0.00000
V1 1.89942017 0.00079756 -0.06539136 0.00016350 0.00000000
PYROPE MG(3)AL(2)SI(3)O(12) Py gPYRO
ST -5938444.00 -6286547.62 266.3591 11.316 1.00000
C1 640.71997 -4542.070 -4701900.000 0. 0.00000
V1 2.25186544 0.00037044 -0.05762090 0.00004420 0.00000000
!PYR3010 SI(3)MG(3)AL(2)O(12) P3010 gPYRo
! ST 0.00 00.00 57.43315 0.000 1.00000
! C1 0.00000 0.000 00.000 0. 0.00000
!COM PYROPE[1]
ALMANDINE SI(3)FE(3)AL(2)O(12) Alm gALMA
ST 0.00 -5267216.00 340.0070 11.511 1.00000
C1 573.96191 -1483.127 -29291968.000 5022076928. 0.00000
V1 1.86125012 0.00074539 -0.05577832 0.00003211 0.00000000
SPESSARTINE SI(3)MN(3)AL(2)O(12) spf SPFA
ST 0.00 -5267216.00 340.0070 11.511 1.00000
C1 573.96191 -1483.127 -29291968.000 5022076928. 0.00000
V1 1.86125012 0.00074539 -0.05577832 0.00003211 0.00000000
! ANDRADITE SI(3)CA(3)FE(2)O(12) Andr ANDR
!! All properties are provisional estimates not included in JPet 1988 paper
! ST 0.00 -5757255.18 318.1403 13.209 3.00000
! C1 596.15930 -2530.370 -13292860.000 1395897310. 0.00000
! V1 2.13334810 0.00047020 -0.06369400 0.00022000 0.00000000
!ANDRADITE CA(3)FE(2)SI(3)O(12) Adr ANDR
! ST -5432285. -5781991.734 293.620 13.185
! C3 475.3358 0.065465 12932606.000
! V1 1.89942017 0.00079756 -0.06539136 0.00016350 0.0
!! ST, C3 from Helgeson et al. (1978)
!! V1 is grossular
!
***** SOLUTION DATA *****
GARNET (MARGULES,SITE) M(3):Ca,Mg,Fe,Mn
GROSSULAR Ca,Ca,Ca
PYROPE Mg,Mg,Mg
ALMANDINE Fe,Fe,Fe
SPESSARTINE Mn,Mn,Mn
***** MARGULES PARAMETER *****
GROSSULAR - PYROPE
112 21560. 18.79 .10
122 69200. 18.79 .10
GROSSULAR - ALMANDINE
112 20320. 5.08 .17
122 2620. 5.08 .09
PYROPE - ALMANDINE
112 230. 0.00 .01
122 3720. 0.00 .06
GROSSULAR - PYROPE - ALMANDINE
123 58825.0 23.87 0.265 0 0
GROSSULAR - PYROPE - SPESSARTINE
123 45424. 18.79 0.1
GROSSULAR - ALMANDINE - SPESSARTINE
123 11470. 5.08 0.13
PYROPE - ALMANDINE - SPESSARTINE
123 1975. 0.0 0.035
!
!========================================================================
! Feldspar
!========================================================================
!
****MINERAL DATA***** from RG Berman (J.Petrology, v.29, 1988) feldspar
ALBITE NA(1)AL(1)SI(3)O(8) Ab 2ALBI
SPC ABSAL
ST -3701276.00 -3921618.20 224.4120 10.083 1.00000
C1 393.63574 -2415.498 -7892826.000 1070636032. 0.00000
V1 2.63072032 0.00032407 -0.19446932 0.00048611 0.00000000
D1 0.00000 0.00000 0.00000 0.00000 0.00000
! Above D1 line needed to compute stable ordering relative to high albite
! using Landau-type function of Salje et al. (1985) in external ALBITE routine
!ALBITE (HIGH) NA(1)AL(1)SI(3)O(8) hAb HALB
! ST -3701276.00 -3921618.20 224.4120 10.083 1.00000
! C1 393.63574 -2415.498 -7892826.000 1070636032. 0.00000
! V1 2.63072032 0.00032407 -0.19446932 0.00048611 0.00000000
!ALBITE (LOW) NA(1)AL(1)SI(3)O(8) lAb LALB
! ST -3711054.00 -3935100.10 207.4430 10.043 1.00000
! C1 393.63574 -2415.498 -7892826.000 1070636032. 0.00000
! V1 2.63072032 0.00032407 -0.19446932 0.00048611 0.00000000
K-FELDSPAR K(1)AL(1)SI(3)O(8) Kfs KSPR
ST -3747058.00 -3970790.78 214.1451 10.869 1.00000
C1 381.37231 -1941.045 -12037252.000 1836425472. 0.00000
V1 1.51450750 0.00054850 -0.18045110 0.00051120 0.00000000
D1 282.98291 -4831.37500 3620706.00000 0.00000 0.00000
D2 -0.15733 0.00003 298.15000 1436.15000 410629.63000
!K-FELDSPAR (HIGH) K(1)AL(1)SI(3)O(8) Sa HKPR
! ST -3740412.00 -3959703.78 229.1571 10.896 1.00000
! C1 381.37231 -1941.045 -12037252.000 1836425472. 0.00000
! V1 1.51450750 0.00054850 -0.18045110 0.00051120 0.00000000
!K-FELDSPAR (LOW) K(1)AL(1)SI(3)O(8) Mcr LKPR
! ST -3747058.00 -3970790.78 214.1451 10.869 1.00000
! C1 381.37231 -1941.045 -12037252.000 1836425472. 0.00000
! V1 1.51450750 0.00054850 -0.18045110 0.00051120 0.00000000
ANORTHITE AL(2)CA(1)SI(2)O(8) An 2ANOR
ST -4002668.00 -4228730.47 200.1861 10.075 1.00000
C1 439.36938 -3734.149 0.000 -317023232. 0.00000
V1 1.09181141 0.00041985 -0.12724268 0.00031762 0.00000000
!
***** SuLUTION DATA ***** Al-avoidence external
FSP (EXT,MARGULES)
ALBITE
K-FELDSPAR
ANORTHITE
***** SOLUTION DATA ***** Al-avoidence with sites
FSP2 (SITE,MARGULES) A(1):Na,Ca,K - T1(2):Al,Si
ALBITE Na - Al,Si
K-FELDSPAR K - Al,Si
ANORTHITE Ca - Si,Si
!
***** MuRGULES PARAMETER *****
! data by Benisek et al. (2010)
ALBITE - K-FELDSPAR
122 22945.000 10.3000 0.327 0
112 17711.000 10.3000 0.461 0
ALBITE - ANORTHITE
122 40000.0000 16.4 0.069 0
112 14000.0000 4.7 -0.049 0
ANORTHITE - K-FELDSPAR
112 90600.0000 0.0 -0.257 0
122 60300.0000 0.0 -0.210 0
ALBITE - ANORTHITE - K-FELDSPAR
123 122778.00 20.85 0.1705 0
123 87300.00 93.90 -0.4670 0
!
***** MuRGULES PARAMETER *****
! data by Benisek et al. (2004)
ALBITE - K-FELDSPAR
122 19550.000 10.5000 0.327 0
112 22820.000 6.3000 0.461 0
ALBITE - ANORTHITE
122 31000.0000 4.5 0.069 0
112 9800.0000 -1.7 -0.049 0
ANORTHITE - K-FELDSPAR
112 90600.0000 29.5 -0.257 0
122 60300.0000 11.2 -0.210 0
ALBITE - ANORTHITE - K-FELDSPAR
123 125035.00 30.15 -0.2965 0
!
***** MARGULES PARAMETER *****
! data by Fuhrman and Lindsley(1988)
ALBITE - K-FELDSPAR
122 18810.000 10.3000 0.394 0
112 27320.000 10.3000 0.394 0
ALBITE - ANORTHITE
122 28226.0000 0 0 0
112 8471.0000 0 0 0
ANORTHITE - K-FELDSPAR
122 52468.0000 0 -0.120 0
112 47396.0000 0 0 0
ALBITE - ANORTHITE - K-FELDSPAR
123 100045.500 10.3 -0.76 0
!
! 123 91345.500 10.3 0.334 0
! 123 8700.000 0.0 -1.094 0
!
***** MuRGULES PARAMETER *****
! data of Ghiorso (1984)
ALBITE - POTASSIUM FELDSPAR
122 30978.336 21.5476 0 0
112 17062.352 0.0000 0 0
12 0.000 0.0000 0.3606608 0
ALBITE - ANORTHITE
122 28225.6824 0 0 0
112 8471.3448 0 0 0
ANORTHITE - POTASSIUM FELDSPAR
122 27983.0104 -20.21022624 0 0
112 67469.0920 11.062496 0 0
ALBITE - ANORTHITE - POTASSIUM FELDSPAR
123 76225.45088 8.43465112 0 0
!
!========================================================================
! Biotite
!========================================================================
!
****MINERAL DATA***** from RG Berman (J.Petrology, v.29, 1988) mica
ANNITE SI(3)FE(3)AL(1)K(1)H(2)O(12) Ann ANNI
ST 0.00 -5142000.00 421.0100 15.483 3.00000
! Properties derived by Berman (Am. Min., 1990), consistent with gt-bio expts;
! modified by McMullin (1990) on basis of SGAM barometry;
! 5142000 to be used with McMullin (1990) biotite mixing properties
!ST 0.0 -5137500.00 420.0000 15.408 3.0
! -5137500 J is value to be used with ideal FeMg mixing in garnet
!ST 0.0 -5160566.00 420.0000 15.477 3.0
! H = -5160566 J is value to be used with Gang & Sax (84) garnet model
C1 727.20800 -4775.040 -13831900.000 2119060000. 0.00000
V1 3.44473262 0.00000000 -0.16969784 0.00000000 0.00000000
MN-BIO SI(3)MN(3)AL(1)K(1)H(2)O(12) mnbi MNBI
ST 0.00 -5142000.00 421.0100 15.483 3.00000
! Diese Daten gleichen denen von Annit (von Thorsten eingefuegt)
C1 727.20800 -4775.040 -13831900.000 2119060000. 0.00000
V1 3.44473262 0.00000000 -0.16969784 0.00000000 0.00000000
PHLOGOPITE K(1)MG(3)AL(1)SI(3)O(12)H(2) Phl PHLG
ST 0.00 -6210391.00 334.3460 14.977 1.00000
C1 610.37988 -2083.781 -21533008.000 2841040896. 0.00000
V1 3.44473262 0.00000000 -0.16969784 0.00000000 0.00000000
!
***** SOLUTION DATA *****
BIOTITE (IDEAL)3 M(3):Fe,Mg,Mn
PHLOGOPITE Mg,Mg,Mg
ANNITE Fe,Fe,Fe
! MN-BIO Mn,Mn,Mn
!
!========================================================================
! White mica
!========================================================================
!
****MINERAL DATA***** from RG Berman (J.Petrology, v.29, 1988) mica
MARGARITE CA(1)AL(4)SI(2)O(12)H(2) Mrg MARG
ST -5850897.00 -6236602.96 265.0841 12.958 1.00000
C1 699.79565 -5587.105 -6807750.000 734321664. 0.00000
V1 2.10185214 0.00124556 -0.11549004 0.00028862 0.00000000
!
****MuNERAL DATA***** from RG Berman (J.Petrology, v.29, 1988)
MUSCOVITE K(1)AL(3)SI(3)O(12)H(2) Ms 1MUSC
ST -5599784.00 -5976740.12 293.1567 14.087 1.00000
C1 651.48926 -3873.229 -18523200.000 2742469376. 0.00000
V1 3.35273302 0.00000000 -0.17169021 0.00042947 0.00000000
!CELADONITE SI(4)AL(1)MG(1)K(1)O(12)H(2) cel 1
! ST 0.0 -5832415.000 288.5270 13.8700
! C3 741.200 -.01874800 -2368800.0 -6616.900
! V1 2.769296315 .000000000 -.113612157 .000000000
CELADONITE SI(4)AL(1)MG(1)K(1)O(12)H(2) cel 1
! Meyre 98 nach Massonne und Szpurka 1997 (Lithos 41, 229 ff.)
ST 0.0 -5832415.000 288.5270 13.8700
C1 645.915 -4129.540000 -13864470. 1978171000.
V1 3.352700000 .000000000 -0.17169000 .000429500
FE-CELADONITE SI(4)AL(1)FE(1)K(1)O(12)H(2) fcel 1
ST 0.0 -5492287.000 303.1480 13.9620
C1 664.755 -4553.1100 -12459890. 1871056000
V1 3.300000000 .000000000 -.170000000 .000400000
!FE-CELADONITE SI(4)AL(1)FE(1)K(1)O(12)H(2) fcel 1
! ST 0.0 -5483730.000 328.0000 14.1500
! C3 756.300 -.01914700 -1586100.0 -6928.700
! V1 2.791519435 .000000000 -.119434629 .000000000
PARAGONITE NA(1)AL(3)SI(3)O(12)H(2) Pg 1PARA
ST -5560976.00 -5944208.14 277.6991 13.216 1.00000
C1 577.56934 -1472.801 -32214368.000 5050077184. 0.00000
V1 3.94240896 0.00059701 -0.19732335 0.00049335 0.00000000
!
***** SuLUTION DATA ****** (Massonne & Szpurka 1997)
PHNG (MARGULES,IDEAL)
MUSCOVITE
CELADONITE
FE-CELADONITE
***** MuRGULES PARAMETER *****
MUSCOVITE - CELADONITE
112 0 15.92 0.187 0 1
122 0 58.598 0.735 0 1
MUSCOVITE - FE-CELADONITE
112 0 24.083 0.3394 0 1
122 0 24.083 0.3394 0 1
!
!***** SOLUTION DATA *****
!WHITE MICA (MARGULES)
! PARAGONITE
! MUSCOVITE
!***** MARGULES PARAMETER *****
!PARAGONITE - MUSCOVITE
!! data of Chatterjee and Froese (1975)
!112 19456.0184 -1.6543536 -0.456056 0
!122 12230.2504 -0.7104432 0.665256 0
!
****MINERAL DATA*****
!
! White mica by Luki Keller
!
MUSCOVITE K(1)AL(3)SI(3)O(12)H(2) Ms
ST -5599784.00 -5976740.12 293.1567 14.087 1.00000
C1 651.48926 -3873.229 -18523200.000 2742469376. 0.00000
V1 3.35273302 0.00000000 -0.17169021 0.00042947 0.00000000
PARAGONITE NA(1)AL(3)SI(3)O(12)H(2) Pg
ST -5560976.00 -5944208.14 277.6991 13.216 1.00000
C1 577.56934 -1472.801 -32214368.000 5050077184. 0.00000
V1 3.94240896 0.00059701 -0.19732335 0.00049335 0.00000000
MCELADONITE K(1)MG(1)AL(1)SI(4)O(12)H(2) MgC
ST 0.00 -5834270.00 288.527 13.870 1.00000
C1 645.915 -4129.54 -13864470.000 1978171000. 0.00000
V1 3.3527 0.00 -0.17169 0.0004295 0.0000000
FCELADONITE K(1)FE(1)AL(1)SI(4)O(12)H(2) FeC
ST 0.00 -5492287.00 303.148 13.962 1.00000
C1 664.755 -4553.11 -12459890.000 1871056000. 0.00000
V1 3.3000 0.00 -0.17000 0.0004000 0.0000000
!
***** SOLUTION DATA ***** White mica by Luki Keller
WHITE MICA (MARGULES,IDEAL) Z(4):Si,Al
MUSCOVITE Si,Si,Si,Al
PARAGONITE Si,Si,Si,Al
MCELADONITE Si,Si,Si,Si
FCELADONITE Si,Si,Si,Si
***** MARGULES PARAMETER *****
MUSCOVITE - PARAGONITE
122 15050.0 0 0.452 0 0
112 6150.0 0 0.452 0 0
MUSCOVITE - MCELADONITE
122 0.0000 58.598 0.735 0 0
112 0.0000 15.920 0.187 0 0
MUSCOVITE - FCELADONITE
112 0.0000 24.083 0.3394 0 0
122 0.0000 24.083 0.3394 0 0
MCELADONITE - PARAGONITE
112 14224 0.00 0.52 0 0
122 14224 0.00 0.52 0 0
FCELADONITE - PARAGONITE
112 16620 0.00 0.34 0 0
122 16620 0.00 0.34 0 0
MUSCOVITE - FCELADONITE - PARAGONITE
123 27220 24.083 1.13 0 0
MUSCOVITE - MCELADONITE - PARAGONITE
123 24824 37.260 1.43 0 0
MCELADONITE - PARAGONITE - FCELADONITE
123 30844 0.00 0.86 0 0
MUSCOVITE - MCELADONITE - FCELADONITE
123 0.00 61.34 0.8 0 0
!
!========================================================================
! Staurolite
!========================================================================
!
****MINERAL DATA***** from RG Berman (J.Petrology, v.29, 1988) staurolite
FE-STAUROLITE SI(7.5)FE(4)AL(18)H(4)O(48) FSt FEST
! Provisional properties not included in JPet 1988 paper; derived Nov., 1989
ST 0.00 -23765364.00 1005.3270 44.676 3.00000
C1 2577.67407 -16265.359 -61185888.000 8692047872. 0.00000
! Cp predicted with equation (11) of Berman and Brown (1985)
V1 1.79363460 0.00090396 -0.08000000 0.00000000 0.00000000
MG-STAUROLITE SI(7.5)AL(18)MG(4)O(48)H(4) MSt MGST
! ST 0.0 -25010130.000 1010.0000 44.2600
ST 0.0 -25112909.000 905.3960 44.2600
C1 2820.500 -24126.000 -13774000.0 0.000
C2 0.000 -0.05936600 0.000 0.000
V1 0.723000452 0.000000000 -0.072300045 0.000000000
! Modified MG-STAU from Holland and Powell
!
***** SOLUTION DATA ***** Nagels solution model for staurolite
STAU (MARGULES,SITE) T(4):Fe,Mg
FE-STAUROLITE Fe,Fe,Fe,Fe
MG-STAUROLITE Mg,Mg,Mg,Mg
***** MARGULES PARAMETER *****
FE-STAUROLITE - MG-STAUROLITE
21 -20000.0 0.0
!
!========================================================================
! Orthopyroxene
!========================================================================
!
****MINERAL DATA***** from RG Berman (J.Petrology, v.29, 1988) pyroxene
ORTHOENSTATITE MG(2)SI(2)O(6) En OREN
ST 0.00 -3091104.50 132.3400 6.266 1.00000
C1 333.15900 -2401.176 -4541120.000 558300672. 0.00000
V1 2.46558172 0.00074670 -0.07493458 0.00004467 0.00000000
!ORTHOENSTATITE MG(1)SI(1)O(3) En OREN
! ST -1459103.00 -1545552.25 66.1700 3.133 1.0
! C1 166.57950 -1200.588 -2270560.000 279150336. 0.0
! V1 2.46558172 0.00074670 -0.07493458 0.00004467 0.0
PROTOENSTATITE MG(1)SI(1)O(3) pEn PREN
ST -1457873.00 -1543958.64 67.4376 3.242 1.00000
C1 166.57950 -1200.588 -2270560.000 279150336. 0.00000
V1 1.68317000 0.00116650 -0.07496000 0.00004480 0.00000000
FERROSILITE SI(2)FE(2)O(6) Fsl FESI
ST 0.00 -2388749.20 191.7644 6.592 1.00000
C1 338.12116 -2385.942 -4194276.000 585065920. 0.00000
V1 3.18076456 0.00075850 -0.09899879 0.00000000 0.00000000
!FERROSILITE SI(1)FE(1)O(3) Fsl FESI
! ST -1118017.00 -1194374.61 95.8822 3.296 1.0
! C1 169.06058 -1192.971 -2097138.000 292532960. 0.0
! V1 3.18076456 0.00075850 -0.09899879 0.0 0.0
!
! Orthopyroxene data by Philipp Hunziker
!
MG.AL-PYROXENE MG(1)AL(2)SI(1)O(6) MgTs
ST 0.00 -3186328.63 129.940000 5.88800
C1 330.418 -2166.493 -6251535.000 715160880.00
V1 2.46558172 0.00074670 -0.07493458 0.00004467 0.00000000
MG.FE-PYROXENE SI(2)FE(1)MG(1)O(6) femg 1
ST 0.00 -2745997.2500 163.8100 6.429 1.00000
C1 335.64008 -2393.559 -4367698.000 571683296. 0.00000
V1 2.82317314 0.00075260 -0.08696670 0.00000000 0.00000000
FE.MG-PYROXENE SI(2)MG(1)FE(1)O(6) mgfe 1
ST 0.00 -2717321.6540 173.9300 6.429 1.00000
C1 335.64008 -2393.559 -4367698.000 571683296. 0.00000
V1 2.82317314 0.00075260 -0.08696670 0.00000000 0.00000000
FE.AL-PYROXENE FE(1)AL(2)SI(1)O(6) FeTs
ST 0.00 -2824362.94 174.563000 6.41800
C1 332.899 -2158.876 -6078113. 728543504.
V1 3.18076456 0.00075850 -0.09899879 0.00000000 0.00000000
!
***** SOLUTION DATA ***** Solution model for Opx (Hunziker)
OPX (IDEAL) M1(1):Mg,Fe - M2(1):Mg,Fe,Al
ORTHOENSTATITE Mg - Mg
FERROSILITE Fe - Fe
MG.FE-PYROXENE Mg - Fe
FE.MG-PYROXENE Fe - Mg
MG.AL-PYROXENE Mg - Al
FE.AL-PYROXENE Fe - Al
!
!========================================================================
! Clinopyroxene
!========================================================================
!
****MINERAL DATA*****
!CLINOENSTATITE MG(1)SI(1)O(3) cEn CLEN
! ST -1458699.00 -1545926.25 66.3251 3.131 1.00000
! C1 139.95824 -497.034 -4400237.000 535708928. 0.00000
! V1 2.19150000 0.00074920 -0.07496000 0.00004480 0.00000000
DIOPSIDE MG(1)CA(1)SI(2)O(6) Di DIOP
ST -3027516.00 -3200583.50 142.5000 6.620 1.00000
C1 305.41333 -1604.931 -7165973.000 921837568. 0.00000
V1 2.77951964 0.00083082 -0.08724773 0.00017069 0.00000000
JADEITE NA(1)AL(1)SI(2)O(6) Jd JADE
ST -2845291.00 -3025118.24 133.5739 6.034 1.00000
C1 311.29297 -2005.101 -5350264.000 662566912. 0.00000
V1 2.31184298 0.00025785 -0.08598017 0.00021488 0.00000000
HEDENBERGITE CA(1)FE(1)SI(2)O(6) Hd HEDE
! All properties are provisional estimates not included in JPet 1988 paper
ST 0.00 -2837144.22 174.3314 6.795 3.00000
C1 296.70260 -1239.150 -7416518.000 893501213. 0.00000
V1 2.37026500 0.00050310 -0.07900000 0.00003923 0.00000000
!HEDENBERGITE CA(1)FE(1)SI(2)O(6) Hed HED
! ST -2676278. -2840727.053 170.403 6.827
! C3 229.4785 0.034206 6284386.800
! V1 2.77951964 0.00083082 -0.08724773 0.00017069 0.0
!! ST, C3 from Helgeson et al. (1978)
!! V1 is diopside
CA-AL PYROXENE CA(1)AL(2)SI(1)O(6) CTs CATS
ST -3120980.00 -3298766.68 140.7508 6.356 1.00000
C1 310.69775 -1671.627 -7455263.000 948781568. 0.00000
V1 2.22504720 0.00052863 -0.08700441 0.00021712 0.00000000
!
***** SOLUTION DATA ***** Omphacite, Meyre et al. ()
OMPH (MARGULES,IDEAL)
DIOPSIDE
JADEITE
HEDENBERGITE
***** MARGULES PARAMETER ******
JADEITE - DIOPSIDE
112 25587.5749 15.9243 0.0000 0.0000 1.0000
122 5920.2458 10.7164 0.0000 0.0000 1.0000
DIOPSIDE - HEDENBERGITE
112 31448.6254 28.8272 0.0000 0.0000 1.0000
122 -65167.1449 -45.4535 0.0000 0.0000 1.0000
JADEITE - HEDENBERGITE
112 3670 9 0 0 1
122 15090 8 0 0 1
! Jd-Hd from Massonne (1992)
!
!========================================================================
! Testing P2/n
!========================================================================
!
****MuNERAL DATA*****
AA MG(1)CA(1)SI(2)O(6) aa aaj
ST 0.00 0.00 0.000 0.0000
C1 0.00000 0.0000 0.000 0.
BB NA(1)AL(1)SI(2)O(6) bb bbj
ST 0.00 0.00 0.000 0.0000
C1 0.00000 0.0000 0.000 0.
***** SuLUTION DATA *****
P2/n (EXT)
AA
BB
***** SaLUTION DATA *****
CPX (SKIP,MARGULES)
CLINOENSTATITE
DIOPSIDE
***** MuRGULES PARAMETER *****
CLINOENSTATITE - DIOPSIDE
112 31216.0 0 -0.0061 0 0
122 25484.0 0 0.0812 0 0
!
!========================================================================
! Chlorite
!========================================================================
!
****MINERAL DATA***** from RG Berman (J.Petrology, v.29, 1988) with
!CLINOCHLORE MG(5)AL(2)SI(3)O(18)H(8) Chl CLIN
! ST -8259134.00 -8909589.84 435.1545 21.147 2.00000
! C1 1214.28442 -11217.129 0.000 -1256253184. 0.00000
! V1 2.64515406 0.00000000 -0.18194707 0.00000000 0.00000000
! MN-CHLOR SI(3)AL(2)MN(5)O(18)H(8) mnch 3
! !from tcds
! ST 0.0 -7149480.000 559.0000 21.3400
! C3 1237.400 .01359400 -3743000.0 -11250.000
! V1 2.764761012 .000000000 -.178069353 .000000000
!
! chlorite data by Philipp Hunziker
!
AMESITE MG(4)AL(4)SI(2)O(18)H(8) Ame 1
ST 0.00 -9022957.58 403.470000 20.52000
C1 1144.45 -8327.2 -20058040. 2820786000.
V1 2.64515406 0.00000000 -0.18194707 0.00000000
PENNINITE MG(5.5)AL(1)SI(3.5)O(18)H(8) Pen 1
ST 0.00 -8822346.74 437.800000 21.37000
C1 1130.32650 -8044.360 -18920800.00 2575230000.00
V1 2.64515406 0.0 -0.18194707 0.0 0.0
FEAMESITE FE(4)AL(4)SI(2)O(18)H(8) FeAm 1
ST 0.00 -7611289.00 505.000000 20.75800
C1 1219.81 -10021.5 -14439720. 2392326000.0
V1 2.64515406 0.0 -0.18194707 0.0
FEPENNINITE FE(5.5)AL(1)SI(3.5)O(18)H(8) FeP 1
ST 00000000.00 -6862527.00 603.8400 21.700
! ST 00000000.00 -6864527.00 603.8400 21.700
C1 1233.9465 -10373.995 -11195610.00 1986097500.00
V1 2.64515406 0.0 -0.18194707 0.0
!
***** SOLUTION DATA ***** Chlorite solution model.
CHL4 (IDEAL)4 T5(4):AL,SI
AMESITE AL,AL,AL,AL
PENNINITE AL,SI,SI,SI
FEAMESITE AL,AL,AL,AL
FEPENNINITE AL,SI,SI,SI
!
!========================================================================
! Cordierite
!========================================================================
!
****MINERAL DATA***** from RG Berman (J.Petrology, v.29, 1988) cordierite
CORDIERITE MG(2)AL(4)SI(5)O(18) Cd CORD
ST -8652708.00 -9158726.56 417.9700 23.311 2.00000
C1 954.38650 -7962.274 -2317258.000 -370214090. 0.00000
V1 0.30028742 0.00018017 -0.11582515 0.00000000 0.00000000
!
HY_CORDIERITE MG(2)AL(4)SI(5)O(18)H(4)O(2) hCd CORD
ST 0.00 -563493.80 156.6650 2.3060 2.00000
C1 270.73860 -9633.768 -8152144.000 0. 0.00000
C2 133752.00 -0.009046 0.000 0. 0.00000
COM CORDIERITE[1]
!
!
FE_CORDIERITE FE(2)AL(4)SI(5)O(18) fCd FECD
! Provisional properties not included in JPet 1988 paper; derived Nov., 1989
ST 0.00 -8469858.00 453.0550 23.707 3.00000
C1 975.30000 -8013.440 -9618520.000 1189278000. 0.00000
V1 1.26000000 0.00000000 -0.11582515 0.00000000 0.00000000
HY_Fe_CORDIERITE FE(2)AL(4)SI(5)O(18)H(4)O(2) hfCd CORD
ST 0.00 -563493.80 156.6650 2.3060 2.00000
C1 270.73860 -9633.768 -8152144.000 0. 0.00000
C2 133752.00 -0.009046 0.000 0. 0.00000
COM FE_CORDIERITE[1]
!
***** SOLUTION DATA *****
CORD (IDEAL)2
CORDIERITE
HY_CORDIERITE
FE_CORDIERITE
HY_Fe_CORDIERITE
!
!========================================================================
! Chloritoid
!========================================================================
!
****MuNERAL DATA***** from RG Berman (J.Petrology, v.29, 1988) chloritoide
MGCTD SI(1.00)AL(2.00)MG(1.00)O(7.00)H(2.00) mctd 1
! ST 0.0 -3559510.000 132.0000 6.8750