source: trunk/src/ua_1_dynspg_flt.h

/***************************************************************************
                 module classe ua_1_dynspg_flt.h  -  description
***************************************************************************/
// Mohamed Berrada
// locean-ipsl.upmc, Paris, April 20, 2009
//===========================================================================
//                          methode forward
forward (YREAL x1,YREAL x2,YREAL x3)
{ /*    1       from    ua_norm_dynzdf          1  i    j    k    t
        2       from    ub                      1  i    j    k    t-1
        3       from    spguv_now_dynspg        1  i    j         t*/


  /* YS1  */
  /*     A 1.982166090909616e-07
         R 1.866715481214184e-07*/
  /* YS1  va clonol*/
  /*     A 1.966848939039096e-07
         R 1.867828604129974e-07*/
  if(Yt==TU)
    YS1=0.;
  else{
    if(Yi==0 || Yi==NX-1 || Yj==0 || Yj==NY-1 || Yk==NZ-1  ){
      YS1=0.;
    }
    else{
      YS1 = ( x2 + r2dt * (x1+x3) ) * umask(Yi,Yj,Yk);
    }
  }
  //
}

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

backward (YREAL x1,YREAL x2,YREAL x3)
{
  YJ1I1=0.;     YJ1I2=0.;     YJ1I3=0.; 
  // if(Yt==6 && Yi==28 && Yj==3)  xtest("ua_1_dynspg_flt",YG1,Yi,Yj,Yk,Yt);
  if(Yt!=TU){  
    if(Yi==0 || Yi==NX-1 || Yj==0 || Yj==NY-1 || Yk==NZ-1  ){
      YJ1I1=0.;     YJ1I2=0.;     YJ1I3=0.; 
    }
    else{
      YJ1I1=r2dt*umask(Yi,Yj,Yk);     
      YJ1I2=umask(Yi,Yj,Yk);
      YJ1I3=YJ1I1; 
    }
  }
  //
  
}

//===========================================================================
//********************  FIN DU MODULE ua_1_dynspg_flt   *********************