pkg/generic_advdiff/gad_grad_x.F

C $Header: /u/gcmpack/models/MITgcmUV/pkg/generic_advdiff/gad_grad_x.F,v 1.1 2001/05/30 19:34:48 adcroft Exp $
C $Name:  $

#include "GAD_OPTIONS.h"

CBOP
C !ROUTINE: GAD_GRAD_X

C !INTERFACE: ==========================================================
      SUBROUTINE GAD_GRAD_X( 
     I           bi,bj,k,
     I           xA,
     I           tracer,
     O           dTdx,
     I           myThid )

C !DESCRIPTION:
C Calculates the zonal gradient of a tracer: $\partial_x \theta$

C !USES: ===============================================================
      IMPLICIT NONE
#include "SIZE.h"
#include "GRID.h"

C !INPUT PARAMETERS: ===================================================
C  bi,bj                :: tile indices
C  k                    :: vertical level
C  tracer               :: tracer field
C  myThid               :: thread number
      INTEGER bi,bj,k
      _RS xA    (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
      _RL tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
      INTEGER myThid

C !OUTPUT PARAMETERS: ==================================================
C  dTdx                 :: zonal gradient
      _RL dTdx  (1-OLx:sNx+OLx,1-OLy:sNy+OLy)

C !LOCAL VARIABLES: ====================================================
C  i,j                  :: loop indices
      INTEGER i,j
CEOP

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.2	C $Header: /u/gcmpack/models/MITgcmUV/pkg/generic_advdiff/gad_grad_x.F,v 1.1 2001/05/30 19:34:48 adcroft Exp $
2			C $Name: $
3	adcroft	1.1
4			#include "GAD_OPTIONS.h"
5
6	adcroft	1.2	CBOP
7			C !ROUTINE: GAD_GRAD_X
8
9			C !INTERFACE: ==========================================================
10	adcroft	1.1	SUBROUTINE GAD_GRAD_X(
11			I bi,bj,k,
12			I xA,
13			I tracer,
14			O dTdx,
15			I myThid )
16	adcroft	1.2
17			C !DESCRIPTION:
18			C Calculates the zonal gradient of a tracer: $\partial_x \theta$
19
20			C !USES: ===============================================================
21	adcroft	1.1	IMPLICIT NONE
22			#include "SIZE.h"
23			#include "GRID.h"
24
25	adcroft	1.2	C !INPUT PARAMETERS: ===================================================
26			C bi,bj :: tile indices
27			C k :: vertical level
28			C tracer :: tracer field
29			C myThid :: thread number
30	adcroft	1.1	INTEGER bi,bj,k
31			_RS xA (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
32			_RL tracer(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
33	adcroft	1.2	INTEGER myThid
34
35			C !OUTPUT PARAMETERS: ==================================================
36			C dTdx :: zonal gradient
37	adcroft	1.1	_RL dTdx (1-OLx:sNx+OLx,1-OLy:sNy+OLy)
38
39	adcroft	1.2	C !LOCAL VARIABLES: ====================================================
40			C i,j :: loop indices
41	adcroft	1.1	INTEGER i,j
42	adcroft	1.2	CEOP
43	adcroft	1.1
44			C o Zonal tracer gradient
45			DO j=1-Oly,sNy+Oly
46			dTdx(1-Olx,j)=0.
47			DO i=1-Olx+1,sNx+Olx
48			dTdx(i,j) = _recip_dxC(i,j,bi,bj)*xA(i,j)
49			& *(tracer(i,j)-tracer(i-1,j))
50			#ifdef COSINEMETH_III
51			& *sqCosFacU(j,bi,bj)
52			#endif
53			ENDDO
54			ENDDO
55
56			RETURN
57			END