16
17
18
19
20
21
22
23
24
25
26
27
28
29
|
# Copyright (c) P.J.Gawthrop 1997.
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.8 1997/05/10 08:12:23 peterg
## Put second argument into _input.
##
## Revision 1.7 1997/05/10 07:01:15 peterg
## _input called in outer loop only.
## time updated in outer loop only.
## Integer time (itime) introduced and updated in outer loop - maybe
|
>
>
>
|
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
|
# Copyright (c) P.J.Gawthrop 1997.
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.9 1997/05/10 09:54:34 peterg
## Moved _input call to after the inner integration loop.
##
## Revision 1.8 1997/05/10 08:12:23 peterg
## Put second argument into _input.
##
## Revision 1.7 1997/05/10 07:01:15 peterg
## _input called in outer loop only.
## time updated in outer loop only.
## Integer time (itime) introduced and updated in outer loop - maybe
|
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
|
write "/* Inner integration loop */"$
write " for (k=1; k<=STEPFACTOR; k++)"$
write " {"$
write " for (i=1; i<=MTTNX; i++)"$
write " x[i] = x[i] + dx[i]*dt;"$
write " $1_ode();"$
write " }"$
write " /* Set up system inputs */"$
write " $1_input(time,itime);"$
|
>
>
>
<
<
|
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
|
write "/* Inner integration loop */"$
write " for (k=1; k<=STEPFACTOR; k++)"$
write " {"$
write " for (i=1; i<=MTTNX; i++)"$
write " x[i] = x[i] + dx[i]*dt;"$
write " /* Set up system inputs */"$
write " $1_input(time,itime);"$
write " $1_ode();"$
write " }"$
write " time = time + DT;"$
write " itime = itime + 1;"$
write " }"$
%Write to files
%Write to output to file
|