forward(YREAL Div_d)
{
  YS1=Div_d + YS1_Div(Yi,Yj,YTemps);
  /*  if (Yi==74 && Yj==35) 
      printf("(%d,%d,%d):%g\n",Yi,Yj,YTemps,YS1);*/
}

backward(YREAL Div_d)
{
  YJ1I1=1;
  /*  if (Yi==74 && Yj==35 && YG1 > 1e-15)
      printf("(%d,%d,%d):%g\n",Yi,Yj,YTemps,YG1);*/
}