Index: mttroot/mtt/bin/trans/make_ode2odes
==================================================================
--- mttroot/mtt/bin/trans/make_ode2odes
+++ mttroot/mtt/bin/trans/make_ode2odes
@@ -7,10 +7,13 @@
 ###############################################################
 ## Version control history
 ###############################################################
 ## $Id$
 ## $Log$
+## Revision 1.88  2004/08/29 13:15:28  geraint
+## Uses sys_sae instead of sys_ae if sorted equations are being used.
+##
 ## Revision 1.87  2004/08/29 01:46:56  geraint
 ## Added rules to create ode2odes for sorted system: sesx and sesy.
 ##
 ## Revision 1.86  2004/08/29 00:19:49  geraint
 ## Defaults to noAlgebraicSolver.
@@ -371,10 +374,11 @@
 Nu=`mtt_getsize $sys u` # Inputs 
 Ny=`mtt_getsize $sys y` # Outputs  
 
 case "$method" in
     "implicit")
+	ae=ae
 	ode=csex
 	odeo=cseo
 	algorithm="mtt_implicit(x,dx,AA,AAx,ddt,$Nx,open_switches)"
 	;;
     "dassl")
@@ -383,14 +387,16 @@
 		ae=sae
 		ode=sesx
 		odeo=sesy
 		;;
 	    "seqn")		# shouldn't be here unless mtt has changed
+		ae=ae
 		ode=sese
 		odeo=sese
 		;;
 	    "no" | *)	
+		ae=ae
 		ode=ode
 		odeo=odeo
 		;;
 	esac
 	algorithm="mtt_dassl(x,u,t,par,dx,ddt,MTTNX,MTTNYZ,open_switches)"
@@ -404,14 +410,16 @@
 		ae=sae
 		ode=sesx
 		odeo=sesy
 		;;
 	    "seqn")		# shouldn't be here unless mtt has changed
+		ae=ae
 		ode=sese
 		odeo=sese
 		;;
 	    "no" | *)	
+		ae=ae
 		ode=ode
 		odeo=odeo
 		;;
 	esac
 	algorithm="mtt_euler(x,dx,ddt,$Nx,open_switches)"