/[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.3 by adcroft, Fri Sep 21 13:11:43 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 area integrated vertical flux due to down-gradient
19    C diffusion of a tracer:
20    C \begin{equation*}
21    C F^r_{diff} = - A^r \kappa_r \frac{1}{\Delta r_c} \delta_k \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  KappaR               :: vertical diffusivity
35    C  tracer               :: tracer field
36    C  myThid               :: thread number
37        INTEGER bi,bj,k        INTEGER bi,bj,k
38        _RL KappaR(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)        _RL KappaR(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
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 dfr   (1-OLx:sNx+OLx,1-OLy:sNy+OLy)  
40        INTEGER myThid        INTEGER myThid
41    
42  C     == Local variables ==  C !OUTPUT PARAMETERS: ==================================================
43    C  dfr                  :: vertical diffusive flux
44          _RL dfr   (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        LOGICAL TOP_LAYER  CEOP
51    
52        km1=max(1,k-1)        km1=max(1,k-1)
53    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22