source: trunk/src/zwt_trazdf_imp.h @ 64

/***************************************************************************
                 module classe zwt_trazdf_imp.h  -  description
***************************************************************************/
// Mohamed Berrada
// locean-ipsl.upmc, Paris, April 17, 2009
//===========================================================================
//                          methode forward
forward (YREAL x1,YREAL x2,YREAL x3,YREAL x4)
{ /*    1       from    zwisd_trazdf_imp        1  i    j    k    t
        2       from    zwisd_trazdf_imp        2  i    j    k-1  t
        3       from    zwisd_trazdf_imp        3  i    j    k    t
        4       from    zwt_trazdf_imp          1  i    j    k-1  t*/
  /* YS1 */
  /*      0
          0*/

  /* if (Yi==30 && Yj==1 && Yt==16)
  {
      printf("Yk=%i,x1=%f,x2=%f,x3=%f,x4=%f\n",Yk,x1,x2,x3,x4);
      }  */
  if(Yt==TU)
    YS1=0.;
  else{
    if(Yi==0 || Yi==NX-1 || Yj==0 || Yj==NY-1){
      YS1=0.;
    }
    else{
      if(Yk==0)
        YS1= x3;
      else
        YS1=x3-x1*x2/x4;
    }
  }
  //
}

//===========================================================================
//                         methode  backward

backward (YREAL x1,YREAL x2,YREAL x3,YREAL x4)
{
  YJ1I1=0.;     YJ1I2=0.;   YJ1I3=0.;     YJ1I4=0.;
  if(Yt!=TU){
    if(Yi==0 || Yi==NX-1 || Yj==0 || Yj==NY-1){
      YJ1I1=0.;     YJ1I2=0.;   YJ1I3=0.;     YJ1I4=0.;
    }
    else{
      if(Yk==0)
        YJ1I3=1.;
      else{
        YJ1I1 = -x2/x4;
        YJ1I2 = -x1/x4;
        YJ1I3 =  1.;
        YJ1I4=   x1*x2/x4/x4;
      }
    }
  }
  //
}

//===========================================================================
//*********************  FIN DU MODULE zwt_trazdf_imp ***********************