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

Diff of /MITgcm/pkg/generic_advdiff/gad_biharm_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_BIHARM_Y
8    
9    C !INTERFACE: ==========================================================
10        SUBROUTINE GAD_BIHARM_Y(        SUBROUTINE GAD_BIHARM_Y(
11       I           bi,bj,k,       I           bi,bj,k,
12       I           yA,del2T,diffK4,       I           yA,del2T,diffK4,
13       O           dfy,       O           dfy,
14       I           myThid )       I           myThid )
       IMPLICIT NONE  
15    
16  C     == GLobal variables ==  C !DESCRIPTION:
17    C Calculates the meridional flux due to bi-harmonic diffusion of a tracer.
18    C Routine takes the laplacian of the tracer as argument and calculates
19    C the meridional gradient:
20    C \begin{equation*}
21    C F^y_{diff} = \kappa_4 \partial_y \nabla^2 \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    
29  C     == Routine arguments ==  C !INPUT PARAMETERS: ===================================================
30    C  bi,bj                :: tile indices
31    C  k                    :: vertical level
32    C  yA                   :: area of face at V points
33    C  del2T                :: Laplacian of tracer
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 del2T(1-OLx:sNx+OLx,1-OLy:sNy+OLy)        _RL del2T(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
38        _RL diffK4        _RL diffK4
       _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  dfy                  :: 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  C     Difference of meridional fluxes ...  C     Difference of meridional fluxes ...
51        DO i=1-Olx,sNx+Olx        DO i=1-Olx,sNx+Olx
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22