Commit c70ea879 authored by Olivier Bonnefon's avatar Olivier Bonnefon
Browse files

add missing files used?

parent 3977e225
load "BinaryIO";
include "../paramsTS.edp";
include "../myMacro.edp";
include "../defGeomGen.edp"
include "../defDofsGen.edp"
include "../buildTobs.edp"
include "defTigerVar.edp"
include "defSimuParamsGen.edp";
real alphaw0,alphaw1;
string SRCDIR="/home/olivierb/aedes/AEDESINTER_AUTOROUTES_ADAPTIV/";
string MODE="DOVISU";
string LOADDIR="/home/olivierb/aedes/AEDESINTER_AUTOROUTES_ADAPTIV/SAVE_NOASSIM/";
//string MODE="COMPUTEBS";
include "defVdepMax.edp"
include "initVdepMax.edp"
include "../postTraitementLoad.edp"
if (ARGV.n>2){
FIRSTYEAR=atoi(ARGV[2]);
}
if (ARGV.n>3){
NBPREDYEAR=atoi(ARGV[3]);
LASTYEAR=FIRSTYEAR+NBPREDYEAR;
}
int curStep=0;
int numRun=0;
int size;
//visu ic
include "defOutputGen.edp";
include "defParamsGen.edp";
include "mat1D2Duseful.edp";
//include "postTraitementLoad.edp";
//include "inputParam.edp";
include "defVarfGen.edp";
FIRSTYEAR=2009;
NBPREDYEAR=1;
cout<<LOADDIR+"/w0s2d0_"+string(numRun)+string(FIRSTYEAR+NBPREDYEAR-2)<<endl;
LoadVec(w0s2d0[],LOADDIR+"/w0s2d0_"+string(numRun)+"_"+string(FIRSTYEAR+NBPREDYEAR-2));
LoadVec(w0s2d1[],LOADDIR+"/w0s2d1_"+string(numRun)+"_"+string(FIRSTYEAR+NBPREDYEAR-2));
LoadVec(w0s2d2[],LOADDIR+"/w0s2d2_"+string(numRun)+"_"+string(FIRSTYEAR+NBPREDYEAR-2));
LoadVec(w0s2d3[],LOADDIR+"/w0s2d3_"+string(numRun)+"_"+string(FIRSTYEAR+NBPREDYEAR-2));
LoadVec(w0s2d4[],LOADDIR+"/w0s2d4_"+string(numRun)+"_"+string(FIRSTYEAR+NBPREDYEAR-2));
LoadVec(w0s2d5[],LOADDIR+"/w0s2d5_"+string(numRun)+"_"+string(FIRSTYEAR+NBPREDYEAR-2));
LoadVec(w0s2d6[],LOADDIR+"/w0s2d6_"+string(numRun)+"_"+string(FIRSTYEAR+NBPREDYEAR-2));
LoadVec(w0s2d7[],LOADDIR+"/w0s2d7_"+string(numRun)+"_"+string(FIRSTYEAR+NBPREDYEAR-2));
LoadVec(w0s2d8[],LOADDIR+"/w0s2d8_"+string(numRun)+"_"+string(FIRSTYEAR+NBPREDYEAR-2));
LoadVec(w0s2d9[],LOADDIR+"/w0s2d9_"+string(numRun)+"_"+string(FIRSTYEAR+NBPREDYEAR-2));
LoadVec(w0s2d10[],LOADDIR+"/w0s2d10_"+string(numRun)+"_"+string(FIRSTYEAR+NBPREDYEAR-2));
LoadVec(w0s2d11[],LOADDIR+"/w0s2d11_"+string(numRun)+"_"+string(FIRSTYEAR+NBPREDYEAR-2));
LoadVec(w0s2d12[],LOADDIR+"/w0s2d12_"+string(numRun)+"_"+string(FIRSTYEAR+NBPREDYEAR-2));
LoadVec(w0s2d13[],LOADDIR+"/w0s2d13_"+string(numRun)+"_"+string(FIRSTYEAR+NBPREDYEAR-2));
LoadVec(w0s2d14[],LOADDIR+"/w0s2d14_"+string(numRun)+"_"+string(FIRSTYEAR+NBPREDYEAR-2));
string lcmm=LOADDIR+"/w0s2d0_"+string(numRun)+"_"+string(FIRSTYEAR+NBPREDYEAR-2);
plot(w0s2d0,w0s2d1,w0s2d2,w0s2d3,w0s2d4,w0s2d5,w0s2d6,w0s2d7,w0s2d8,w0s2d9,w0s2d10,w0s2d11,w0s2d12,w0s2d13,w0s2d14,cmm=lcmm,ps=lcmm+".eps",value=true,fill=true,wait=true);
LOADDIR="/home/olivierb/aedes/AEDESINTER_AUTOROUTES_ADAPTIV/SAVE/";
cout<<LOADDIR+"/w0s2d0_"+string(numRun)+"_"+string(NBPREDYEAR-1)+"_"+string(FIRSTYEAR)<<endl;
LoadVec(w0s2d0[],LOADDIR+"/w0s2d0_"+string(numRun)+"_"+string(NBPREDYEAR-1)+"_"+string(FIRSTYEAR));cout<<"done"<<endl;
LoadVec(w0s2d1[],LOADDIR+"/w0s2d1_"+string(numRun)+"_"+string(NBPREDYEAR-1)+"_"+string(FIRSTYEAR));
LoadVec(w0s2d2[],LOADDIR+"/w0s2d2_"+string(numRun)+"_"+string(NBPREDYEAR-1)+"_"+string(FIRSTYEAR));
LoadVec(w0s2d3[],LOADDIR+"/w0s2d3_"+string(numRun)+"_"+string(NBPREDYEAR-1)+"_"+string(FIRSTYEAR));
LoadVec(w0s2d4[],LOADDIR+"/w0s2d4_"+string(numRun)+"_"+string(NBPREDYEAR-1)+"_"+string(FIRSTYEAR));
LoadVec(w0s2d5[],LOADDIR+"/w0s2d5_"+string(numRun)+"_"+string(NBPREDYEAR-1)+"_"+string(FIRSTYEAR));
LoadVec(w0s2d6[],LOADDIR+"/w0s2d6_"+string(numRun)+"_"+string(NBPREDYEAR-1)+"_"+string(FIRSTYEAR));
LoadVec(w0s2d7[],LOADDIR+"/w0s2d7_"+string(numRun)+"_"+string(NBPREDYEAR-1)+"_"+string(FIRSTYEAR));
LoadVec(w0s2d8[],LOADDIR+"/w0s2d8_"+string(numRun)+"_"+string(NBPREDYEAR-1)+"_"+string(FIRSTYEAR));
LoadVec(w0s2d9[],LOADDIR+"/w0s2d9_"+string(numRun)+"_"+string(NBPREDYEAR-1)+"_"+string(FIRSTYEAR));
LoadVec(w0s2d10[],LOADDIR+"/w0s2d10_"+string(numRun)+"_"+string(NBPREDYEAR-1)+"_"+string(FIRSTYEAR));
LoadVec(w0s2d11[],LOADDIR+"/w0s2d11_"+string(numRun)+"_"+string(NBPREDYEAR-1)+"_"+string(FIRSTYEAR));
LoadVec(w0s2d12[],LOADDIR+"/w0s2d12_"+string(numRun)+"_"+string(NBPREDYEAR-1)+"_"+string(FIRSTYEAR));
LoadVec(w0s2d13[],LOADDIR+"/w0s2d13_"+string(numRun)+"_"+string(NBPREDYEAR-1)+"_"+string(FIRSTYEAR));
LoadVec(w0s2d14[],LOADDIR+"/w0s2d14_"+string(numRun)+"_"+string(NBPREDYEAR-1)+"_"+string(FIRSTYEAR));
lcmm=LOADDIR+"/w0s2d14_"+string(numRun)+"_"+string(NBPREDYEAR-1)+"_"+string(FIRSTYEAR);
plot(w0s2d0,w0s2d1,w0s2d2,w0s2d3,w0s2d4,w0s2d5,w0s2d6,w0s2d7,w0s2d8,w0s2d9,w0s2d10,w0s2d11,w0s2d12,w0s2d13,w0s2d14,cmm=lcmm,ps=lcmm+".eps",value=true,fill=true,wait=true);
FreeFem++-nw loopOnSimuNoAssim.edp 2009 0
FreeFem++-nw loopOnSimuNoAssim.edp 2010 0
FreeFem++-nw loopOnSimuNoAssim.edp 2011 0
FreeFem++-nw loopOnSimuNoAssim.edp 2012 0
FreeFem++-nw loopOnSimuNoAssim.edp 2013 0
FreeFem++-nw loopOnSimuNoAssim.edp 2014 0
FreeFem++-nw loopOnSimuNoAssim.edp 2015 0
FreeFem++-nw loopOnSimuNoAssim.edp 2016 0
FreeFem++-nw loopOnSimuNoAssim.edp 2017 0
FreeFem++-nw loopOnSimuNoAssim.edp 2018 0
FreeFem++-nw loopOnSimuNoAssim.edp 2019 0
load "BinaryIO";
include "../paramsTS.edp";
include "../myMacro.edp";
include "../defGeomGen.edp"
include "../defDofsGen.edp"
include "../buildTobs.edp"
string SRCDIR="/home/olivierb/aedes/AEDESINTER_AUTOROUTES_ADAPTIV/";
if (ARGV.n>2){
FIRSTYEAR=atoi(ARGV[2]);
}
if (ARGV.n>3){
NBPREDYEAR=atoi(ARGV[3]);
LASTYEAR=FIRSTYEAR+NBPREDYEAR;
}
include "defVdepMax.edp"
include "initVdepMax.edp"
include "../postTraitementLoad.edp"
int curStep=0;
int numRun=0;
int size;
real BSScore;
int indexYear;
int cmpData;
real thetaAssim=1;
{
CURYEAR=FIRSTYEAR;
include "initVdepMax.edp";
include "updateVdepFromTobs.edp";
include "updateVdepMax.edp";
for (int i=0;i<114;i++)
cout<<" vdepMax[i]="<< vdepMax[i]<<endl;
}
CURYEAR=LASTYEAR;
real coefP=1;
include "defP.edp";
include "computeBS.edp";
cout<<BSScore<<endl;
int[int] useDep(114);
for (int i=0;i<114;i++)
useDep[i]=1;
useDep[3]=0;
useDep[23]=0;
useDep[25]=0;
useDep[26]=0;
useDep[28]=0;
useDep[36]=0;
useDep[46]=0;
useDep[48]=0;
useDep[52]=0;
useDep[53]=0;
useDep[54]=0;
useDep[61]=0;
useDep[63]=0;
useDep[71]=0;
useDep[74]=0;
useDep[92]=0;
useDep[105]=0;
useDep[107]=0;
useDep[108]=0;
useDep[109]=0;
useDep[110]=0;
useDep[111]=0;
useDep[113]=0;
string Kstring="";
if (NBPREDYEAR>0)
Kstring=string(NBPREDYEAR)+"_";
//string SAVEDIR="../SAVE/";
cout<<SAVEDIR+"u2d0_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun)<<endl;
GetSizeVec(size,SAVEDIR+"u2d0_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u0(size);
//cout<<"size="<<size<<" "<<Vh0.ndof<<endl;
nbDays=size/Vh0.ndof;
LoadVec(u0,SAVEDIR+"u2d0_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,SAVEDIR+"u2d1_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u1(size);
LoadVec(u1,SAVEDIR+"u2d1_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,SAVEDIR+"u2d2_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u2(size);
LoadVec(u2,SAVEDIR+"u2d2_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,SAVEDIR+"u2d3_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u3(size);
LoadVec(u3,SAVEDIR+"u2d3_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,SAVEDIR+"u2d4_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u4(size);
LoadVec(u4,SAVEDIR+"u2d4_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,SAVEDIR+"u2d5_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u5(size);
LoadVec(u5,SAVEDIR+"u2d5_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,SAVEDIR+"u2d6_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u6(size);
LoadVec(u6,SAVEDIR+"u2d6_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,SAVEDIR+"u2d7_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u7(size);
LoadVec(u7,SAVEDIR+"u2d7_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,SAVEDIR+"u2d8_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u8(size);
LoadVec(u8,SAVEDIR+"u2d8_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,SAVEDIR+"u2d9_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u9(size);
LoadVec(u9,SAVEDIR+"u2d9_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,SAVEDIR+"u2d10_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u10(size);
LoadVec(u10,SAVEDIR+"u2d10_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,SAVEDIR+"u2d11_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u11(size);
LoadVec(u11,SAVEDIR+"u2d11_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,SAVEDIR+"u2d12_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u12(size);
LoadVec(u12,SAVEDIR+"u2d12_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,SAVEDIR+"u2d13_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u13(size);
LoadVec(u13,SAVEDIR+"u2d13_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,SAVEDIR+"u2d14_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u14(size);
LoadVec(u14,SAVEDIR+"u2d14_"+Kstring+string(CURYEAR)+"_rep_"+string(numRun));
load "BinaryIO";
include "../paramsTS.edp";
include "../myMacro.edp";
include "../defGeomGen.edp"
include "../defDofsGen.edp"
include "../buildTobs.edp"
include "defTigerVar.edp"
include "defSimuParamsGen.edp";
real alphaw0,alphaw1;
string SRCDIR="/home/olivierb/aedes/AEDESINTER_AUTOROUTES_ADAPTIV/";
//string MODE="DOVISU";
string MODE="COMPUTEBS";
string SAVEDIR="../SAVE/";
include "defVdepMax.edp"
include "initVdepMax.edp"
include "../postTraitementLoad.edp"
if (ARGV.n>2){
FIRSTYEAR=atoi(ARGV[2]);
}
if (ARGV.n>3){
NBPREDYEAR=atoi(ARGV[3]);
LASTYEAR=FIRSTYEAR+NBPREDYEAR-1;
}
CURYEAR=FIRSTYEAR;
int curStep=0;
int numRun=0;
int size;
int nbDays;
include "loadState.edp";
int nbCallSave=0;
int firstStep=364-nbDays;
for ( curStep=firstStep;curStep<364;curStep++){
u2d0[]=u0(nbCallSave*u2d0.n:(nbCallSave+1)*u2d0.n);
u2d1[]=u1(nbCallSave*u2d1.n:(nbCallSave+1)*u2d1.n);
u2d2[]=u2(nbCallSave*u2d2.n:(nbCallSave+1)*u2d2.n);
u2d3[]=u3(nbCallSave*u2d3.n:(nbCallSave+1)*u2d3.n);
u2d4[]=u4(nbCallSave*u2d4.n:(nbCallSave+1)*u2d4.n);
u2d5[]=u5(nbCallSave*u2d5.n:(nbCallSave+1)*u2d5.n);
u2d6[]=u6(nbCallSave*u2d6.n:(nbCallSave+1)*u2d6.n);
u2d7[]=u7(nbCallSave*u2d7.n:(nbCallSave+1)*u2d7.n);
u2d8[]=u8(nbCallSave*u2d8.n:(nbCallSave+1)*u2d8.n);
u2d9[]=u9(nbCallSave*u2d9.n:(nbCallSave+1)*u2d9.n);
u2d10[]=u10(nbCallSave*u2d10.n:(nbCallSave+1)*u2d10.n);
u2d11[]=u11(nbCallSave*u2d11.n:(nbCallSave+1)*u2d11.n);
u2d12[]=u12(nbCallSave*u2d12.n:(nbCallSave+1)*u2d12.n);
u2d13[]=u13(nbCallSave*u2d13.n:(nbCallSave+1)*u2d13.n);
u2d14[]=u14(nbCallSave*u2d14.n:(nbCallSave+1)*u2d14.n);
if (MODE=="COMPUTEBS"){
include "stateToDep.edp";
include "updateVdepMax.edp";
}
if (MODE=="DOVISU")
if (!(nbCallSave % 10))
plot(u2d0,u2d1,u2d2,u2d3,u2d4,u2d5,u2d6,u2d7,u2d8,u2d9,u2d10,u2d11,u2d12,u2d13,u2d14,
cmm="rep "+string(numRun)+" year="+string(CURYEAR)+" step="+string(curStep),value=true,fill=true,wait=false);
nbCallSave++;
}
ofstream bsFile("resBS.txt",append);
if (MODE=="COMPUTEBS"){
real BSScore=0;
int indexYear;
int cmpData=0;
real coefP=2;
include "defP.edp";
include "computeBS.edp";
cout<<BSScore<<endl;
bsFile<<FIRSTYEAR<<" "<<NBPREDYEAR<<" "<<BSScore<<" ";
BSScore=0;
include "initVdepMax.edp";
curStep=0;
numRun=0;
real thetaAssim=1;
{
CURYEAR=FIRSTYEAR-1;
include "initVdepMax.edp";
include "updateVdepFromTobs.edp";
include "updateVdepMax.edp";
for (int i=0;i<114;i++)
cout<<" vdepMax[i]="<< vdepMax[i]<<endl;
CURYEAR=LASTYEAR;
coefP=1;
include "defP.edp";
include "computeBS.edp";
bsFile<<BSScore<<endl;
}
}
load "BinaryIO";
include "../paramsTS.edp";
include "../myMacro.edp";
include "../defGeomGen.edp"
include "../defDofsGen.edp"
include "../buildTobs.edp"
include "defTigerVar.edp"
include "defSimuParamsGen.edp";
real alphaw0,alphaw1;
string SRCDIR="/home/olivierb/aedes/AEDESINTER_AUTOROUTES_ADAPTIV/";
//string MODE="DOVISU";
string MODE="COMPUTEBS";
string SAVEDIR="../SAVE_NOASSIM/";
include "defVdepMax.edp"
include "initVdepMax.edp"
include "../postTraitementLoad.edp"
if (ARGV.n>2){
FIRSTYEAR=atoi(ARGV[2]);
}
if (ARGV.n>3){
NBPREDYEAR=atoi(ARGV[3]);
LASTYEAR=FIRSTYEAR+NBPREDYEAR;
}
CURYEAR=FIRSTYEAR;
int curStep=0;
int numRun=0;
int size;
int nbDays;
include "loadState.edp";
int nbCallSave=0;
int firstStep=364-nbDays;
for ( curStep=firstStep;curStep<364;curStep++){
u2d0[]=u0(nbCallSave*u2d0.n:(nbCallSave+1)*u2d0.n);
u2d1[]=u1(nbCallSave*u2d1.n:(nbCallSave+1)*u2d1.n);
u2d2[]=u2(nbCallSave*u2d2.n:(nbCallSave+1)*u2d2.n);
u2d3[]=u3(nbCallSave*u2d3.n:(nbCallSave+1)*u2d3.n);
u2d4[]=u4(nbCallSave*u2d4.n:(nbCallSave+1)*u2d4.n);
u2d5[]=u5(nbCallSave*u2d5.n:(nbCallSave+1)*u2d5.n);
u2d6[]=u6(nbCallSave*u2d6.n:(nbCallSave+1)*u2d6.n);
u2d7[]=u7(nbCallSave*u2d7.n:(nbCallSave+1)*u2d7.n);
u2d8[]=u8(nbCallSave*u2d8.n:(nbCallSave+1)*u2d8.n);
u2d9[]=u9(nbCallSave*u2d9.n:(nbCallSave+1)*u2d9.n);
u2d10[]=u10(nbCallSave*u2d10.n:(nbCallSave+1)*u2d10.n);
u2d11[]=u11(nbCallSave*u2d11.n:(nbCallSave+1)*u2d11.n);
u2d12[]=u12(nbCallSave*u2d12.n:(nbCallSave+1)*u2d12.n);
u2d13[]=u13(nbCallSave*u2d13.n:(nbCallSave+1)*u2d13.n);
u2d14[]=u14(nbCallSave*u2d14.n:(nbCallSave+1)*u2d14.n);
if (MODE=="COMPUTEBS"){
include "stateToDep.edp";
include "updateVdepMax.edp";
}
if (MODE=="DOVISU")
if (!(nbCallSave % 10))
plot(u2d0,u2d1,u2d2,u2d3,u2d4,u2d5,u2d6,u2d7,u2d8,u2d9,u2d10,u2d11,u2d12,u2d13,u2d14,
cmm="rep "+string(numRun)+" year="+string(CURYEAR)+" step="+string(curStep),value=true,fill=true,wait=false);
nbCallSave++;
}
ofstream bsFile("resBSNoAssim.txt",append);
if (MODE=="COMPUTEBS"){
real BSScore=0;
int indexYear=LASTYEAR-2003;
int cmpData=0;
real coefP=2;
include "defP.edp";
include "computeBS.edp";
cout<<BSScore<<endl;
bsFile<<FIRSTYEAR<<" "<<NBPREDYEAR<<" "<<BSScore<<" ";
BSScore=0;
include "initVdepMax.edp";
curStep=0;
numRun=0;
real thetaAssim=1;
{
CURYEAR=FIRSTYEAR;
include "initVdepMax.edp";
include "updateVdepFromTobs.edp";
include "updateVdepMax.edp";
for (int i=0;i<114;i++)
cout<<" vdepMax[i]="<< vdepMax[i]<<endl;
CURYEAR=LASTYEAR;
coefP=1;
include "defP.edp";
include "computeBS.edp";
bsFile<<BSScore<<" ";
}
{
CURYEAR=2008;
include "initVdepMax.edp";
include "updateVdepFromTobs.edp";
include "updateVdepMax.edp";
for (int i=0;i<114;i++)
cout<<" vdepMax[i]="<< vdepMax[i]<<endl;
CURYEAR=LASTYEAR;
coefP=1;
include "defP.edp";
include "computeBS.edp";
bsFile<<BSScore<<endl;
}
}
FreeFem++-nw loopOnSimuNK.edp 2019 1
FreeFem++-nw loopOnSimuNK.edp 2018 2
FreeFem++-nw loopOnSimuNK.edp 2017 3
FreeFem++-nw loopOnSimuNK.edp 2016 4
FreeFem++-nw loopOnSimuNK.edp 2015 5
FreeFem++-nw loopOnSimuNK.edp 2014 6
FreeFem++-nw loopOnSimuNK.edp 2013 7
FreeFem++-nw loopOnSimuNK.edp 2012 8
FreeFem++-nw loopOnSimuNK.edp 2011 9
FreeFem++-nw loopOnSimuNK.edp 2010 10
FreeFem++-nw loopOnSimuNK.edp 2009 11
load "BinaryIO";
include "../paramsTS.edp";
include "../myMacro.edp";
include "../defGeomGen.edp"
include "../defDofsGen.edp"
include "../buildTobs.edp"
include "defTigerVar.edp"
include "defSimuParamsGen.edp";
real alphaw0,alphaw1;
string SRCDIR="/home/olivierb/aedes/AEDESINTER_AUTOROUTES_ADAPTIV/";
string MODE="DOVISU";
//string MODE="COMPUTEBS";
include "defVdepMax.edp"
include "initVdepMax.edp"
include "../postTraitementLoad.edp"
if (ARGV.n>2){
FIRSTYEAR=atoi(ARGV[2]);
}
if (ARGV.n>3){
NBPREDYEAR=atoi(ARGV[3]);
LASTYEAR=FIRSTYEAR+NBPREDYEAR;
}
int curStep=0;
int numRun=0;
int size;
//visu ic
include "defOutputGen.edp";
include "defParamsGen.edp";
include "mat1D2Duseful.edp";
//include "postTraitementLoad.edp";
//include "inputParam.edp";
include "defVarfGen.edp";
include "defIC.edp";
include "defIC1dNull.edp";
load "BinaryIO";
include "./paramsTS.edp";
include "./myMacro.edp";
include "./defGeomGen.edp"
include "./defDofsGen.edp"
include "./buildTobs.edp"
include "defTigerVar.edp"
include "defSimuParamsGen.edp";
real alphaw0,alphaw1;
string SRCDIR="/home/olivierb/aedes/AEDESINTER_AUTOROUTES_ADAPTIV/";
string MODE="DOVISU";
string LOADDIR="/home/olivierb/aedes/AEDESINTER_AUTOROUTES_ADAPTIV/SAVE/";
//string MODE="COMPUTEBS";
include "./postTraitementLoad.edp"
if (ARGV.n>2){
FIRSTYEAR=atoi(ARGV[2]);
}
if (ARGV.n>3){
NBPREDYEAR=atoi(ARGV[3]);
LASTYEAR=FIRSTYEAR+NBPREDYEAR;
}
int curStep=0;
int numRun=0;
int size;
//visu ic
include "defOutputGen.edp";
include "defParamsGen.edp";
include "mat1D2Duseful.edp";
//include "postTraitementLoad.edp";
//include "inputParam.edp";
include "defVarfGen.edp";
LoadVec(w0s2d0[],LOADDIR+"/w0s2d0_"+string(numRun)+"_2008");
LoadVec(w0s2d1[],LOADDIR+"/w0s2d1_"+string(numRun)+"_2008");
LoadVec(w0s2d2[],LOADDIR+"/w0s2d2_"+string(numRun)+"_2008");
LoadVec(w0s2d3[],LOADDIR+"/w0s2d3_"+string(numRun)+"_2008");
LoadVec(w0s2d4[],LOADDIR+"/w0s2d4_"+string(numRun)+"_2008");
LoadVec(w0s2d5[],LOADDIR+"/w0s2d5_"+string(numRun)+"_2008");
LoadVec(w0s2d6[],LOADDIR+"/w0s2d6_"+string(numRun)+"_2008");
LoadVec(w0s2d7[],LOADDIR+"/w0s2d7_"+string(numRun)+"_2008");
LoadVec(w0s2d8[],LOADDIR+"/w0s2d8_"+string(numRun)+"_2008");
LoadVec(w0s2d9[],LOADDIR+"/w0s2d9_"+string(numRun)+"_2008");
LoadVec(w0s2d10[],LOADDIR+"/w0s2d10_"+string(numRun)+"_2008");
LoadVec(w0s2d11[],LOADDIR+"/w0s2d11_"+string(numRun)+"_2008");
LoadVec(w0s2d12[],LOADDIR+"/w0s2d12_"+string(numRun)+"_2008");
LoadVec(w0s2d13[],LOADDIR+"/w0s2d13_"+string(numRun)+"_2008");
LoadVec(w0s2d14[],LOADDIR+"/w0s2d14_"+string(numRun)+"_2008");
string lcmm=LOADDIR+"/w0s2d0_"+string(numRun)+"_2008";
plot(w0s2d0,w0s2d1,w0s2d2,w0s2d3,w0s2d4,w0s2d5,w0s2d6,w0s2d7,w0s2d8,w0s2d9,w0s2d10,w0s2d11,w0s2d12,w0s2d13,w0s2d14,cmm=lcmm,ps=lcmm+".eps",value=true,fill=true,wait=true);
FreeFem++-nw cas1.edp
FreeFem++-nw cas1.edp 2009 1
FreeFem++-nw cas1.edp 2010 1
FreeFem++-nw cas1.edp 2011 1
FreeFem++-nw cas1.edp 2012 1
FreeFem++-nw cas1.edp 2013 1
FreeFem++-nw cas1.edp 2014 1
FreeFem++-nw cas1.edp 2015 1
FreeFem++-nw cas1.edp 2016 1
FreeFem++-nw cas1.edp 2017 1
FreeFem++-nw cas1.edp 2018 1
FreeFem++-nw cas1.edp 2019 1
FreeFem++-nw cas1.edp 2009 2
FreeFem++-nw cas1.edp 2010 2
FreeFem++-nw cas1.edp 2011 2
FreeFem++-nw cas1.edp 2012 2
FreeFem++-nw cas1.edp 2013 2
FreeFem++-nw cas1.edp 2014 2
FreeFem++-nw cas1.edp 2015 2