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

Diff of /MITgcm/pkg/generic_advdiff/gad_diff_y.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_Y
8    
9    C !INTERFACE: ==========================================================
10        SUBROUTINE GAD_DIFF_Y(        SUBROUTINE GAD_DIFF_Y(
11       I           bi,bj,k,       I           bi,bj,k,
12       I           yA, diffKh,       I           yA, diffKh,
13       I           tracer,       I           tracer,
14       O           dfy,       O           dfy,
15       I           myThid )       I           myThid )
       IMPLICIT NONE  
16    
17  C     == GLobal variables ==  C !DESCRIPTION:
18    C Calculates the meridional flux due to down-gradient diffusion of a tracer:
19    C \begin{equation*}
20    C F^y_{diff} = - \kappa_h \partial_y \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    
28  C     == Routine arguments ==  C !INPUT PARAMETERS: ===================================================
29    C  bi,bj                :: tile indices
30    C  k                    :: vertical level
31    C  yA                   :: area of face at V points
32    C  diffKh               :: horizontal diffusivity
33    C  tracer               :: tracer field
34    C  myThid               :: thread number
35        INTEGER bi,bj,k        INTEGER bi,bj,k
36        _RS yA    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RS yA    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
37        _RL diffKh        _RL diffKh
38        _RL tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
       _RL dfy   (1-OLx:sNx+OLx,1-OLy:sNy+OLy)  
39        INTEGER myThid        INTEGER myThid
40    
41  C     == Local variables ==  C !OUTPUT PARAMETERS: ==================================================
42    C  dfx                  :: meridional diffusive flux
43          _RL dfy   (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
44    
45    C !LOCAL VARIABLES: ====================================================
46    C  i,j                  :: loop indices
47        INTEGER i,j        INTEGER i,j
48    CEOP
49    
50        DO i=1-Olx,sNx+Olx        DO i=1-Olx,sNx+Olx
51         dfy(i,1-Oly)=0.         dfy(i,1-Oly)=0.
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22