view.edp 12.5 KB
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
load "BinaryIO";
load "iovtk";

include "declareEspace.edp";
vh0 T0;
vh1 T1;
vh2 T2;
vh3 T3;
vh4 T4;
vh5 T5;
vh6 T6;
vh7 T7;
vh8 T8;
vh9 T9;
vh10 T10;
vh11 T11;
vh12 T12;
vh13 T13;
vh14 T14;
string DIRNAME;
DIRNAME="/home/olivierb/solvers/trunk/SandBox/FF/AedesAlbopictus/intraAnnuel/donnees_SAFRAN/AUTOROUTES/DATA_FF/";
string FIELD;
//FIELD = "T_BY_HEURE/";
//FIELD = "JANUARY_AVERAGE/";
//FIELD = "PHOTO_PERIOD_BY_DAY/";
//FIELD = "END_DIAPOSE_BY_DAY/";
//FIELD = "BEGIN_DIAPOSE_BY_DAY/";
//FIELD = "GDDTHRESHOLD_BY_DAY/";
29
30
FIELD = "GT_BY_DAY/";
//FIELD = "GP_BY_DAY/";
31
32
33
34
35
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
63
64
65
66
int step=1;
// LoadVec(T[],"T_BY_HEURE/meshT0");
// plot(T,wait=true,cmm="T0");
// LoadVec(T[],"T_BY_HEURE/meshT1");
// plot(T,wait=true,cmm="T1");
// LoadVec(T[],"T_BY_HEURE/meshT2");
// plot(T,wait=true,cmm="T2");
// LoadVec(T[],"T_BY_HEURE/meshT3");
// plot(T,wait=true,cmm="T3");
int firstY=2011;
int lastY=2016;
int numA=2011;
if (FIELD == "PHOTO_PERIOD_BY_DAY/"){
   for (int j=0;j<365;j++){
         LoadVec(T0[],DIRNAME +FIELD+"j0s"+string(j));
         LoadVec(T1[],DIRNAME +FIELD+"j1s"+string(j));
         LoadVec(T2[],DIRNAME +FIELD+"j2s"+string(j));
         LoadVec(T3[],DIRNAME +FIELD+"j3s"+string(j));
         LoadVec(T4[],DIRNAME +FIELD+"j4s"+string(j));
         LoadVec(T5[],DIRNAME +FIELD+"j5s"+string(j));
         LoadVec(T6[],DIRNAME +FIELD+"j6s"+string(j));
         LoadVec(T7[],DIRNAME +FIELD+"j7s"+string(j));
         LoadVec(T8[],DIRNAME +FIELD+"j8s"+string(j));
         LoadVec(T9[],DIRNAME +FIELD+"j9s"+string(j));
         LoadVec(T10[],DIRNAME +FIELD+"j10s"+string(j));
         LoadVec(T11[],DIRNAME +FIELD+"j11s"+string(j));
         LoadVec(T12[],DIRNAME +FIELD+"j12s"+string(j));
         LoadVec(T13[],DIRNAME +FIELD+"j13s"+string(j));
         LoadVec(T14[],DIRNAME +FIELD+"j14s"+string(j));
         plot(T0,T1,T2,T3,T4,T5,T6,T7,T8,T9,T10,T11,T12,T13,T14,cmm="photo period j "+string(j),value=true,fill=true,wait=false);
   }
   
}


if (FIELD == "GP_BY_DAY/"){
67
   for (int j=1;j<365;j++){
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
         LoadVec(T0[],DIRNAME +FIELD+"/gp0s"+string(j));
         LoadVec(T1[],DIRNAME +FIELD+"/gp1s"+string(j));
         LoadVec(T2[],DIRNAME +FIELD+"/gp2s"+string(j));
         LoadVec(T3[],DIRNAME +FIELD+"/gp3s"+string(j));
         LoadVec(T4[],DIRNAME +FIELD+"/gp4s"+string(j));
         LoadVec(T5[],DIRNAME +FIELD+"/gp5s"+string(j));
         LoadVec(T6[],DIRNAME +FIELD+"/gp6s"+string(j));
         LoadVec(T7[],DIRNAME +FIELD+"/gp7s"+string(j));
         LoadVec(T8[],DIRNAME +FIELD+"/gp8s"+string(j));
         LoadVec(T9[],DIRNAME +FIELD+"/gp9s"+string(j));
         LoadVec(T10[],DIRNAME +FIELD+"/gp10s"+string(j));
         LoadVec(T11[],DIRNAME +FIELD+"/gp11s"+string(j));
         LoadVec(T12[],DIRNAME +FIELD+"/gp12s"+string(j));
         LoadVec(T13[],DIRNAME +FIELD+"/gp13s"+string(j));
         LoadVec(T14[],DIRNAME +FIELD+"/gp14s"+string(j));
         plot(T0,T1,T2,T3,T4,T5,T6,T7,T8,T9,T10,T11,T12,T13,T14,cmm="GP j "+string(j),value=true,fill=true,wait=true);
   }
   
}

if (FIELD == "GT_BY_DAY/"){
   for (int j=180;j<365;j++){
         LoadVec(T0[],DIRNAME +FIELD+string(numA)+"/mesh0s"+string(j));
         LoadVec(T1[],DIRNAME +FIELD+string(numA)+"/mesh1s"+string(j));
         LoadVec(T2[],DIRNAME +FIELD+string(numA)+"/mesh2s"+string(j));
         LoadVec(T3[],DIRNAME +FIELD+string(numA)+"/mesh3s"+string(j));
         LoadVec(T4[],DIRNAME +FIELD+string(numA)+"/mesh4s"+string(j));
         LoadVec(T5[],DIRNAME +FIELD+string(numA)+"/mesh5s"+string(j));
         LoadVec(T6[],DIRNAME +FIELD+string(numA)+"/mesh6s"+string(j));
         LoadVec(T7[],DIRNAME +FIELD+string(numA)+"/mesh7s"+string(j));
         LoadVec(T8[],DIRNAME +FIELD+string(numA)+"/mesh8s"+string(j));
         LoadVec(T9[],DIRNAME +FIELD+string(numA)+"/mesh9s"+string(j));
         LoadVec(T10[],DIRNAME +FIELD+string(numA)+"/mesh10s"+string(j));
         LoadVec(T11[],DIRNAME +FIELD+string(numA)+"/mesh11s"+string(j));
         LoadVec(T12[],DIRNAME +FIELD+string(numA)+"/mesh12s"+string(j));
         LoadVec(T13[],DIRNAME +FIELD+string(numA)+"/mesh13s"+string(j));
         LoadVec(T14[],DIRNAME +FIELD+string(numA)+"/mesh14s"+string(j));
         plot(T0,T1,T2,T3,T4,T5,T6,T7,T8,T9,T10,T11,T12,T13,T14,cmm="GT j "+string(j),value=true,fill=true,wait=true);
   }
   
}

110
if (FIELD == "GDDTHRESHOLD_BY_DAY/"){
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
   for (int j=1;j<365-59;j++){
         LoadVec(T0[],DIRNAME +FIELD+string(numA)+"/mesh0s"+string(j));
         LoadVec(T1[],DIRNAME +FIELD+string(numA)+"/mesh1s"+string(j));
         LoadVec(T2[],DIRNAME +FIELD+string(numA)+"/mesh2s"+string(j));
         LoadVec(T3[],DIRNAME +FIELD+string(numA)+"/mesh3s"+string(j));
         LoadVec(T4[],DIRNAME +FIELD+string(numA)+"/mesh4s"+string(j));
         LoadVec(T5[],DIRNAME +FIELD+string(numA)+"/mesh5s"+string(j));
         LoadVec(T6[],DIRNAME +FIELD+string(numA)+"/mesh6s"+string(j));
         LoadVec(T7[],DIRNAME +FIELD+string(numA)+"/mesh7s"+string(j));
         LoadVec(T8[],DIRNAME +FIELD+string(numA)+"/mesh8s"+string(j));
         LoadVec(T9[],DIRNAME +FIELD+string(numA)+"/mesh9s"+string(j));
         LoadVec(T10[],DIRNAME +FIELD+string(numA)+"/mesh10s"+string(j));
         LoadVec(T11[],DIRNAME +FIELD+string(numA)+"/mesh11s"+string(j));
         LoadVec(T12[],DIRNAME +FIELD+string(numA)+"/mesh12s"+string(j));
         LoadVec(T13[],DIRNAME +FIELD+string(numA)+"/mesh13s"+string(j));
         LoadVec(T14[],DIRNAME +FIELD+string(numA)+"/mesh14s"+string(j));
127
         plot(T0,T1,T2,T3,T4,T5,T6,T7,T8,T9,T10,T11,T12,T13,T14,cmm="GDDTHRESHOLD_BY_DAY j "+string(j),value=true,fill=true,wait=false);
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
   }
   
}


if (FIELD == "END_DIAPOSE_BY_DAY/"){
   for (int j=0;j<365;j++){
         LoadVec(T0[],DIRNAME +FIELD+"end0s"+string(j));
         LoadVec(T1[],DIRNAME +FIELD+"end1s"+string(j));
         LoadVec(T2[],DIRNAME +FIELD+"end2s"+string(j));
         LoadVec(T3[],DIRNAME +FIELD+"end3s"+string(j));
         LoadVec(T4[],DIRNAME +FIELD+"end4s"+string(j));
         LoadVec(T5[],DIRNAME +FIELD+"end5s"+string(j));
         LoadVec(T6[],DIRNAME +FIELD+"end6s"+string(j));
         LoadVec(T7[],DIRNAME +FIELD+"end7s"+string(j));
         LoadVec(T8[],DIRNAME +FIELD+"end8s"+string(j));
         LoadVec(T9[],DIRNAME +FIELD+"end9s"+string(j));
         LoadVec(T10[],DIRNAME +FIELD+"end10s"+string(j));
         LoadVec(T11[],DIRNAME +FIELD+"end11s"+string(j));
         LoadVec(T12[],DIRNAME +FIELD+"end12s"+string(j));
         LoadVec(T13[],DIRNAME +FIELD+"end13s"+string(j));
         LoadVec(T14[],DIRNAME +FIELD+"end14s"+string(j));
         plot(T0,T1,T2,T3,T4,T5,T6,T7,T8,T9,T10,T11,T12,T13,T14,cmm="end diapos j "+string(j),value=true,fill=true,wait=true);
   }
   
}
if (FIELD == "BEGIN_DIAPOSE_BY_DAY/"){
   for (int j=170;j<365;j+=10){
         LoadVec(T0[],DIRNAME +FIELD+"begin0s"+string(j));
         LoadVec(T1[],DIRNAME +FIELD+"begin1s"+string(j));
         LoadVec(T2[],DIRNAME +FIELD+"begin2s"+string(j));
         LoadVec(T3[],DIRNAME +FIELD+"begin3s"+string(j));
         LoadVec(T4[],DIRNAME +FIELD+"begin4s"+string(j));
         LoadVec(T5[],DIRNAME +FIELD+"begin5s"+string(j));
         LoadVec(T6[],DIRNAME +FIELD+"begin6s"+string(j));
         LoadVec(T7[],DIRNAME +FIELD+"begin7s"+string(j));
         LoadVec(T8[],DIRNAME +FIELD+"begin8s"+string(j));
         LoadVec(T9[],DIRNAME +FIELD+"begin9s"+string(j));
         LoadVec(T10[],DIRNAME +FIELD+"begin10s"+string(j));
         LoadVec(T11[],DIRNAME +FIELD+"begin11s"+string(j));
         LoadVec(T12[],DIRNAME +FIELD+"begin12s"+string(j));
         LoadVec(T13[],DIRNAME +FIELD+"begin13s"+string(j));
         LoadVec(T14[],DIRNAME +FIELD+"begin14s"+string(j));
         plot(T0,T1,T2,T3,T4,T5,T6,T7,T8,T9,T10,T11,T12,T13,T14,cmm="begin diapos j "+string(j),value=true,fill=true,wait=true);
   }
   
}

if (FIELD != "JANUARY_AVERAGE/")
  for (int h=1;h<=365*24;h=h+step){
    LoadVec(T0[],DIRNAME +FIELD+string(numA)+"/mesh0s"+string(h));
    LoadVec(T1[],DIRNAME +FIELD+string(numA)+"/mesh1s"+string(h));
    LoadVec(T2[],DIRNAME +FIELD+string(numA)+"/mesh2s"+string(h));
    LoadVec(T3[],DIRNAME +FIELD+string(numA)+"/mesh3s"+string(h));
    LoadVec(T4[],DIRNAME +FIELD+string(numA)+"/mesh4s"+string(h));
    LoadVec(T5[],DIRNAME +FIELD+string(numA)+"/mesh5s"+string(h));
    LoadVec(T6[],DIRNAME +FIELD+string(numA)+"/mesh6s"+string(h));
    LoadVec(T7[],DIRNAME +FIELD+string(numA)+"/mesh7s"+string(h));
    LoadVec(T8[],DIRNAME +FIELD+string(numA)+"/mesh8s"+string(h));
    LoadVec(T9[],DIRNAME +FIELD+string(numA)+"/mesh9s"+string(h));
    LoadVec(T10[],DIRNAME +FIELD+string(numA)+"/mesh10s"+string(h));
    LoadVec(T11[],DIRNAME +FIELD+string(numA)+"/mesh11s"+string(h));
    LoadVec(T12[],DIRNAME +FIELD+string(numA)+"/mesh12s"+string(h));
    LoadVec(T13[],DIRNAME +FIELD+string(numA)+"/mesh13s"+string(h));
    LoadVec(T14[],DIRNAME +FIELD+string(numA)+"/mesh14s"+string(h));
    plot(T0,T1,T2,T3,T4,T5,T6,T7,T8,T9,T10,T11,T12,T13,T14,cmm="annee "+string(numA)+" "+string(h/24)+" jours "+string(h%24)+ " h ",value=true,fill=true,wait=false);
//  cout<<"T min="<<T[].min<<" Tmax="<<T[].max<<endl;
//  savevtk("T_BY_HEURE/meshT"+string(h)+".vtk",ThF,T,dataname="UU");
  }
for (numA=firstY; numA<lastY+1;numA++){
if (FIELD == "JANUARY_AVERAGE/"){
  cout<<DIRNAME +FIELD+string(numA)+"/mesh0"<<endl;
    LoadVec(T0[],DIRNAME +FIELD+string(numA)+"/mesh0");
    LoadVec(T1[],DIRNAME +FIELD+string(numA)+"/mesh1");
    LoadVec(T2[],DIRNAME +FIELD+string(numA)+"/mesh2");
    LoadVec(T3[],DIRNAME +FIELD+string(numA)+"/mesh3");
    LoadVec(T4[],DIRNAME +FIELD+string(numA)+"/mesh4");
    LoadVec(T5[],DIRNAME +FIELD+string(numA)+"/mesh5");
    LoadVec(T6[],DIRNAME +FIELD+string(numA)+"/mesh6");
    LoadVec(T7[],DIRNAME +FIELD+string(numA)+"/mesh7");
    LoadVec(T8[],DIRNAME +FIELD+string(numA)+"/mesh8");
    LoadVec(T9[],DIRNAME +FIELD+string(numA)+"/mesh9");
    LoadVec(T10[],DIRNAME +FIELD+string(numA)+"/mesh10");
    LoadVec(T11[],DIRNAME +FIELD+string(numA)+"/mesh11");
    LoadVec(T12[],DIRNAME +FIELD+string(numA)+"/mesh12");
    LoadVec(T13[],DIRNAME +FIELD+string(numA)+"/mesh13");
    LoadVec(T14[],DIRNAME +FIELD+string(numA)+"/mesh14");
    plot(T0,T1,T2,T3,T4,T5,T6,T7,T8,T9,T10,T11,T12,T13,T14,cmm="annee "+string(numA),value=true,fill=true,wait=false,ps="average_january_"+string(numA)+".eps");
}
}
exit(0);
int cmpJ=152;
int h;
if (1)
for ( h=8760-3666;h<8760;h++){
  if (h%24 == 0)
    cmpJ++;
  LoadVec(T0[],"DATA_2002/T_BY_HEURE/meshT0s"+string(h));
  LoadVec(T1[],"DATA_2002/T_BY_HEURE/meshT1s"+string(h));
  LoadVec(T2[],"DATA_2002/T_BY_HEURE/meshT2s"+string(h));
  LoadVec(T3[],"DATA_2002/T_BY_HEURE/meshT3s"+string(h));
  LoadVec(T4[],"DATA_2002/T_BY_HEURE/meshT4s"+string(h));
  LoadVec(T5[],"DATA_2002/T_BY_HEURE/meshT5s"+string(h));
  LoadVec(T6[],"DATA_2002/T_BY_HEURE/meshT6s"+string(h));
  LoadVec(T7[],"DATA_2002/T_BY_HEURE/meshT7s"+string(h));
  LoadVec(T8[],"DATA_2002/T_BY_HEURE/meshT8s"+string(h));
  LoadVec(T9[],"DATA_2002/T_BY_HEURE/meshT9s"+string(h));
  LoadVec(T10[],"DATA_2002/T_BY_HEURE/meshT10s"+string(h));
  LoadVec(T11[],"DATA_2002/T_BY_HEURE/meshT11s"+string(h));
  LoadVec(T12[],"DATA_2002/T_BY_HEURE/meshT12s"+string(h));
  LoadVec(T13[],"DATA_2002/T_BY_HEURE/meshT13s"+string(h));
  LoadVec(T14[],"DATA_2002/T_BY_HEURE/meshT14s"+string(h));
  string leg="jour="+string(cmpJ)+" h="+string(h%24);
  plot(T0,T1,T2,T3,T4,T5,T6,T7,T8,T9,T10,T11,T12,T13,T14,cmm=leg,value=true,fill=true,wait=true);
//  cout<<"T min="<<T[].min<<" Tmax="<<T[].max<<endl;
//  savevtk("T_BY_HEURE/meshT"+string(h)+".vtk",ThF,T,dataname="UU");
}
exit(0);
cmpJ=-1;
for (h=(11*30+15)*24;h<8760;h++){
  if (h%24 == 0)
    cmpJ++;
  LoadVec(T0[],"DATA_2015/T_BY_HEURE/meshT0s"+string(h));
  LoadVec(T1[],"DATA_2015/T_BY_HEURE/meshT1s"+string(h));
  LoadVec(T2[],"DATA_2015/T_BY_HEURE/meshT2s"+string(h));
  LoadVec(T3[],"DATA_2015/T_BY_HEURE/meshT3s"+string(h));
  LoadVec(T4[],"DATA_2015/T_BY_HEURE/meshT4s"+string(h));
  LoadVec(T5[],"DATA_2015/T_BY_HEURE/meshT5s"+string(h));
  LoadVec(T6[],"DATA_2015/T_BY_HEURE/meshT6s"+string(h));
  LoadVec(T7[],"DATA_2015/T_BY_HEURE/meshT7s"+string(h));
  LoadVec(T8[],"DATA_2015/T_BY_HEURE/meshT8s"+string(h));
  LoadVec(T9[],"DATA_2015/T_BY_HEURE/meshT9s"+string(h));
  LoadVec(T10[],"DATA_2015/T_BY_HEURE/meshT10s"+string(h));
  LoadVec(T11[],"DATA_2015/T_BY_HEURE/meshT11s"+string(h));
  LoadVec(T12[],"DATA_2015/T_BY_HEURE/meshT12s"+string(h));
  LoadVec(T13[],"DATA_2015/T_BY_HEURE/meshT13s"+string(h));
  LoadVec(T14[],"DATA_2015/T_BY_HEURE/meshT14s"+string(h));
  string leg="jour="+string(cmpJ)+" h="+string(h%24);
  plot(T0,T1,T2,T3,T4,T5,T6,T7,T8,T9,T10,T11,T12,T13,T14,cmm=leg,value=true,fill=true,wait=true);
//  cout<<"T min="<<T[].min<<" Tmax="<<T[].max<<endl;
//  savevtk("T_BY_HEURE/meshT"+string(h)+".vtk",ThF,T,dataname="UU");
}