source: trunk/src/zhpi_dyn_hpg_zco.h @ 64

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

  /* YS1 */
  /*     A 
         R */
  if(Yt==TU)
    YS1=0.;
  else{
    double zcoef0 = - grav * 0.5;
    double zcoef1 = zcoef0 * fse3w(Yi,Yj,Yk);
    if(Yi==0 || Yi==NX-1 || Yj==0 || Yj==NY-1 || Yk==NZ-1){
      YS1=0.;
    }
    else{
      if(Yk==0)
        YS1=zcoef1 * ( x1 - x2)/ e1u(Yi,Yj);
      else{
        YS1=x5+zcoef1*((x1+x3)-(x2+x4)) / e1u(Yi,Yj);
      }
    }
  }
  //
}

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

backward (YREAL x1,YREAL x2,YREAL x3,YREAL x4,YREAL x5)
{
  
}

  //===========================================================================
//*************************  FIN DU MODULE zhpi_dyn_hpg_zco   *************************