source: tags/version-1.1/trunk/src/sum_spgv_dynspg_flt.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.3 KB
Line 
1
2/***************************************************************************
3                 module classe sum_spgv_dynspg_flt.h  -  description
4***************************************************************************/
5// Mohamed Berrada
6// locean-ipsl.upmc, Paris, April 20, 2009
7//===========================================================================
8//                          methode forward
9forward ()
10{ /*    1..NZ   from    va_1_dynspg_flt         1  i    j   1..NZ */
11
12  /* YS1  */
13  /*     A
14         R */
15  if(Yt==TU)
16    YS1=0.;
17  else{
18    if(Yi==0 || Yi==NX-1 || Yj==0 || Yj==NY-1  ){
19      YS1=0.;
20    }
21    else{
22      double sum=0.;
23      for(int k=0;k<NZ-1;k++)
24        sum=sum+ fse3v(Yi,Yj,k) * Yting[k];
25      YS1=sum * e1v(Yi,Yj); 
26    }
27  }
28  //
29}
30
31//===========================================================================
32//                         methode  backward
33
34backward ()
35{
36  for(int k=0;k<NZ;k++)  Yjac[0][k]    =0.;
37  if(Yt!=TU){ 
38    if(Yi==0 || Yi==NX-1 || Yj==0 || Yj==NY-1  ){
39      for(int k=0;k<NZ;k++)  Yjac[0][k]    =0.;
40    }
41    else{
42      for(int k=0;k<NZ-1;k++)  Yjac[0][k]    =fse3v(Yi,Yj,k)*e1v(Yi,Yj);
43    }
44  }
45  //
46 
47}
48
49//===========================================================================
50//******************  FIN DU MODULE sum_spgv_dynspg_flt   *******************
51
Note: See TracBrowser for help on using the repository browser.