% ==================================================================
% Scriptfile for .... calcs 
%
% in the prefsfile, these scripts are needed
%           dataset   62          
%           calcmode 1 
%           scriptfile gridproject  % if this file keeps its name
% 
% Use "%" for commenting (nothing is read between % and the line ending).
% No lines are read after the "*" character.
% ==================================================================

axfile mp50MnNCKFMASHTO  % axfile, tc-<system>.txt, becomes axfile <system>

pseudosection no         % this is not a pseudosection calc (so no bulk comp is needed)

autoexit

% =======================
% phases
% =======================

with g chl st bi ky

% or:
% with someof g chl st bi ky ctd
                        % list of phases in calculation (if not enough for 
                        % univariant or divariant equilibrium, you will be asked to 
                        % specify compositional isopleths)

inexcess  mu q H2O      % list of phases to be included in all calculations

% =======================
% conditions
% =======================

diagramPT  4 14  400 700     % P in kbar; T in degC 

calcP  6 10 1                % Plow Phigh Pinc
calcT  500 700 50            % Tlow Thigh Tinc

calctatp                     % do calculate T at P (so, info in calcT ignored)

% =======================
% what output?
% =======================

drawpd yes              % print output for plotting with drawpd software

% =======================
% set up x-eos
% =======================

% samecoding mu pa      % duplicating x-eos: make paragonite from muscovite

% !!!! 
% Add your own xyzguess blocks for all phases (including any made by samecoding)
% !!!! 

% ---------------------------------------------------------- #16
xyzguess x(g)          0.887427
xyzguess z(g)                 0  % boiler
xyzguess m(g)                 0  % boiler
xyzguess f(g)                 0  % boiler
% -----------------------------
xyzguess x(bi)         0.674182
xyzguess m(bi)                0  % boiler
xyzguess y(bi)         0.260333
xyzguess f(bi)                0  % boiler
xyzguess t(bi)                0  % boiler
xyzguess Q(bi)        0.0969282
% -----------------------------
xyzguess x(mu)         0.548754
xyzguess y(mu)         0.909621
xyzguess f(mu)                0  % boiler
xyzguess n(mu)                0  % boiler
xyzguess c(mu)                0  % boiler
% -----------------------------
xyzguess x(st)         0.858563
xyzguess m(st)                0  % boiler
xyzguess f(st)                0  % boiler
xyzguess t(st)                0  % boiler
% ---------------------------------------------------------- #21
xyzguess x(chl)        0.465451
xyzguess y(chl)        0.805095
xyzguess f(chl)               0  % boiler
xyzguess m(chl)               0  % boiler
xyzguess QAl(chl)      0.194402  range -1.000 1.000
xyzguess Q1(chl)      0.0352754  range -1.000 1.000
xyzguess Q4(chl)      0.0888073  range -1.000 1.000
% -----------------------------
xyzguess x(ctd)        0.772861
xyzguess m(ctd)               0  % boiler
xyzguess f(ctd)               0  % boiler
% -----------------------------

* 

with the univariant assemblage, THERMOCALC produces:

<==========================================================>
phases: g bi st chl (mu, q, fluid) 

------------------------------------------------------------
 P(kbar)     T(°C)      x(g)     x(bi)     y(bi)     Q(bi)     x(mu)     y(mu) ...
   6.000    549.65    0.9596    0.8719    0.2224   0.04184    0.7447    0.9289 ...
    29g + 48mu + 18chl = 48bi + 10st + 72q + 53H2O
   7.000    571.17    0.9109    0.7241    0.2277   0.08844    0.5813    0.9109 ...
    27g + 50mu + 20chl = 50bi + 10st + 75q + 61H2O
   8.000    589.88    0.8641    0.6008    0.2278    0.1194    0.4842    0.8947 ...
    25g + 52mu + 22chl = 52bi + 10st + 78q + 69H2O
   9.000    606.20    0.8178    0.5010    0.2236    0.1341    0.4155    0.8793 ...
    23g + 54mu + 24chl = 54bi + 10st + 79q + 77H2O
  10.000    620.52    0.7713    0.4203    0.2159    0.1362    0.3619    0.8640 ...
    21g + 56mu + 26chl = 56bi + 10st + 81q + 85H2O


