source: trunk/src/gcx_dynspg_flt.h @ 49

/***************************************************************************
                 module classe gcx_dynspg_flt.h  -  description
***************************************************************************/
// Mohamed Berrada
// locean-ipsl.upmc, Paris, April 20, 2009
//===========================================================================
//                          methode forward
forward (YREAL x1,YREAL x2)
{ /*    1       from    gcx2            1  i    j         t-1
        2       from    gcx2            1  i    j         t-2*/

  /* YS1  */
  /*     A 
         R */ 

  if(Yt<=TU){
    YS1=0.;
  }
  else{
    if(Yi==0 || Yi==NX-1 || Yj==0 || Yj==NY-1  ){
      YS1=0.;
    }
    else{
      //YS1 = 2. * x1   - x2;
      YS1 = 2. * x1   - x2;
    }
  }
  //if(Yi==28 && Yj==3)  xtest("gcx_dynspg_flt",YS1,Yi,Yj,Yk,Yt);

  //
}

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

backward (YREAL x1,YREAL x2)
{
  YJ1I1=0.;     YJ1I2=0.;
  if(Yt>TU){  
    if(Yi==0 || Yi==NX-1 || Yj==0 || Yj==NY-1  ){
      YJ1I1=0.;     YJ1I2=0.;
    }
    else{
      YJ1I1=2.;     YJ1I2=-1.;
    }
  }
  //  
}

//===========================================================================
//*********************  FIN DU MODULE gcx_dynspg_flt   *********************