29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
|
# PARAMS="$PARAMS ;"
$MATRIX << EOF > sm2ir_m.log 2>mtt_error.txt
% Read in the numeric system parameters
mttpar = $1_numpar;
[nx,ny,nu,nz,nyz] = $1_def;
t=0; %Just in case it appears in the parameter list.
$PARAMS
if exist('T')==0
T=[0:0.1:10]
end;
|
>
>
>
<
<
|
<
|
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
|
# PARAMS="$PARAMS ;"
$MATRIX << EOF > sm2ir_m.log 2>mtt_error.txt
% Read in the numeric system parameters
mttpar = $1_numpar;
% Read in the simulation parameters
simpar = $1_simpar;
[nx,ny,nu,nz,nyz] = $1_def;
t=0; %Just in case it appears in the parameter list.
%Defaults
T=[0:simpar.dt:simpar.last]';
if exist('u0')==0
u0 = ones(nu,1);
end;
if exist('x0')==0
x0 = zeros(nx,1);
|