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

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

Import initial

  • Property svn:eol-style set to native
File size: 1.4 KB
Line 
1
2/***************************************************************************
3                 module classe zcent_ver_traadv_cen2.h  -  description
4***************************************************************************/
5// Mohamed Berrada
6// locean-ipsl.upmc, Paris, October 20, 2009
7//===========================================================================
8//                          methode forward
9forward (YREAL x1,YREAL x2,YREAL x3)
10{ /*
11        1       from    wa                      1  i    j    k    t-1
12        2       from    ta                      1  i    j    k    t-1
13        3       from    ta                      1  i    j    k-1  t-1
14  */
15
16  if(Yt==TU)
17    YS1=0.;
18  else{
19    if(Yi==0 || Yi==NX-1 || Yj==0 || Yj==NY-1 || Yk==0 || Yk==NZ-1){
20      YS1=0.;
21      if(Yk==0) YS1= x1*x2;
22    }
23    else{
24      YS1=0.5*x1*(x2+x3);
25    }
26  }
27  //
28}
29
30//===========================================================================
31//                         methode  backward
32
33backward (YREAL x1,YREAL x2,YREAL x3)
34{
35  YJ1I1=0.;     YJ1I2=0.;     YJ1I3=0.; 
36
37  if(Yt!=TU){ 
38    if(Yi==0 || Yi==NX-1 || Yj==0 || Yj==NY-1 || Yk==0 || Yk==NZ-1){
39      YJ1I1=0.;     YJ1I2=0.;     YJ1I3=0.; 
40      if(Yk==0){
41        YJ1I1=x2;     YJ1I2=x1;
42      }
43    }
44    else{
45      YJ1I1= 0.5*(x2+x3); 
46      YJ1I2= 0.5*x1;
47      YJ1I3= 0.5*x1;
48    }
49  }
50  //
51 
52 
53}
54
55  //===========================================================================
56//*************************  FIN DU MODULE zcent_ver_traadv_cen2   *************************
57
Note: See TracBrowser for help on using the repository browser.