@@ LA-SR-MN-CR-OXIDE, Povoden-Karadeniz, version 1.0, 2012-04-12
@@ In the following, @@ denotes a commentary line, all others are command lines
@@
@@ This oxide database is designed for the use in LSCM for solid oxide fuel cells. 
@@ 
@@ This ThermoCalc macrofile (.tcm) can be easily transformed into .tdb format 
@@ (thermodynamic databases) by using Thermocalc --> Gibbs module (Command: go Gibbs)
@@ --> command list-data --> file --> give a name, e.g. LSCM_v1.0 --> N
@@
 SET-LOG LSCM_v1.0.log,,
@@ sets a log file of the setup file
 GO G
@@ go to Gibbs module of ThermoCalc
 ENTER-ELEMENT LA SR MN CR O VA
@@ definition of elements
@@ELEMENT NAME    REF. STATE             ATOMIC MASS  H0          S0
 AMEND-ELEMENT LA DOUBLE_HCP(ABAC)       1.3891E+02   6.6651E+03  5.6902E+01,,
 AMEND-ELEMENT SR SR_FCC_A1              8.7620E+01   6.5680E+03  5.5694E+01,,
 AMEND-ELEMENT MN CBCC_A12               5.4938E+01   4.9960E+03  3.2008E+01,,
 AMEND-ELEMENT CR BCC_A2                 5.1996E+01   4.0500E+03  2.35429E+01,,
 AMEND-ELEMENT O 1/2_MOLE_O2(G)          1.5999E+01   4.3410E+03  1.0252E+02,,
 AMEND-ELEMENT VA VACUUM                 0            0           0,,
@@ Species
 ENTER-SPECIES  LA+2                               LA/+2
 ENTER-SPECIES  LA+3                               LA/+3
 ENTER-SPECIES  SR2                                SR2
 ENTER-SPECIES  SR+2                               SR/+2
 ENTER-SPECIES  MN+2                               MN/+2
 ENTER-SPECIES  MN+3                               MN/+3
 ENTER-SPECIES  MN+4                               MN/+4
 ENTER-SPECIES  CR2                                CR2
 ENTER-SPECIES  CR+2                               CR/+2
 ENTER-SPECIES  CR+3                               CR/+3
 ENTER-SPECIES  CR+4                               CR/+4
 ENTER-SPECIES  CR+5                               CR/+5
 ENTER-SPECIES  CR+6                               CR/+6
 ENTER-SPECIES  O1                                 O1
 ENTER-SPECIES  O2                                 O2
 ENTER-SPECIES  O3                                 O3
 ENTER-SPECIES  O-2                                O/-2
 ENTER-SPECIES  O2-2                               O2/-2
 ENTER-SPECIES  LA1O1                              LA1O1
 ENTER-SPECIES  LA2O1                              LA2O1
 ENTER-SPECIES  LA2O2                              LA2O2
 ENTER-SPECIES  SRO                                SRO
 ENTER-SPECIES  SRO2                               SRO2
 ENTER-SPECIES  SR2O                               SR2O
 ENTER-SPECIES  MN1O1                              MN1O1
 ENTER-SPECIES  MN1O2                              MN1O2
 ENTER-SPECIES  MN2O3                              MN2O3
 ENTER-SPECIES  MN2O7                              MN2O7
 ENTER-SPECIES  MN3O4                              MN3O4 
 ENTER-SPECIES  CR1O1                              CR1O1
 ENTER-SPECIES  CR1O2                              CR1O2
 ENTER-SPECIES  CR1O3                              CR1O3
 ENTER-SPECIES  CR2O2                              CR2O2
 ENTER-SPECIES  CR2O3                              CR2O3
 ENTER-SPECIES  CR3O4                              CR3O4
 ENTER-SPECIES  LAMNO3                             LAMNO3
 ENTER-SPECIES  LACRO3                             LACRO3
 ENTER-SPECIES  SRCRO3                             SRCRO3
 ENTER-SPECIES  SRCRO4                             SRCRO4
 ENTER-SPECIES  SRMNO3                             SRMNO3
@@ ----------------------------------------------------------
@@ Gibbs energy functions
@@ ----------------------------------------------------------
@@ SER Lattice stabilities, Dinsdale 1991 [9]
@@ La, double hcp
 ENTER-SYMBOL FUNCTION GHSERLA 298.15 
	-7968.403+120.284604*T-26.34*T*LN(T)-.001295165*T**2;
 550	Y
	-3381.413+59.06113*T-17.1659411*T*LN(T)-.008371705*T**2
	+6.8932E-07*T**3-399448*T**(-1); 
 2000	Y
	-15608.882+181.390071*T-34.3088*T*LN(T);  
 4000	N 
@@ Gibbs energy function name is GHSERLA, low temperature limit is 298.15 K. 
@@ First segment of Gibbs energy function goes up to 550 K, Y means that a further segment follows.
@@ Last segment of Gibbs energy ends at 4000 K, N stands for no further segments. 
@@ Sr, fcc
 ENTER-SYMBOL FUNCTION GHSERSR 298.15  
	-7532.367+107.183879*T-23.905*T*LN(T)-4.61225E-3*T**2
	-1.67477E-07*T**3-2055*T**(-1);
 820	Y  
	-13380.102+153.196104*T-30.0905432*T*LN(T)-3.251266E-3*T**2
	+1.84189E-07*T**3+850134*T**(-1);
 3000	N  
@@ Mn, cbcca12
 ENTER-SYMBOL FUNCTION GHSERMN 298.15
	-8115.27966+130.059572*T-23.4582*T*LN(T)-0.00734768*T**2
	+69827.1*T**(-1);
 1519	Y
	-28733.41+312.2648*T-48*T*LN(T)+1.656847E30*T**(-9);
 2000	N
@@ Cr, bcc
 ENTER-SYMBOL FUNCTION GHSERCR 298.15
	-8856.94+157.48*T-26.908*t*LN(T)+0.00189435*T**2
	-1.47721E-06*T**3+139250*T**(-1);
 2180	Y
	-34869.344+344.18*T-50.0*T*LN(T)-2.88526E+32*T**(-9);
 6000	N
@@ O1, (1/2 O2)
 ENTER-SYMBOL FUNCTION GHSEROO 298.15
	-3480.872255-25.5028601*T-11.1355068*T*LN(T)-0.005098873*T**2
	+6.6184604E-07*T**3-38364.8742*T**(-1);
 1000	Y  
	-6568.76015+12.66000166*T-16.8138015*T*LN(T)-5.9579637E-04*T**2
	+6.78055555E-09*T**3+262904.778*T**(-1);
 3300	Y  
	-13986.728+31.259625*T-18.9536*T*LN(T)-4.25243E-04*T**2
	+1.0721E-08*T**3+4383200*T**(-1);
 6000	N
@@-------------------------------------------------------------------
@@ Solid metals, Dinsdale 1991 [9]
@@ La, bcc
 ENTER-SYMBOL FUNCTION GLABCC 298.15
	-3952.161+88.072353*T-21.7919*T*LN(T)-0.004045175*T**2
	-5.25865E-07*T**3; 
 800	Y
	+321682.673-3565.08252*T+513.440708*T*LN(T)-0.387295093*T**2
	+4.9547989E-05*T**3-36581228*T**(-1); 
 1134	Y
	-16377.894+218.492988*T-39.5388*T*LN(T);
 1193	Y
	-136609.91+1123.34397*T-163.413074*T*LN(T)+0.053968535*T**2
	-4.056395E-06*T**3+21167204*T**(-1); 
 2000	Y
	-8205.988+174.836315*T-34.3088*T*LN(T); 
 4000	N
@@ La, fcc
 ENTER-SYMBOL FUNCTION GLAFCC 298.15
	-6109.797+89.878761*T-21.7919*T*LN(T)-0.004045175*T**2
	-5.25865E-07*T**3;
 1134	Y
	-124598.976+955.878375*T-139.346741*T*LN(T)+0.042032405*T**2
	-3.066199E-06*T**3+20994153*T**(-1); 
 2000	Y
	-12599.386+178.54399*T-34.3088*T*LN(T); 
 4000	N
@@ Sr, bcc
ENTER-SYMBOL FUNCTION GSRBCC 298.15  
	-6779.234+116.583654*T-25.6708365*T*LN(T)-3.126762E-3*T**2
	+2.2965E-7*T**3+27649*T**(-1);
 820	Y  
	-6970.594+122.067301*T-26.57*T*LN(T)-1.9493E-3*T**2
	-1.7895E-8*T**3+16495*T**(-1);
 1050	Y  
	+8168.357+0.423037*T-9.7788593*T*LN(T)-9.539908E-3*T**2
	+5.20221E-7*T**3-2414794*T**(-1);  
 3000	N
@@ Mn, bcc
 ENTER-SYMBOL FUNCTION GMNBCC 298.15
	-3235.3+127.85*T-23.7*T*LN(T)-0.00744271*T**2
	+60000*T**(-1);
 1519	Y
	-23188.83+307.7043*T-48*T*LN(T)+1.265152E30*T**(-9);
 3000	N
@@ Mn, fcc
 ENTER-SYMBOL FUNCTION GMNFCC 298.15
	-3439.3+131.884*T-24.5177*T*LN(T)-0.006*T**2
	+69600*T**(-1);
 1519	Y
	-26070.1+309.6664*T-48*T*LN(T)-3.86196E30*T**(-9);
 3000	N
@@ Mn, cuba13
 ENTER-SYMBOL FUNCTION GMNCUB 298.15  
	-5800.4+135.995*T-24.8785*T*LN(T)-0.00583359*T**2
	+70269*T**(-1);
 1519	Y 
	-28290.76+311.2933*T-48*T*LN(T)+3.96757E30*T**(-9);
 3000	N
@@ Cr, fcc based on GHSERCR
 ENTER-SYMBOL FUNCTION GCRFCC 298.15
	+GHSERCR+7284+0.163*T;
 1519	N
@@ Cr, hcpa3
 ENTER-SYMBOL FUNCTION GCRHCP 298.15
	-4418.94+157.48*T-26.908*T*LN(T)+0.00189435*T**2
	-1.47721E-06*T**3+139250*T**(-1);
 2180	Y
	-30431.344+344.18*T-50*T*LN(T)-2.88526E+32*T**(-9);
 6000	N 
@@ Cr, cuba13
 ENTER-SYMBOL FUNCTION GCRCUB  298.15
	+7042.06+158.1076*T-26.908*T*LN(T)+0.00189435*T**2
	-1.47721E-06*T**3+139250*T**(-1);
 2180	Y
	-18970.344+344.8076*T-50*T*LN(T)-2.88526E+32*T**(-9);
 6000	N 
@@ Cr, bcca12
 ENTER-SYMBOL FUNCTION GCRBCCA12 298.15
	+2230.06+160.1996*T-26.908*T*LN(T)+0.00189435*T**2
	-1.47721E-06*T**3+139250*T**(-1);
 2180	Y
	-23782.344+346.8996*T-50*T*LN(T)-2.88526E+32*T**(-9);
 6000	N 
@@------------------------------------------------------------
@@ Binary oxides, optimized
@@ La-oxides, Grundy et al. 2004 [17]
 ENTER-SYMBOL FUNCTION GLAO 298.15
	+0.3333333*GLA2O3D+0.3333333*GHSERLA+62000;,,,
@@ Zinkevich et al. 2006 [18]
 ENTER-SYMBOL FUNCTION GLA2O3D 298.15  
	-1833257+692.9664*T-120.629*T*LN(T)-0.006854*T**2
      +808000*T**(-1)-1E7*T**(-2);,,,
 ENTER-SYMBOL FUNCTION GLA2O3H 298.15  
	32350-13.986*T+GLA2O3D;,,,
 ENTER-SYMBOL FUNCTION GLA2O3X 298.15  
	43192-18.555*T+GLA2O3D;,,,
@@
@@ SrO, Risold et al. 1996 [50]
 ENTER-SYMBOL FUNCTION GSROSOL 298.15  
	-607870+268.9*T-47.56*T*LN(T)-0.00307*T**2
      +190000*T**(-1);,,,
@@ SrO2
 ENTER-SYMBOL FUNCTION GSRO2SOL 298.15  
	+GSROSOL+GHSEROO-43740+70*T;,,,
@@ MANGANOSITE, MNO, Grundy et al. 2005 [21]
 ENTER-SYMBOL FUNCTION GMN1O1 298.15
	-4.02477557E+05+2.59355626E+02*T-4.68352649E+01*T*LN(T)
	-3.85001409E-03*T**2+2.12922234E+05*T**(-1);,,,
@@ PYROLYSITE, MN1O2, Grundy et al. 2005 [22]
 ENTER-SYMBOL FUNCTION GMN1O2 298.15
	-5.45091278E+05+3.95379396E+02*T-6.52766201E+01*T*LN(T)
	-7.80284521E-03*T**2+6.64955386E+05*T**(-1);,,,
@@
@@ MN2O3-FUNCTION, MODIFIED, Grundy et al. 2005 [22]
 ENTER-SYMBOL FUNCTION GMN2O3 298.15
	-9.96393E+05+5.6846E+02*T-9.911E+01*T*LN(T)-2.056E-02*T**2
	+6.0822E+05*T**(-1);,,,
@@ ALPHA-HAUSMANNITE, ALPHA-MN3O4 (DISTORTED), Grundy et al. 2005 [22]
 ENTER-SYMBOL FUNCTION GTMN3O4 298.15
	-1.43703676E+06+8.89567858E+02*T-1.54747566E+02*T*LN(T)
	-1.74079033E-02*T**2+9.86138663E+05*T**(-1);,,,
@@ BETA-HAUSMANNITE, BETA-MN3O4 (CUBIC), Grundy et al. 2005 [22]
 ENTER-SYMBOL FUNCTION GCMN3O4 298.15
	-1.41618912E+06+8.75120338E+02*T-1.54747566E+02*T*LN(T)
	-1.74079033E-02*T**2+9.86138663E+05*T**(-1);,,,
@@
@@ METASTABLE CRO Povoden et al. 2006 [19]
 ENTER-SYMBOL FUNCTION GCR1O1 298.15
        +0.5*GCR2O3-0.5*GHSEROO+255269-53.82*T;,,,
@@ ESKOLAITE, Cr2O3, Povoden et al. 2006 [19]
 ENTER-SYMBOL FUNCTION GCR2O3 298.15 
	-1.164542E+06+7.2856E+02*T-119.8*T*LN(T)-4.97E-03*T**2
	+1.05E+06*T**(-1);,,,
@@ reduced neutral endmember of CR2O3, Povoden et al. 2006 [19]
 ENTER-SYMBOL FUNCTION GCRO0 298.15 
 	+108305+GCR2O3+0.66666666667*GHSERCR;,,,
@@ CR-SPINEL CR3O4, Povoden et al. 2006 [19]
 ENTER-SYMBOL FUNCTION GCR3O4 298.15
	+1.5*GCR2O3-0.5*GHSEROO+280045-93.76*T;,,,
@@ --------------------------------------------------------------
@@ Ternary oxides, optimized (except perovskite functions)
@@
@@ Functions of the La-Sr-O system, Grundy et al. 2005 [25]
 ENTER-SYMBOL FUNCTION SR_ALPHA 298.15
	+2*GSROSOL+2.5E+04;,,,
 ENTER-SYMBOL FUNCTION SRH_ALPH 298.15
	+2*GSROSOL+2.5E+04;,,,
 ENTER-SYMBOL FUNCTION SRX_ALPH 298.15
	+2*GSROSOL+2.5E+04;,,,
 ENTER-SYMBOL FUNCTION LA_BETA 298.15
	+GLA2O3D+2.158E+04;,,,
 ENTER-SYMBOL FUNCTION RE_ALPHA 298.15
	0;,,,
 ENTER-SYMBOL FUNCTION RE_BETA 298.15
	0;,,,
@@
@@ Functions of the La-Mn-O system, Grundy et al. 2005 [23]
 ENTER-SYMBOL FUNCTION GL2MNO4 298.15 
	+GLA2O3D+GMN1O1+6.26029731E+04-3.71704891E+01*T;,,,
@@
@@ Functions of the La-Cr-O system, Povoden et al. 2009 [20]
@@ LA2CRO6 ([], slightly modified) 
 ENTER-SYMBOL FUNCTION GLA2CRO6 298.15 
 GLA2O3D+0.5*GCR2O3+1.5*GHSEROO-73045-4.14*T;,,,
@@ LA2CR3O12 (GLA2CR3O12 in ref [])
 ENTER-SYMBOL FUNCTION GLA2CR3 298.15  
 GLA2O3D+1.5*GCR2O3+4.5*GHSEROO-371557+205*T;,,,
@@
@@ Functions of the Sr-Mn-O system, Grundy et al. 2005 [26]
@@ HEX Phase
 ENTER-SYMBOL FUNCTION GSM4_HEX 298.15 
	+GSROSOL+GMN1O2-1.11300000E+05;,,,
 ENTER-SYMBOL FUNCTION GSM3_HEX 298.15 
	+GSROSOL+0.5*GMN2O3-7.73000000E+03
	-1.70000000E+01*T;,,,
@@ SrMn3Oz as SrMnO3_Mn2O3
 ENTER-SYMBOL FUNCTION GSM4OZ 298.15 
	+GMN2O3+GSM4_HEX-8.79100000E+03;,,,
 ENTER-SYMBOL FUNCTION GSM3OZ 298.15 
	+GMN2O3+GSM3_HEX-2.19200000E+04;,,,
@@ RP1 (Ruddlesden-Popper Phase 1)
 ENTER-SYMBOL FUNCTION GS4O_RP1 298.15 
	+2*GSROSOL+GMN1O2-1.32830000E+05;,,,
 ENTER-SYMBOL FUNCTION GL3O_RP1 298.15 
	+GSROSOL+0.5*GLA2O3D+0.5*GMN2O3-68300;,,,
@@ RP2
 ENTER-SYMBOL FUNCTION GS4O_RP2 298.15 
	+3*GSROSOL+2*GMN1O2-8.99100000E+04-90*T;,,,
 ENTER-SYMBOL FUNCTION GL3O_RP2 298.15 
	+GSROSOL+GLA2O3D+GMN2O3-137400;,,,
@@ RP3 
 ENTER-SYMBOL FUNCTION GSM4_RP3 298.15 
	+4*GSROSOL+3*GMN1O2-3.78500000E+05;,,,
@@ Sr7Mn4O1
 ENTER-SYMBOL FUNCTION GS7M4 298.15 
	+7*GSROSOL+4*GMN1O2-6.12450000E+05+50*T;,,,
@@
@@ Functions of the Sr-Cr-O system, This work
@@ SRCR2O4
 ENTER-SYMBOL FUNCTION GSC2O4 
    298.15 +GSROSOL+GCR2O3+20000-28*T;,,,
@@ SR2CRO4
 ENTER-SYMBOL FUNCTION GS2CO4 
    298.15 2*GSROSOL+0.5*GCR2O3+0.5*GHSEROO-145000+50*T;,,,
@@ SR3CR2O8
 ENTER-SYMBOL FUNCTION GS3C2O8 298.15  
	+2.66666667*GSROSOL+GCR2O3+2.333333333333*GHSEROO
	-508507+219*T;,,,
@@ SRCRO4
 ENTER-SYMBOL FUNCTION GSCO4 298.15 
	+GSROSOL+0.5*GCR2O3+1.5*GHSEROO-273771
      +131.6*T;,,,
@@
@@ Functions of the Mn-Cr-O system, Povoden et al. 2006 [21] 
@@ CUBIC SPINEL
 ENTER-SYMBOL FUNCTION GSPINEL 298.15 
      0.666666667*GCR3O4+.33333333334*GCMN3O4-210795.3+61.69*T;,,,
@@ TETRAGONALLY DISTORTED SPINEL
 ENTER-SYMBOL FUNCTION GTSPINEL 298.15 
      0.666666667*GCR3O4+.33333333334*GTMN3O4-200942+75.1*T;,,,
@@ ----------------------------------------------------------------
@@ Perovskite functions, Grundy et al. 2005 [24]
 ENTER-SYMBOL FUNCTION GL3O 298.15 
	+0.5*GLA2O3D+0.5*GMN2O3-63367+51.77*T-7.19*T*LN(T)
	+232934*T**(-1);,,,
@@ GL3OL slightly modified from ref [24]
 ENTER-SYMBOL FUNCTION GL3OL 298.15 
	+0.5*GLA2O3D+0.5*GMN2O3-63367+51.77*T-7.19*T*LN(T)
	+232934*T**(-1)-3429+4.72*t;,,,
@@ GL3OR slightly modified from ref [24]
 ENTER-SYMBOL FUNCTION GL3OR 298.15 
	+0.5*GLA2O3D+0.5*GMN2O3-63367+51.77*T-7.19*T*LN(T)
	+232934*T**(-1)+400-0.4*t;,,,
 ENTER-SYMBOL FUNCTION GL2O 298.15 
	+0.5*GLA2O3D+GMN1O1+27672;,,,
 ENTER-SYMBOL FUNCTION GL4O 298.15 
	+0.5*GLA2O3D+0.75*GMN1O2-91857+20.31*T;,,,
 ENTER-SYMBOL FUNCTION GV4O 298.15 
	+0.333333*GLA2O3D+GMN1O2-53760;,,,
 ENTER-SYMBOL FUNCTION GVVV 298.15 
	+6*GL2O+4*GL4O+3*GV4O-12*GL3O-254212;,,,
 ENTER-SYMBOL  FUNCTION GMS3O 298.15 
	+GSROSOL+0.5*GMN2O3-7.73000000E+03-1.44550000E+04
	-1.70000000E+01*T;,,,
 ENTER-SYMBOL  FUNCTION GMS4O 298.15 
	+GSROSOL+GMN1O2-1.11300000E+05+2.26500000E+04
	-7.69000000E+00*T;,,,
 ENTER-SYMBOL FUNCTION ANTI 298.15 
	+547422;,,,
@@ Reciprocals: all 0!
@@
@@ Povoden et al. 2009 [20]
@@ LaCrO3-PEROVSKITE
 ENTER-SYMBOL FUNCTION GLACRO3 298.15 
 	0.5*GLA2O3D+0.5*GCR2O3-73591+2.38*T-0.68*T*LN(T);,,,
 ENTER-SYMBOL FUNCTION GALACRO3 298.15 
 	+GLACRO3-340+0.63*t;,,,
@@ (LaSr)CrO3+/-delta-Perovskite
@@ Reference SrCrVa3
 ENTER-SYMBOL FUNCTION GS4V 298.15
	GSROSOL+0.5*GCR2O3-2.5*GHSEROO;,,,
@@ FUNCTION FOR NEUTRAL SRCRO3 endmember, modified from ref [20]
 ENTER-SYMBOL FUNCTION GS4O 298.15
	+GSROSOL+0.5*GCR2O3+0.5*GHSEROO+30800-40*T;,,,
@@ FUNCTION FOR NEUTRAL ENDMEMBER SR(CR+3,VA)O3 (modified GS3O or GS3V in ref [20])
 ENTER-SYMBOL FUNCTION GS3V 298.15
	+GSROSOL+0.5*GCR2O3+11.2386*T+140300-67*T;,,,
 ent-sym fun GVCR4O 298.15 
      +0.5*GCR2O3+0.5*GHSEROO-279000-265*T;,,,
 ent-sym fun GLACR4O 298.15   
      +0.33333*GLA2O3D+.5*GCR2O3+0.5*GHSEROO-200000;,,,
@@ ----------------------------------------------------------------------
@@ Liquid metal functions, Dinsdale 1991 [9]
@@ La
 ENTER-SYMBOL FUNCTION GLALIQ  298.15
	+5332.653+18.23012*T-11.0188191*T*LN(T)-0.020171603*T**2
	+2.93775E-06*T**3-133541*T**(-1);       
 1134	Y
	-3942.004+171.018431*T-34.3088*T*LN(T); 
 4000	N
@@ Sr
ENTER-SYMBOL FUNCTION GSRLIQ 298.15  
	+2194.997-10.118994*T-5.0668978*T*LN(T)-3.1840595E-2*T**2
	+4.981237E-06*T**3-265559*T**(-1);
 1050	Y  
	-10855.29+213.406219*T-39.463*T*LN(T);  
 3000	N
@@ Mn
 ENTER-SYMBOL FUNCTION GMN_L 298.15
	+GHSERMN+17859.91-12.6208*T-4.41929E-21*T**7;
 1519	Y
	+GHSERMN+18739.51-13.2288*T-1.656847E30*T**(-9);
 3000	N
@@ Cr
 ENTER-SYMBOL FUNCTION GCR_L 298.15
	+15483.015+146.059775*T-26.908*T*LN(T)+.00189435*T**2
	-1.47721E-06*T**3+139250*T**(-1)+2.37615E-21*T**7;   
 2180	Y
	-16459.984+335.616317*T-50*T*LN(T);   
 6000	N
@@ ---------------------------------------------------------------------
@@ Liquid oxide functions, optimized
@@ liquid La2O3, Zinkevich et al. 2006 [18]
 ENTER-SYMBOL FUNCTION GLA2O3LIQ  298.15  
     -1833257+692.9664*T-120.629*T*LN(T)-0.006854*T**2
	+808000*T**(-1)-1E7*T**(-2)+141329-56.6220*T;,,,
@@ liquid SrO, Risold et al. 1996 [50]
 ENTER-SYMBOL FUNCTION GSROLIQ  
    298.15  -566346+449*T-73.1*T*LN(T);,,,
@@ liquid Mn oxides, Grundy et al. 2004 [22] 
 ENTER-SYMBOL FUNCTION GMN1O1_L 298.15
   GMN1O1+4.39465890E+04-2.06284295E+01*T;,,,
 ENTER-SYMBOL FUNCTION GMN2O3_L 298.15
          +2*GMN1O1+GHSEROO-6.49525609E+04+4.31437957E+01*T;,,,
@@ liquid Cr oxides, Povoden et al. 2006 [19]
 ENTER-SYMBOL FUNCTION GCR1O1_L  298.15
   0.5*GCR2O3-0.5*GHSEROO+339673-121.4*T;,,,
 ENTER-SYMBOL FUNCTION GCR2O3_L  298.15
  GCR2O3+439078-169*T;,,,
@@------------------------------------------------------------
@@ Gas functions from SGTE 1996 [29]
@@ La gas
@@ La(g)
 ENTER-SYMBOL FUNCTION F12026T  298.15 
	+422273.955-30.3347881*T-22.06299*T*LN(T)-0.005444405*T**2
	+4.71447833E-07*T**3+102710.1*T**(-1); 
 600	Y
	+426628.905-85.4786162*T-13.83676*T*LN(T)-0.011938995*T**2
	+1.33826017E-06*T**3-312130.2*T**(-1); 
 1300	Y
	+404460.17+114.016725*T-42.00406*T*LN(T)+0.0037094435*T**2
	-2.70261E-07*T**3+2891891*T**(-1); 
 3200	Y
	+497751.747-246.085237*T+2.791973*T*LN(T)-0.006002155*T**2
	+1.30043383E-07*T**3-34158815*T**(-1);
 8200	Y
	-92343.0441+773.338363*T-111.0188*T*LN(T)+0.0037862445*T**2
	-2.82257667E-08*T**3+5.418475E+08*T**(-1); 
 10000	N
@@ Sr gas
@@ Sr(g)
 ENTER-SYMBOL FUNCTION F15323T 298.15  
	+154227.522-24.1431703*T-20.98549*T*LN(T)+1.951298E-4*T**2
	-3.09095833E-8*T**3+4675.2365*T**(-1);
 1800	Y  
	+111247.483+242.365806*T-56.52776*T*LN(T)+0.0133862*T**2
	-9.57800833E-7*T**3+9843260*T**(-1);
 3300	Y  
	+770872.513-2114.76782*T+233.253*T*LN(T)-0.04337796*T**2
	+1.134592E-6*T**3-2.7250735E+8*T**(-1);
 4900	Y  
	-196742.694+263.327068*T-44.45892*T*LN(T)-0.008078665*T**2
	+2.96671167E-7*T**3+3.57637E+8*T**(-1);
 6200	Y  
	-949056.902+1952.13337*T-239.3059*T*LN(T)+0.01421437*T**2
	-1.79062E-7*T**3+8.9842E+8*T**(-1);
 9600	Y  
	+34305.7758+474.957384*T-77.25547*T*LN(T)+0.00232914*T**2
	-1.54504333E-8*T**3-2.2245325E+8*T**(-1);
 10000	N
@@ Sr2(g)
 ENTER-SYMBOL FUNCTION F15338T 298.15  
	+295010.66+61.845039*T-54.13634*T*LN(T)+0.040485225*T**2
	-9.264165E-6*T**3-70453.75*T**(-1);
 500	Y  
	+307156.188-147.411671*T-20.95926*T*LN(T)+1.012636E-4*T**2
	-8.03856667E-9*T**3-905190.5*T**(-1);
 3000	N
@@ Cr gas
@@ Cr (g)
 ENTER-SYMBOL FUNCTION F7491T 298.15 
	+390765.331-31.5192158*T-21.36083*T*LN(T)+7.253215E-04*T**2
	-1.588679E-07*T**3+10285.15*T**(-1); 
 1100	Y
	+393886.928-44.1074654*T-19.96003*T*LN(T)+.001513089*T**2
	-4.23648333E-07*T**3-722515*T**(-1); 
 6000	N
@@ Cr2 (g)
 ENTER-SYMBOL FUNCTION F7763T 298.15 
	+598511.403+41.5353212*T-40.56798*T*LN(T)+.004961847*T**2
	-1.61216717E-06*T**3+154422.85*T**(-1);
 800	Y
	+613345.232-104.207991*T-19.7643*T*LN(T)-.007085085*T**2
	-4.69883E-07*T**3-1738066.5*T**(-1);  
 1400	Y
	+642608.843-369.28626*T+17.64743*T*LN(T)-.02767321*T**2
	+1.605906E-06*T**3-5831655*T**(-1);  
 2300	Y
	+553119.895+159.188555*T-52.07969*T*LN(T)-.004229401*T**2
	+1.5939925E-07*T**3+14793625*T**(-1);  
 3900	Y
	+347492.34+623.137623*T-105.0428*T*LN(T)+3.9699545E-04*T**2
	+1.51783483E-07*T**3+1.4843765E+08*T**(-1);  
 5800	Y
	-484185.055+2598.25559*T-334.7145*T*LN(T)+.028597625*T**2
	-4.97520167E-07*T**3+7.135805E+08*T**(-1);  
 6000	N
@@ oxygen gas
@@ O(g)
 ENTER-SYMBOL FUNCTION F7397T 298.15  
	+243206.529-42897.0876*T**(-1)-0.7513421*T-21.0155542*T*LN(T)
	+1.26870532E-04*T**2-1.23131285E-08*T**3;   
 2950	Y  
	+252301.473-3973170.33*T**(-1)-51.974853*T-17.2118798*T*LN(T)
	-5.41356254E-04*T**2+7.64520703E-09*T**3;  
 6000	N
@@ O2(g)
 ENTER-SYMBOL FUNCTION F7531T 298.15  
	-6961.74451-76729.7484*T**(-1)-51.0057202*T-22.2710136*T*LN(T)
	-.0101977469*T**2+1.32369208E-06*T**3;  
 1000	Y  
	-13137.5203+525809.556*T**(-1)+25.3200332*T-33.627603*T*LN(T)
	-.00119159274*T**2+1.35611111E-08*T**3;  
 3300	Y  
	-27973.4908+8766421.4*T**(-1)+62.5195726*T-37.9072074*T*LN(T)
	-8.50483772E-04*T**2+2.14409777E-08*T**3;  
 6000	N
@@ O3(g)
 ENTER-SYMBOL FUNCTION F7683T 298.15  
	+133772.042-11328.9959*T**(-1)-84.8602165*T-19.8314069*T*LN(T)
	-.0392015696*T**2+7.90727187E-06*T**3;  
 600	Y  
	+120765.524+997137.156*T**(-1)+120.113376*T-51.8410152*T*LN(T)
	-.00353983136*T**2+3.20640143E-07*T**3;  
 1500	Y  
	+115412.196+1878139.02*T**(-1)+164.679664*T-58.069736*T*LN(T)
	-2.84399032E-04*T**2+5.95650279E-10*T**3;  
 6000	N
@@ La oxide gas
@@ LaO(g)
 ENTER-SYMBOL FUNCTION F12049T 298.15 
	-131496.968-24.5469483*T-31.53764*T*LN(T)-0.0051956*T**2
	+7.60442333E-07*T**3+103677.85*T**(-1);
 900	Y
	-133112.849+7.93847638*T-36.65559*T*LN(T)+2.4937065E-04*T**2
	-2.05688333E-07*T**3+108868.35*T**(-1);
 2500	Y
	-137735.323-23.9414477*T-31.58251*T*LN(T)-0.003177688*T**2
	+6.84986667E-08*T**3+5676870*T**(-1);
 5400	Y
	-39118.6731-213.786313*T-10.21743*T*LN(T)-0.005021225*T**2
	+9.162985E-08*T**3-74562000*T**(-1);
 10000	N
@@ La2O(g)
 ENTER-SYMBOL FUNCTION F12085T 298.15  
	-69316.3279+46.9195461*T-51.12563*T*LN(T)-0.005701935*T**2
	+8.637425E-07*T**3+212452.95*T**(-1);
 1000	Y
	-73167.5826+93.9280358*T-58.13034*T*LN(T)-1.332372E-05*T**2
	+4.41584333E-10*T**3+616730*T**(-1); 
 6000	N 
@@ La2O2(g)
 ENTER-SYMBOL FUNCTION F12089T 298.15 
	-642532.065+54.8487791*T-51.72813*T*LN(T)-0.028452875*T**2
	+4.99643833E-06*T**3+271002.95*T**(-1);
 700	Y
	-657497.385+256.452173*T-82.32033*T*LN(T)-1.8245965E-04*T**2
	+6.891315E-09*T**3+1664162*T**(-1);
 5100	Y
	-672122.137+293.83814*T-86.72291*T*LN(T)+4.319301E-04*T**2
	-9.75906E-09*T**3+11187120*T**(-1); 
 6000	N 
@@ Sr oxide Gas
@@ SrO(g)
 ENTER-SYMBOL FUNCTION F13511T 298.15  
	-25476.9742+3.04351985*T-34.37623*T*LN(T)-0.0026980695*T**2
	+3.78874167E-7*T**3+120146.05*T**(-1);
 900	Y  
	-44602.142+205.651627*T-63.83687*T*LN(T)+0.017645965*T**2
	-2.284235E-6*T**3+2463047*T**(-1);
 1800	Y  
	+243278.077-1500.21201*T+161.9497*T*LN(T)-0.0612273*T**2
	+2.896125E-6*T**3-66468000*T**(-1);
 2900	Y  
	-571113.316+1685.71589*T-234.6556*T*LN(T)+0.024571595*T**2
	-5.82819833E-7*T**3+2.468897E+8*T**(-1);
 4500	Y  
	-14433.8514+256.066959*T-66.76292*T*LN(T)+0.002226246*T**2
	-2.98498E-8*T**3-97083400*T**(-1);
 8800	Y  
	+52967.3441+134.904343*T-53.17021*T*LN(T)+0.001008387*T**2
	-9.46948833E-9*T**3-1.6008755E+8*T**(-1);
 10000	N 
@@ Mn oxide gas
@@ MN1O1(g)
 ENTER-SYMBOL FUNCTION F12454T 298.15
	+151356.092-24.1662496*T-31.02229*T*LN(T)-.005519525*T**2
	+8.47063667E-07*T**3+97921.9*T**(-1);
 900	Y
	+146503.882+29.4389965*T-38.90434*T*LN(T)+4.8722855E-04*T**2
	-5.46958333E-08*T**3+687003.5*T**(-1);
 4100	Y
	+163891.797+1.24225715*T-35.94693*T*LN(T)+6.220655E-04*T**2
	-7.949125E-08*T**3-12461235*T**(-1);
 6000	N
@@ MN1O2(g) 
 ENTER-SYMBOL FUNCTION F12461T 298.15
	+9095.11807-9.18675027*T-36.92694*T*LN(T)-.01944391*T**2
	+3.43552E-06*T**3+201368.45*T**(-1);
 700	Y 
	-4061.05468+167.414546*T-63.68895*T*LN(T)+.004964408*T**2
	-7.80481667E-07*T**3+1423101*T**(-1);
 2000	Y
	+74138.9718-312.782575*T+.1646534*T*LN(T)-.01801756*T**2
	+7.59669E-07*T**3-16355180*T**(-1);
 3900	Y 
	-269662.307+860.503282*T-143.1466*T*LN(T)+.008591745*T**2
	-1.67854167E-07*T**3+1.388169E+08*T**(-1);
 6000	N
@@ Cr-oxide Gas
@@ Cr1O1(g) 
 ENTER-SYMBOL FUNCTION F7705T 298.15
	+176483.869-31.9513659*T-30.2897*T*LN(T)-.00607059*T**2
	+9.229905E-07*T**3+35263.135*T**(-1); 
 900	Y
	+170853.62+32.1684007*T-39.74749*T*LN(T)+.00119977*T**2
	-1.52515733E-07*T**3+682877*T**(-1); 
 4000	Y 
	+307209.502-414.237405*T+14.48744*T*LN(T)-.008463125*T**2
	+1.722975E-07*T**3-64209900*T**(-1); 
 8400	Y
	-403765.708+805.224944*T-121.5329*T*LN(T)+.003139382*T**2
	-1.36845867E-08*T**3+6.35563E+08*T**(-1); 
 10000	N
 @ CR1O2(g)
 ENTER-SYMBOL FUNCTION F7713T 298.15
	-121743.409-9.3513573*T-36.29523*T*LN(T)-.01993456*T**2
	+3.5091E-06*T**3+193490.75*T**(-1); 
 700	Y
	-130807.622+121.470522*T-56.35306*T*LN(T)-3.3428045E-04*T**2
	-6.06589333E-08*T**3+956146*T**(-1); 
 3400	Y 
	-110536.941+26.3566562*T-44.22742*T*LN(T)-.003455125*T**2
	+8.58097667E-08*T**3-5175510*T**(-1); 
 6000	N
 @ CR1O3(g)  
 ENTER-SYMBOL FUNCTION F7726T 298.15
	-342387.72+103.718924*T-52.84504*T*LN(T)-.02759153*T**2
	+4.86524667E-06*T**3+358400*T**(-1); 
 700	Y
	-356869.109+299.291672*T-82.53502*T*LN(T)-7.979305E-05*T**2
	-1.88415667E-09*T**3-2.06280000E+01*T**(-1); 
 4000	Y 
	-362423.325+334.459538*T-87.10177*T*LN(T)+.001170557*T**2
	-5.86174E-08*T**3+1764429*T**(-1); 
 6000	N

--------------------------------------------------------------------
@@ Phases
@@ -------------------------------------------------------------------
@@ Binary oxides
@@
@@ Sr oxides, Risold et al. 1996 [50]
@@
 ENTER-PHASE SRO2,, 1 SRO2;,,,
 ENTER-PAR G(SRO2,SRO2;0) 298.15  +GSRO2SOL;,,,
@@ 
@@ Mn oxides, Grundy et al. 2005 [22]
@@
@@ STOICHIOMETRIC PYROLYSITE, MN1O2
 ENTER-PHASE MN1O2,, 2 1 2 MN; O;,,,
@@ Phase with 2 sublattices
@@ 1 atom of species Mn in the first sublattice and 2 atoms of species O in the second sublattice
 ENTER-PAR G(MN1O2,MN:O;0) 298.15  +GMN1O2;,,,
@@
@@ ----------------------------------------------------------------------
@@ Pseudoternary oxides
@@
@@ LA-SR-O, Grundy et al. 2004 [25] 
@@
 ENTER-PHASE LA2O3SS,, 2 2 3 LA+2 LA+3 SR+2; O-2 VA;,,,
 ENTER-PAR G(LA2O3SS,LA+2:O-2;0) 298.15  +GLAO;,,,
 ENTER-PAR G(LA2O3SS,LA+2:VA;0),,  +GLAO-GHSEROO;,,,
 ENTER-PAR G(LA2O3SS,LA+3:O-2;0) 298.15  +GLA2O3D;,,,
 ENTER-PAR G(LA2O3SS,LA+3:VA;0),,  +GLA2O3D-3*GHSEROO;,,,
 ENTER-PAR G(LA2O3SS,SR+2:O-2;0),,  +SR_ALPHA+GHSEROO
    +15.87691*T;,,,
 ENTER-PAR G(LA2O3SS,SR+2:VA;0),,  +SR_ALPHA-2*GHSEROO
    +15.87691*T;,,,
 ENTER-PAR L(LA2O3SS,LA+3,SR+2:O-2;0) 298.15  +2.149E+05-7.81E+01*T;,,,
 ENTER-PAR L(LA2O3SS,LA+3,SR+2:VA;0) 298.15  +2.149E+05-7.81E+01*T;,,,
@@ Interaction energies of zeroth order for La+3 and Sr+2 in the first sublattice
@@
 ENTER-PHASE LA2O3_HEXSS,, 2 2 3 LA+3 SR+2; O-2 VA;,,,
 ENTER-PAR G(LA2O3_HEXSS,LA+3:O-2;0),,  +GLA2O3H;,,,
 ENTER-PAR G(LA2O3_HEXSS,LA+3:VA;0),,  +GLA2O3H-3*GHSEROO;,,,
 ENTER-PAR G(LA2O3_HEXSS,SR+2:O-2;0),,  +SRH_ALPH+GHSEROO+15.87691*T;,,,
 ENTER-PAR G(LA2O3_HEXSS,SR+2:VA;0),,  +SRH_ALPH-2*GHSEROO+15.87691*T;,,,
 ENTER-PAR L(LA2O3_HEXSS,LA+3,SR+2:O-2;0) 298.15  +193600-78.1*T;,,,
 ENTER-PAR L(LA2O3_HEXSS,LA+3,SR+2:VA;0) 298.15  +193600-78.1*T;,,,
@@
 ENTER-PHASE LA2O3_CUBSS,, 2 2 3 LA+3 SR+2; O-2 VA;,,,
 ENTER-PAR G(LA2O3_CUBSS,LA+3:O-2;0),,  +GLA2O3X;,,,
 ENTER-PAR G(LA2O3_CUBSS,LA+3:VA;0),,  +GLA2O3X-3*GHSEROO;,,,
 ENTER-PAR G(LA2O3_CUBSS,SR+2:O-2;0),,  +SRX_ALPH+GHSEROO+15.87691*T;,,,
 ENTER-PAR G(LA2O3_CUBSS,SR+2:VA;0),,  +SRX_ALPH-2*GHSEROO+15.87691*T;,,,
 ENTER-PAR L(LA2O3_CUBSS,LA+3,SR+2:O-2;0) 298.15  +168700-78.1*T;,,,
 ENTER-PAR L(LA2O3_CUBSS,LA+3,SR+2:O-2;1) 298.15  -20000;,,,
@@ Interaction energy of first order for La+3 and Sr+2 in the first sublattice
 ENTER-PAR L(LA2O3_CUBSS,LA+3,SR+2:VA;0) 298.15  +168700-78.1*T;,,,
 ENTER-PAR L(LA2O3_CUBSS,LA+3,SR+2:VA;1) 298.15  -20000;,,,
@@ La4SrO7 as BETA phase
 ENTER-PHASE BETA,, 2 2 3 LA+3 SR+2; O-2 VA;,,,
 ENTER-PAR G(BETA,LA+3:O-2;0),,  +LA_BETA;,,,
 ENTER-PAR G(BETA,LA+3:VA;0),,  +LA_BETA-3*GHSEROO;,,,
 ENTER-PAR G(BETA,SR+2:O-2;0),,  +SR_ALPHA+416100+GHSEROO
    -0.33333333*RE_BETA+15.87691*T;,,,
 ENTER-PAR G(BETA,SR+2:VA;0),,  +SR_ALPHA+416100-2*GHSEROO
    +0.66666667*RE_BETA+15.87691*T;,,,
 ENTER-PAR L(BETA,LA+3,SR+2:O-2;0) 298.15  -121000-237.8*T;,,,
 ENTER-PAR L(BETA,LA+3,SR+2:VA;0) 298.15  -121000-237.8*T;,,,
@@ La4Sr3O9, Stoichiometric
 ENTER-PHASE LA4SR3O9,, 3 4 3 9 LA+3; SR+2; O-2;,,,
 ENTER-PAR  G(LA4SR3O9,LA+3:SR+2:O-2;0),,  +2*GLA2O3D+3*GSROSOL+229800
    -136.75*T;,,,
@@ SrO Solid Solution 
 ENTER-PHASE SRO,, 2 1 1 LA+3 SR+2 VA; O-2;,,,
 ENTER-PAR G(SRO,LA+3:O-2;0),,  +0.5*GLA2O3D+113700;,,,
 ENTER-PAR G(SRO,SR+2:O-2;0),,  +GSROSOL;,,,
 ENTER-PAR G(SRO,VA:O-2;0),,  0;,,,
@@
@@ LA-MN-O, Grundy et al. 2005 [23] 
@@
 ENTER-PHASE L2MNO4,, 3 2 1 4 LA+3; MN+2; O-2;,,,
 ENTER-PAR G(L2MNO4,LA+3:MN+2:O-2;0) 298.15  +GL2MNO4;,,,
@@
 ENTER-PHASE LMN2O5,, 4 1 1 1 5 LA+3; MN+3; MN+4; O-2;,,,
 ENTER-PAR G(LMN2O5,LA+3:MN+3:MN+4:O-2;0) 298.15  +GLMN2O5;,,,
@@
@@ LA-CR-O, Povoden et al. 2009 [20]
@@
@@ STOICHIOMETRIC LA2CRO6
 ENTER-PHASE LA2CRO6,, 3 2 1 6 LA+3; CR+6; O-2;,,,
 ENTER-PAR G(LA2CRO6,LA+3:CR+6:O-2;0) 298.15  +GLA2CRO6;,,,
@@
@@ STOICHIOMETRIC LA2CR3O12 
 ENTER-PHASE LA2CR3,, 3 2 3 12 LA+3; CR+6; O-2;,,,
 ENTER-PAR G(LA2CR3,LA+3:CR+6:O-2;0) 298.15  +GLA2CR3;,,,
@@
@@ SR-MN-O, Grundy et al. 2005 [26]
@@
 ENTER-PHASE SR7MN4O15,, 3 7 4 15 SR+2; MN+4; O-2;,,,
 ENTER-PAR G(SR7MN4O15,SR+2:MN+4:O-2;0),,  +GS7M4;,,,
@@
 ENTER-PHASE SRMN3O6,, 5 1 2 3 1 3 SR+2; MN+3; O-2; MN+3 MN+4; O-2 VA;,,,
 ENTER-PAR G(SRMN3O6,SR+2:MN+3:O-2:MN+3:O-2;0),,  +GSM3OZ+0.5*GHSEROO
    +11.23859*T;,,,
 ENTER-PAR G(SRMN3O6,SR+2:MN+3:O-2:MN+3:VA;0),,  +GSM3OZ-2.5*GHSEROO
    +11.23859*T;,,,
 ENTER-PAR G(SRMN3O6,SR+2:MN+3:O-2:MN+4:O-2;0),,  +GSM4OZ;,,,
 ENTER-PAR G(SRMN3O6,SR+2:MN+3:O-2:MN+4:VA;0),,  +GSM4OZ-3*GHSEROO;,,,
@@
 ENTER-PHASE SR4MN3O10,, 3 4 3 10 SR+2; MN+4; O-2;,,,
 ENTER-PAR G(SR4MN3O10,SR+2:MN+4:O-2;0),,  +GSM4_RP3;,,,
@@
 ENTER-PHASE SRMNO3_HEX,, 3 1 1 3 SR+2; MN+3 MN+4; O-2 VA;,,, 
 ENTER-PAR G(SRMNO3_HEX,SR+2:MN+3:O-2;0),,  +GSM3_HEX+0.5*GHSEROO+11.23859*T;,,,
 ENTER-PAR G(SRMNO3_HEX,SR+2:MN+3:VA;0),,  +GSM3_HEX-2.5*GHSEROO+11.23859*T;,,,
 ENTER-PAR G(SRMNO3_HEX,SR+2:MN+4:O-2;0),,  +GSM4_HEX;,,,
 ENTER-PAR G(SRMNO3_HEX,SR+2:MN+4:VA;0),,  +GSM4_HEX-3*GHSEROO;,,,
@@
@@ RP2, Grundy et al. 2005 [26]
 ENTER-PHASE RP2,, 4 1 2 2 7 SR+2; SR+2; MN+3 MN+4; O-2;,,,
 ENTER-PAR G(RP2,SR+2:SR+2:MN+3:O-2;0) 298.15  +3*GSROSOL+GMN2O3+GHSEROO;,,,
 ENTER-PAR G(RP2,SR+2:SR+2:MN+4:O-2;0) 298.15  +GS4O_RP2;,,,
@@
@@ SR-CR-O
@@
@@ SRCR2O4
 ENTER-PHASE SC2O4,, 3 1 2 4 SR+2; CR+3; O-2;,,,
 ENTER-PAR G(SC2O4,SR+2:CR+3:O-2;0) 298.15  +GSC2O4;,,,
@@ SR2CRO4
 ENTER-PHASE S2CO4,, 3 2 1 4 SR+2; CR+4; O-2;,,,
 ENTER-PAR G(S2CO4,SR+2:CR+4:O-2;0) 298.15  +GS2CO4;,,,
@@
@@ SR2.67CR2O8
 ENTER-PHASE S3C2N,, 3 2.666667 2 8 SR; CR; O;,,,
 ENTER-PAR G(S3C2N,SR:CR:O;0) 298.15  +GS3C2O8;,,,
@@ SRCRO4
 ENTER-PHASE SCO4,, 3 1 1 4 SR+2; CR+6; O-2;,,,
 ENTER-PAR G(SCO4,SR+2:CR+6:O-2;0) 298.15  +GSCO4;,,,
@@
@@ CR-MN-O, [21,22] 
@@
@@ NONSTOICHIOMETRIC MANGANOSITE (MNO) SOLID SOLUTION [19,22]
 ENTER-PHASE MNO_HALIT,, 2 1 1 MN+2 MN+3 CR+3 VA; O-2;,,,
 ENTER-PAR G(MNO_HALIT,MN+2:O-2;0) 298.15  +GMN1O1;,,,
 ENTER-PAR G(MNO_HALIT,MN+3:O-2;0) 298.15  +GMN1O1-21883.5213
    -22.1853365*T;,,,
 ENTER-PAR G(MNO_HALIT,CR+3:O-2;0) 298.15  +0.5*GCR2O3-7.93845*T
    +71549.3;,,,
 ENTER-PAR G(MNO_HALIT,VA:O-2;0) 298.15  0;,,,
 ENTER-PAR L(MNO_HALIT,MN+2,MN+3:O-2;0) 298.15  -4.21048766E+04;,,,
 ENTER-PAR L(MNO_HALIT,MN+2,MN+3:O-2;1) 298.15  +4.65131533E+04;,,,
@@
@@ Mn2O3 [19,22]
 ENTER-PHASE MN2O3,, 3 2 3 1 MN+3 CR+3; O-2 VA; O-2 VA;,,,
 AMEND-PHASE MN2O3 MAGN,, .28
 ENTER-PAR TC(MN2O3,MN+3:O-2:VA;0) 298.15  +309;,,,
 ENTER-PAR TC(MN2O3,CR+3:O-2:VA;0) 298.15  +308.6;,,,
 ENTER-PAR BMAGN(MN2O3,MN+3:O-2:VA;0) 298.15  +0.59;,,,
 ENTER-PAR BMAGN(MN2O3,CR+3:O-2:VA;0) 298.15  +3;,,,
 ENTER-PAR G(MN2O3,MN+3:O-2:VA;0) 298.15  +GMN2O3;,,,
 ENTER-PAR G(MN2O3,MN+3:O-2:O-2;0) 298.15  +GMN2O3+GHSEROO
    +100000+15.87691*T;,,,
 ENTER-PAR G(MN2O3,MN+3:VA:O-2;0) 298.15  +GMN2O3-2*GHSEROO
    +100000+15.87691*T;,,,
 ENTER-PAR G(MN2O3,MN+3:VA:VA;0) 298.15  +GMN2O3-3*GHSEROO;,,,
 ENTER-PAR G(MN2O3,CR+3:O-2:VA;0) 298.15  +GCR2O3+3459;,,,
 ENTER-PAR G(MN2O3,CR+3:O-2:O-2;0) 298.15  +GCR2O3+GHSEROO
    +100000+15.87691*T;,,,
 ENTER-PAR G(MN2O3,CR+3:VA:O-2;0) 298.15  +GCR2O3-2*GHSEROO
    +100000+15.87691*T;,,,
 ENTER-PAR G(MN2O3,CR+3:VA:VA;0) 298.15  +GCR2O3-3*GHSEROO;,,,
@@ 
@@ ESKOLAITE, NONSTOICHIOMETRIC CR2O3 SOLID SOLUTION [19,22]
 ENTER-PHASE CR2O3,, 3 2 1 3 MN+3 CR+2 CR+3; CR+3 VA; O-2;,,,
 AMEND-PHASE CR2O3 MAGN,, .28
 ENTER-PAR TC(CR2O3,MN+3:CR+3:O-2;0) 298.15  +309;,,,
 ENTER-PAR TC(CR2O3,MN+3:VA:O-2;0) 298.15  +309;,,,
 ENTER-PAR TC(CR2O3,CR+2:CR+3:O-2;0) 298.15  +308.6;,,,
 ENTER-PAR TC(CR2O3,CR+2:VA:O-2;0) 298.15  +308.6;,,,
 ENTER-PAR TC(CR2O3,CR+3:CR+3:O-2;0) 298.15  +308.6;,,,
 ENTER-PAR TC(CR2O3,CR+3:VA:O-2;0) 298.15  +308.6;,,,
 ENTER-PAR BMAGN(CR2O3,MN+3:CR+3:O-2;0) 298.15  +0.59;,,, 
 ENTER-PAR BMAGN(CR2O3,MN+3:VA:O-2;0) 298.15  +0.59;,,,
 ENTER-PAR BMAGN(CR2O3,CR+2:CR+3:O-2;0) 298.15  +3;,,,
 ENTER-PAR BMAGN(CR2O3,CR+2:VA:O-2;0) 298.15  +3;,,,
 ENTER-PAR BMAGN(CR2O3,CR+3:CR+3:O-2;0) 298.15  +3;,,,
 ENTER-PAR BMAGN(CR2O3,CR+3:VA:O-2;0) 298.15  +3;,,,
 ENTER-PAR G(CR2O3,MN+3:CR+3:O-2;0) 298.15  +GMN2O3+GSERCR+39503;,,,
 ENTER-PAR G(CR2O3,MN+3:VA:O-2;0) 298.15  +GMN2O3+39503;,,,
 ENTER-PAR G(CR2O3,CR+2:CR+3:O-2;0) 298.15  +GCRO0
    +.33333333334*GHSERCR-5.2923*T;,,,
 ENTER-PAR G(CR2O3,CR+2:VA:O-2;0) 298.15  +GCRO0
    -0.666666666667*GHSERCR-5.2923*T;,,,
 ENTER-PAR G(CR2O3,CR+3:CR+3:O-2;0) 298.15  +GCR2O3+GHSERCR;,,,
 ENTER-PAR G(CR2O3,CR+3:VA:O-2;0) 298.15  +GCR2O3;,,,
@@
@@ CUBIC SPINEL [21,22]
 ENTER-PHASE CMNCR2O4,, 3 1 2 4 MN+2 CR+2; MN+3 CR+3; O-2;,,,
 ENTER-PAR G(CMNCR2O4,MN+2:CR+3:O-2;0) 298.15  +GSPINEL;,,,
 ENTER-PAR G(CMNCR2O4,CR+2:CR+3:O-2;0) 298.15  +GCR3O4;,,,
 ENTER-PAR G(CMNCR2O4,MN+2:MN+3:O-2;0) 298.15  +GCMN3O4;,,,
 ENTER-PAR G(CMNCR2O4,CR+2:MN+3:O-2;0) 298.15  +GCR3O4+GCMN3O4-GSPINEL;,,,
@@ DISTORTED_SPINEL, Grundy, Povoden [21,22]
 ENTER-PHASE TMNCR2O4,, 3 1 2 4 MN+2; MN+3 CR+3; O-2;,,,
 ENTER-PAR G(TMNCR2O4,MN+2:MN+3:O-2;0) 298.15  +GTMN3O4;,,,
 ENTER-PAR G(TMNCR2O4,MN+2:CR+3:O-2;0) 298.15  +GTSPINEL;,,,
@@ ----------------------------------------------------------------------------
@@ Pseudoquaternary oxides
@@
@@ La-Sr-Mn-O, Grundy et al. 2005 [27]
@@ 
 ENTER-PHASE LS3MN2O7,, 4 1 2 2 7 SR+2; LA+3 SR+2; MN+3 MN+4; O-2;,,,
 ENTER-PAR G(LS3MN2O7,SR+2:LA+3:MN+3:O-2;0) 298.15  +GL3O_RP2;,,,
 ENTER-PAR G(LS3MN2O7,SR+2:LA+3:MN+4:O-2;0) 298.15  +GL3O_RP2+GS4O_RP2-3*GSROSOL
    -GMN2O3-GHSEROO-R_RP2;,,,
 ENTER-PAR G(LS3MN2O7,SR+2:SR+2:MN+3:O-2;0) 298.15  +3*GSROSOL+GMN2O3+GHSEROO;,,,
 ENTER-PAR G(LS3MN2O7,SR+2:SR+2:MN+4:O-2;0) 298.15  +GS4O_RP2;,,,
@@ ---------------------------------------------------------------------------
@@ Pseudoquinary phases, ideal extensions from lower-order (pseudoternary and pseudoquaternary) systems
@@
@@ High Temperature RHOMBOHEDRAL perovskite
 ENTER-PHASE RPEROV,, 3 1 1 3 LA+3 SR+2 VA; MN+2 MN+3 MN+4 CR+3 CR+4 VA; 
    O-2 VA;,,,
@@
 ENTER-PAR G(RPEROV,SR+2:MN+2:O-2;0),,  +GMS3O+GL2O-GL3OR+GHSEROO+22.47717*T;,,,
 ENTER-PAR G(RPEROV,SR+2:MN+3:O-2;0),,	  +GMS3O+0.5*GHSEROO+11.23859*T;,,,
 ENTER-PAR G(RPEROV,SR+2:MN+4:O-2;0),,	 +GMS4O;,,,
 ENTER-PAR G(RPEROV,SR+2:MN+2:VA;0),,  +GMS3O+GL2O-GL3OR-2*GHSEROO+22.47717*T;,,,
 ENTER-PAR G(RPEROV,SR+2:MN+3:VA;0),,  +GMS3O-2.5*GHSEROO+11.23859*T;,,,
 ENTER-PAR G(RPEROV,SR+2:MN+4:VA;0),,  +GMS4O-3*GHSEROO;,,,
 ENTER-PAR G(RPEROV,SR+2:CR+3:O-2;0),,	  +GS3V+0.1666667*GS4O-0.16666667*GS4V;,,,
 ENTER-PAR G(RPEROV,SR+2:CR+3:VA;0),,  +GS3V-0.8333333*GS4O+0.8333333*GS4V;,,,
 ENTER-PAR G(RPEROV,SR+2:VA:O-2;0),,  +GMS3O-GL3OR+2*GL4O-1.5*GV4O+0.5*GVVV
    +2*GHSEROO+12.62121*T;,,,
 ENTER-PAR G(RPEROV,SR+2:VA:VA;0),,  +GMS3O+2*GL4O-1.5*GV4O+0.5*GVVV-GL3OR
    -GHSEROO+12.62121*T;,,,
 ENTER-PAR G(RPEROV,LA+3:MN+2:O-2;0),,  +GL2O+0.5*GHSEROO+11.2386*T;,,,
 ENTER-PAR G(RPEROV,LA+3:MN+2:VA;0),,  +GL2O-2.5*GHSEROO+11.2386*T;,,,
 ENTER-PAR G(RPEROV,LA+3:MN+3:O-2;0),,  +GL3OR;,,,
 ENTER-PAR G(RPEROV,LA+3:MN+3:VA;0),,  +GL3OR-3*GHSEROO;,,,
 ENTER-PAR G(RPEROV,LA+3:MN+4:O-2;0),,  +0.66667*GL4O+0.5*GV4O-0.166667*GVVV
    -0.5*GHSEROO+5.76318*T;,,,
 ENTER-PAR G(RPEROV,LA+3:MN+4:VA;0),,  +0.66667*GL4O+0.5*GV4O-0.166667*GVVV
    -3.5*GHSEROO+5.76318*T;,,,
 ENTER-PAR G(RPEROV,LA+3:VA:O-2;0),,  +2*GL4O-1.5*GV4O+0.5*GVVV+1.5*GHSEROO
    +1.41263*T;,,,
 ENTER-PAR G(RPEROV,LA+3:VA:VA;0),,  +2*GL4O+0.5*GVVV-1.5*GV4O-1.5*GHSEROO
    +1.41263*T;,,,
 ENTER-PAR G(RPEROV,VA:MN+3:O-2;0),,  +GL3OR+1.5*GV4O+0.5*GVVV-2*GL4O
    +1.5*GHSEROO-1.41263*T;,,,
 ENTER-PAR G(RPEROV,VA:MN+4:O-2;0),,  +2*GV4O+0.33333*GVVV-1.33333*GL4O
    +GHSEROO+4.35056*T;,,,
 ENTER-PAR G(RPEROV,VA:MN+2:VA;0),,  +GL2O+1.5*GV4O+0.5*GVVV-2*GL4O
    -GHSEROO+9.82596*T;,,,
 ENTER-PAR G(RPEROV,VA:MN+3:VA;0),,  +GL3OR+1.5*GV4O+0.5*GVVV-2*GL4O
    -1.5*GHSEROO-1.41263*T;,,,
 ENTER-PAR G(RPEROV,VA:MN+4:VA;0),,  +2*GV4O+0.3333*GVVV-1.333*GL4O
    -2*GHSEROO+4.35057*T;,,,
 ENTER-PAR G(RPEROV,VA:VA:O-2;0),,  +GVVV+3*GHSEROO;,,,
 ENTER-PAR G(RPEROV,VA:VA:VA;0),,  +GVVV;,,,
 ENTER-PAR G(RPEROV,LA+3:CR+3:O-2;0),,  +GLACRO3;,,,
 ENTER-PAR G(RPEROV,LA+3:CR+3:VA;0),,  +GLACRO3-3*GHSEROO;,,,
 ENTER-PAR G(RPEROV,VA:CR+3:O-2;0),,  +GLACRO3+1.5*GVCR4O+0.5*GVVV-2*GLACR4O
    +1.5*GHSEROO-1.41263*T;,,,
 ENTER-PAR G(RPEROV,VA:CR+3:VA;0),,  +GLACRO3+1.5*GVCR4O+0.5*GVVV-2*GLACR4O
    -1.5*GHSEROO-1.41263*T;,,,
 ENTER-PAR G(RPEROV,LA+3:CR+4:VA;0),,  +GS4O-GS3V-0.1666667*GS4O
    +0.16666667*GS4V+GLACRO3-3*GHSEROO;,,,
 ENTER-PAR G(RPEROV,LA+3:CR+4:O-2;0),,  +GS4O+GLACRO3-GS3V-0.1666667*GS4O
    +0.16666667*GS4V;,,,
 ENTER-PAR G(RPEROV,SR+2:CR+4:O-2;0),,  +GS4O;,,,
 ENTER-PAR G(RPEROV,SR+2:CR+4:VA;0),,  +GS4V;,,,
 ENTER-PAR G(RPEROV,VA:CR+4:O-2;0),,  +2*GVCR4O+0.33333*GVVV-1.33333*GLACR4O
    +GHSEROO+4.35056*T;,,,
 ENTER-PAR G(RPEROV,VA:CR+4:VA;0),,  +2*GVCR4O+0.3333*GVVV-1.333*GLACR4O
    -2*GHSEROO+4.35057*T;,,,
@@ Optimized interactions
 ENTER-PAR L(RPEROV,LA+3,SR+2:MN+2:VA;1),,  -136600;,,,
 ENTER-PAR L(RPEROV,LA+3,SR+2:MN+2:O-2;1),,  -136600;,,,
 ENTER-PAR L(RPEROV,LA+3,SR+2:MN+4:VA;1),,  -117000;,,,
 ENTER-PAR L(RPEROV,LA+3,SR+2:MN+4:O-2;1),,  -117000;,,,
 ENTER-PAR L(RPEROV,LA+3:CR+3,MN+3:O-2;0),,  +8500;,,,
 ENTER-PAR L(RPEROV,LA+3,SR+2:MN+3:O-2;0) 298.15  -1.5*T;,,,
 ENTER-PAR L(RPEROV,LA+3:CR+4,VA:O-2;0),,  +250000;,,,
 ENTER-PAR L(RPEROV,LA+3:CR+3,VA:O-2;0),,  +250000;,,,
 ENTER-PAR L(RPEROV,LA+3:MN+4,MN+3:O-2;1) 298.15  +185;,,,
 ENTER-PAR L(RPEROV,LA+3:MN+4,VA:O-2;0),,  -2;,,,
 ENTER-PAR L(RPEROV,LA+3,VA:MN+4:O-2;0),,  +20;,,,
 ENTER-PAR L(RPEROV,LA+3:CR+3,MN+3:O-2;0) 298.15  +3766;,,,
 ENTER-PAR L(RPEROV,LA+3:CR+3,MN+3:O-2;1) 298.15  -1297;,,,
 ENTER-PAR L(RPEROV,LA+3,SR+2:CR+3:O-2;0),,  -39.5*T;,,,
 ENTER-PAR L(RPEROV,LA+3,SR+2:CR+4:O-2;0),,  -39.5*T;,,,
 ENTER-PAR L(RPEROV,SR+2:CR+4,MN+3:O-2;0) 298.15  0;,,,
 ENTER-PAR L(RPEROV,SR+2:CR+4,MN+4:O-2;0) 298.15  0;,,,
@@
--------------------------------------------------------------------------
@@ Liquid, ideal extension from lower-order (binary and pseudoternary) systems
 ENTER-PHASE IONIC_LIQUID Y  LA+3 SR+2 MN+2 MN+3 CR+2 CR+3; O-2 VA
   SRCRO3 SRCRO4;,,,
 AMEND-PHASE IONIC_LIQUID COMP 2,,,,,,,,
@@ Consider miscibility gaps
 ENTER-PAR G(ION,LA+3:O-2;0) 298.15  +GLA2O3LIQ;,,,
 ENTER-PAR G(ION,LA+3:VA;0) 298.15  +GLALIQ;,,,
 ENTER-PAR G(ION,SR+2:O-2;0) 298.15  +2*GSROLIQ;,,,
 ENTER-PAR G(ION,SR+2:VA;0) 298.15  +GSRLIQ;,,,
 ENTER-PAR G(ION,MN+2:O-2;0) 298.15  +2*GMN1O1_L;,,,
 ENTER-PAR G(ION,MN+2:VA;0) 298.15  +GMN_L;,,,
 ENTER-PAR G(ION,MN+3:O-2;0) 298.15  +GMN2O3_L;,,,
 ENTER-PAR G(ION,MN+3:VA;0) 298.15  +2*GMN_L+GMN2O3_L-3*GMN1O1_L;,,,
 ENTER-PAR G(ION,CR+2:O-2;0) 298.15  +2*GCR1O1_L;,,,
 ENTER-PAR G(ION,CR+2:VA;0) 298.15  +GCR_L;,,,
 ENTER-PAR G(ION,CR+3:O-2;0) 298.15  +GCR2O3_L;,,,
 ENTER-PAR G(ION,CR+3:VA;0) 298.15  +2*GCR_L+GCR2O3_L-3*GCR1O1_L;,,,
@@ Interaction parameters from binaries
 ENTER-PAR  L(ION,MN+2:O-2,VA;0)  298.15  +1.29519000E+05;,,,
 ENTER-PAR  L(ION,MN+2:O-2,VA;1)  298.15  -4.54590000E+04;,,,
 ENTER-PAR  L(ION,MN+2,MN+3:O-2;0)  298.15  -3.38590000E+04;,,,
 ENTER-PAR L(ION,CR+2:O-2,VA;0) 298.15  +121000;,,,
 ENTER-PAR L(ION,CR+3:O-2,VA;0) 298.15  +121000;,,,
@@ Interaction parameters from (pseudo)ternaries
 ENTER-PAR L(ION,LA+3,MN+2:O-2;0) 298.15  -119062;,,,
 ENTER-PAR L(ION,LA+3,MN+2:VA;0)  298.15  +11368;,,,
 ENTER-PAR L(ION,LA+3,MN+2:VA;1)  298.15  -11316.4;,,,
 ENTER-PAR L(ION,LA+3,MN+2:VA;4)  298.15  -9111;,,,
 ENTER-PAR L(ION,LA+3,MN+3:O-2;0) 298.15  -119062;,,,
 ENTER-PAR L(ION,SR+2,MN+2:O-2;0) 298.15  -176300;,,,
 ENTER-PAR L(ION,SR+2,MN+2:VA;0) 298.15  +46000;,,,
 ENTER-PAR L(ION,SR+2,MN+3:O-2;0) 298.15  -176300;,,,
 ENTER-PAR L(ION,MN+2,CR+2:VA;0) 298.15  -15009+13.6587*T;,,,
 ENTER-PAR L(ION,MN+2,CR+2:VA;1) 298.15  +504+0.9479*T;,,,
 ENTER-PAR L(ION,MN+3,CR+3:O-2;0) 298.15  -188487.7;,,,
 ENTER-PAR L(ION,LA+3,CR+2:O-2;0) 298.15  -101850;,,,
 ENTER-PAR L(ION,LA+3,CR+3:O-2;0) 298.15  -101850;,,,
 ENTER-PAR L(ION,LA+3,CR+2:O-2;1) 298.15  -39016;,,,
 ENTER-PAR L(ION,LA+3,CR+3:O-2;1) 298.15  -39016;,,,
 ENTER-PAR L(ION,LA+3,CR+2:VA;0) 298.15  +60713-5.49*T;,,,
 ENTER-PAR L(ION,LA+3,CR+2:VA;1) 298.15  +64573-23*T;,,,
 ENTER-PAR L(ION,SR+2,CR+2:VA;0) 298.15  +200000;,,,
 ENTER-PAR L(ION,SR+2,CR+3:O-2;0) 298.15  -275000;,,,
 ENTER-PAR L(ION,SR+2,CR+3:O-2;1) 298.15  -100000;,,,
 ENTER-PAR L(ION,SR+2,CR+2:O-2;0) 298.15  -275000;,,,
 ENTER-PAR L(ION,SR+2,CR+2:O-2;1) 298.15  -100000;,,,
 ENTER-PAR L(ION,SR+2,CR+3:O-2;2) 298.15  -250000;,,,
 ENTER-PAR L(ION,SR+2,CR+2:O-2;2) 298.15  -250000;,,,
 ENTER-PAR G(ION,SRCRO3;0) 298.15  +GS4O+29200-10*T;,,,
 ENTER-PAR G(ION,SRCRO4;0) 298.15  +GSCO4+50000-31*T;,,,
 ENTER-PAR G(ION,SR+2:SRCRO3,O-2;0) 298.15  -50000;,,,
 ENTER-PAR G(ION,SR+2:SRCRO4,O-2;0) 298.15  -100000;,,,
@@------------------------------------------------------------
@@ Gas
 ENTER-PHASE GAS G 1 LA SR SR2  MN CR CR2 O1 O2 O3 LA1O1 
                     LA2O1 LA2O2 SRO SR2O MN1O1 MN1O2  
                     CR1O1 CR1O2 CR1O3;,,,
 ENTER-PAR G(GAS,LA;0) 298.15  +F12026T+RTLNP;,,,
 ENTER-PAR G(GAS,SR;0) 298.15  +F15323T+RTLNP;,,,
 ENTER-PAR G(GAS,SR2;0) 298.15  +F15338T+RTLNP;,,,
 ENTER-PAR G(GAS,SR;0) 298.15  +SRGAS+RTLNP;,,,
 ENTER-PAR G(GAS,SR2;0) 298.15  +SR2GAS+RTLNP;,,,
 ENTER-PAR G(GAS,MN;0) 298.15  +F12439T+RTLNP;,,,   
 ENTER-PAR G(GAS,CR;0) 298.15  +F7491T+RTLNP;,,,,
 ENTER-PAR G(GAS,CR2;0) 298.15  +F7763T+RTLNP;,,,
 ENTER-PAR G(GAS,O1;0) 298.15  +F7397T+RTLNP;,,,
 ENTER-PAR G(GAS,O2;0) 298.15  +2*GHSEROO+RTLNP;,,,
 ENTER-PAR G(GAS,O3;0) 298.15  +F7683T+RTLNP;,,,
 ENTER-PAR G(GAS,LA1O1;0) 298.15  +F12049T+RTLNP;,,,    
 ENTER-PAR G(GAS,LA2O1;0) 298.15  +F12085T+RTLNP;,,,  
 ENTER-PAR G(GAS,LA2O2;0) 298.15  +F12089T+RTLNP;,,,
 ENTER-PAR G(GAS,SRO;0) 298.15  +F13511T+RTLNP;,,,
 ENTER-PAR G(GAS,SRO;0) 298.15  +SROGAS+RTLNP;,,,
 ENTER-PAR G(GAS,SR2O;0) 298.15  +SR2OGAS+RTLNP;,,,
 ENTER-PAR G(GAS,MN1O1;0) 298.15  +F12454T+RTLNP;,,,   
 ENTER-PAR G(GAS,MN1O2;0) 298.15  +F12461T+RTLNP;,,,
 ENTER-PAR G(GAS,CR1O1;0) 298.15  +F7705T+RTLNP;,,,
 ENTER-PAR G(GAS,CR1O2;0) 298.15  +F7713T+RTLNP;,,,
 ENTER-PAR G(GAS,CR1O3;0) 298.15  +F7726T+RTLNP;,,,
@@
@@ O2 gas
 ENTER-PHASE O2GAS,, 1 O2;,,,
 ENTER-PAR G(O2GAS,O2;0) 298.15  +2*GHSEROO+RTLNP;,,,
@@
 GO PAR
@@ Go to PARROT module of ThermoCalc
@@
 SET-OUT-LEVEL,,,,, N,,
 CREATE LSCM_v1.0
@@ Creates PAR-file 
 set-interactive