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

Timestamp:

2020-07-16T19:08:23+02:00 (4 years ago)

Author:

dford

Message:

Allow nitrogen balancing when assimilating 3D chlorophyll data. See Met Office utils ticket 346.

File:

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

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

@@ -46 +46 @@
                                    ln_foam_medusa,                         &
                                    pgrow_avg, ploss_avg, phyt_avg,         &
+                                   pgrow_avg_3d, ploss_avg_3d, phyt_avg_3d, &
                                    xkphd, xkphn, xkzme, xkzmi,             &
                                    xmetapd, xmetapn, xmetazme, xmetazmi,   &
@@ -229 +230 @@
                                      ((zphn(ji,jj) + zphd(ji,jj)) *       &
                                       fse3t(ji,jj,jk) * fq0)
+                  !!
+                  pgrow_avg_3d(ji,jj,jk) = (fprn(ji,jj) * zphn(ji,jj)) +  &
+                                           (fprd(ji,jj) * zphd(ji,jj))
+                  ploss_avg_3d(ji,jj,jk) = fgmepd(ji,jj) + fdpd(ji,jj) +  &
+                                           fdpd2(ji,jj)                +  &
+                                           fgmepn(ji,jj) + fdpn(ji,jj) +  &
+                                           fdpn2(ji,jj)  + fgmipn(ji,jj)
+                  phyt_avg_3d(ji,jj,jk)  = zphn(ji,jj) + zphd(ji,jj)
                ENDIF
             ENDDO