Changeset 8495 for branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc_v2/NEMOGCM/NEMO/TOP_SRC/MEDUSA/plankton.F90

Timestamp:

2017-09-04T18:50:10+02:00 (7 years ago)

Author:

dford

Message:

Merge in changes from dev_r5518_GO6_package_asm_surf_bgc, and adapt to the updated MEDUSA structure.

File:

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

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

@@ -36 +36 @@
                                    fdpn, fdpn2, fdzme, fdzme2,             &
                                    fdzmi, fdzmi2, fsdiss, fsin,            &
+# if defined key_foam_medusa
+                                   fdep1, fprn, fprd,                      &
+                                   fgmepd, fgmepn, fgmipn,                 &
+# endif
                                    zphd, zphn, zpds, zzme, zzmi
+# if defined key_foam_medusa
+      USE dom_oce,           ONLY: e3t_0, e3t_n, gdepw_0, gdepw_n, tmask
+      USE par_kind,          ONLY: wp
+# else
       USE dom_oce,           ONLY: tmask
+# endif
       USE par_oce,           ONLY: jpim1, jpjm1
       USE phytoplankton_mod, ONLY: phytoplankton
       USE sms_medusa,        ONLY: jmpd, jmpn, jmzme, jmzmi,               &
+# if defined key_foam_medusa
+                                   pgrow_avg, ploss_avg, phyt_avg,         &
+# endif
                                    xkphd, xkphn, xkzme, xkzmi,             &
                                    xmetapd, xmetapn, xmetazme, xmetazmi,   &
                                    xmpd, xmpn, xmzme, xmzmi, xsdiss
+# if defined key_foam_medusa
+      USE zdfmxl,            ONLY: hmld
+# endif
       USE zooplankton_mod,   ONLY: zooplankton
+
+# if defined key_foam_medusa
+   !!* Substitution
+#  include "domzgr_substitute.h90"
+# endif
 
       !! Level
@@ -50 +70 @@
 
       INTEGER :: ji, jj
+
+# if defined key_foam_medusa
+      REAL(wp) :: fq0
+# endif
 
       !!-------------------------------------------------------------------
@@ -188 +212 @@
       ENDDO
 
+# if defined key_foam_medusa
+      !! Mixed layer averages for ocean colour assimilation
+      !!
+      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
+
    END SUBROUTINE plankton