Index: mttroot/mtt/bin/trans/dm2ir_m2dat ================================================================== --- mttroot/mtt/bin/trans/dm2ir_m2dat +++ mttroot/mtt/bin/trans/dm2ir_m2dat @@ -14,30 +14,46 @@ ############################################################### ## Version control history ############################################################### ## $Id$ ## $Log$ +## Revision 1.2 1996/08/11 10:40:00 peter +## Now handles e and NaN numbers. +## ## Revision 1.1 1996/08/10 14:15:03 peter ## Initial revision ## ############################################################### echo Creating $1_ir.dat + +if [ "$2" = "" ]; +then + PARAMS='t=[0:0.1:10]' + echo Using default parameter $PARAMS +else + PARAMS=$2; +fi + $MATRIX << EOF | awk '$1 !~ /[A-MO-Zb-cf-z]/ && NF>0 {print} ' >$1_ir.dat - tmax=2; tstep=tmax/100; - T = [0:tstep:tmax]'; + $PARAMS; + [n,m]=size(t); + if m>n + t=t'; + end; + [A,B,C,D,E] = $1_dm($1_numpar); - y = dm2ir(A,B,C,D,E,T); + y = dm2ir(A,B,C,D,E,t); [Ny,Nu]=size(D); NN = Ny*Nu; - xy = [T y(:,1)]; + xy = [t y(:,1)]; for i=2:NN - xy = [xy T y(:,i)]; + xy = [xy t y(:,i)]; end; xy EOF Index: mttroot/mtt/bin/trans/sm2sr_m ================================================================== --- mttroot/mtt/bin/trans/sm2sr_m +++ mttroot/mtt/bin/trans/sm2sr_m @@ -14,33 +14,51 @@ ############################################################### ## Version control history ############################################################### ## $Id$ ## $Log$ +## Revision 1.3 1996/08/11 10:38:30 peter +## Minor debugging changes. +## ## Revision 1.2 1996/08/11 10:06:50 peter ## Now allows e format and NaN numbers. ## ## Revision 1.1 1996/08/11 09:45:46 peter ## Initial revision ## ############################################################### echo Creating $1_sr.dat + + +if [ "$2" = "" ]; +then + PARAMS='t=[0:0.1:10]' + echo Using default parameter $PARAMS +else + PARAMS=$2; +fi + $MATRIX << EOF | awk '$1 !~ /[A-MO-Zb-cf-z]/ && NF>0 {print} ' >$1_sr.dat - tmax=10.0; tstep=tmax/100; - T = [0:tstep:tmax]'; + $PARAMS; + [n,m]=size(t); + if m>n + t=t'; + end; + + [A,B,C,D,E] = $1_dm($1_numpar); - y = dm2sr(A,B,C,D,E,T); + y = dm2sr(A,B,C,D,E,t); [Ny,Nu]=size(D); NN = Ny*Nu; - xy = [T y(:,1)]; + xy = [t y(:,1)]; for i=2:NN - xy = [xy T y(:,i)]; + xy = [xy t y(:,i)]; end; xy EOF