Index: mttroot/mtt/bin/trans/ode2odes_r2c ================================================================== --- mttroot/mtt/bin/trans/ode2odes_r2c +++ mttroot/mtt/bin/trans/ode2odes_r2c @@ -18,10 +18,13 @@ ############################################################### ## Version control history ############################################################### ## $Id$ ## $Log$ +## Revision 1.13 1998/02/24 13:34:45 peterg +## Back under RCS +## # Revision 1.12 1997/05/15 08:39:56 peterg # Don't initialise states - now done in numpar file. # # Revision 1.11 1997/05/12 16:00:54 peterg # Removed itime again., @@ -77,11 +80,11 @@ %Set up the number of argument variables to zero in case the user has forgotten MTTNVar := 0; %Read the symbolic parameters file -in "$1_sympar.r"; +%%in "$1_sympar.r"; ON BigFloat, NumVal; PRECISION 16; %Compatible with Matlab %OFF Nat; @@ -111,14 +114,15 @@ write "#define MTTNX2 ", MTTNx+2 $ write "#define MTTNY2 ", MTTNy+2 $ write "#define MTTNU2 ", MTTNu+2 $ write "#include "$ -write "#include ""$1_odes.h"" "$ +write "#include ""$1_simpar.h"" "$ write "#include ""$1_ode.c"" "$ write "#include ""$1_input.c"" "$ write "#include ""$1_numpar.c"" "$ +write "#include ""$1_state.c"" "$ %External (global) variable list write "#include ""$1_sympar.h"" "$ @@ -154,10 +158,14 @@ %Set up user-defined constants write "/* Set up user-defined constants */"$ write " $1_numpar();"$ +%Set up initial state +write "/* Set up initial state */"$ +write " $1_state();"$ + %Initialise main (Euler) integration loop write "/* Initialise main (Euler) integration loop */"$ write " time = 0.0;"$ write " dt = DT/STEPFACTOR;"$ %% write " for (i=1; i<=MTTNX; i++)"$