Add scavenging rate diagnostic

# Conflicts: # pkg/darwin/darwin_forcing.F # Conflicts: # pkg/darwin/darwin_forcing.F # Conflicts: # pkg/darwin/darwin_forcing.F # Conflicts: # pkg/darwin/darwin_forcing.F
darwinproject · Jan 3, 2023 · 89a8b1e · 89a8b1e
1 parent 0446592
commit 89a8b1e
Show file tree

Hide file tree

Showing 2 changed files with 10 additions and 0 deletions.
diff --git a/pkg/darwin/darwin_diagnostics_init.F b/pkg/darwin/darwin_diagnostics_init.F
@@ -541,6 +541,13 @@ SUBROUTINE DARWIN_DIAGNOSTICS_INIT( myThid )
       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
      I           diagName, diagCode, diagUnits, diagTitle, 0, myThid )
 
+      WRITE(diagName,'(A)')'scavRate'
+      WRITE(diagTitle,'(A)')'Iron scavenging rate'
+      diagUnits = '1/s             '
+      diagCode  = 'SMRP    MR      '
+      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
+     I           diagName, diagCode, diagUnits, diagTitle, 0, myThid )
+
       WRITE(diagName,'(A)')'sedFe'
       WRITE(diagTitle,'(A)')'Iron input from sediment'
       diagUnits = 'mmol Fe/m3/s    '

diff --git a/pkg/darwin/darwin_forcing.F b/pkg/darwin/darwin_forcing.F
@@ -93,6 +93,7 @@ SUBROUTINE DARWIN_FORCING( Ptrdummy,
       _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 scavRate(1-OLx:sNx+OLx, 1-OLy:sNy+OLy, Nr)
       _RL scavLoss(1-OLx:sNx+OLx, 1-OLy:sNy+OLy, Nr)
       _RL scv,scav_pom
       _RL flx, POCl
@@ -391,6 +392,7 @@ SUBROUTINE DARWIN_FORCING( Ptrdummy,
 #else
         scv = scav
 #endif
+        scavRate(i,j,k) = scv
         scavLoss(i,j,k) = scv*MAX(0 _d 0, freefe(i,j,k))
         gPtr(i,j,k,iFeT) = gPtr(i,j,k,iFeT) - scavLoss(i,j,k)
 #ifdef DARWIN_ALLOW_CONS
@@ -895,6 +897,7 @@ SUBROUTINE DARWIN_FORCING( Ptrdummy,
        ENDDO
 #endif
 #endif
+       CALL DIAGNOSTICS_FILL(scavRate,'scavRate',0,Nr,2,bi,bj,myThid)
        CALL DIAGNOSTICS_FILL(scavLoss,'scvLosFe',0,Nr,2,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)