stockEggs.edp 3.71 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
29
30
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
67
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
110
111
112
113
114
115
116
117
118
119
120
121
122
load "BinaryIO";
include "declareEspace.edp";
cout<<"rep00 \n"<<endl;
vh0 u2d0;
vh1 u2d1;
vh2 u2d2;
vh3 u2d3;
vh4 u2d4;
vh5 u2d5;
vh6 u2d6;
vh7 u2d7;
vh8 u2d8;
vh9 u2d9;
vh10 u2d10;
vh11 u2d11;
vh12 u2d12;
vh13 u2d13;
vh14 u2d14;

LoadVec(u2d0[],"./END/u2d0");
LoadVec(u2d1[],"./END/u2d1");
LoadVec(u2d2[],"./END/u2d2");
LoadVec(u2d3[],"./END/u2d3");
LoadVec(u2d4[],"./END/u2d4");
LoadVec(u2d5[],"./END/u2d5");
LoadVec(u2d6[],"./END/u2d6");
LoadVec(u2d7[],"./END/u2d7");
LoadVec(u2d8[],"./END/u2d8");
LoadVec(u2d9[],"./END/u2d9");
LoadVec(u2d10[],"./END/u2d10");
LoadVec(u2d11[],"./END/u2d11");
LoadVec(u2d12[],"./END/u2d12");
LoadVec(u2d13[],"./END/u2d13");
LoadVec(u2d14[],"./END/u2d14");

plot(u2d0,u2d1,u2d2,u2d3,u2d4,u2d5,u2d6,u2d7,u2d8,u2d9,u2d10,u2d11,u2d12,u2d13,u2d14,cmm="u end",value=true,fill=true,wait=false);

vh0 january2d0;
vh1 january2d1;
vh2 january2d2;
vh3 january2d3;
vh4 january2d4;
vh5 january2d5;
vh6 january2d6;
vh7 january2d7;
vh8 january2d8;
vh9 january2d9;
vh10 january2d10;
vh11 january2d11;
vh12 january2d12;
vh13 january2d13;
vh14 january2d14;

vh0 egg2d0;
vh1 egg2d1;
vh2 egg2d2;
vh3 egg2d3;
vh4 egg2d4;
vh5 egg2d5;
vh6 egg2d6;
vh7 egg2d7;
vh8 egg2d8;
vh9 egg2d9;
vh10 egg2d10;
vh11 egg2d11;
vh12 egg2d12;
vh13 egg2d13;
vh14 egg2d14;

string DATASAFRAN;
DATASAFRAN="/home/olivierb/solvers/trunk/SandBox/FF/AedesAlbopictus/intraAnnuel/donnees_SAFRAN/AUTOROUTES/DATA_FF/";
int numA=2011;
LoadVec(january2d0[],DATASAFRAN+"JANUARY_AVERAGE/"+string(numA)+"/mesh"+"0");
LoadVec(january2d1[],DATASAFRAN+"JANUARY_AVERAGE/"+string(numA)+"/mesh"+"1");
LoadVec(january2d2[],DATASAFRAN+"JANUARY_AVERAGE/"+string(numA)+"/mesh"+"2");
LoadVec(january2d3[],DATASAFRAN+"JANUARY_AVERAGE/"+string(numA)+"/mesh"+"3");
LoadVec(january2d4[],DATASAFRAN+"JANUARY_AVERAGE/"+string(numA)+"/mesh"+"4");
LoadVec(january2d5[],DATASAFRAN+"JANUARY_AVERAGE/"+string(numA)+"/mesh"+"5");
LoadVec(january2d6[],DATASAFRAN+"JANUARY_AVERAGE/"+string(numA)+"/mesh"+"6");
LoadVec(january2d7[],DATASAFRAN+"JANUARY_AVERAGE/"+string(numA)+"/mesh"+"7");
LoadVec(january2d8[],DATASAFRAN+"JANUARY_AVERAGE/"+string(numA)+"/mesh"+"8");
LoadVec(january2d9[],DATASAFRAN+"JANUARY_AVERAGE/"+string(numA)+"/mesh"+"9");
LoadVec(january2d10[],DATASAFRAN+"JANUARY_AVERAGE/"+string(numA)+"/mesh"+"10");
LoadVec(january2d11[],DATASAFRAN+"JANUARY_AVERAGE/"+string(numA)+"/mesh"+"11");
LoadVec(january2d12[],DATASAFRAN+"JANUARY_AVERAGE/"+string(numA)+"/mesh"+"12");
LoadVec(january2d13[],DATASAFRAN+"JANUARY_AVERAGE/"+string(numA)+"/mesh"+"13");
LoadVec(january2d14[],DATASAFRAN+"JANUARY_AVERAGE/"+string(numA)+"/mesh"+"14");

plot(january2d0,january2d1,january2d2,january2d3,january2d4,january2d5,january2d6,january2d7,january2d8,january2d9,january2d10,january2d11,january2d12,january2d13,january2d14,cmm="T",value=true,fill=true,wait=false);




macro trunckEggs(VVHH,EGGINPUT,EGGOUTPUT,TJ){
  int i;
  int n=VVHH.ndof;
  for (i=0;i<n;i++){
    if (TJ[][i]>5){
      EGGOUTPUT[][i]=EGGINPUT[][i];
    }else{
      EGGOUTPUT[][i]=0;
    }
  }
}//

trunckEggs(vh0,u2d0,egg2d0,january2d0);
trunckEggs(vh1,u2d1,egg2d1,january2d1);
trunckEggs(vh2,u2d2,egg2d2,january2d2);
trunckEggs(vh3,u2d3,egg2d3,january2d3);
trunckEggs(vh4,u2d4,egg2d4,january2d4);
trunckEggs(vh5,u2d5,egg2d5,january2d5);
trunckEggs(vh6,u2d6,egg2d6,january2d6);
trunckEggs(vh7,u2d7,egg2d7,january2d7);
trunckEggs(vh8,u2d8,egg2d8,january2d8);
trunckEggs(vh9,u2d9,egg2d9,january2d9);
trunckEggs(vh10,u2d10,egg2d10,january2d10);
trunckEggs(vh11,u2d11,egg2d11,january2d11);
trunckEggs(vh12,u2d12,egg2d12,january2d12);
trunckEggs(vh13,u2d13,egg2d13,january2d13);
trunckEggs(vh14,u2d14,egg2d14,january2d14);

plot(egg2d0,egg2d1,egg2d2,egg2d3,egg2d4,egg2d5,egg2d6,egg2d7,egg2d8,egg2d9,egg2d10,egg2d11,egg2d12,egg2d13,egg2d14,cmm="egg",value=true,fill=true,wait=false);