pkg/darwin/dic_diags_init.F

#include "CPP_OPTIONS.h"
#include "PTRACERS_OPTIONS.h"
#include "DARWIN_OPTIONS.h"

#ifdef ALLOW_PTRACERS
#ifdef ALLOW_DARWIN

#ifdef ALLOW_CARBON

CBOP
C !ROUTINE: DIC_DIAGS_INIT

C !INTERFACE: ==========================================================
      SUBROUTINE DIC_DIAGS_INIT( myThid )

C !DESCRIPTION:
C  DIC diagnostic initialization

C !USES: ===============================================================
      IMPLICIT NONE
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "DARWIN_FLUX.h"

c!INPUT PARAMETERS: ===================================================
C  myThid               :: thread number
      INTEGER myThid


C !LOCAL VARIABLES: ====================================================
      INTEGER bi, bj, k
      CHARACTER*(MAX_LEN_MBUF) suff
CEOP

#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(SURave,   1,  bi, bj, myThid)
          CALL TIMEAVE_RESET(SURCave,   1,  bi, bj, myThid)
          CALL TIMEAVE_RESET(SUROave,   1,  bi, bj, myThid)
#ifdef pH_3D
          CALL TIMEAVE_RESET(pCO2ave,   Nr,  bi, bj, myThid)
          CALL TIMEAVE_RESET(pHave,   Nr,  bi, bj, myThid)
#else
          CALL TIMEAVE_RESET(pCO2ave,   1,  bi, bj, myThid)
          CALL TIMEAVE_RESET(pHave,   1,  bi, bj, myThid)
#endif
          print*,'QQ start timeave'
          do k=1,Nr
            DIC_TimeAve(bi,bj,k)=0.d0
          enddo
         ENDDO
        ENDDO
#endif /* ALLOW_TIMEAVE */

        RETURN
        END
#endif  /*ALLOW_CARBON*/

#endif  /*DARWIN*/
#endif  /*ALLOW_PTRACERS*/
c ==================================================================

1	jahn	1.1	#include "CPP_OPTIONS.h"
2			#include "PTRACERS_OPTIONS.h"
3			#include "DARWIN_OPTIONS.h"
4
5			#ifdef ALLOW_PTRACERS
6			#ifdef ALLOW_DARWIN
7
8			#ifdef ALLOW_CARBON
9
10			CBOP
11			C !ROUTINE: DIC_DIAGS_INIT
12
13			C !INTERFACE: ==========================================================
14			SUBROUTINE DIC_DIAGS_INIT( myThid )
15
16			C !DESCRIPTION:
17			C DIC diagnostic initialization
18
19			C !USES: ===============================================================
20			IMPLICIT NONE
21			#include "SIZE.h"
22			#include "EEPARAMS.h"
23			#include "PARAMS.h"
24			#include "DARWIN_FLUX.h"
25
26			c!INPUT PARAMETERS: ===================================================
27			C myThid :: thread number
28			INTEGER myThid
29
30
31			C !LOCAL VARIABLES: ====================================================
32			INTEGER bi, bj, k
33			CHARACTER*(MAX_LEN_MBUF) suff
34			CEOP
35
36			#ifdef ALLOW_TIMEAVE
37			c set arrays to zero if first timestep
38			DO bj = myByLo(myThid), myByHi(myThid)
39			DO bi = myBxLo(myThid), myBxHi(myThid)
40			CALL TIMEAVE_RESET(SURave, 1, bi, bj, myThid)
41			CALL TIMEAVE_RESET(SURCave, 1, bi, bj, myThid)
42			CALL TIMEAVE_RESET(SUROave, 1, bi, bj, myThid)
43	stephd	1.2	#ifdef pH_3D
44			CALL TIMEAVE_RESET(pCO2ave, Nr, bi, bj, myThid)
45			CALL TIMEAVE_RESET(pHave, Nr, bi, bj, myThid)
46			#else
47	jahn	1.1	CALL TIMEAVE_RESET(pCO2ave, 1, bi, bj, myThid)
48			CALL TIMEAVE_RESET(pHave, 1, bi, bj, myThid)
49	stephd	1.2	#endif
50	jahn	1.1	print*,'QQ start timeave'
51			do k=1,Nr
52			DIC_TimeAve(bi,bj,k)=0.d0
53			enddo
54			ENDDO
55			ENDDO
56			#endif /* ALLOW_TIMEAVE */
57
58			RETURN
59			END
60			#endif /ALLOW_CARBON/
61
62			#endif /DARWIN/
63			#endif /ALLOW_PTRACERS/
64			c ==================================================================
65