/***************************************************************************
                 module classe rhd.h  -  description
***************************************************************************/
// Mohamed Berrada
// locean-ipsl.upmc, Paris, April 17, 2009
//===========================================================================
//                          methode forward
forward (YREAL x1,YREAL x2)
{ /*	1	from	tb			1  i 	j    k    t
	2	from	sb			1  i 	j    k    t*/
  /* YS1  */
  /*     7.835114356366673e-16
	 2.331434299545338e-16*/
  if(Yk==NZ-1){
    YS1=0.;
  }
  else{
    YS1 = (rbeta* x2 - ralpha * x1 ) * tmask(Yi,Yj,Yk);
  }
  //
}

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

backward (YREAL x1,YREAL x2)
{
  YJ1I1=0.;  YJ1I2=0.;
  if(Yk==NZ-1){
    YJ1I1=0.;  YJ1I2=0.;
  }
  else{
    YJ1I1 =  - ralpha * tmask(Yi,Yj,Yk);
    YJ1I2 =    rbeta * tmask(Yi,Yj,Yk);
  }
  //
  
}

//===========================================================================
//**************************  FIN DU MODULE rhd   ***************************