Add 2d sediment flux diagnostic

# Conflicts:
#	pkg/darwin/darwin_cons.F
#	pkg/darwin/darwin_diagnostics_init.F

# Conflicts:
#	pkg/darwin/darwin_cons.F
#	pkg/darwin/darwin_forcing.F

Author: jahn

pkg/darwin/darwin_diagnostics_init.F

@@ -548,6 +548,13 @@ SUBROUTINE DARWIN_DIAGNOSTICS_INIT( myThid )
       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
      I           diagName, diagCode, diagUnits, diagTitle, 0, myThid )
 
+      WRITE(diagName,'(A)')'sedFlxFe'
+      WRITE(diagTitle,'(A)')'Iron sediment flux'
+      diagUnits = 'mmol Fe/m2/s    '
+      diagCode  = 'SM P    U1      '
+      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
+     I           diagName, diagCode, diagUnits, diagTitle, 0, myThid )
+
 #ifdef DARWIN_ALLOW_CARBON
       WRITE(diagName,'(A)')'pH'
       WRITE(diagTitle,'(A)')'pH'

pkg/darwin/darwin_forcing.F

@@ -91,6 +91,7 @@ SUBROUTINE DARWIN_FORCING( Ptrdummy,
       _RL gO2(1-OLx:sNx+OLx, 1-OLy:sNy+OLy)
       _RL surfFe(1-OLx:sNx+OLx, 1-OLy:sNy+OLy)
       _RL freeFe(1-OLx:sNx+OLx, 1-OLy:sNy+OLy, Nr)
+      _RL sedFlxFe(1-OLx:sNx+OLx, 1-OLy:sNy+OLy)
       _RL tmp3d(1-OLx:sNx+OLx, 1-OLy:sNy+OLy, Nr)
       _RL scv,scav_poc
       _RL flx, POCl
@@ -448,10 +449,13 @@ SUBROUTINE DARWIN_FORCING( Ptrdummy,
 #else
           flx = fesedflux
 #endif
+          sedFlxFe(i,j) = flx
           sedFe(i,j,k) = flx/(drF(k)*hFacC(i,j,k,bi,bj))
 #ifdef DARWIN_ALLOW_CONS
           ironSedFlx(i,j,bi,bj) = ironSedFlx(i,j,bi,bj) + dTsub(k)*flx
 #endif
+        ELSE
+         sedFlxFe(i,j) = 0 _d 0
         ENDIF
        ENDDO
       ENDDO
@@ -887,7 +891,8 @@ SUBROUTINE DARWIN_FORCING( Ptrdummy,
        ENDDO
 #endif
 #endif
-       CALL DIAGNOSTICS_FILL(sedFe, 'sedFe   ',0,Nr,1,bi,bj,myThid)
+       CALL DIAGNOSTICS_FILL(sedFe,   'sedFe   ',0,Nr,1,bi,bj,myThid)
+       CALL DIAGNOSTICS_FILL(sedFlxFe,'sedFlxFe',0,1,2,bi,bj,myThid)
        IF (DIAGNOSTICS_IS_ON('sfcSolFe', myThid)) THEN
         DO j=jMin,jMax
          DO i=iMin,iMax