pkg/darwin/darwin_cons_reset.F

5e411acc9e Oliv*0001 #include "DARWIN_OPTIONS.h"
                0002
                0003 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0004 CBOP
                0005 C     !ROUTINE: DARWIN_CONS_RESET
                0006
                0007 C     !INTERFACE:
                0008       SUBROUTINE DARWIN_CONS_RESET( myThid )
                0009
                0010 C     !DESCRIPTION:
                0011 C     Reset conservation correction accumulator arrays
                0012
                0013 C     !USES:
                0014       IMPLICIT NONE
                0015 #include "EEPARAMS.h"
                0016 #include "SIZE.h"
                0017 #include "DARWIN_SIZE.h"
                0018 #include "DARWIN_FIELDS.h"
                0019
                0020 C     !INPUT PARAMETERS:
                0021       INTEGER myThid
                0022 CEOP
                0023
                0024 #ifdef ALLOW_DARWIN
                0025 #ifdef DARWIN_ALLOW_CONS
                0026
                0027 C     !LOCAL VARIABLES:
                0028       INTEGER i,j,k,bi,bj
                0029
                0030       DO bj=myByLo(myThid),myByHi(myThid)
                0031        DO bi=myBxLo(myThid),myBxHi(myThid)
                0032         DO k=1,Nr
                0033          DO j=1,sNy
                0034           DO i=1,sNx
                0035            DARWIN_Nfix(i,j,k,bi,bj) = 0 _d 0
                0036            DARWIN_Ndenit(i,j,k,bi,bj) = 0 _d 0
                0037            DARWIN_partScav(i,j,k,bi,bj) = 0 _d 0
                0038            DARWIN_minFeLoss(i,j,k,bi,bj) = 0 _d 0
                0039            DARWIN_O2prod(i,j,k,bi,bj) = 0 _d 0
                0040            DARWIN_O2cons(i,j,k,bi,bj) = 0 _d 0
                0041            DARWIN_AlkSrc(i,j,k,bi,bj) = 0 _d 0
                0042           ENDDO
                0043          ENDDO
                0044         ENDDO
                0045         DO j=1,sNy
                0046          DO i=1,sNx
                0047           ironSedFlx(i,j,bi,bj) = 0 _d 0
2c93eb88ef Oliv*0048 #ifdef DARWIN_ALLOW_HYDROTHERMAL_VENTS
                0049           ironVentFlx(i,j,bi,bj) = 0 _d 0
                0050 #endif
5e411acc9e Oliv*0051           carbSfcFlx(i,j,bi,bj) = 0 _d 0
                0052           carbVirFlx(i,j,bi,bj) = 0 _d 0
                0053           oxySfcFlx(i,j,bi,bj) = 0 _d 0
                0054           alkVirFlx(i,j,bi,bj) = 0 _d 0
248275f1c4 Oliv*0055 #ifdef DARWIN_BOTTOM_SINK
                0056           botSnkC(i,j,bi,bj) = 0 _d 0
                0057           botSnkN(i,j,bi,bj) = 0 _d 0
                0058           botSnkP(i,j,bi,bj) = 0 _d 0
                0059           botSnkFe(i,j,bi,bj) = 0 _d 0
                0060           botSnkSi(i,j,bi,bj) = 0 _d 0
2a97638e3d Oliv*0061 #endif
8a10163480 Oliv*0062 #ifdef DARWIN_ALLOW_RADIv1
2a97638e3d Oliv*0063           radiFluxC(i,j,bi,bj) = 0 _d 0
                0064           radiFluxA(i,j,bi,bj) = 0 _d 0
                0065           radiFluxO(i,j,bi,bj) = 0 _d 0
8a10163480 Oliv*0066 #endif
                0067 #ifdef DARWIN_ALLOW_RADIv2
                0068           radiFluxC(i,j,bi,bj) = 0 _d 0
                0069           radiFluxA(i,j,bi,bj) = 0 _d 0
                0070           radiFluxO(i,j,bi,bj) = 0 _d 0
                0071           radiFluxP(i,j,bi,bj) = 0 _d 0
                0072           radiFluxN(i,j,bi,bj) = 0 _d 0
248275f1c4 Oliv*0073 #endif
5e411acc9e Oliv*0074          ENDDO
                0075         ENDDO
                0076        ENDDO
                0077       ENDDO
                0078
                0079 #endif /* DARWIN_ALLOW_CONS */
                0080 #endif /* ALLOW_DARWIN */
                0081
                0082       RETURN
                0083       END
                0084