source: tags/version-1.1/trunk/src/ta_trasbc.h @ 24

/***************************************************************************
                 module classe ta_traadv_musc.h  -  description
***************************************************************************/
// Mohamed Berrada
// locean-ipsl.upmc, Paris, Mars 16, 2009
//===========================================================================
//                          methode forward
forward (YREAL x1,YREAL x2,YREAL x3,YREAL x4,YREAL x5,YREAL x6,YREAL x7)
{/*    1        from    ta_trasbc               1  i    j    k    t
       2        from    musf_zt1_traadv_mus     1  i    j    k    t
       3        from    musf_zt1_traadv_mus     1  i-1  j    k    t
       4        from    musf_zt2_traadv_mus     1  i    j    k    t 
       5        from    musf_zt2_traadv_mus     1  i    j-1  k    t
       6        from    vaf_zt1_traadv_musc     1  i    j    k    t
       7        from    vaf_zt1_traadv_musc     1  i    j    k+1  t*/

  if(Yt==TU)
    YS1=0.;
  else{
    double zbtr1,zta1,zbtr2,zta2;
    if(Yk==NZ-1 || Yj==0 ||  Yj==NY-1 || Yi==0 ||  Yi==NX-1)
      YS1=0.;
    else{
#if defined key_zco
      zbtr1 = 1. / ( e1t(Yi,Yj)*e2t(Yi,Yj) );
#else
      zbtr1 = 1. / ( e1t(Yi,Yj)*e2t(Yi,Yj)*fse3t(Yi,Yj,Yk) );
#endif
      // horizontal advective trends
      zta1 = - zbtr1 * (x2 - x3 + x4-x5);
      
      zbtr2 = 1. / fse3t(Yi,Yj,Yk);
      //! horizontal advective trends
      zta2 = - zbtr2 * ( x6-x7);
      //    ! add it to the general tracer trends
      YS1=x1+zta1+zta2;
    }
  }
  //  
}

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

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

  //===========================================================================
//********************* FIN DU MODULE sa_traadv_musc **********************