Index: mttroot/mtt/bin/trans/make_ode2odes ================================================================== --- mttroot/mtt/bin/trans/make_ode2odes +++ mttroot/mtt/bin/trans/make_ode2odes @@ -1,10 +1,18 @@ #! /bin/sh ###################################### ##### Model Transformation Tools ##### ###################################### + +############################################################### +## Version control history +############################################################### +## $Id$ +## $Log$ +############################################################### + # Bourne shell script: make_ode2odes # Copyright (c) P.J.Gawthrop July 1998. @@ -40,17 +48,15 @@ mtt_write(t,MTTx,MTTy,$Nx,$Ny); # And write them #Integration loop for it = 1:iLast for i = 1:STEPFACTOR - MTTu = $1_input(MTTx,t);# Input - MTTx = $1_switch(MTTx); # Switches - MTTdx = $1_ode(MTTx,MTTu,t); - for j = 1:$Nx - MTTx(j) = MTTx(j) + MTTdx(j)*DDT; - end; - MTTy = $1_odeo(MTTx,MTTu,t); + [MTTu] = $1_input(MTTx,t);# Input + [MTTx] = $1_switch(MTTx); # Switches + [MTTdx] = $1_ode(MTTx,MTTu,t); + [MTTx] = mtt_update(MTTdx,MTTx,DDT,$Nx,METHOD); + [MTTy] = $1_odeo(MTTx,MTTu,t); t = t + DDT; end; mtt_write(t,MTTx,MTTy,$Nx,$Ny); end;