/[MITgcm]/MITgcm/pkg/generic_advdiff/gad_c2_adv_r.F
ViewVC logotype

Diff of /MITgcm/pkg/generic_advdiff/gad_c2_adv_r.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.2 by jmc, Thu Jul 12 00:30:27 2001 UTC revision 1.3 by adcroft, Thu Sep 20 20:07:06 2001 UTC
# Line 3  C $Name$ Line 3  C $Name$
3    
4  #include "GAD_OPTIONS.h"  #include "GAD_OPTIONS.h"
5    
6    CBOP
7    C !ROUTINE: GAD_C2_ADV_R
8    
9    C !INTERFACE: ==========================================================
10        SUBROUTINE GAD_C2_ADV_R(        SUBROUTINE GAD_C2_ADV_R(
11       I           bi,bj,k,       I           bi,bj,k,
12       I           rTrans,       I           rTrans,
13       I           tracer,       I           tracer,
14       O           wT,       O           wT,
15       I           myThid )       I           myThid )
 C     /==========================================================\  
 C     | SUBROUTINE GAD_C2_ADV_R                                  |  
 C     | o Compute vertical advective Flux of Tracer using        |  
 C     |   2nd Order Centered Scheme                              |  
 C     |==========================================================|  
       IMPLICIT NONE  
16    
17  C     == GLobal variables ==  C !DESCRIPTION:
18    C Calculates the vertical flux due to advection of a tracer using
19    C centered second-order interpolation:
20    C \begin{equation*}
21    C F^r_{adv} = w \theta
22    C \end{equation*}
23    
24    C !USES: ===============================================================
25          IMPLICIT NONE
26  #include "SIZE.h"  #include "SIZE.h"
27  #include "GRID.h"  #include "GRID.h"
28  #include "EEPARAMS.h"  #include "EEPARAMS.h"
29  #include "PARAMS.h"  #include "PARAMS.h"
30    
31  C     == Routine arguments ==  C !INPUT PARAMETERS: ===================================================
32    C  bi,bj                :: tile indices
33    C  k                    :: vertical level
34    C  rTrans               :: vertical volume transport
35    C  tracer               :: tracer field
36    C  myThid               :: thread number
37        INTEGER bi,bj,k        INTEGER bi,bj,k
38        _RL rTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL rTrans(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
39        _RL tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)        _RL tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
       _RL wT    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)  
40        INTEGER myThid        INTEGER myThid
41    
42  C     == Local variables ==  C !OUTPUT PARAMETERS: ==================================================
43    C  wT                   :: vertical advective flux
44          _RL wT    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
45    
46    C !LOCAL VARIABLES: ====================================================
47    C  i,j                  :: loop indices
48    C  km1                  :: =k+1 for k<Nr, =Nr for k>=Nr
49        INTEGER i,j,km1        INTEGER i,j,km1
50    CEOP
51    
52        km1=max(1,k-1)        km1=max(1,k-1)
53    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22