pkg/generic_advdiff/gad_grad_x.F

C $Header: $
C $Name: $

#include "GAD_OPTIONS.h"

      SUBROUTINE GAD_GRAD_X( 
     I           bi,bj,k,
     I           xA,
     I           tracer,
     O           dTdx,
     I           myThid )
      IMPLICIT NONE

C     == GLobal variables ==
#include "SIZE.h"
#include "GRID.h"

C     == Routine arguments ==
      INTEGER bi,bj,k
      _RS xA    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
      _RL tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
      _RL dTdx  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
      INTEGER myThid

C     == Local variables ==
      INTEGER i,j

C     o Zonal tracer gradient
      DO j=1-Oly,sNy+Oly
       dTdx(1-Olx,j)=0.
       DO i=1-Olx+1,sNx+Olx
        dTdx(i,j) = _recip_dxC(i,j,bi,bj)*xA(i,j)
     &               *(tracer(i,j)-tracer(i-1,j))
#ifdef COSINEMETH_III
     &               *sqCosFacU(j,bi,bj)
#endif
       ENDDO
      ENDDO

      RETURN
      END
1	adcroft	1.1	C $Header: $
2			C $Name: $
3
4			#include "GAD_OPTIONS.h"
5
6			SUBROUTINE GAD_GRAD_X(
7			I bi,bj,k,
8			I xA,
9			I tracer,
10			O dTdx,
11			I myThid )
12			IMPLICIT NONE
13
14			C == GLobal variables ==
15			#include "SIZE.h"
16			#include "GRID.h"
17
18			C == Routine arguments ==
19			INTEGER bi,bj,k
20			_RS xA (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
21			_RL tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
22			_RL dTdx (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
23			INTEGER myThid
24
25			C == Local variables ==
26			INTEGER i,j
27
28			C o Zonal tracer gradient
29			DO j=1-Oly,sNy+Oly
30			dTdx(1-Olx,j)=0.
31			DO i=1-Olx+1,sNx+Olx
32			dTdx(i,j) = _recip_dxC(i,j,bi,bj)*xA(i,j)
33			& *(tracer(i,j)-tracer(i-1,j))
34			#ifdef COSINEMETH_III
35			& *sqCosFacU(j,bi,bj)
36			#endif
37			ENDDO
38			ENDDO
39
40			RETURN
41			END