source: trunk/src/gcb.h @ 1

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

Import initial
  • Property svn:eol-style set to native
File size: 1.2 KB
Line 
1
2/***************************************************************************
3                 module classe gcb.h  -  description
4***************************************************************************/
5// Mohamed Berrada
6// locean-ipsl.upmc, Paris, April 20, 2009
7//===========================================================================
8//                          methode forward
9forward ()
10{ /*    1..NZ   from    ua_1_dynspg_flt         1  i    j   1..NZ */
11
12  /* YS1  */
13  /*     A 2.118821770102610e-01
14         R 1.866715481224456e-07*/ 
15  /* YS1  v clonol*/
16  /*     A 2.102448613964121e-01
17         R 1.867828604113514e-07*/
18
19  if(Yt==TU){
20    YS1=0.;
21  }
22  else{
23    if(Yi==0 || Yi==NX-1 || Yj==0 || Yj==NY-1  ){
24      YS1=0.;//lbc_lnk boundaries
25    }
26    else{
27      double sum=0.;
28      for(int k=0;k<NZ-1;k++)
29        sum=sum+ fse3u(Yi,Yj,k) * Yting[k];
30      YS1=sum * e2u(Yi,Yj); 
31    }
32  }
33  //
34}
35
36//===========================================================================
37//                         methode  backward
38
39backward ()
40{
41 
42}
43
44  //===========================================================================
45//*************************  FIN DU MODULE gcb   *************************
46
Note: See TracBrowser for help on using the repository browser.