Changeset 8012 for branches/UKMO/dev_r5518_medusa_chg_trc_bio_medusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/plankton.F90

Timestamp:

2017-05-10T09:47:35+02:00 (7 years ago)

Author:

marc

Message:

Pulled further code from trcbio_medusa.F90 into other files

File:

• branches/UKMO/dev_r5518_medusa_chg_trc_bio_medusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/plankton.F90 (modified) (3 diffs)

branches/UKMO/dev_r5518_medusa_chg_trc_bio_medusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/plankton.F90

@@ -35 +35 @@
       USE bio_medusa_mod,    ONLY: fdpd, fdpd2, fdpds, fdpds2,             &
                                    fdpn, fdpn2, fdzme, fdzme2,             &
-                                   fdzmi, fdzmi2, fsin,                    &
+                                   fdzmi, fdzmi2, fsdiss, fsin,            &
                                    zphd, zphn, zpds, zzme, zzmi
       USE dom_oce,           ONLY: tmask
@@ -43 +43 @@
                                    xkphd, xkphn, xkzme, xkzmi,             &
                                    xmetapd, xmetapn, xmetazme, xmetazmi,   &
-                                   xmpd, xmpn, xmzme, xmzmi
+                                   xmpd, xmpn, xmzme, xmzmi, xsdiss
       USE zooplankton_mod,   ONLY: zooplankton
 
@@ -178 +178 @@
       ENDDO
 
+      !! diatom frustule dissolution. This section is moved from just
+      !! below CALL to iron_chem_scav in trcbio_medusa.F90 - marc 9/5/17
+      DO jj = 2,jpjm1
+         DO ji = 2,jpim1
+            IF (tmask(ji,jj,jk) == 1) THEN
+               fsdiss(ji,jj)  = xsdiss * zpds(ji,jj)
+            ENDIF
+         ENDDO
+      ENDDO
+
    END SUBROUTINE plankton