Changeset 8495 for branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc_v2/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_fin.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/bio_medusa_fin.F90 (modified) (3 diffs)

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

@@ -45 +45 @@
                                    f3_co3, f3_h2co3, f3_hco3,           &
                                    f3_omarg, f3_omcal, f3_pH,           &
+# if defined key_foam_medusa
+                                   mld_max, pgrow_avg,                  &
+                                   ploss_avg, phyt_avg,                 &
+# endif
                                    za_sed_c, za_sed_ca, za_sed_fe,      &
                                    za_sed_n, za_sed_si,                 &
@@ -53 +57 @@
       USE trc,               ONLY: med_diag, nittrc000 
       USE trcnam_trp,        ONLY: ln_trcadv_cen2, ln_trcadv_tvd
+# if defined key_foam_medusa
+      USE zdfmxl,            ONLY: hmld
+# endif
  
       !! time (integer timestep)
@@ -218 +225 @@
       endif
 
+# if defined key_foam_medusa
+!!----------------------------------------------------------------------
+!! Diagnostics required for ocean colour assimilation:
+!! Mixed layer average phytoplankton growth, loss and concentration
+!! Maximum mixed layer depth
+!!----------------------------------------------------------------------
+!!
+      DO jj = 2,jpjm1
+         DO ji = 2,jpim1
+            pgrow_avg(ji,jj) = pgrow_avg(ji,jj) / hmld(ji,jj)
+            ploss_avg(ji,jj) = ploss_avg(ji,jj) / hmld(ji,jj)
+            phyt_avg(ji,jj)  = phyt_avg(ji,jj)  / hmld(ji,jj)
+            IF ( hmld(ji,jj) .GT. mld_max(ji,jj) ) THEN
+               mld_max(ji,jj) = hmld(ji,jj)
+            ENDIF
+         END DO
+      END DO
+# endif
+
 #  if defined key_debug_medusa
          !! AXY (12/07/17)