EvanescentBoundaryJ
Dans cette fonction, dans la première boucle:
for (int i=0; i<nr; i++){
real omega = sqrt(G*1.0/Rmed[i]/Rmed[i]/Rmed[i]);
for (h=0; h<ni; h++){
int k = h+i*ni;
if(Stellar || Radiative || grid.physic().adiabaticEoS ){
int j = 0;
int l = j +h*ns + i*ni*ns;
dens[l] = Densj0[k];
vrad[l] = Vradj0[k];
vphi[l] = Vphij0[k];
vtheta[l] = Vthetaj0[k];
j=ns-1;
l = j +h*ns + i*ni*ns;
dens[l] = Densjs[k];
vrad[l] = Vradjs[k];
vphi[l] = Vphijs[k];
vtheta[l] = Vthetajs[k];
} else {
int j=0;
int l = j +h*ns + i*ni*ns;
dens[l] = SigmaMed[k];
vtheta[l]= omega*Rmed[i]*sin(phiMed[h])*sqrt(1.0-(pow(aspectRatio[i],2.0)*(2.+SIGMASLOPE)));
vtheta[l] -= Rmed[i]*sin(phiMed[h])*OmegaFrame;
if (ViscosityAlpha) {
vrad[l] = -3.0*visc[k]/Rmed[i]*(-SIGMASLOPE+2.0*grid.physic().flaringIndex+1.0);
} else {
vrad[l] = -3.0*visc[k]/Rmed[i]*(-SIGMASLOPE+.5);
}
vphi[l] = 0.;
}
dens[l+ns-1]= dens[l]; **// PROBLÈME ICI **
vtheta[l+ns-1]=vtheta[l];
vrad[l+ns-1] = vrad[l] ;
vphi[l+ns-1] = 0.;
}
}
l
n'est pas clairement définie à la fin