pkg/gmredi/gmredi_init.F

C $Header: /u/gcmpack/models/MITgcmUV/pkg/gmredi/gmredi_init.F,v 1.4 2001/02/04 14:38:49 cnh Exp $
C $Name:  $

#include "GMREDI_OPTIONS.h"

      SUBROUTINE GMREDI_INIT( myThid )
C     /==========================================================\
C     | SUBROUTINE GMREDI_INIT                                   |
C     | o Routine to initialize GM/Redi variables and constants. |
C     |==========================================================|
C     | Initialize GM/Redi parameters, read in data.gmredi       |
C     \==========================================================/
      IMPLICIT NONE

C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "GMREDI.h"
#include "GMREDI_DIAGS.h"

C     === Routine arguments ===
C     myThid -  Number of this instance of KPP_INIT
      INTEGER myThid

#ifdef ALLOW_GMREDI

C     === Local variables ===
C     msgBuf      - Informational/error meesage buffer
C     iUnit       - Work variable for IO unit number
      CHARACTER*(MAX_LEN_MBUF) msgBuf
      INTEGER iUnit
      INTEGER k,bi,bj

#ifdef ALLOW_TIMEAVE
C     Initialize averages to zero
      DO bj = myByLo(myThid), myByHi(myThid)
       DO bi = myBxLo(myThid), myBxHi(myThid)
        CALL TIMEAVE_RESET(GM_Kwx_T,Nr, bi,bj,myThid)
        CALL TIMEAVE_RESET(GM_Kwy_T,Nr, bi,bj,myThid)
        CALL TIMEAVE_RESET(GM_Kwz_T,Nr, bi,bj,myThid)
        DO k=1,Nr
         GM_TimeAve(k,bi,bj)=0.
        ENDDO
#ifdef GM_VISBECK_VARIABLE_K
        CALL TIMEAVE_RESET(Visbeck_K_T, 1, bi,bj,myThid)
#endif
       ENDDO
      ENDDO
#endif /* ALLOW_TIMEAVE */

#endif /* ALLOW_GMREDI */

      return
      end
1	jmc	1.5	C $Header: /u/gcmpack/models/MITgcmUV/pkg/gmredi/gmredi_init.F,v 1.4 2001/02/04 14:38:49 cnh Exp $
2			C $Name: $
3	adcroft	1.1
4			#include "GMREDI_OPTIONS.h"
5
6			SUBROUTINE GMREDI_INIT( myThid )
7			C /==========================================================\
8			C \| SUBROUTINE GMREDI_INIT \|
9			C \| o Routine to initialize GM/Redi variables and constants. \|
10			C \|==========================================================\|
11			C \| Initialize GM/Redi parameters, read in data.gmredi \|
12			C \==========================================================/
13			IMPLICIT NONE
14
15			C === Global variables ===
16			#include "SIZE.h"
17			#include "EEPARAMS.h"
18			#include "PARAMS.h"
19			#include "GRID.h"
20			#include "GMREDI.h"
21			#include "GMREDI_DIAGS.h"
22
23			C === Routine arguments ===
24			C myThid - Number of this instance of KPP_INIT
25			INTEGER myThid
26
27			#ifdef ALLOW_GMREDI
28
29			C === Local variables ===
30			C msgBuf - Informational/error meesage buffer
31			C iUnit - Work variable for IO unit number
32			CHARACTER*(MAX_LEN_MBUF) msgBuf
33			INTEGER iUnit
34			INTEGER k,bi,bj
35
36	jmc	1.5	#ifdef ALLOW_TIMEAVE
37	adcroft	1.1	C Initialize averages to zero
38			DO bj = myByLo(myThid), myByHi(myThid)
39			DO bi = myBxLo(myThid), myBxHi(myThid)
40	jmc	1.5	CALL TIMEAVE_RESET(GM_Kwx_T,Nr, bi,bj,myThid)
41			CALL TIMEAVE_RESET(GM_Kwy_T,Nr, bi,bj,myThid)
42			CALL TIMEAVE_RESET(GM_Kwz_T,Nr, bi,bj,myThid)
43	adcroft	1.1	DO k=1,Nr
44			GM_TimeAve(k,bi,bj)=0.
45			ENDDO
46			#ifdef GM_VISBECK_VARIABLE_K
47	jmc	1.5	CALL TIMEAVE_RESET(Visbeck_K_T, 1, bi,bj,myThid)
48	adcroft	1.1	#endif
49			ENDDO
50			ENDDO
51	jmc	1.5	#endif /* ALLOW_TIMEAVE */
52	adcroft	1.1
53			#endif /* ALLOW_GMREDI */
54
55			return
56			end