[display/print with fixed width font (eg Monaco)] calculation type :phase diagram calculations (from prefs) datafile being read: "tcdp1.txt" reading phase info from this file ... chl bi cd st ctd g mu opx liq ksp ky sill and q H2O log output is in "tclog.txt" main output is in the file, "tcop1.txt" other (eg drawpd) output is in the file, "tcdr.txt" THERMOCALC 3.26 running at 18.29 on Fri 14 Jul,2006, with tcds55.txt produced at 19:29:59 on 22 Nov 2003 (with sigma fit = 1.067) phases ignored this run: cd ctd opx liq ksp sill and (from script) choose from: chl bi st g mu ky q H2O mu q H2O - automatically included (from script) which phases : (nothing input) fluid in excess (from script) phases in excess: mu q (from script) variance of required equilibria : 0 = invariant 1 = univariant 2 = divariant ... n = n-variant variance : 1 specification of PT window: P range over which T of reactions to be calculated P window: P low,high : (nothing input) T window within which reactions expected to lie T window: T low,high : (nothing input) P window :0.1 <-> 15 kbar :P interval : 2 -------------------------------------- chl site fractions x(Fe,M23) = x x(Mg,M23) = 1 - x x(Al,M1) = y - Q x(Fe,M1) = x (1 - y + Q) x(Mg,M1) = (1 - x) (1 - y + Q) x(Al,M4) = y + Q x(Fe,M4) = x (1 - y - Q) x(Mg,M4) = (1 - x) (1 - y - Q) x(Al,T2) = y x(Si,T2) = 1 - y proportions afchl = 1 - y - Q clin = 2Q + (-2/5 x) (3 - y) daph = 2/5 x (3 - y) ames = y - Q ideal mixing actvities afchl = x(Mg,M23)^4 x(Mg,M1) x(Mg,M4) x(Si,T2)^2 clin = 4 x(Mg,M23)^4 x(Mg,M1) x(Al,M4) x(Al,T2) x(Si,T2) daph = 4 x(Fe,M23)^4 x(Fe,M1) x(Al,M4) x(Al,T2) x(Si,T2) ames = x(Mg,M23)^4 x(Al,M1) x(Al,M4) x(Al,T2)^2 non-ideality by symmetric formalism W(afchl,clin) = 18.0 W(afchl,daph) = 14.5 W(afchl,ames) = 20.0 W(clin,daph) = 2.5 W(clin,ames) = 18.0 W(daph,ames) = 13.5 -------------------------------------- bi site fractions x(Al,M1) = y x(Fe,M1) = x (1 - y) + 2/3 Q x(Mg,M1) = (1 - x) (1 - y) + (-2/3 Q) x(Fe,M2) = x - 1/3 Q x(Mg,M2) = 1 - x + 1/3 Q x(Al,T1) = 1/2 + 1/2 y x(Si,T1) = 1/2 - 1/2 y proportions phl = (1 - x) (1 - y) + (-2/3 Q) ann = x - 1/3 Q east = y obi = (-x) y + Q ideal mixing actvities phl = 4 x(Mg,M1) x(Mg,M2)^2 x(Al,T1) x(Si,T1) ann = 4 x(Fe,M1) x(Fe,M2)^2 x(Al,T1) x(Si,T1) east = x(Al,M1) x(Mg,M2)^2 x(Al,T1)^2 obi = 4 x(Fe,M1) x(Mg,M2)^2 x(Al,T1) x(Si,T1) non-ideality by symmetric formalism W(phl,ann) = 9.0 W(phl,east) = 10.0 W(phl,obi) = 3.0 W(ann,east) = -1.0 W(ann,obi) = 6.0 W(east,obi) = 10.0 dependent ("make") end-members obi = 2/3 phl + 1/3 ann DQF increments obi = -10.73 -------------------------------------- st site fractions x(Mg) = 1 - x x(Fe) = x proportions mst = 1 - x fst = x ideal mixing actvities mst = x(Mg)^4 fst = x(Fe)^4 non-ideality by symmetric formalism W(mst,fst) = -8.0 -------------------------------------- g site fractions x(Mg) = 1 - x x(Fe) = x proportions py = 1 - x alm = x ideal mixing actvities py = x(Mg)^3 alm = x(Fe)^3 non-ideality by symmetric formalism W(py,alm) = 2.5 -------------------------------------- mu site fractions x(Al,M2A) = y x(Mg,M2A) = (1 - x) (1 - y) x(Fe,M2A) = x (1 - y) x(Al,T1) = 1/2 y x(Si,T1) = 1 - 1/2 y proportions mu = y cel = (1 - x) (1 - y) fcel = x (1 - y) activities mu = 4 x(Al,M2A) x(Al,T1) x(Si,T1) cel = x(Mg,M2A) x(Si,T1)^2 fcel = x(Fe,M2A) x(Si,T1)^2 <==================================================> phases : chl, bi, st, g, (mu, q, fluid) or [ky] P(kbar) T(¡C) x(chl) y(chl) Q(chl) x(bi) y(bi) Q(bi) x(st) x(g) x(mu) y(mu) 221 6.10 562.1 0.6782 0.6111 0.3888 0.7509 0.4313 0.1573 0.9413 0.9371 0.6660 0.9150 24chl + 29g + 65mu = 65bi + 10st + 113q + 75H2O 8.10 591.7 0.5058 0.5787 0.4213 0.5505 0.2867 0.3438 0.8822 0.8725 0.5084 0.8843 28chl + 20g + 61mu = 61bi + 10st + 92q + 94H2O 10.10 614.1 0.3901 0.5619 0.4381 0.4293 0.2224 0.4299 0.8195 0.8027 0.4018 0.8526 29chl + 19g + 61mu = 61bi + 10st + 88q + 96H2O 12.10 634.1 0.2973 0.5507 0.4492 0.3398 0.1907 0.4410 0.7461 0.7191 0.3136 0.8189 29chl + 20g + 61mu = 61bi + 10st + 91q + 95H2O 14.10 653.0 0.2187 0.5427 0.4573 0.2617 0.1717 0.3984 0.6576 0.6167 0.2362 0.7837 29chl + 21g + 62mu = 62bi + 10st + 95q + 94H2O 15.00 661.2 0.1871 0.5397 0.4602 0.2280 0.1650 0.3653 0.6116 0.5635 0.2041 0.7677 28chl + 21g + 63mu = 63bi + 10st + 97q + 94H2O <==================================================> phases : chl, bi, st, ky, (mu, q, fluid) or [g] P(kbar) T(¡C) x(chl) y(chl) Q(chl) x(bi) y(bi) Q(bi) x(st) x(mu) y(mu) 0.10 427.8 0.3513 0.6083 0.3917 0.4310 0.4549 0.3407 0.8864 0.3046 0.9805 47chl + 10st + 105mu = 105bi + 202ky + 53q + 208H2O 2.10 538.5 0.3516 0.6575 0.3424 0.4306 0.5011 0.2736 0.8370 0.3444 0.9626 43chl + 10st + 98mu = 98bi + 192ky + 52q + 192H2O 4.10 568.6 0.3419 0.6199 0.3801 0.4094 0.4236 0.3194 0.8177 0.3456 0.9431 50chl + 10st + 110mu = 110bi + 209ky + 56q + 219H2O 6.10 590.6 0.3333 0.5943 0.4057 0.3915 0.3506 0.3620 0.8008 0.3431 0.9188 57chl + 10st + 122mu = 122bi + 227ky + 60q + 247H2O 8.10 607.8 0.3250 0.5758 0.4242 0.3762 0.2868 0.3975 0.7852 0.3383 0.8897 64chl + 10st + 134mu = 134bi + 244ky + 65q + 275H2O 10.10 621.6 0.3167 0.5617 0.4383 0.3625 0.2330 0.4252 0.7701 0.3319 0.8560 70chl + 10st + 146mu = 146bi + 259ky + 72q + 300H2O 12.10 632.9 0.3080 0.5507 0.4493 0.3495 0.1884 0.4454 0.7550 0.3239 0.8182 76chl + 10st + 157mu = 157bi + 273ky + 81q + 325H2O 14.10 642.2 0.2987 0.5418 0.4581 0.3366 0.1520 0.4589 0.7397 0.3146 0.7770 82chl + 10st + 169mu = 169bi + 285ky + 93q + 348H2O 15.00 646.0 0.2942 0.5384 0.4615 0.3306 0.1378 0.4629 0.7325 0.3099 0.7576 85chl + 10st + 175mu = 175bi + 290ky + 100q + 359H2O <==================================================> phases : chl, bi, g, ky, (mu, q, fluid) or [st] P(kbar) T(¡C) x(chl) y(chl) Q(chl) x(bi) y(bi) Q(bi) x(g) x(mu) y(mu) 11 4.10 549.5 0.7482 0.7422 0.2577 0.8369 0.6149 0.06546 0.9538 0.7256 0.9515 10chl + 12g + 35mu = 35bi + 38ky + 46q + 40H2O 6.10 582.0 0.5585 0.6313 0.3687 0.6316 0.4383 0.2194 0.8948 0.5552 0.9226 22chl + 10g + 54mu = 54bi + 64ky + 54q + 88H2O 8.10 603.5 0.4490 0.5886 0.4114 0.4970 0.3097 0.3534 0.8408 0.4571 0.8914 43chl + 10g + 92mu = 92bi + 114ky + 76q + 171H2O 10.10 619.8 0.3678 0.5654 0.4345 0.4095 0.2355 0.4229 0.7847 0.3816 0.8566 53chl + 10g + 110mu = 110bi + 139ky + 86q + 213H2O 12.10 633.2 0.3004 0.5503 0.4497 0.3427 0.1886 0.4435 0.7227 0.3165 0.8181 53chl + 10g + 109mu = 109bi + 137ky + 86q + 210H2O 14.10 644.7 0.2419 0.5395 0.4605 0.2847 0.1554 0.4280 0.6524 0.2581 0.7763 49chl + 10g + 104mu = 104bi + 128ky + 86q + 198H2O 15.00 649.4 0.2178 0.5356 0.4643 0.2599 0.1431 0.4113 0.6176 0.2336 0.7566 49chl + 10g + 103mu = 103bi + 125ky + 87q + 194H2O <==================================================> phases : chl, st, g, ky, (mu, q, fluid) or [bi] P(kbar) T(¡C) x(chl) y(chl) Q(chl) x(st) x(g) x(mu) y(mu) 0.10 504.8 0.1873 0.5310 0.4690 0.7422 0.6666 0.1879 0.9393 10chl + 78st + 340q = 121g + 594ky + 196H2O 2.10 652.7 0.1918 0.5697 0.4303 0.6559 0.5772 0.2235 0.9175 10chl + 59st + 263q = 96g + 451ky + 159H2O 4.10 668.6 0.2112 0.5712 0.4287 0.6662 0.5983 0.2451 0.9005 10chl + 69st + 304q = 108g + 526ky + 178H2O 6.10 669.3 0.2325 0.5681 0.4318 0.6846 0.6278 0.2653 0.8823 10chl + 85st + 367q = 129g + 643ky + 209H2O 8.10 661.6 0.2558 0.5631 0.4368 0.7069 0.6610 0.2851 0.8629 10chl + 108st + 464q = 160g + 821ky + 255H2O 10.10 647.7 0.2824 0.5573 0.4427 0.7325 0.6977 0.3057 0.8423 10chl + 146st + 622q = 211g + 1112ky + 331H2O 12.10 628.4 0.3140 0.5512 0.4487 0.7614 0.7375 0.3285 0.8206 10chl + 217st + 920q = 306g + 1660ky + 474H2O 14.10 604.2 0.3529 0.5452 0.4547 0.7937 0.7803 0.3552 0.7981 10chl + 391st + 1645q = 538g + 2995ky + 823H2O 15.00 591.6 0.3738 0.5427 0.4573 0.8094 0.8004 0.3692 0.7879 10chl + 582st + 2442q = 793g + 4460ky + 1205H2O <==================================================> phases : bi, st, g, ky, (mu, q, fluid) or [chl] P(kbar) T(¡C) x(bi) y(bi) Q(bi) x(st) x(g) x(mu) y(mu) 11 4.10 687.5 0.9300 0.7012 0.01886 0.9739 0.9726 0.8836 0.9474 10bi + 169st + 724q = 233g + 10mu + 1285ky + 338H2O 6.10 682.6 0.6793 0.5492 0.1373 0.8916 0.8825 0.6142 0.9193 10bi + 132st + 568q = 184g + 10mu + 1000ky + 264H2O 8.10 669.8 0.5107 0.3921 0.2805 0.8349 0.8190 0.4789 0.8898 10bi + 121st + 520q = 169g + 10mu + 913ky + 242H2O 10.10 651.2 0.4089 0.2702 0.3880 0.7903 0.7689 0.3874 0.8565 10bi + 109st + 468q = 153g + 10mu + 818ky + 217H2O 12.10 627.9 0.3435 0.1838 0.4492 0.7521 0.7261 0.3160 0.8179 10bi + 97st + 417q = 137g + 10mu + 725ky + 193H2O 14.10 600.6 0.2956 0.1220 0.4752 0.7179 0.6879 0.2562 0.7728 10bi + 86st + 373q = 124g + 10mu + 644ky + 172H2O 15.00 587.1 0.2772 0.1002 0.4781 0.7034 0.6717 0.2321 0.7499 10bi + 82st + 355q = 118g + 10mu + 611k