Commit 2ee985eb authored by Olivier Bonnefon's avatar Olivier Bonnefon
Browse files

add BS and set data until 2019

parent 717a9e53
real BSScore=0;
int indexYear=CURYEAR-2003;
for (int i=0;i<114;i++){
real lBS=0;
//premier terme
lBS= (1-vdepP[i])^2;
if (Tobs[indexYear*114+i] == 0)
lBS+=-1;
BSScore+=lBS*lBS;
//second terme
lBS= (1-vdepP[i])*vdepP[i];
if (Tobs[indexYear*114+i] == 1)
lBS+= -1;
BSScore+=lBS*lBS;
//troisieme terme
lBS= (vdepP[i])^2;
if (Tobs[indexYear*114+i] == 2)
lBS+= -1;
BSScore+=lBS*lBS;
}
real[int] vdepP(114);
real coefP=1;
for (int i=0;i<114;i++){
vdepP[i]=1;
if (coefP*vdepMax[i]<1)
vdepP[i]=coefP*vdepMax[i];
}
includepath += "../../Tools1D/"
includepath +="../"
includepath += "../../Recorder/"
includepath += "../../examples1D/"
loadpath += "/opt/freefem++-3.58/examples++-load/"
load "BinaryIO";
include "../paramsTS.edp";
include "../myMacro.edp";
include "../defGeomGen.edp"
include "../defDofsGen.edp"
include "../buildTobs.edp"
//string MODE="DOVISU";
string MODE="COMPUTEBS";
include "defVdepMax.edp"
include "initVdepMax.edp"
include "../postTraitementLoad.edp"
int curStep=0;
int numRun=0;
int size;
CURYEAR=2019;
{
cout<<"../SAVE/u2d0_"+string(CURYEAR)+"_rep_"+string(numRun)<<endl;
GetSizeVec(size,"../SAVE/u2d0_"+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u0(size);
//cout<<"size="<<size<<" "<<Vh0.ndof<<endl;
int nbDays=size/Vh0.ndof;
LoadVec(u0,"../SAVE/u2d0_"+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,"../SAVE/u2d1_"+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u1(size);
LoadVec(u1,"../SAVE/u2d1_"+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,"../SAVE/u2d2_"+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u2(size);
LoadVec(u2,"../SAVE/u2d2_"+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,"../SAVE/u2d3_"+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u3(size);
LoadVec(u3,"../SAVE/u2d3_"+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,"../SAVE/u2d4_"+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u4(size);
LoadVec(u4,"../SAVE/u2d4_"+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,"../SAVE/u2d5_"+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u5(size);
LoadVec(u5,"../SAVE/u2d5_"+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,"../SAVE/u2d6_"+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u6(size);
LoadVec(u6,"../SAVE/u2d6_"+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,"../SAVE/u2d7_"+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u7(size);
LoadVec(u7,"../SAVE/u2d7_"+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,"../SAVE/u2d8_"+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u8(size);
LoadVec(u8,"../SAVE/u2d8_"+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,"../SAVE/u2d9_"+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u9(size);
LoadVec(u9,"../SAVE/u2d9_"+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,"../SAVE/u2d10_"+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u10(size);
LoadVec(u10,"../SAVE/u2d10_"+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,"../SAVE/u2d11_"+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u11(size);
LoadVec(u11,"../SAVE/u2d11_"+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,"../SAVE/u2d12_"+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u12(size);
LoadVec(u12,"../SAVE/u2d12_"+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,"../SAVE/u2d13_"+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u13(size);
LoadVec(u13,"../SAVE/u2d13_"+string(CURYEAR)+"_rep_"+string(numRun));
GetSizeVec(size,"../SAVE/u2d14_"+string(CURYEAR)+"_rep_"+string(numRun));
real[int] u14(size);
LoadVec(u14,"../SAVE/u2d14_"+string(CURYEAR)+"_rep_"+string(numRun));
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++;
}
if (MODE=="COMPUTEBS"){
include "defP.edp";
include "computeBS.edp";
cout<<BSScore<<endl;
}
}
include "stateToFrance.edp";
vdep0[]=IVDep0*vall[];
vdep1[]=IVDep1*vall[];
vdep2[]=IVDep2*vall[];
vdep3[]=IVDep3*vall[];
vdep4[]=IVDep4*vall[];
vdep5[]=IVDep5*vall[];
vdep6[]=IVDep6*vall[];
vdep7[]=IVDep7*vall[];
vdep8[]=IVDep8*vall[];
vdep9[]=IVDep9*vall[];
vdep10[]=IVDep10*vall[];
vdep11[]=IVDep11*vall[];
vdep12[]=IVDep12*vall[];
vdep13[]=IVDep13*vall[];
vdep14[]=IVDep14*vall[];
vdep15[]=IVDep15*vall[];
vdep16[]=IVDep16*vall[];
vdep17[]=IVDep17*vall[];
vdep18[]=IVDep18*vall[];
vdep19[]=IVDep19*vall[];
vdep20[]=IVDep20*vall[];
vdep21[]=IVDep21*vall[];
vdep22[]=IVDep22*vall[];
vdep23[]=IVDep23*vall[];
vdep24[]=IVDep24*vall[];
vdep25[]=IVDep25*vall[];
vdep26[]=IVDep26*vall[];
vdep27[]=IVDep27*vall[];
vdep28[]=IVDep28*vall[];
vdep29[]=IVDep29*vall[];
vdep30[]=IVDep30*vall[];
vdep31[]=IVDep31*vall[];
vdep32[]=IVDep32*vall[];
vdep33[]=IVDep33*vall[];
vdep34[]=IVDep34*vall[];
vdep35[]=IVDep35*vall[];
vdep36[]=IVDep36*vall[];
vdep37[]=IVDep37*vall[];
vdep38[]=IVDep38*vall[];
vdep39[]=IVDep39*vall[];
vdep40[]=IVDep40*vall[];
vdep41[]=IVDep41*vall[];
vdep42[]=IVDep42*vall[];
vdep43[]=IVDep43*vall[];
vdep44[]=IVDep44*vall[];
vdep45[]=IVDep45*vall[];
vdep46[]=IVDep46*vall[];
vdep47[]=IVDep47*vall[];
vdep48[]=IVDep48*vall[];
vdep49[]=IVDep49*vall[];
vdep50[]=IVDep50*vall[];
vdep51[]=IVDep51*vall[];
vdep52[]=IVDep52*vall[];
vdep53[]=IVDep53*vall[];
vdep54[]=IVDep54*vall[];
vdep55[]=IVDep55*vall[];
vdep56[]=IVDep56*vall[];
vdep57[]=IVDep57*vall[];
vdep58[]=IVDep58*vall[];
vdep59[]=IVDep59*vall[];
vdep60[]=IVDep60*vall[];
vdep61[]=IVDep61*vall[];
vdep62[]=IVDep62*vall[];
vdep63[]=IVDep63*vall[];
vdep64[]=IVDep64*vall[];
vdep65[]=IVDep65*vall[];
vdep66[]=IVDep66*vall[];
vdep67[]=IVDep67*vall[];
vdep68[]=IVDep68*vall[];
vdep69[]=IVDep69*vall[];
vdep70[]=IVDep70*vall[];
vdep71[]=IVDep71*vall[];
vdep72[]=IVDep72*vall[];
vdep73[]=IVDep73*vall[];
vdep74[]=IVDep74*vall[];
vdep75[]=IVDep75*vall[];
vdep76[]=IVDep76*vall[];
vdep77[]=IVDep77*vall[];
vdep78[]=IVDep78*vall[];
vdep79[]=IVDep79*vall[];
vdep80[]=IVDep80*vall[];
vdep81[]=IVDep81*vall[];
vdep82[]=IVDep82*vall[];
vdep83[]=IVDep83*vall[];
vdep84[]=IVDep84*vall[];
vdep85[]=IVDep85*vall[];
vdep86[]=IVDep86*vall[];
vdep87[]=IVDep87*vall[];
vdep88[]=IVDep88*vall[];
vdep89[]=IVDep89*vall[];
vdep90[]=IVDep90*vall[];
vdep91[]=IVDep91*vall[];
vdep92[]=IVDep92*vall[];
vdep93[]=IVDep93*vall[];
vdep94[]=IVDep94*vall[];
vdep95[]=IVDep95*vall[];
vdep96[]=IVDep96*vall[];
vdep97[]=IVDep97*vall[];
vdep98[]=IVDep98*vall[];
vdep99[]=IVDep99*vall[];
vdep100[]=IVDep100*vall[];
vdep101[]=IVDep101*vall[];
vdep102[]=IVDep102*vall[];
vdep103[]=IVDep103*vall[];
vdep104[]=IVDep104*vall[];
vdep105[]=IVDep105*vall[];
vdep106[]=IVDep106*vall[];
vdep107[]=IVDep107*vall[];
vdep108[]=IVDep108*vall[];
vdep109[]=IVDep109*vall[];
vdep110[]=IVDep110*vall[];
vdep111[]=IVDep111*vall[];
vdep112[]=IVDep112*vall[];
vdep113[]=IVDep113*vall[];
vall=0;
vabuf[]=IVA0*u2d0[];
adddofs(vall,vabuf);
vabuf[]=IVA1*u2d1[];
adddofs(vall,vabuf);
vabuf[]=IVA2*u2d2[];
adddofs(vall,vabuf);
vabuf[]=IVA3*u2d3[];
adddofs(vall,vabuf);
vabuf[]=IVA4*u2d4[];
adddofs(vall,vabuf);
vabuf[]=IVA5*u2d5[];
adddofs(vall,vabuf);
vabuf[]=IVA6*u2d6[];
adddofs(vall,vabuf);
vabuf[]=IVA7*u2d7[];
adddofs(vall,vabuf);
vabuf[]=IVA8*u2d8[];
adddofs(vall,vabuf);
vabuf[]=IVA9*u2d9[];
adddofs(vall,vabuf);
vabuf[]=IVA10*u2d10[];
adddofs(vall,vabuf);
vabuf[]=IVA11*u2d11[];
adddofs(vall,vabuf);
vabuf[]=IVA12*u2d12[];
adddofs(vall,vabuf);
vabuf[]=IVA13*u2d13[];
adddofs(vall,vabuf);
vabuf[]=IVA14*u2d14[];
adddofs(vall,vabuf);
real aux;
aux=vdep0[].max;
if (aux>vdepMax[0])
vdepMax[0]=aux;
aux=vdep1[].max;
if (aux>vdepMax[1])
vdepMax[1]=aux;
aux=vdep2[].max;
if (aux>vdepMax[2])
vdepMax[2]=aux;
aux=vdep3[].max;
if (aux>vdepMax[3])
vdepMax[3]=aux;
aux=vdep4[].max;
if (aux>vdepMax[4])
vdepMax[4]=aux;
aux=vdep5[].max;
if (aux>vdepMax[5])
vdepMax[5]=aux;
aux=vdep6[].max;
if (aux>vdepMax[6])
vdepMax[6]=aux;
aux=vdep7[].max;
if (aux>vdepMax[7])
vdepMax[7]=aux;
aux=vdep8[].max;
if (aux>vdepMax[8])
vdepMax[8]=aux;
aux=vdep9[].max;
if (aux>vdepMax[9])
vdepMax[9]=aux;
aux=vdep10[].max;
if (aux>vdepMax[10])
vdepMax[10]=aux;
aux=vdep11[].max;
if (aux>vdepMax[11])
vdepMax[11]=aux;
aux=vdep12[].max;
if (aux>vdepMax[12])
vdepMax[12]=aux;
aux=vdep13[].max;
if (aux>vdepMax[13])
vdepMax[13]=aux;
aux=vdep14[].max;
if (aux>vdepMax[14])
vdepMax[14]=aux;
aux=vdep15[].max;
if (aux>vdepMax[15])
vdepMax[15]=aux;
aux=vdep16[].max;
if (aux>vdepMax[16])
vdepMax[16]=aux;
aux=vdep17[].max;
if (aux>vdepMax[17])
vdepMax[17]=aux;
aux=vdep18[].max;
if (aux>vdepMax[18])
vdepMax[18]=aux;
aux=vdep19[].max;
if (aux>vdepMax[19])
vdepMax[19]=aux;
aux=vdep20[].max;
if (aux>vdepMax[20])
vdepMax[20]=aux;
aux=vdep21[].max;
if (aux>vdepMax[21])
vdepMax[21]=aux;
aux=vdep22[].max;
if (aux>vdepMax[22])
vdepMax[22]=aux;
aux=vdep23[].max;
if (aux>vdepMax[23])
vdepMax[23]=aux;
aux=vdep24[].max;
if (aux>vdepMax[24])
vdepMax[24]=aux;
aux=vdep25[].max;
if (aux>vdepMax[25])
vdepMax[25]=aux;
aux=vdep26[].max;
if (aux>vdepMax[26])
vdepMax[26]=aux;
aux=vdep27[].max;
if (aux>vdepMax[27])
vdepMax[27]=aux;
aux=vdep28[].max;
if (aux>vdepMax[28])
vdepMax[28]=aux;
aux=vdep29[].max;
if (aux>vdepMax[29])
vdepMax[29]=aux;
aux=vdep30[].max;
if (aux>vdepMax[30])
vdepMax[30]=aux;
aux=vdep31[].max;
if (aux>vdepMax[31])
vdepMax[31]=aux;
aux=vdep32[].max;
if (aux>vdepMax[32])
vdepMax[32]=aux;
aux=vdep33[].max;
if (aux>vdepMax[33])
vdepMax[33]=aux;
aux=vdep34[].max;
if (aux>vdepMax[34])
vdepMax[34]=aux;
aux=vdep35[].max;
if (aux>vdepMax[35])
vdepMax[35]=aux;
aux=vdep36[].max;
if (aux>vdepMax[36])
vdepMax[36]=aux;
aux=vdep37[].max;
if (aux>vdepMax[37])
vdepMax[37]=aux;
aux=vdep38[].max;
if (aux>vdepMax[38])
vdepMax[38]=aux;
aux=vdep39[].max;
if (aux>vdepMax[39])
vdepMax[39]=aux;
aux=vdep40[].max;
if (aux>vdepMax[40])
vdepMax[40]=aux;
aux=vdep41[].max;
if (aux>vdepMax[41])
vdepMax[41]=aux;
aux=vdep42[].max;
if (aux>vdepMax[42])
vdepMax[42]=aux;
aux=vdep43[].max;
if (aux>vdepMax[43])
vdepMax[43]=aux;
aux=vdep44[].max;
if (aux>vdepMax[44])
vdepMax[44]=aux;
aux=vdep45[].max;
if (aux>vdepMax[45])
vdepMax[45]=aux;
aux=vdep46[].max;
if (aux>vdepMax[46])
vdepMax[46]=aux;
aux=vdep47[].max;
if (aux>vdepMax[47])
vdepMax[47]=aux;
aux=vdep48[].max;
if (aux>vdepMax[48])
vdepMax[48]=aux;
aux=vdep49[].max;
if (aux>vdepMax[49])
vdepMax[49]=aux;
aux=vdep50[].max;
if (aux>vdepMax[50])
vdepMax[50]=aux;
aux=vdep51[].max;
if (aux>vdepMax[51])
vdepMax[51]=aux;
aux=vdep52[].max;
if (aux>vdepMax[52])
vdepMax[52]=aux;
aux=vdep53[].max;
if (aux>vdepMax[53])
vdepMax[53]=aux;
aux=vdep54[].max;
if (aux>vdepMax[54])
vdepMax[54]=aux;
aux=vdep55[].max;
if (aux>vdepMax[55])
vdepMax[55]=aux;
aux=vdep56[].max;
if (aux>vdepMax[56])
vdepMax[56]=aux;
aux=vdep57[].max;
if (aux>vdepMax[57])
vdepMax[57]=aux;
aux=vdep58[].max;
if (aux>vdepMax[58])
vdepMax[58]=aux;
aux=vdep59[].max;
if (aux>vdepMax[59])
vdepMax[59]=aux;
aux=vdep60[].max;
if (aux>vdepMax[60])
vdepMax[60]=aux;
aux=vdep61[].max;
if (aux>vdepMax[61])
vdepMax[61]=aux;
aux=vdep62[].max;
if (aux>vdepMax[62])
vdepMax[62]=aux;
aux=vdep63[].max;
if (aux>vdepMax[63])
vdepMax[63]=aux;
aux=vdep64[].max;
if (aux>vdepMax[64])
vdepMax[64]=aux;
aux=vdep65[].max;
if (aux>vdepMax[65])
vdepMax[65]=aux;
aux=vdep66[].max;
if (aux>vdepMax[66])
vdepMax[66]=aux;
aux=vdep67[].max;
if (aux>vdepMax[67])
vdepMax[67]=aux;
aux=vdep68[].max;
if (aux>vdepMax[68])
vdepMax[68]=aux;
aux=vdep69[].max;
if (aux>vdepMax[69])
vdepMax[69]=aux;
aux=vdep70[].max;
if (aux>vdepMax[70])
vdepMax[70]=aux;
aux=vdep71[].max;
if (aux>vdepMax[71])
vdepMax[71]=aux;
aux=vdep72[].max;
if (aux>vdepMax[72])
vdepMax[72]=aux;
aux=vdep73[].max;
if (aux>vdepMax[73])
vdepMax[73]=aux;
aux=vdep74[].max;
if (aux>vdepMax[74])
vdepMax[74]=aux;
aux=vdep75[].max;
if (aux>vdepMax[75])
vdepMax[75]=aux;
aux=vdep76[].max;
if (aux>vdepMax[76])
vdepMax[76]=aux;
aux=vdep77[].max;
if (aux>vdepMax[77])
vdepMax[77]=aux;
aux=vdep78[].max;
if (aux>vdepMax[78])
vdepMax[78]=aux;