pkg/dic/dic_biotic_init.F

C $Header: /u/gcmpack/MITgcm/pkg/dic/dic_biotic_diags.F,v 1.2 2003/07/09 19:59:18 stephd Exp $
C $Name:  $

#include "CPP_OPTIONS.h"
#include "GCHEM_OPTIONS.h"

CStartOfInterface
      SUBROUTINE DIC_BIOTIC_INIT( myThid )
C     /==========================================================\
C     | SUBROUTINE DIC_BIOTIC_INIT                              |
C     | o Do DIC diagnostic initialization                   |
C     |==========================================================|
C     \==========================================================/
      IMPLICIT NONE

C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#ifdef DIC_BIOTIC
#include "DIC_BIOTIC.h"
#endif

C     == Routine arguments ==
C     myThid -  Number of this instance of INI_FORCING
      INTEGER myThid
CEndOfInterface


C     == Local variables ==
      LOGICAL  DIFFERENT_MULTIPLE
      EXTERNAL DIFFERENT_MULTIPLE
      INTEGER bi, bj, k
      CHARACTER*(MAX_LEN_MBUF) suff

#ifdef DIC_BIOTIC


#ifdef ALLOW_TIMEAVE

c set arrays to zero if first timestep
        DO bj = myByLo(myThid), myByHi(myThid)
         DO bi = myBxLo(myThid), myBxHi(myThid)
          CALL TIMEAVE_RESET(BIOave,   Nr,  bi, bj, myThid)
          CALL TIMEAVE_RESET(CARave,   Nr,  bi, bj, myThid)
          CALL TIMEAVE_RESET(SURave,   1,  bi, bj, myThid)
          CALL TIMEAVE_RESET(SUROave,   1,  bi, bj, myThid)
          CALL TIMEAVE_RESET(pCO2ave,   1,  bi, bj, myThid)
          CALL TIMEAVE_RESET(pHave,   1,  bi, bj, myThid)
          CALL TIMEAVE_RESET(fluxCO2ave,   1,  bi, bj, myThid)
          print*,'QQ start timeave'
          do k=1,Nr
            DIC_TimeAve(bi,bj,k)=0.d0
          enddo
         ENDDO
        ENDDO

#endif /* ALLOW_TIMEAVE */
      
#endif /* DIC_BIOTIC   */
      
      RETURN
      END
1	stephd	1.1	C $Header: /u/gcmpack/MITgcm/pkg/dic/dic_biotic_diags.F,v 1.2 2003/07/09 19:59:18 stephd Exp $
2			C $Name: $
3
4			#include "CPP_OPTIONS.h"
5			#include "GCHEM_OPTIONS.h"
6
7			CStartOfInterface
8			SUBROUTINE DIC_BIOTIC_INIT( myThid )
9			C /==========================================================\
10			C \| SUBROUTINE DIC_BIOTIC_INIT \|
11			C \| o Do DIC diagnostic initialization \|
12			C \|==========================================================\|
13			C \==========================================================/
14			IMPLICIT NONE
15
16			C === Global variables ===
17			#include "SIZE.h"
18			#include "EEPARAMS.h"
19			#include "PARAMS.h"
20			#ifdef DIC_BIOTIC
21			#include "DIC_BIOTIC.h"
22			#endif
23
24			C == Routine arguments ==
25			C myThid - Number of this instance of INI_FORCING
26			INTEGER myThid
27			CEndOfInterface
28
29
30			C == Local variables ==
31			LOGICAL DIFFERENT_MULTIPLE
32			EXTERNAL DIFFERENT_MULTIPLE
33			INTEGER bi, bj, k
34			CHARACTER*(MAX_LEN_MBUF) suff
35
36			#ifdef DIC_BIOTIC
37
38
39			#ifdef ALLOW_TIMEAVE
40
41			c set arrays to zero if first timestep
42			DO bj = myByLo(myThid), myByHi(myThid)
43			DO bi = myBxLo(myThid), myBxHi(myThid)
44			CALL TIMEAVE_RESET(BIOave, Nr, bi, bj, myThid)
45			CALL TIMEAVE_RESET(CARave, Nr, bi, bj, myThid)
46			CALL TIMEAVE_RESET(SURave, 1, bi, bj, myThid)
47			CALL TIMEAVE_RESET(SUROave, 1, bi, bj, myThid)
48			CALL TIMEAVE_RESET(pCO2ave, 1, bi, bj, myThid)
49			CALL TIMEAVE_RESET(pHave, 1, bi, bj, myThid)
50			CALL TIMEAVE_RESET(fluxCO2ave, 1, bi, bj, myThid)
51			print*,'QQ start timeave'
52			do k=1,Nr
53			DIC_TimeAve(bi,bj,k)=0.d0
54			enddo
55			ENDDO
56			ENDDO
57
58			#endif /* ALLOW_TIMEAVE */
59
60			#endif /* DIC_BIOTIC */
61
62			RETURN
63			END