pkg/autodiff/ini_autodiff.F

C $Header: /u/gcmpack/MITgcm/pkg/autodiff/ini_autodiff.F,v 1.7 2004/09/17 23:02:01 heimbach Exp $
C $Name:  $

#include "AUTODIFF_OPTIONS.h"

CBOP
C     !ROUTINE: INI_AUTODIFF
C     !INTERFACE:
      SUBROUTINE INI_AUTODIFF( myThid )
C     !DESCRIPTION: \bv
C     *==========================================================*
C     | SUBROUTINE INI_AUTODIFF                                    
C     | o Initialise to zero some active arrays  
C     |   These ini                
C     *==========================================================*
C     \ev

C     !USES:
      IMPLICIT NONE
C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "FFIELDS.h"
#ifdef EXACT_CONSERV
# include "SURFACE.h"
#endif

C     !INPUT/OUTPUT PARAMETERS:
C     == Routine arguments ==
C     myThid -  Number of this instance of INI_UVEL 
      INTEGER myThid

#ifdef ALLOW_AUTODIFF_TAMC

C     !LOCAL VARIABLES:
C     == Local variables ==
C     bi,bj  - Loop counters
C     I,J,K
      INTEGER bi, bj
      INTEGER  I, J, K
CEOP

C--   Over all tiles
      DO bj = myByLo(myThid), myByHi(myThid)
       DO bi = myBxLo(myThid), myBxHi(myThid)

C-      3D arrays
        DO K=1,Nr
         DO J=1-Oly,sNy+Oly
          DO I=1-Olx,sNx+Olx
#ifdef ALLOW_EP_FLUX
           EfluxY(i,j,k,bi,bj)=0.
           EfluxP(i,j,k,bi,bj)=0.
#endif
          ENDDO
         ENDDO
        ENDDO

C-      2D arrays
        DO J=1-Oly,sNy+Oly
         DO I=1-Olx,sNx+Olx
#ifdef EXACT_CONSERV
         dEtaHdt(i,j,bi,bj) = 0. _d 0
         PmEpR(i,j,bi,bj) = 0. _d 0
cph          hDivFlow(I,J,bi,bj)=0. _d 0
#endif
#ifdef ALLOW_EXF
          EmPmR(I,J,bi,bj) = 0. _d 0
#endif
         ENDDO
        ENDDO

       ENDDO
      ENDDO

#endif /* ALLOW_AUTODIFF_TAMC */

      END
1	heimbach	1.8	C $Header: /u/gcmpack/MITgcm/pkg/autodiff/ini_autodiff.F,v 1.7 2004/09/17 23:02:01 heimbach Exp $
2	edhill	1.5	C $Name: $
3	heimbach	1.1
4	edhill	1.4	#include "AUTODIFF_OPTIONS.h"
5	heimbach	1.1
6			CBOP
7			C !ROUTINE: INI_AUTODIFF
8			C !INTERFACE:
9			SUBROUTINE INI_AUTODIFF( myThid )
10			C !DESCRIPTION: \bv
11			C ==========================================================
12			C \| SUBROUTINE INI_AUTODIFF
13	heimbach	1.3	C \| o Initialise to zero some active arrays
14			C \| These ini
15	heimbach	1.1	C ==========================================================
16			C \ev
17
18			C !USES:
19			IMPLICIT NONE
20			C === Global variables ===
21			#include "SIZE.h"
22			#include "EEPARAMS.h"
23			#include "PARAMS.h"
24	heimbach	1.8	#include "FFIELDS.h"
25	heimbach	1.3	#ifdef EXACT_CONSERV
26			# include "SURFACE.h"
27			#endif
28	heimbach	1.1
29			C !INPUT/OUTPUT PARAMETERS:
30			C == Routine arguments ==
31			C myThid - Number of this instance of INI_UVEL
32			INTEGER myThid
33
34			#ifdef ALLOW_AUTODIFF_TAMC
35
36			C !LOCAL VARIABLES:
37			C == Local variables ==
38			C bi,bj - Loop counters
39			C I,J,K
40			INTEGER bi, bj
41	heimbach	1.8	INTEGER I, J, K
42	heimbach	1.1	CEOP
43
44			C-- Over all tiles
45			DO bj = myByLo(myThid), myByHi(myThid)
46			DO bi = myBxLo(myThid), myBxHi(myThid)
47
48			C- 3D arrays
49	heimbach	1.8	DO K=1,Nr
50			DO J=1-Oly,sNy+Oly
51			DO I=1-Olx,sNx+Olx
52			#ifdef ALLOW_EP_FLUX
53			EfluxY(i,j,k,bi,bj)=0.
54			EfluxP(i,j,k,bi,bj)=0.
55			#endif
56			ENDDO
57			ENDDO
58			ENDDO
59	heimbach	1.1
60			C- 2D arrays
61			DO J=1-Oly,sNy+Oly
62			DO I=1-Olx,sNx+Olx
63			#ifdef EXACT_CONSERV
64	heimbach	1.6	dEtaHdt(i,j,bi,bj) = 0. _d 0
65	heimbach	1.7	PmEpR(i,j,bi,bj) = 0. _d 0
66	heimbach	1.6	cph hDivFlow(I,J,bi,bj)=0. _d 0
67	heimbach	1.3	#endif
68	edhill	1.5	#ifdef ALLOW_EXF
69	heimbach	1.3	EmPmR(I,J,bi,bj) = 0. _d 0
70	heimbach	1.1	#endif
71			ENDDO
72			ENDDO
73
74			ENDDO
75			ENDDO
76
77			#endif /* ALLOW_AUTODIFF_TAMC */
78
79			END