Index: mttroot/mtt/bin/trans/ode2odes_r2c
==================================================================
--- mttroot/mtt/bin/trans/ode2odes_r2c
+++ mttroot/mtt/bin/trans/ode2odes_r2c
@@ -18,10 +18,14 @@
 ###############################################################
 ## Version control history
 ###############################################################
 ## $Id$
 ## $Log$
+## Revision 1.15  1998/05/14 08:20:29  peterg
+## Start time at DT - the result at time zero is computed outside the
+## loop
+##
 ## Revision 1.14  1998/05/13 08:57:27  peterg
 ## Now uses simpar.h in place of odes.h
 ##
 ## Revision 1.13  1998/02/24 13:34:45  peterg
 ## Back under RCS
@@ -167,11 +171,11 @@
 write "/* Set up initial state */"$
 write "  $1_state();"$
 
 %Initialise main (Euler) integration loop
 write "/* Initialise main (Euler) integration loop */"$
-write "  time = DT;"$
+write "  time = 0;"$
 write "  dt = DT/STEPFACTOR;"$
 %% write "  for (i=1; i<=MTTNX; i++)"$
 %% write "      x[i] = 0.0;"$
 
 %Set up system inputs
@@ -191,11 +195,11 @@
 write "  $1_ode(y,dx,x,u);"$
 
 %Main (Euler) integration loop
 write "/* Main (Euler) integration loop */"$
 
-write "  while (time<LAST)"$
+write "  while (time<(LAST-DT))"$
 write "  {"$
 
 %Write to output to file
 write "/* Write to output file */"$
 write "    fprintf(fpso, ""%5.4g "",time);"$
@@ -208,10 +212,11 @@
 write "    fprintf(fps, ""%5.4g "",time);"$
 write "    for (i=1; i<=MTTNX; i++)"$
 write "      fprintf(fps, ""%5.4g "", x[i]);"$
 write "    fprintf(fps, ""\n"");"$
 
+write "      time = time + DT;"$
 
 write "/* Inner integration loop */"$
 write "    for (k=1; k<=STEPFACTOR; k++)"$
 write "    {"$
 write "      for (i=1; i<=MTTNX; i++)"$
@@ -220,12 +225,10 @@
 write "      /* Set up system inputs */"$
 write "      $1_input(time);"$
 write "      $1_ode();"$
 write "    }"$
 
-
-write "      time = time + DT;"$
 write "  }"$
 
 %Write to files
 %Write to output to file
 write "/* Write to output file */"$