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

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

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

revision 1.1 by adcroft, Wed May 30 19:34:48 2001 UTC revision 1.2 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_DIFF_R
8    
9    C !INTERFACE: ==========================================================
10        SUBROUTINE GAD_DIFF_R(        SUBROUTINE GAD_DIFF_R(
11       I           bi,bj,k,       I           bi,bj,k,
12       I           KappaR,       I           KappaR,
13       I           tracer,       I           tracer,
14       O           dfr,       O           dfr,
15       I           myThid )       I           myThid )
       IMPLICIT NONE  
16    
17  C     == GLobal variables ==  C !DESCRIPTION:
18    C Calculates the vertical flux due to down-gradient diffusion of a tracer:
19    C \begin{equation*}
20    C F^r_{diff} = - \kappa_r \partial_r \theta
21    C \end{equation*}
22    
23    C !USES: ===============================================================
24          IMPLICIT NONE
25  #include "SIZE.h"  #include "SIZE.h"
26  #include "GRID.h"  #include "GRID.h"
27  #include "EEPARAMS.h"  #include "EEPARAMS.h"
28  #include "PARAMS.h"  #include "PARAMS.h"
29    
30  C     == Routine arguments ==  C !INPUT PARAMETERS: ===================================================
31    C  bi,bj                :: tile indices
32    C  k                    :: vertical level
33    C  KappaR               :: vertical diffusivity
34    C  tracer               :: tracer field
35    C  myThid               :: thread number
36        INTEGER bi,bj,k        INTEGER bi,bj,k
37        _RL KappaR(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL KappaR(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
38        _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 dfr   (1-OLx:sNx+OLx,1-OLy:sNy+OLy)  
39        INTEGER myThid        INTEGER myThid
40    
41  C     == Local variables ==  C !OUTPUT PARAMETERS: ==================================================
42    C  dfr                  :: vertical diffusive flux
43          _RL dfr   (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
44    
45    C !LOCAL VARIABLES: ====================================================
46    C  i,j                  :: loop indices
47    C  km1                  :: =k+1 for k<Nr, =Nr for k>-Nr
48        INTEGER i,j,km1        INTEGER i,j,km1
       LOGICAL TOP_LAYER  
49    
50        km1=max(1,k-1)        km1=max(1,k-1)
51    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22