Commit 8802bdf5 authored by Olivier Bonnefon's avatar Olivier Bonnefon
Browse files

remove unused files

parent aff9e448
int ndofs=0;
mesh Th0=readmesh("th0.msh");
fespace vh0(Th0,P1);
vh0 XX0=x;
vh0 YY0=y;
//SaveVec(XX0[],"X");
//SaveVec(YY0[],"Y");
ndofs+=vh0.ndof;
mesh Th1=readmesh("th1.msh");
fespace vh1(Th1,P1);
vh1 XX1=x;
vh1 YY1=y;
//SaveVec(XX1[],"X");
//SaveVec(YY1[],"Y");
ndofs+=vh1.ndof;
mesh Th2=readmesh("th2.msh");
fespace vh2(Th2,P1);
vh2 XX2=x;
vh2 YY2=y;
//SaveVec(XX2[],"X");
//SaveVec(YY2[],"Y");
ndofs+=vh2.ndof;
mesh Th3=readmesh("th3.msh");
fespace vh3(Th3,P1);
vh3 XX3=x;
vh3 YY3=y;
//SaveVec(XX3[],"X");
//SaveVec(YY3[],"Y");
ndofs+=vh3.ndof;
mesh Th4=readmesh("th4.msh");
fespace vh4(Th4,P1);
vh4 XX4=x;
vh4 YY4=y;
//SaveVec(XX4[],"X");
//SaveVec(YY4[],"Y");
ndofs+=vh4.ndof;
mesh Th5=readmesh("th5.msh");
fespace vh5(Th5,P1);
vh5 XX5=x;
vh5 YY5=y;
//SaveVec(XX5[],"X");
//SaveVec(YY5[],"Y");
ndofs+=vh5.ndof;
mesh Th6=readmesh("th6.msh");
fespace vh6(Th6,P1);
vh6 XX6=x;
vh6 YY6=y;
//SaveVec(XX6[],"X");
//SaveVec(YY6[],"Y");
ndofs+=vh6.ndof;
mesh Th7=readmesh("th7.msh");
fespace vh7(Th7,P1);
vh7 XX7=x;
vh7 YY7=y;
//SaveVec(XX7[],"X");
//SaveVec(YY7[],"Y");
ndofs+=vh7.ndof;
mesh Th8=readmesh("th8.msh");
fespace vh8(Th8,P1);
vh8 XX8=x;
vh8 YY8=y;
//SaveVec(XX8[],"X");
//SaveVec(YY8[],"Y");
ndofs+=vh8.ndof;
mesh Th9=readmesh("th9.msh");
fespace vh9(Th9,P1);
vh9 XX9=x;
vh9 YY9=y;
//SaveVec(XX9[],"X");
//SaveVec(YY9[],"Y");
ndofs+=vh9.ndof;
mesh Th10=readmesh("th10.msh");
fespace vh10(Th10,P1);
vh10 XX10=x;
vh10 YY10=y;
//SaveVec(XX10[],"X");
//SaveVec(YY10[],"Y");
ndofs+=vh10.ndof;
mesh Th11=readmesh("th11.msh");
fespace vh11(Th11,P1);
vh11 XX11=x;
vh11 YY11=y;
//SaveVec(XX11[],"X");
//SaveVec(YY11[],"Y");
ndofs+=vh11.ndof;
mesh Th12=readmesh("th12.msh");
fespace vh12(Th12,P1);
vh12 XX12=x;
vh12 YY12=y;
//SaveVec(XX12[],"X");
//SaveVec(YY12[],"Y");
ndofs+=vh12.ndof;
mesh Th13=readmesh("th13.msh");
fespace vh13(Th13,P1);
vh13 XX13=x;
vh13 YY13=y;
//SaveVec(XX13[],"X");
//SaveVec(YY13[],"Y");
ndofs+=vh13.ndof;
mesh Th14=readmesh("th14.msh");
fespace vh14(Th14,P1);
vh14 XX14=x;
vh14 YY14=y;
//SaveVec(XX14[],"X");
//SaveVec(YY14[],"Y");
ndofs+=vh14.ndof;
cout<<"ndofs total="<<ndofs<<endl;;
includepath += "../Tools1D/"
includepath += "../Recorder/"
includepath += "../examples1D/"
loadpath += "/home/olivierb/SOFT/freefem++-3.50/examples++-load/"
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);
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;
cout<<"rep0 \n"<<endl;
real[int] XX(20000);
cout<<"u2d0 loading\n"<<endl;
LoadVec(XX[],"./END/u2d0");
cout<<"u2d0 loaded\n"<<endl;
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+string(numA)+"/mesh"+"0");
LoadVec(january2d1[],DATASAFRAN+string(numA)+"/mesh"+"1");
LoadVec(january2d2[],DATASAFRAN+string(numA)+"/mesh"+"2");
LoadVec(january2d3[],DATASAFRAN+string(numA)+"/mesh"+"3");
LoadVec(january2d4[],DATASAFRAN+string(numA)+"/mesh"+"4");
LoadVec(january2d5[],DATASAFRAN+string(numA)+"/mesh"+"5");
LoadVec(january2d6[],DATASAFRAN+string(numA)+"/mesh"+"6");
LoadVec(january2d7[],DATASAFRAN+string(numA)+"/mesh"+"7");
LoadVec(january2d8[],DATASAFRAN+string(numA)+"/mesh"+"8");
LoadVec(january2d9[],DATASAFRAN+string(numA)+"/mesh"+"9");
LoadVec(january2d10[],DATASAFRAN+string(numA)+"/mesh"+"10");
LoadVec(january2d11[],DATASAFRAN+string(numA)+"/mesh"+"11");
LoadVec(january2d12[],DATASAFRAN+string(numA)+"/mesh"+"12");
LoadVec(january2d13[],DATASAFRAN+string(numA)+"/mesh"+"13");
LoadVec(january2d14[],DATASAFRAN+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);
This diff is collapsed.
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment