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

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

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

revision 1.5 by jmc, Wed Jun 22 00:27:47 2005 UTC revision 1.6 by jmc, Thu Oct 20 17:03:08 2005 UTC
# Line 4  C $Name$ Line 4  C $Name$
4  #include "GAD_OPTIONS.h"  #include "GAD_OPTIONS.h"
5    
6  CBOP  CBOP
7  C     !ROUTINE: GAD_FLUXLIMIT_IMPL_R  C     !ROUTINE: GAD_U3C4_IMPL_R
8  C     !INTERFACE:  C     !INTERFACE:
9        SUBROUTINE GAD_U3C4_IMPL_R(        SUBROUTINE GAD_U3C4_IMPL_R(
10       I           bi,bj,k, iMin,iMax,jMin,jMax,       I           bi,bj,k, iMin,iMax,jMin,jMax,
# Line 13  C     !INTERFACE: Line 13  C     !INTERFACE:
13       I           myThid )       I           myThid )
14    
15  C     !DESCRIPTION:  C     !DESCRIPTION:
16    C     Compute matrix element to solve vertical advection implicitly
17  C     Compute matrix element to solve vertical advection  C      using 3rd order upwind advection scheme,
18  C     \begin{enumerate}  C         or 4th order Centered advection scheme.
19  C     \item implicitly using 3rd order upwind, or  C     Method:
20  C     \item 4th order Centered advection schemes.  C      contribution of vertical transport at interface k is added
21  C     \end{enumerate}  C      to matrix lines k and k-1
 C     Also, the contribution of vertical transport at interface k  
 C     is added to matrix lines k and k-1  
22    
23  C     !USES:  C     !USES:
24        IMPLICIT NONE        IMPLICIT NONE
# Line 32  C     == Global variables === Line 30  C     == Global variables ===
30  #include "PARAMS.h"  #include "PARAMS.h"
31  #include "GAD.h"  #include "GAD.h"
32    
33  C !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
34  C == Routine Arguments ==  C     == Routine Arguments ==
35  C bi,bj           :: tile indices  C     bi,bj           :: tile indices
36  C k               :: vertical level  C     k               :: vertical level
37  C iMin,iMax       :: computation domain  C     iMin,iMax       :: computation domain
38  C jMin,jMax       :: computation domain  C     jMin,jMax       :: computation domain
39  C advectionScheme :: advection scheme to use  C     advectionScheme :: advection scheme to use
40  C deltaTarg       :: time step  C     deltaTarg       :: time step
41  C rTrans          :: vertical volume transport  C     rTrans          :: vertical volume transport
42  C tFld            :: tracer field  C     tFld            :: tracer field
43  C a5d             :: 2nd  lower diag of pentadiagonal matrix  C     a5d             :: 2nd  lower diag of pentadiagonal matrix
44  C b5d             :: 1rst lower diag of pentadiagonal matrix  C     b5d             :: 1rst lower diag of pentadiagonal matrix
45  C c5d             :: main diag       of pentadiagonal matrix  C     c5d             :: main diag       of pentadiagonal matrix
46  C d5d             :: 1rst upper diag of pentadiagonal matrix  C     d5d             :: 1rst upper diag of pentadiagonal matrix
47  C e5d             :: 2nd  upper diag of pentadiagonal matrix  C     e5d             :: 2nd  upper diag of pentadiagonal matrix
48  C myThid          :: thread number  C     myThid          :: thread number
49        INTEGER bi,bj,k        INTEGER bi,bj,k
50        INTEGER iMin,iMax,jMin,jMax        INTEGER iMin,iMax,jMin,jMax
51        INTEGER advectionScheme        INTEGER advectionScheme
# Line 60  C myThid          :: thread number Line 58  C myThid          :: thread number
58        _RL e5d   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL e5d   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
59        INTEGER myThid        INTEGER myThid
60    
61  C == Local Variables ==  C     == Local Variables ==
62  C i,j             :: loop indices  C     i,j             :: loop indices
63  C kp1             :: =min( k+1 , Nr )  C     kp1             :: =min( k+1 , Nr )
64  C km2             :: =max( k-2 , 1 )  C     km2             :: =max( k-2 , 1 )
65  C rCenter         :: centered contribution  C     rCenter         :: centered contribution
66  C rUpwind         :: upwind   contribution  C     rUpwind         :: upwind   contribution
67        LOGICAL flagC4        LOGICAL flagC4
68        INTEGER i,j,kp1,km2        INTEGER i,j,kp1,km2
69        _RL rCenter, rUpwind        _RL rCenter, rUpwind
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22