10
11
12
13
14
15
16
17
18
19
20
21
22
23
|
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
|
+
+
+
|
# Copyright (c) P.J.Gawthrop 1999
###############################################################
## Version control history
###############################################################
## $Id$
## $Log$
## Revision 1.8 2000/09/21 10:07:06 peterg
## Now correctely writes out the state into _odes.dat rep.
##
## Revision 1.7 2000/04/06 09:02:46 peterg
## Removed bash let commands as unsupported by sh
## Do arithmetic within awk - yuk
##
## Revision 1.6 2000/04/06 08:36:15 peterg
## *** empty log message ***
##
|
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
|
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
|
-
+
-
+
-
+
-
+
|
# Find Nx and Ny
Nx=`mtt_getsize $1 x` # States
Ny=`mtt_getsize $1 y` # Outputs
case $2 in
odeso)
declare -i First=1; declare -i Last=$First+$Ny
First=1; Last=`expr $First + $Ny`
;;
odes)
declare -i First=$Ny+2; declare -i Last=$First+$Nx
First=`expr $Ny + 2`; Last=`expr $First + $Nx`
;;
*)
declare -i First=1; declare -i Last=$Ny+1
First=1; Last=`expr $Ny + 1`
;;
esac
echo $Nx $Ny $First $Last
echo Creating $1_$2.dat
awk '{
## Extract the data
if (index($1,"#") != 1){
for (i=First;i<=Last;i++)
printf("%s ", $i);
printf("\n");
}
}' First=$First Last=$Last <$1_odes.dat2 > $1_$2.dat
|