Changeset 10302 for branches/UKMO/dev_r5518_GO6_package_FOAMv14/NEMOGCM/NEMO/TOP_SRC/MEDUSA/plankton.F90

Timestamp:

2018-11-13T18:21:16+01:00 (5 years ago)

Author:

dford

Message:

Merge in revisions 8447:10159 of dev_r5518_GO6_package.

File:

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

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

@@ -36 +36 @@
                                    fdpn, fdpn2, fdzme, fdzme2,             &
                                    fdzmi, fdzmi2, fsdiss, fsin,            &
+                                   fdep1, fprn, fprd,                      &
+                                   fgmepd, fgmepn, fgmipn,                 &
                                    zphd, zphn, zpds, zzme, zzmi
-      USE dom_oce,           ONLY: tmask
+      USE dom_oce,           ONLY: e3t_0, e3t_n, gdepw_0, gdepw_n, tmask
+      USE par_kind,          ONLY: wp
       USE par_oce,           ONLY: jpim1, jpjm1
       USE phytoplankton_mod, ONLY: phytoplankton
       USE sms_medusa,        ONLY: jmpd, jmpn, jmzme, jmzmi,               &
+                                   ln_foam_medusa,                         &
+                                   pgrow_avg, ploss_avg, phyt_avg,         &
                                    xkphd, xkphn, xkzme, xkzmi,             &
                                    xmetapd, xmetapn, xmetazme, xmetazmi,   &
                                    xmpd, xmpn, xmzme, xmzmi, xsdiss
+      USE zdfmxl,            ONLY: hmld
       USE zooplankton_mod,   ONLY: zooplankton
+
+   !!* Substitution
+#  include "domzgr_substitute.h90"
 
       !! Level
@@ -50 +59 @@
 
       INTEGER :: ji, jj
+
+      REAL(wp) :: fq0
 
       !!-------------------------------------------------------------------
@@ -188 +199 @@
       ENDDO
 
+      IF ( ln_foam_medusa ) THEN
+         !! Mixed layer averages for ocean colour assimilation
+         !!
+         DO jj = 2,jpjm1
+            DO ji = 2,jpim1
+               IF (tmask(ji,jj,jk) == 1) THEN
+                  if (fdep1(ji,jj).le.hmld(ji,jj)) then
+                     !! this level is entirely in the mixed layer
+                     fq0 = 1.0
+                  elseif (fsdepw(ji,jj,jk).ge.hmld(ji,jj)) then
+                     !! this level is entirely below the mixed layer
+                     fq0 = 0.0
+                  else
+                     !! this level straddles the mixed layer
+                     fq0 = (hmld(ji,jj) - fsdepw(ji,jj,jk)) / fse3t(ji,jj,jk)
+                  endif
+                  !!
+                  pgrow_avg(ji,jj) = pgrow_avg(ji,jj) +                   &
+                                     (((fprn(ji,jj) * zphn(ji,jj)) +      &
+                                       (fprd(ji,jj) * zphd(ji,jj))  ) *   &
+                                      fse3t(ji,jj,jk) * fq0)
+                  ploss_avg(ji,jj) = ploss_avg(ji,jj) +                   &
+                                     ((fgmepd(ji,jj) + fdpd(ji,jj) +      &
+                                       fdpd2(ji,jj)                +      &
+                                       fgmepn(ji,jj) + fdpn(ji,jj) +      &
+                                       fdpn2(ji,jj)  + fgmipn(ji,jj) ) *  &
+                                      fse3t(ji,jj,jk) * fq0)
+                  phyt_avg(ji,jj)  = phyt_avg(ji,jj)  +                   &
+                                     ((zphn(ji,jj) + zphd(ji,jj)) *       &
+                                      fse3t(ji,jj,jk) * fq0)
+               ENDIF
+            ENDDO
+         ENDDO
+      ENDIF
+
    END SUBROUTINE plankton