Commit 15a2153b authored by Olivier Bonnefon's avatar Olivier Bonnefon
Browse files

mv INIT_W1_2011

parent 8802bdf5
This source diff could not be displayed because it is too large. You can view the blob instead.
load "BinaryIO";
include "../myMacro.edp";
include "../defGeomGen.edp";
include "../defParamsGen.edp";
include "../defTigerVar.edp";
mesh ThALL=readmesh("../france.msh");
plot(ThALL,wait=true);
fespace VALL(ThALL,P1);
VALL vall;
VALL vabuf;
matrix IVA0=interpolate(VALL,Vh0,inside=1);
matrix IVA1=interpolate(VALL,Vh1,inside=1);
matrix IVA2=interpolate(VALL,Vh2,inside=1);
matrix IVA3=interpolate(VALL,Vh3,inside=1);
matrix IVA4=interpolate(VALL,Vh4,inside=1);
matrix IVA5=interpolate(VALL,Vh5,inside=1);
matrix IVA6=interpolate(VALL,Vh6,inside=1);
matrix IVA7=interpolate(VALL,Vh7,inside=1);
matrix IVA8=interpolate(VALL,Vh8,inside=1);
matrix IVA9=interpolate(VALL,Vh9,inside=1);
matrix IVA10=interpolate(VALL,Vh10,inside=1);
matrix IVA11=interpolate(VALL,Vh11,inside=1);
matrix IVA12=interpolate(VALL,Vh12,inside=1);
matrix IVA13=interpolate(VALL,Vh13,inside=1);
matrix IVA14=interpolate(VALL,Vh14,inside=1);
matrix I0VA=interpolate(Vh0,VALL);
matrix I1VA=interpolate(Vh1,VALL);
matrix I2VA=interpolate(Vh2,VALL);
matrix I3VA=interpolate(Vh3,VALL);
matrix I4VA=interpolate(Vh4,VALL);
matrix I5VA=interpolate(Vh5,VALL);
matrix I6VA=interpolate(Vh6,VALL);
matrix I7VA=interpolate(Vh7,VALL);
matrix I8VA=interpolate(Vh8,VALL);
matrix I9VA=interpolate(Vh9,VALL);
matrix I10VA=interpolate(Vh10,VALL);
matrix I11VA=interpolate(Vh11,VALL);
matrix I12VA=interpolate(Vh12,VALL);
matrix I13VA=interpolate(Vh13,VALL);
matrix I14VA=interpolate(Vh14,VALL);
vall=0;
int nnaux=vall.n;
//it strats to 1 because of R num
mesh ThDep1=readmesh("../dep0.msh");
fespace VDep1(ThDep1,P1);
VDep1 vdep1;
matrix IVDep1=interpolate(VDep1,VALL);
mesh ThDep2=readmesh("../dep1.msh");
fespace VDep2(ThDep2,P1);
VDep2 vdep2;
matrix IVDep2=interpolate(VDep2,VALL);
mesh ThDep3=readmesh("../dep2.msh");
fespace VDep3(ThDep3,P1);
VDep3 vdep3;
matrix IVDep3=interpolate(VDep3,VALL);
mesh ThDep4=readmesh("../dep3.msh");
fespace VDep4(ThDep4,P1);
VDep4 vdep4;
matrix IVDep4=interpolate(VDep4,VALL);
mesh ThDep5=readmesh("../dep4.msh");
fespace VDep5(ThDep5,P1);
VDep5 vdep5;
matrix IVDep5=interpolate(VDep5,VALL);
mesh ThDep6=readmesh("../dep5.msh");
fespace VDep6(ThDep6,P1);
VDep6 vdep6;
matrix IVDep6=interpolate(VDep6,VALL);
mesh ThDep7=readmesh("../dep6.msh");
fespace VDep7(ThDep7,P1);
VDep7 vdep7;
matrix IVDep7=interpolate(VDep7,VALL);
mesh ThDep8=readmesh("../dep7.msh");
fespace VDep8(ThDep8,P1);
VDep8 vdep8;
matrix IVDep8=interpolate(VDep8,VALL);
mesh ThDep9=readmesh("../dep8.msh");
fespace VDep9(ThDep9,P1);
VDep9 vdep9;
matrix IVDep9=interpolate(VDep9,VALL);
mesh ThDep10=readmesh("../dep9.msh");
fespace VDep10(ThDep10,P1);
VDep10 vdep10;
matrix IVDep10=interpolate(VDep10,VALL);
mesh ThDep11=readmesh("../dep10.msh");
fespace VDep11(ThDep11,P1);
VDep11 vdep11;
matrix IVDep11=interpolate(VDep11,VALL);
mesh ThDep12=readmesh("../dep11.msh");
fespace VDep12(ThDep12,P1);
VDep12 vdep12;
matrix IVDep12=interpolate(VDep12,VALL);
mesh ThDep13=readmesh("../dep12.msh");
fespace VDep13(ThDep13,P1);
VDep13 vdep13;
matrix IVDep13=interpolate(VDep13,VALL);
mesh ThDep14=readmesh("../dep13.msh");
fespace VDep14(ThDep14,P1);
VDep14 vdep14;
matrix IVDep14=interpolate(VDep14,VALL);
mesh ThDep15=readmesh("../dep14.msh");
fespace VDep15(ThDep15,P1);
VDep15 vdep15;
matrix IVDep15=interpolate(VDep15,VALL);
mesh ThDep16=readmesh("../dep15.msh");
fespace VDep16(ThDep16,P1);
VDep16 vdep16;
matrix IVDep16=interpolate(VDep16,VALL);
mesh ThDep17=readmesh("../dep16.msh");
fespace VDep17(ThDep17,P1);
VDep17 vdep17;
matrix IVDep17=interpolate(VDep17,VALL);
mesh ThDep18=readmesh("../dep17.msh");
fespace VDep18(ThDep18,P1);
VDep18 vdep18;
matrix IVDep18=interpolate(VDep18,VALL);
mesh ThDep19=readmesh("../dep18.msh");
fespace VDep19(ThDep19,P1);
VDep19 vdep19;
matrix IVDep19=interpolate(VDep19,VALL);
mesh ThDep20=readmesh("../dep19.msh");
fespace VDep20(ThDep20,P1);
VDep20 vdep20;
matrix IVDep20=interpolate(VDep20,VALL);
mesh ThDep21=readmesh("../dep20.msh");
fespace VDep21(ThDep21,P1);
VDep21 vdep21;
matrix IVDep21=interpolate(VDep21,VALL);
mesh ThDep22=readmesh("../dep21.msh");
fespace VDep22(ThDep22,P1);
VDep22 vdep22;
matrix IVDep22=interpolate(VDep22,VALL);
mesh ThDep23=readmesh("../dep22.msh");
fespace VDep23(ThDep23,P1);
VDep23 vdep23;
matrix IVDep23=interpolate(VDep23,VALL);
mesh ThDep24=readmesh("../dep23.msh");
fespace VDep24(ThDep24,P1);
VDep24 vdep24;
matrix IVDep24=interpolate(VDep24,VALL);
mesh ThDep25=readmesh("../dep24.msh");
fespace VDep25(ThDep25,P1);
VDep25 vdep25;
matrix IVDep25=interpolate(VDep25,VALL);
mesh ThDep26=readmesh("../dep25.msh");
fespace VDep26(ThDep26,P1);
VDep26 vdep26;
matrix IVDep26=interpolate(VDep26,VALL);
mesh ThDep27=readmesh("../dep26.msh");
fespace VDep27(ThDep27,P1);
VDep27 vdep27;
matrix IVDep27=interpolate(VDep27,VALL);
mesh ThDep28=readmesh("../dep27.msh");
fespace VDep28(ThDep28,P1);
VDep28 vdep28;
matrix IVDep28=interpolate(VDep28,VALL);
mesh ThDep29=readmesh("../dep28.msh");
fespace VDep29(ThDep29,P1);
VDep29 vdep29;
matrix IVDep29=interpolate(VDep29,VALL);
mesh ThDep30=readmesh("../dep29.msh");
fespace VDep30(ThDep30,P1);
VDep30 vdep30;
matrix IVDep30=interpolate(VDep30,VALL);
mesh ThDep31=readmesh("../dep30.msh");
fespace VDep31(ThDep31,P1);
VDep31 vdep31;
matrix IVDep31=interpolate(VDep31,VALL);
mesh ThDep32=readmesh("../dep31.msh");
fespace VDep32(ThDep32,P1);
VDep32 vdep32;
matrix IVDep32=interpolate(VDep32,VALL);
mesh ThDep33=readmesh("../dep32.msh");
fespace VDep33(ThDep33,P1);
VDep33 vdep33;
matrix IVDep33=interpolate(VDep33,VALL);
mesh ThDep34=readmesh("../dep33.msh");
fespace VDep34(ThDep34,P1);
VDep34 vdep34;
matrix IVDep34=interpolate(VDep34,VALL);
mesh ThDep35=readmesh("../dep34.msh");
fespace VDep35(ThDep35,P1);
VDep35 vdep35;
matrix IVDep35=interpolate(VDep35,VALL);
mesh ThDep36=readmesh("../dep35.msh");
fespace VDep36(ThDep36,P1);
VDep36 vdep36;
matrix IVDep36=interpolate(VDep36,VALL);
mesh ThDep37=readmesh("../dep36.msh");
fespace VDep37(ThDep37,P1);
VDep37 vdep37;
matrix IVDep37=interpolate(VDep37,VALL);
mesh ThDep38=readmesh("../dep37.msh");
fespace VDep38(ThDep38,P1);
VDep38 vdep38;
matrix IVDep38=interpolate(VDep38,VALL);
mesh ThDep39=readmesh("../dep38.msh");
fespace VDep39(ThDep39,P1);
VDep39 vdep39;
matrix IVDep39=interpolate(VDep39,VALL);
mesh ThDep40=readmesh("../dep39.msh");
fespace VDep40(ThDep40,P1);
VDep40 vdep40;
matrix IVDep40=interpolate(VDep40,VALL);
mesh ThDep41=readmesh("../dep40.msh");
fespace VDep41(ThDep41,P1);
VDep41 vdep41;
matrix IVDep41=interpolate(VDep41,VALL);
mesh ThDep42=readmesh("../dep41.msh");
fespace VDep42(ThDep42,P1);
VDep42 vdep42;
matrix IVDep42=interpolate(VDep42,VALL);
mesh ThDep43=readmesh("../dep42.msh");
fespace VDep43(ThDep43,P1);
VDep43 vdep43;
matrix IVDep43=interpolate(VDep43,VALL);
mesh ThDep44=readmesh("../dep43.msh");
fespace VDep44(ThDep44,P1);
VDep44 vdep44;
matrix IVDep44=interpolate(VDep44,VALL);
mesh ThDep45=readmesh("../dep44.msh");
fespace VDep45(ThDep45,P1);
VDep45 vdep45;
matrix IVDep45=interpolate(VDep45,VALL);
mesh ThDep46=readmesh("../dep45.msh");
fespace VDep46(ThDep46,P1);
VDep46 vdep46;
matrix IVDep46=interpolate(VDep46,VALL);
mesh ThDep47=readmesh("../dep46.msh");
fespace VDep47(ThDep47,P1);
VDep47 vdep47;
matrix IVDep47=interpolate(VDep47,VALL);
mesh ThDep48=readmesh("../dep47.msh");
fespace VDep48(ThDep48,P1);
VDep48 vdep48;
matrix IVDep48=interpolate(VDep48,VALL);
mesh ThDep49=readmesh("../dep48.msh");
fespace VDep49(ThDep49,P1);
VDep49 vdep49;
matrix IVDep49=interpolate(VDep49,VALL);
mesh ThDep50=readmesh("../dep49.msh");
fespace VDep50(ThDep50,P1);
VDep50 vdep50;
matrix IVDep50=interpolate(VDep50,VALL);
mesh ThDep51=readmesh("../dep50.msh");
fespace VDep51(ThDep51,P1);
VDep51 vdep51;
matrix IVDep51=interpolate(VDep51,VALL);
mesh ThDep52=readmesh("../dep51.msh");
fespace VDep52(ThDep52,P1);
VDep52 vdep52;
matrix IVDep52=interpolate(VDep52,VALL);
mesh ThDep53=readmesh("../dep52.msh");
fespace VDep53(ThDep53,P1);
VDep53 vdep53;
matrix IVDep53=interpolate(VDep53,VALL);
mesh ThDep54=readmesh("../dep53.msh");
fespace VDep54(ThDep54,P1);
VDep54 vdep54;
matrix IVDep54=interpolate(VDep54,VALL);
mesh ThDep55=readmesh("../dep54.msh");
fespace VDep55(ThDep55,P1);
VDep55 vdep55;
matrix IVDep55=interpolate(VDep55,VALL);
mesh ThDep56=readmesh("../dep55.msh");
fespace VDep56(ThDep56,P1);
VDep56 vdep56;
matrix IVDep56=interpolate(VDep56,VALL);
mesh ThDep57=readmesh("../dep56.msh");
fespace VDep57(ThDep57,P1);
VDep57 vdep57;
matrix IVDep57=interpolate(VDep57,VALL);
mesh ThDep58=readmesh("../dep57.msh");
fespace VDep58(ThDep58,P1);
VDep58 vdep58;
matrix IVDep58=interpolate(VDep58,VALL);
mesh ThDep59=readmesh("../dep58.msh");
fespace VDep59(ThDep59,P1);
VDep59 vdep59;
matrix IVDep59=interpolate(VDep59,VALL);
mesh ThDep60=readmesh("../dep59.msh");
fespace VDep60(ThDep60,P1);
VDep60 vdep60;
matrix IVDep60=interpolate(VDep60,VALL);
mesh ThDep61=readmesh("../dep60.msh");
fespace VDep61(ThDep61,P1);
VDep61 vdep61;
matrix IVDep61=interpolate(VDep61,VALL);
mesh ThDep62=readmesh("../dep61.msh");
fespace VDep62(ThDep62,P1);
VDep62 vdep62;
matrix IVDep62=interpolate(VDep62,VALL);
mesh ThDep63=readmesh("../dep62.msh");
fespace VDep63(ThDep63,P1);
VDep63 vdep63;
matrix IVDep63=interpolate(VDep63,VALL);
mesh ThDep64=readmesh("../dep63.msh");
fespace VDep64(ThDep64,P1);
VDep64 vdep64;
matrix IVDep64=interpolate(VDep64,VALL);
mesh ThDep65=readmesh("../dep64.msh");
fespace VDep65(ThDep65,P1);
VDep65 vdep65;
matrix IVDep65=interpolate(VDep65,VALL);
mesh ThDep66=readmesh("../dep65.msh");
fespace VDep66(ThDep66,P1);
VDep66 vdep66;
matrix IVDep66=interpolate(VDep66,VALL);
mesh ThDep67=readmesh("../dep66.msh");
fespace VDep67(ThDep67,P1);
VDep67 vdep67;
matrix IVDep67=interpolate(VDep67,VALL);
mesh ThDep68=readmesh("../dep67.msh");
fespace VDep68(ThDep68,P1);
VDep68 vdep68;
matrix IVDep68=interpolate(VDep68,VALL);
mesh ThDep69=readmesh("../dep68.msh");
fespace VDep69(ThDep69,P1);
VDep69 vdep69;
matrix IVDep69=interpolate(VDep69,VALL);
mesh ThDep70=readmesh("../dep69.msh");
fespace VDep70(ThDep70,P1);
VDep70 vdep70;
matrix IVDep70=interpolate(VDep70,VALL);
mesh ThDep71=readmesh("../dep70.msh");
fespace VDep71(ThDep71,P1);
VDep71 vdep71;
matrix IVDep71=interpolate(VDep71,VALL);
mesh ThDep72=readmesh("../dep71.msh");
fespace VDep72(ThDep72,P1);
VDep72 vdep72;
matrix IVDep72=interpolate(VDep72,VALL);
mesh ThDep73=readmesh("../dep72.msh");
fespace VDep73(ThDep73,P1);
VDep73 vdep73;
matrix IVDep73=interpolate(VDep73,VALL);
mesh ThDep74=readmesh("../dep73.msh");
fespace VDep74(ThDep74,P1);
VDep74 vdep74;
matrix IVDep74=interpolate(VDep74,VALL);
mesh ThDep75=readmesh("../dep74.msh");
fespace VDep75(ThDep75,P1);
VDep75 vdep75;
matrix IVDep75=interpolate(VDep75,VALL);
mesh ThDep76=readmesh("../dep75.msh");
fespace VDep76(ThDep76,P1);
VDep76 vdep76;
matrix IVDep76=interpolate(VDep76,VALL);
mesh ThDep77=readmesh("../dep76.msh");
fespace VDep77(ThDep77,P1);
VDep77 vdep77;
matrix IVDep77=interpolate(VDep77,VALL);
mesh ThDep78=readmesh("../dep77.msh");
fespace VDep78(ThDep78,P1);
VDep78 vdep78;
matrix IVDep78=interpolate(VDep78,VALL);
mesh ThDep79=readmesh("../dep78.msh");
fespace VDep79(ThDep79,P1);
VDep79 vdep79;
matrix IVDep79=interpolate(VDep79,VALL);
mesh ThDep80=readmesh("../dep79.msh");
fespace VDep80(ThDep80,P1);
VDep80 vdep80;
matrix IVDep80=interpolate(VDep80,VALL);
mesh ThDep81=readmesh("../dep80.msh");
fespace VDep81(ThDep81,P1);
VDep81 vdep81;
matrix IVDep81=interpolate(VDep81,VALL);
mesh ThDep82=readmesh("../dep81.msh");
fespace VDep82(ThDep82,P1);
VDep82 vdep82;
matrix IVDep82=interpolate(VDep82,VALL);
mesh ThDep83=readmesh("../dep82.msh");
fespace VDep83(ThDep83,P1);
VDep83 vdep83;
matrix IVDep83=interpolate(VDep83,VALL);
mesh ThDep84=readmesh("../dep83.msh");
fespace VDep84(ThDep84,P1);
VDep84 vdep84;
matrix IVDep84=interpolate(VDep84,VALL);
mesh ThDep85=readmesh("../dep84.msh");
fespace VDep85(ThDep85,P1);
VDep85 vdep85;
matrix IVDep85=interpolate(VDep85,VALL);
mesh ThDep86=readmesh("../dep85.msh");
fespace VDep86(ThDep86,P1);
VDep86 vdep86;
matrix IVDep86=interpolate(VDep86,VALL);
mesh ThDep87=readmesh("../dep86.msh");
fespace VDep87(ThDep87,P1);
VDep87 vdep87;
matrix IVDep87=interpolate(VDep87,VALL);
mesh ThDep88=readmesh("../dep87.msh");
fespace VDep88(ThDep88,P1);
VDep88 vdep88;
matrix IVDep88=interpolate(VDep88,VALL);
mesh ThDep89=readmesh("../dep88.msh");
fespace VDep89(ThDep89,P1);
VDep89 vdep89;
matrix IVDep89=interpolate(VDep89,VALL);
mesh ThDep90=readmesh("../dep89.msh");
fespace VDep90(ThDep90,P1);
VDep90 vdep90;
matrix IVDep90=interpolate(VDep90,VALL);
mesh ThDep91=readmesh("../dep90.msh");
fespace VDep91(ThDep91,P1);
VDep91 vdep91;
matrix IVDep91=interpolate(VDep91,VALL);
mesh ThDep92=readmesh("../dep91.msh");
fespace VDep92(ThDep92,P1);
VDep92 vdep92;
matrix IVDep92=interpolate(VDep92,VALL);
mesh ThDep93=readmesh("../dep92.msh");
fespace VDep93(ThDep93,P1);
VDep93 vdep93;
matrix IVDep93=interpolate(VDep93,VALL);
mesh ThDep94=readmesh("../dep93.msh");
fespace VDep94(ThDep94,P1);
VDep94 vdep94;
matrix IVDep94=interpolate(VDep94,VALL);
mesh ThDep95=readmesh("../dep94.msh");
fespace VDep95(ThDep95,P1);
VDep95 vdep95;
matrix IVDep95=interpolate(VDep95,VALL);
mesh ThDep96=readmesh("../dep95.msh");
fespace VDep96(ThDep96,P1);
VDep96 vdep96;
matrix IVDep96=interpolate(VDep96,VALL);
mesh ThDep97=readmesh("../dep96.msh");
fespace VDep97(ThDep97,P1);
VDep97 vdep97;
matrix IVDep97=interpolate(VDep97,VALL);
mesh ThDep98=readmesh("../dep97.msh");
fespace VDep98(ThDep98,P1);
VDep98 vdep98;
matrix IVDep98=interpolate(VDep98,VALL);
mesh ThDep99=readmesh("../dep98.msh");
fespace VDep99(ThDep99,P1);
VDep99 vdep99;
matrix IVDep99=interpolate(VDep99,VALL);
mesh ThDep100=readmesh("../dep99.msh");
fespace VDep100(ThDep100,P1);
VDep100 vdep100;
matrix IVDep100=interpolate(VDep100,VALL);
mesh ThDep101=readmesh("../dep100.msh");
fespace VDep101(ThDep101,P1);
VDep101 vdep101;
matrix IVDep101=interpolate(VDep101,VALL);
mesh ThDep102=readmesh("../dep101.msh");
fespace VDep102(ThDep102,P1);
VDep102 vdep102;
matrix IVDep102=interpolate(VDep102,VALL);
mesh ThDep103=readmesh("../dep102.msh");
fespace VDep103(ThDep103,P1);
VDep103 vdep103;
matrix IVDep103=interpolate(VDep103,VALL);
mesh ThDep104=readmesh("../dep103.msh");
fespace VDep104(ThDep104,P1);
VDep104 vdep104;
matrix IVDep104=interpolate(VDep104,VALL);
mesh ThDep105=readmesh("../dep104.msh");
fespace VDep105(ThDep105,P1);
VDep105 vdep105;
matrix IVDep105=interpolate(VDep105,VALL);
mesh ThDep106=readmesh("../dep105.msh");
fespace VDep106(ThDep106,P1);
VDep106 vdep106;
matrix IVDep106=interpolate(VDep106,VALL);
mesh ThDep107=readmesh("../dep106.msh");
fespace VDep107(ThDep107,P1);
VDep107 vdep107;
matrix IVDep107=interpolate(VDep107,VALL);
mesh ThDep108=readmesh("../dep107.msh");
fespace VDep108(ThDep108,P1);
VDep108 vdep108;
matrix IVDep108=interpolate(VDep108,VALL);
mesh ThDep109=readmesh("../dep108.msh");
fespace VDep109(ThDep109,P1);
VDep109 vdep109;
matrix IVDep109=interpolate(VDep109,VALL);
mesh ThDep110=readmesh("../dep109.msh");
fespace VDep110(ThDep110,P1);
VDep110 vdep110;
matrix IVDep110=interpolate(VDep110,VALL);
mesh ThDep111=readmesh("../dep110.msh");
fespace VDep111(ThDep111,P1);
VDep111 vdep111;
matrix IVDep111=interpolate(VDep111,VALL);
mesh ThDep112=readmesh("../dep111.msh");
fespace VDep112(ThDep112,P1);
VDep112 vdep112;
matrix IVDep112=interpolate(VDep112,VALL);
mesh ThDep113=readmesh("../dep112.msh");
fespace VDep113(ThDep113,P1);
VDep113 vdep113;
matrix IVDep113=interpolate(VDep113,VALL);
mesh ThDep114=readmesh("../dep113.msh");
fespace VDep114(ThDep114,P1);
VDep114 vdep114;
matrix IVDep114=interpolate(VDep114,VALL);
matrix IVbisDep1=interpolate(VALL,VDep1,inside=1);
matrix IVbisDep2=interpolate(VALL,VDep2,inside=1);
matrix IVbisDep3=interpolate(VALL,VDep3,inside=1);
matrix IVbisDep4=interpolate(VALL,VDep4,inside=1);
matrix IVbisDep5=interpolate(VALL,VDep5,inside=1);
matrix IVbisDep6=interpolate(VALL,VDep6,inside=1);
matrix IVbisDep7=interpolate(VALL,VDep7,inside=1);
matrix IVbisDep8=interpolate(VALL,VDep8,inside=1);
matrix IVbisDep9=interpolate(VALL,VDep9,inside=1);
matrix IVbisDep10=interpolate(VALL,VDep10,inside=1);
matrix IVbisDep11=interpolate(VALL,VDep11,inside=1);
matrix IVbisDep12=interpolate(VALL,VDep12,inside=1);
matrix IVbisDep13=interpolate(VALL,VDep13,inside=1);
matrix IVbisDep14=interpolate(VALL,VDep14,inside=1);
matrix IVbisDep15=interpolate(VALL,VDep15,inside=1);
matrix IVbisDep16=interpolate(VALL,VDep16,inside=1);
matrix IVbisDep17=interpolate(VALL,VDep17,inside=1);
matrix IVbisDep18=interpolate(VALL,VDep18,inside=1);
matrix IVbisDep19=interpolate(VALL,VDep19,inside=1);
matrix IVbisDep20=interpolate(VALL,VDep20,inside=1);
matrix IVbisDep21=interpolate(VALL,VDep21,inside=1);
matrix IVbisDep22=interpolate(VALL,VDep22,inside=1);
matrix IVbisDep23=interpolate(VALL,VDep23,inside=1);
matrix IVbisDep24=interpolate(VALL,VDep24,inside=1);
matrix IVbisDep25=interpolate(VALL,VDep25,inside=1);
matrix IVbisDep26=interpolate(VALL,VDep26,inside=1);
matrix IVbisDep27=interpolate(VALL,VDep27,inside=1);
matrix IVbisDep28=interpolate(VALL,VDep28,inside=1);
matrix IVbisDep29=interpolate(VALL,VDep29,inside=1);
matrix IVbisDep30=interpolate(VALL,VDep30,inside=1);
matrix IVbisDep31=interpolate(VALL,VDep31,inside=1);
matrix IVbisDep32=interpolate(VALL,VDep32,inside=1);
matrix IVbisDep33=interpolate(VALL,VDep33,inside=1);
matrix IVbisDep34=interpolate(VALL,VDep34,inside=1);
matrix IVbisDep35=interpolate(VALL,VDep35,inside=1);
matrix IVbisDep36=interpolate(VALL,VDep36,inside=1);
matrix IVbisDep37=interpolate(VALL,VDep37,inside=1);
matrix IVbisDep38=interpolate(VALL,VDep38,inside=1);
matrix IVbisDep39=interpolate(VALL,VDep39,inside=1);
matrix IVbisDep40=interpolate(VALL,VDep40,inside=1);
matrix IVbisDep41=interpolate(VALL,VDep41,inside=1);
matrix IVbisDep42=interpolate(VALL,VDep42,inside=1);