Index: mttroot/mtt/lib/comp/simple/ES_eqn.m ================================================================== --- mttroot/mtt/lib/comp/simple/ES_eqn.m +++ mttroot/mtt/lib/comp/simple/ES_eqn.m @@ -15,10 +15,13 @@ % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % %% Version control history % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % %% $Id$ % %% $Log$ +% %% Revision 1.2 1998/02/25 16:59:56 peterg +% %% Corrected enthalpy-entropy conversion -- use temp on the entropy side. +% %% % %% Revision 1.1 1998/02/25 15:03:51 peterg % %% Initial revision % %% % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% @@ -41,15 +44,17 @@ e_bond = bonds(1,:); s_bond = bonds(2,:); % The efforts are the same - but the enthalpy side has relative T if e_bond(1)==1 %effort imposed on the enthalpy bond - fprintf(eqnfile,'%s := %s + t_0;\n' , ... +% fprintf(eqnfile,'%s := %s + t_0;\n' , ... + fprintf(eqnfile,'%s := %s;\n' , ... varname(name,bond_number(2), 1), ... varname(name,bond_number(1), 1) ); else %effort imposed on the entropy bond - fprintf(eqnfile,'%s := %s - t_0;\n' , ... +% fprintf(eqnfile,'%s := %s - t_0;\n' , ... + fprintf(eqnfile,'%s := %s;\n' , ... varname(name,bond_number(1), 1), ... varname(name,bond_number(2), 1) ); end; % The flows need to be converted - use the absolute temp on the