pkg/autodiff/ini_autodiff.F

C $Header: /u/gcmpack/MITgcm/pkg/autodiff/ini_autodiff.F,v 1.5 2003/11/01 04:50:02 edhill 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"
#ifdef EXACT_CONSERV
# include "SURFACE.h"
#endif
#ifdef ALLOW_EXF
# include "FFIELDS.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

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
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.6	C $Header: /u/gcmpack/MITgcm/pkg/autodiff/ini_autodiff.F,v 1.5 2003/11/01 04:50:02 edhill 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.3	#ifdef EXACT_CONSERV
25			# include "SURFACE.h"
26			#endif
27	edhill	1.5	#ifdef ALLOW_EXF
28	heimbach	1.3	# include "FFIELDS.h"
29			#endif
30	heimbach	1.1
31			C !INPUT/OUTPUT PARAMETERS:
32			C == Routine arguments ==
33			C myThid - Number of this instance of INI_UVEL
34			INTEGER myThid
35
36			#ifdef ALLOW_AUTODIFF_TAMC
37
38			C !LOCAL VARIABLES:
39			C == Local variables ==
40			C bi,bj - Loop counters
41			C I,J,K
42			INTEGER bi, bj
43			INTEGER I, J, K
44			CEOP
45
46			C-- Over all tiles
47			DO bj = myByLo(myThid), myByHi(myThid)
48			DO bi = myBxLo(myThid), myBxHi(myThid)
49
50			C- 3D arrays
51
52			C- 2D arrays
53			DO J=1-Oly,sNy+Oly
54			DO I=1-Olx,sNx+Olx
55			#ifdef EXACT_CONSERV
56	heimbach	1.6	dEtaHdt(i,j,bi,bj) = 0. _d 0
57			cph hDivFlow(I,J,bi,bj)=0. _d 0
58	heimbach	1.3	#endif
59	edhill	1.5	#ifdef ALLOW_EXF
60	heimbach	1.3	EmPmR(I,J,bi,bj) = 0. _d 0
61	heimbach	1.1	#endif
62			ENDDO
63			ENDDO
64
65			ENDDO
66			ENDDO
67
68			#endif /* ALLOW_AUTODIFF_TAMC */
69
70			END