source: trunk/src/zsuv_traldf_lap.h @ 65

Last change on this file since 65 was 1, checked in by jbrlod, 13 years ago

Import initial
  • Property svn:eol-style set to native
File size: 1.5 KB
RevLine 
[1]1/***************************************************************************
2                 module classe zsuv_traldf_lap.h  -  description
3***************************************************************************/
4// Mohamed Berrada
5// locean-ipsl.upmc, Paris, April 2, 2009
6//===========================================================================
7//                          methode forward
8forward (YREAL x1, YREAL x2, YREAL x3)
9{ /*    1       from    sa              1  i    j    k    t-2
10        2       from    sa              1  i+1  j    k    t-2
11        3       from    sa              1  i    j+1  k    t-2*/
12  if(Yt==TU)
13    YS1=0.;
14  else{
15    if(Yi==NX-1 || Yj==NY-1 || Yk== NZ-1)
16      YS1=0.;
17    else{
18#if defined key_zco
19      double  zabe1 = fsahtu(Yi,Yj,Yk) * umask(Yi,Yj,Yk) * ze1ur(Yi,Yj);
20      double  zabe2 = fsahtv(Yi,Yj,Yk) * vmask(Yi,Yj,Yk) * ze2vr(Yi,Yj);
21#else
22      double  zabe1 = fsahtu(Yi,Yj,Yk) * umask(Yi,Yj,Yk) * ze1ur(Yi,Yj) * fse3u(Yi,Yj,Yk);
23      double  zabe2 = fsahtv(Yi,Yj,Yk) * vmask(Yi,Yj,Yk) * ze2vr(Yi,Yj) * fse3v(Yi,Yj,Yk);
24#endif
25      YS1 = zabe1*(x2-x1);
26      //ztu(ji,jj,jk) = zabe1 * ( tb(ji+1,jj  ,jk) - tb(ji,jj,jk) );
27      YS2 = zabe2*(x3-x1);
28      //  ztv(ji,jj,jk) = zabe2 * ( tb(ji  ,jj+1,jk) - tb(ji,jj,jk) );
29    }
30  }
31  //
32}
33//===========================================================================
34//                         methode  backward
35
36backward (YREAL x1, YREAL x2, YREAL x3)
37{
38 
39}
40
41  //===========================================================================
42//*************************  FIN DU MODULE zsuv_traldf_lap   *************************
43
Note: See TracBrowser for help on using the repository browser.