% --------------------------------------------------------
  Samecoding

  The samecoding script lines below should be copied into
  a THERMOCALC script file to tell it to make additional 
  copies of an x-eos with their own names and composition/order
  starting guesses. For example, the script
       samecoding cpx pig
  says "make a copy of cpx, called pig" (multiple copies 
  can be made).
  To use the samecoding script, there must be a xyzguess 
  starting guess block, as shown below, for the duplicate phase.
% --------------------------------------------------------

samecoding cpx pig    % clinopyroxene, pigeonite
samecoding spn cm mt  % Al-spinel, Cr-spinel, magnetite




% --------------------------------------------------------
  Starting guesses for composition/order variables

  Some possible starting guesses are given here that MIGHT help 
  (but make sure they represent a plausible composition for 
  the phases in your calcs). 
% --------------------------------------------------------

% silicate melt, hydrous (NCKFMASHTOCr)
xyzguess wo(L)         0.107216
xyzguess sl(L)        0.0873506
xyzguess fo(L)        0.0142023
xyzguess fa(L)        0.0109878
xyzguess jd(L)         0.151907
xyzguess hm(L)       0.00154624
xyzguess ek(L)       9.81910e-6
xyzguess ti(L)       0.00413534
xyzguess kj(L)        0.0368694
xyzguess yct(L)       0.0885231
xyzguess h2o(L)        0.219727

% silicate melt, dry (NCKFMASTOCr)
xyzguess wo(L)         0.228033
xyzguess sl(L)         0.123908
xyzguess fo(L)         0.383275
xyzguess fa(L)        0.0784050
xyzguess jd(L)        0.0408634
xyzguess hm(L)        0.0156940
xyzguess ek(L)        0.0178215
xyzguess ti(L)       0.00817432
xyzguess kj(L)       0.00163650
xyzguess yct(L)     0.000695507

% aqueous fluid with dissolved silicates
xyzguess wo(fl)     0.000103865
xyzguess sl(fl)      5.87807e-5
xyzguess fo(fl)      1.41812e-8
xyzguess fa(fl)      3.11791e-8
xyzguess jd(fl)     0.000862050
xyzguess hm(fl)     0.000159068
xyzguess ek(fl)      3.29610e-7
xyzguess ti(fl)     0.000162396
xyzguess kj(fl)      0.00562408
xyzguess h2o(fl)       0.985936

% plagioclase
xyzguess ca(pl4tr)        0.8
xyzguess k(pl4tr)        0.0116338


% olivine
xyzguess x(ol)         0.102434
xyzguess c(ol)       0.00275804
xyzguess Q(ol)      0.000156353

% orthopyroxene
xyzguess x(opx)       0.0956777
xyzguess y(opx)        0.182193  range  0.000 2.000
xyzguess c(opx)       0.0489271
xyzguess Q(opx)      -0.0475598  range -1.000 1.000
xyzguess f(opx)       0.0177886
xyzguess t(opx)      0.00823030
xyzguess cr(opx)      0.0205952
xyzguess j(opx)      0.00906326

% clinopyroxene
xyzguess x(cpx)       0.0964208
xyzguess y(cpx)       0.0825558  range  0.000 2.000
xyzguess o(cpx)        0.219892
xyzguess n(cpx)        0.118338
xyzguess Q(cpx)      -0.0282419  range -1.000 1.000
xyzguess f(cpx)       0.0209603
xyzguess cr(cpx)      0.0154784
xyzguess t(cpx)       0.0127859
xyzguess k(cpx)      0.00527429

% pigeonite
xyzguess x(pig)        0.124
xyzguess y(pig)        0.112  range  0.000 2.000
xyzguess o(pig)        0.888
xyzguess n(pig)        0.028
xyzguess Q(pig)       -0.0115  range -1.000 1.000
xyzguess f(pig)        0.004
xyzguess cr(pig)       0.001
xyzguess t(pig)        0.148
xyzguess k(pig)        0.001

% spinel
xyzguess x(spn)        0.144990
xyzguess y(spn)       0.0139965
xyzguess c(spn)       0.0451070
xyzguess t(spn)       0.0110937
xyzguess Q1(spn)       0.503697
xyzguess Q2(spn)      0.0807768
xyzguess Q3(spn)      0.0218537

% K-feldspar
xyzguess na(k4tr)       0.1
xyzguess ca(k4tr)       0.004

% muscovite
xyzguess x(mu)          0.275719
xyzguess y(mu)          0.334730
xyzguess f(mu)          0.0164856
xyzguess n(mu)          0.00226622
xyzguess c(mu)          3.36760e-5

% biotite
xyzguess x(bi)         0.457528
xyzguess y(bi)        0.0772161
xyzguess f(bi)        0.0158764
xyzguess t(bi)         0.119336
xyzguess Q(bi)         0.228603

% garnet
xyzguess x(g)          0.455851
xyzguess c(g)          0.286243
xyzguess f(g)        0.00253871
xyzguess cr(g)      0.000841005
xyzguess t(g)         0.0113394

% epidote
xyzguess f(ep)          0.189427
xyzguess Q(ep)          0.182107  range  0.000 0.500

% cordierite
xyzguess x(cd)            0.3
xyzguess h(cd)            0.7

% Al-spinel
xyzguess x(spn)            0.14499
xyzguess y(spn)            0.01398
xyzguess c(spn)            0.04504
xyzguess t(spn)            0.01114
xyzguess Q1(spn)           0.50372 range -1 1
xyzguess Q2(spn)           0.08077 range -1 1
xyzguess Q3(spn)           0.02183 range -1 1

% Cr-spinel
xyzguess x(cm)            0.2
xyzguess y(cm)            0.1
xyzguess c(cm)            0.8
xyzguess t(cm)            0.05
xyzguess Q1(cm)           0.05 range -1 1
xyzguess Q2(cm)           0.05 range -1 1
xyzguess Q3(cm)           0.05 range -1 1

% magnetite
xyzguess x(mt)            0.8
xyzguess y(mt)            0.55
xyzguess c(mt)            0.2
xyzguess t(mt)            0.01
xyzguess Q1(mt)           0.33 range -1 1
xyzguess Q2(mt)           0.14 range -1 1
xyzguess Q3(mt)           0.25 range -1 1

% hornblende
xyzguess x(hb)         0.355436
xyzguess y(hb)         0.393748
xyzguess z(hb)        0.0255955
xyzguess a(hb)         0.476446
xyzguess k(hb)        0.0647451
xyzguess c(hb)         0.924162
xyzguess f(hb)         0.103286
xyzguess t(hb)        0.0360617
xyzguess Q1(hb)      -0.0186478  range -1.000 1.000
xyzguess Q2(hb)        0.111970  range -1.000 1.000

% ilmenite
xyzguess x(ilm)        0.913364
xyzguess Q(ilm)        0.627896  range -0.990 0.990


% --------------------------------------------------------