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.1  1998/02/25 15:03:51  peterg
+% %% Initial revision
+% %%
 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 
 
 % Copyright (c) P.J. Gawthrop, 1997.
 
@@ -36,11 +39,11 @@
 
 % There are 2 ports; extract the information
 e_bond = bonds(1,:);
 s_bond = bonds(2,:);
 
-% The efforts are the same
+% 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' , ...
     varname(name,bond_number(2), 1), ...
     varname(name,bond_number(1), 1) );
 else %effort imposed on the entropy bond 
@@ -47,17 +50,18 @@
   fprintf(eqnfile,'%s := %s - t_0;\n' , ...
     varname(name,bond_number(1), 1), ...
     varname(name,bond_number(2), 1) );
 end;
 
-% The flows need to be converted
+% The flows need to be converted - use the absolute temp on the
+% entropy side
 if e_bond(2)==-1 %flow imposed on the enthalpy bond 
-  fprintf(eqnfile,'%s := %s/(%s + t_0);\n' , ...
+  fprintf(eqnfile,'%s := %s/%s;\n' , ...
     varname(name,bond_number(2), -1), ...
     varname(name,bond_number(1), -1), ...
     varname(name,bond_number(2), 1) );
 else % flow imposed on the entropy bond 
-    fprintf(eqnfile,'%s := %s*(%s - t_0);\n' , ...
+    fprintf(eqnfile,'%s := %s*%s;\n' , ...
     varname(name,bond_number(1), -1), ...
     varname(name,bond_number(2), -1), ...
     varname(name,bond_number(2), 1) );
 end;