visuSimu.edp 5.06 KB
Newer Older
Olivier Bonnefon's avatar
Olivier Bonnefon committed
1
2
3
4
5
6
load "BinaryIO";
include "paramsTS.edp";
include "myMacro.edp";
include "defGeomGen.edp"
include "defDofsGen.edp"
int curStep=0;
7
int numRun=99;
Olivier Bonnefon's avatar
Olivier Bonnefon committed
8
9
10
11
int size;

{
  GetSizeVec(size,"../SAVE/u2d0__rep_"+string(numRun));
12
13
14
15
16
17
18
19
20
  cout<<size<<endl;
//  cout<<Vh0.ndof<<endl;
//  cout<<((364-58)*Vh0.ndof)<<endl;
//  cout<<"utm12d0.n "<<utm12d0.n<<endl;
  int saux=(364-58)*7*utm12d0.n;
  cout<<"saux="<<saux<<endl;
  // SetSizeVec(saux,"../SAVE/u2d0_"+"_rep_"+string(numRun));
//  GetSizeVec(size,"../SAVE/u2d0__rep_"+string(numRun));
//  cout<<size<<endl;
Olivier Bonnefon's avatar
Olivier Bonnefon committed
21
22
23
  real[int] u0(size);
  LoadVec(u0,"../SAVE/u2d0__rep_"+string(numRun));
  
24
25
  saux=(364-58)*7*utm12d1.n;
  GetSizeVec(size,"../SAVE/u2d1__rep_"+string(numRun));cout<<size<<" "<<saux<<endl;
Olivier Bonnefon's avatar
Olivier Bonnefon committed
26
27
28
  real[int] u1(size);
  LoadVec(u1,"../SAVE/u2d1__rep_"+string(numRun));
  
29
30
  saux=(364-58)*7*utm12d2.n;
  GetSizeVec(size,"../SAVE/u2d2__rep_"+string(numRun));cout<<size<<" "<<saux<<endl;
Olivier Bonnefon's avatar
Olivier Bonnefon committed
31
32
33
  real[int] u2(size);
  LoadVec(u2,"../SAVE/u2d2__rep_"+string(numRun));
  
34
35
  saux=(364-58)*7*utm12d3.n;
  GetSizeVec(size,"../SAVE/u2d3__rep_"+string(numRun));cout<<size<<" "<<saux<<endl;
Olivier Bonnefon's avatar
Olivier Bonnefon committed
36
37
38
  real[int] u3(size);
  LoadVec(u3,"../SAVE/u2d3__rep_"+string(numRun));
  
39
40
  saux=(364-58)*7*utm12d4.n;
  SetSizeVec(saux,"../SAVE/u2d4_"+"_rep_"+string(numRun));
Olivier Bonnefon's avatar
Olivier Bonnefon committed
41
42
43
44
  GetSizeVec(size,"../SAVE/u2d4__rep_"+string(numRun));
  real[int] u4(size);
  LoadVec(u4,"../SAVE/u2d4__rep_"+string(numRun));
  
45
46
  saux=(364-58)*7*utm12d5.n;
  SetSizeVec(saux,"../SAVE/u2d5_"+"_rep_"+string(numRun));
Olivier Bonnefon's avatar
Olivier Bonnefon committed
47
48
49
50
  GetSizeVec(size,"../SAVE/u2d5__rep_"+string(numRun));
  real[int] u5(size);
  LoadVec(u5,"../SAVE/u2d5__rep_"+string(numRun));
  
51
52
53
   saux=(364-58)*7*utm12d6.n;
  SetSizeVec(saux,"../SAVE/u2d6_"+"_rep_"+string(numRun));
 GetSizeVec(size,"../SAVE/u2d6__rep_"+string(numRun));
Olivier Bonnefon's avatar
Olivier Bonnefon committed
54
55
56
  real[int] u6(size);
  LoadVec(u6,"../SAVE/u2d6__rep_"+string(numRun));
  
57
58
  saux=(364-58)*7*utm12d7.n;
  SetSizeVec(saux,"../SAVE/u2d7_"+"_rep_"+string(numRun));
Olivier Bonnefon's avatar
Olivier Bonnefon committed
59
60
61
62
  GetSizeVec(size,"../SAVE/u2d7__rep_"+string(numRun));
  real[int] u7(size);
  LoadVec(u7,"../SAVE/u2d7__rep_"+string(numRun));
  
63
64
  saux=(364-58)*7*utm12d8.n;
  SetSizeVec(saux,"../SAVE/u2d8_"+"_rep_"+string(numRun));
Olivier Bonnefon's avatar
Olivier Bonnefon committed
65
66
67
68
  GetSizeVec(size,"../SAVE/u2d8__rep_"+string(numRun));
  real[int] u8(size);
  LoadVec(u8,"../SAVE/u2d8__rep_"+string(numRun));
  
69
70
  saux=(364-58)*7*utm12d9.n;
  SetSizeVec(saux,"../SAVE/u2d9_"+"_rep_"+string(numRun));
Olivier Bonnefon's avatar
Olivier Bonnefon committed
71
72
73
74
  GetSizeVec(size,"../SAVE/u2d9__rep_"+string(numRun));
  real[int] u9(size);
  LoadVec(u9,"../SAVE/u2d9__rep_"+string(numRun));
  
75
76
  saux=(364-58)*7*utm12d10.n;
  SetSizeVec(saux,"../SAVE/u2d10_"+"_rep_"+string(numRun));
Olivier Bonnefon's avatar
Olivier Bonnefon committed
77
78
79
80
  GetSizeVec(size,"../SAVE/u2d10__rep_"+string(numRun));
  real[int] u10(size);
  LoadVec(u10,"../SAVE/u2d10__rep_"+string(numRun));
  
81
82
83
   saux=(364-58)*7*utm12d11.n;
  SetSizeVec(saux,"../SAVE/u2d11_"+"_rep_"+string(numRun));
 GetSizeVec(size,"../SAVE/u2d11__rep_"+string(numRun));
Olivier Bonnefon's avatar
Olivier Bonnefon committed
84
85
86
  real[int] u11(size);
  LoadVec(u11,"../SAVE/u2d11__rep_"+string(numRun));
  
87
88
  saux=(364-58)*7*utm12d12.n;
  SetSizeVec(saux,"../SAVE/u2d12_"+"_rep_"+string(numRun));
Olivier Bonnefon's avatar
Olivier Bonnefon committed
89
90
91
92
  GetSizeVec(size,"../SAVE/u2d12__rep_"+string(numRun));
  real[int] u12(size);
  LoadVec(u12,"../SAVE/u2d12__rep_"+string(numRun));
  
93
94
  saux=(364-58)*7*utm12d13.n;
  SetSizeVec(saux,"../SAVE/u2d13_"+"_rep_"+string(numRun));
Olivier Bonnefon's avatar
Olivier Bonnefon committed
95
96
97
98
  GetSizeVec(size,"../SAVE/u2d13__rep_"+string(numRun));
  real[int] u13(size);
  LoadVec(u13,"../SAVE/u2d13__rep_"+string(numRun));
  
99
100
101
   saux=(364-58)*7*utm12d14.n;
  SetSizeVec(saux,"../SAVE/u2d14_"+"_rep_"+string(numRun));
 GetSizeVec(size,"../SAVE/u2d14__rep_"+string(numRun));
Olivier Bonnefon's avatar
Olivier Bonnefon committed
102
103
104
105
  real[int] u14(size);
  LoadVec(u14,"../SAVE/u2d14__rep_"+string(numRun));
  int nbCallSave=0;
  for (CURYEAR=FIRSTYEAR; CURYEAR<=LASTYEAR; CURYEAR++){
106
    int firstStep=59;
Olivier Bonnefon's avatar
Olivier Bonnefon committed
107
108
    if (CURYEAR==2011)
      firstStep=362-122;
109
    for ( curStep=firstStep;curStep<365;curStep++){
Olivier Bonnefon's avatar
Olivier Bonnefon committed
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
      utm12d0[]=u0(nbCallSave*utm12d0.n:(nbCallSave+1)*utm12d0.n);
      utm12d1[]=u1(nbCallSave*utm12d1.n:(nbCallSave+1)*utm12d1.n);
      utm12d2[]=u2(nbCallSave*utm12d2.n:(nbCallSave+1)*utm12d2.n);
      utm12d3[]=u3(nbCallSave*utm12d3.n:(nbCallSave+1)*utm12d3.n);
      utm12d4[]=u4(nbCallSave*utm12d4.n:(nbCallSave+1)*utm12d4.n);
      utm12d5[]=u5(nbCallSave*utm12d5.n:(nbCallSave+1)*utm12d5.n);
      utm12d6[]=u6(nbCallSave*utm12d6.n:(nbCallSave+1)*utm12d6.n);
      utm12d7[]=u7(nbCallSave*utm12d7.n:(nbCallSave+1)*utm12d7.n);
      utm12d8[]=u8(nbCallSave*utm12d8.n:(nbCallSave+1)*utm12d8.n);
      utm12d9[]=u9(nbCallSave*utm12d9.n:(nbCallSave+1)*utm12d9.n);
      utm12d10[]=u10(nbCallSave*utm12d10.n:(nbCallSave+1)*utm12d10.n);
      utm12d11[]=u11(nbCallSave*utm12d11.n:(nbCallSave+1)*utm12d11.n);
      utm12d12[]=u12(nbCallSave*utm12d12.n:(nbCallSave+1)*utm12d12.n);
      utm12d13[]=u13(nbCallSave*utm12d13.n:(nbCallSave+1)*utm12d13.n);
      utm12d14[]=u14(nbCallSave*utm12d14.n:(nbCallSave+1)*utm12d14.n);
      if (!(nbCallSave % 10))
        plot(utm12d0,utm12d1,utm12d2,utm12d3,utm12d4,utm12d5,utm12d6,utm12d7,utm12d8,utm12d9,utm12d10,utm12d11,utm12d12,utm12d13,utm12d14,
             cmm="rep "+string(numRun)+" year="+string(CURYEAR)+" step="+string(curStep),value=true,fill=true,wait=false);
      nbCallSave++;
    }
  }
}