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

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

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

revision 1.2 by adcroft, Wed Mar 24 15:20:35 2004 UTC revision 1.3 by edhill, Mon Mar 29 03:33:51 2004 UTC
# Line 12  C     !INTERFACE: Line 12  C     !INTERFACE:
12       O           a3d, b3d, c3d,       O           a3d, b3d, c3d,
13       I           myThid )       I           myThid )
14    
15  C     !DESCRIPTION: \bv  C     !DESCRIPTION:
16  C     *==========================================================*  C     Compute matrix element to solve vertical advection implicitly
17  C     | S/R GAD_FLUXLIMIT_IMPL_R  C     using flux--limiter advection scheme.  The contribution of
18  C     | o Compute matrix element to solve vertical advection  C     vertical transport at interface k is added to matrix lines k and
19  C     |   implicitly using flux-limiter advection scheme  C     k-1.
 C     *==========================================================*  
 C     | o contribution of vertical transport at interface k  
 C     |   is added to matrix lines k & k-1  
 C     *==========================================================*  
 C     \ev  
20    
21  C     !USES:  C     !USES:
22        IMPLICIT NONE        IMPLICIT NONE
# Line 34  C     == Global variables === Line 29  C     == Global variables ===
29    
30  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
31  C     == Routine Arguments ==  C     == Routine Arguments ==
32  C     bi,bj                :: tile indices  C     bi,bj        :: tile indices
33  C     k                    :: vertical level  C     k            :: vertical level
34  C     iMin,iMax,jMin,jMax  :: computation domain  C     iMin,iMax    :: computation domain
35  C     deltaTarg            :: time step  C     jMin,jMax    :: computation domain
36  C     rTrans               :: vertical volume transport  C     deltaTarg    :: time step
37  C     tFld                 :: tracer field  C     rTrans       :: vertical volume transport
38  C     a3d                  :: lower diagonal of the tridiagonal matrix  C     tFld         :: tracer field
39  C     b3d                  :: main  diagonal of the tridiagonal matrix  C     a3d          :: lower diagonal of the tridiagonal matrix
40  C     c3d                  :: upper diagonal of the tridiagonal matrix  C     b3d          :: main  diagonal of the tridiagonal matrix
41  C     myThid               :: thread number  C     c3d          :: upper diagonal of the tridiagonal matrix
42    C     myThid       :: thread number
43        INTEGER bi,bj,k        INTEGER bi,bj,k
44        INTEGER iMin,iMax,jMin,jMax        INTEGER iMin,iMax,jMin,jMax
45        _RL deltaTarg        _RL deltaTarg
# Line 55  C     myThid               :: thread num Line 51  C     myThid               :: thread num
51        INTEGER myThid        INTEGER myThid
52    
53  C     == Local Variables ==  C     == Local Variables ==
54  C     i,j                  :: loop indices  C     i,j          :: loop indices
55  C     kp1                  :: =min( k+1 , Nr )  C     kp1          :: =min( k+1 , Nr )
56  C     km1                  :: =max( k-1 , 1 )  C     km1          :: =max( k-1 , 1 )
57  C     km2                  :: =max( k-2 , 1 )  C     km2          :: =max( k-2 , 1 )
58  C     Cr                   :: slope ratio  C     Cr           :: slope ratio
59  C     Rjm,Rj,Rjp           :: differences at i-1,i,i+1  C     Rjm,Rj,Rjp   :: differences at i-1,i,i+1
60  C     w_CFL                :: Courant-Friedrich-Levy number  C     w_CFL        :: Courant-Friedrich-Levy number
61  C     upwindFac            :: upwind factor  C     upwindFac    :: upwind factor
62  C     rCenter              :: centered contribution  C     rCenter      :: centered contribution
63  C     rUpwind              :: upwind   contribution  C     rUpwind      :: upwind   contribution
64        INTEGER i,j,kp1,km1,km2        INTEGER i,j,kp1,km1,km2
65        _RL Cr,Rjm,Rj,Rjp, w_CFL        _RL Cr,Rjm,Rj,Rjp, w_CFL
66        _RL upwindFac(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL upwindFac(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22