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	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