%MTDATA script - driving force
%24 April 2009, by Mathew Peet
% This is a quick example script to show calculation of
% driving force for formation of ferrite from austenite.
mu
define system ÔFe,Cr,Ni,Mo,N,C,MnÕ source tcfe !
% fe Cr Ni Mo N
set w(1) undef w(2) 23.00 w(3) 9.0 w(4) 3.0 w(5) 0.15 !
% C Mn
set w(6) 0.02 w(7) 0.8 !
set w 100 !
classify absent p(*) !
step temp 973 1273 10 !
%
% calculating for FCC
%
classify normal p(FCC_A1) !
compute print brief print mol !
compute print graphics_output !
units temperature celsius !
ordinate gibbs system !
plot tabulate spreadsheet !
%
% calculating for BCC + FCC
%
classify normal p(BCC_A2) !
compute print brief print mol !
compute print graphics_output !
units temperature celsius !
ordinate gibbs system !
plot tabulate spreadsheet !

% This gives the answer in Joules per 100 kg.
% For joules per kg set w=1 as follows (replace earlier lines with these)
%set w(1) undef w(2) 0.23 w(3) 0.09 w(4) 0.03 w(5) 0.0015 !
%set w(6) 0.0002 w(7) 0.008 !
%set w 1 !