Changeset 8436 for branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90

Timestamp:

2017-08-14T15:22:09+02:00 (7 years ago)

Author:

dford

Message:

Implement initial version of surface chlorophyll assimilation for MEDUSA.

File:

• branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90 (modified) (7 diffs)

branches/UKMO/dev_r5518_GO6_package_asm_surf_bgc/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90

@@ -80 +80 @@
 #  else
       USE trcco2_medusa
+#  if defined key_foam_medusa
+      USE mocsy_mainmod
+#  endif
 #  endif
       USE trcoxy_medusa
@@ -581 +584 @@
       fslownflux(:,:) = 0.0
       fslowcflux(:,:) = 0.0
+
+# if defined key_foam_medusa
+      pgrow_avg(:,:) = 0.0
+      ploss_avg(:,:) = 0.0
+      phyt_avg(:,:)  = 0.0
+      IF( kt == nittrc000 ) THEN
+         mld_max(:,:) = 0.0
+      ENDIF
+# endif
 
       !!
@@ -1347 +1359 @@
       !! We want this to be start of month or if starting afresh from 
       !! climatology - marc 20/6/17
+#if defined key_foam_medusa
+      !! DAF (Aug 2017): For FOAM we want to run daily
+      If ( (kt == nittrc000 .AND. .NOT.ln_rsttr) .OR.                        &
+           (mod(kt*rdt,86400.) == rdt) ) THEN
+#else
       If ( (kt == nittrc000 .AND. .NOT.ln_rsttr) .OR.                        &
            ((86400*mod(nn_date0,100) + mod(kt*rdt,2592000.)) == rdt) ) THEN
+#endif
          !!----------------------------------------------------------------------
          !! Calculate the carbonate chemistry for the whole ocean on the first
@@ -1795 +1813 @@
                      iters, ' AT (', ji, ', ', jj, ', ', jk, ') AT ', kt
                   endif
+#     if defined key_foam_medusa
+                  !! DAF (Aug 2017): calculate fCO2 for observation operator
+                  CALL p2fCO2( f_pco2w, ztmp, f_pp0, 0.0, 1, f_fco2w )
+#     endif
 #   endif
 #  else
@@ -1932 +1954 @@
                      fgco2(ji,jj) = f_co2flux * fthk * CO2flux_conv  !! mmol-C/m3/d -> kg-CO2/m2/s
                   !! ENDIF
+#   if defined key_foam_medusa
+                  !! DAF (Aug 2017): Save pCO2 and fCO2 for observation operator
+                  f2_pco2w(ji,jj) = f_pco2w
+                  f2_fco2w(ji,jj) = f_pco2w
+#   endif
                   IF ( lk_iomput ) THEN
                       IF( med_diag%ATM_PCO2%dgsave ) THEN
@@ -3603 +3630 @@
                CALL flush(numout)
 
+# if defined key_foam_medusa
+               !!----------------------------------------------------------------------
+               !! Mixed layer averages for ocean colour assimilation
+               !!----------------------------------------------------------------------
+               !!
+               if (fdep1.le.hmld(ji,jj)) then
+                  !! this level is entirely in the mixed layer
+                  fq0 = 1.0
+               elseif (fdep.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) - fdep) / fthk
+               endif
+               !!
+               pgrow_avg(ji,jj) = pgrow_avg(ji,jj) + ( &
+                  ( (fprn * zphn) + (fprd * zphd) ) * fthk * fq0)
+               ploss_avg(ji,jj) = ploss_avg(ji,jj) + ( &
+                  ( fgmipn + fgmepn + fdpn + fdpn2 + fgmepd + fdpd + fdpd2 ) * fthk * fq0 )
+               phyt_avg(ji,jj)  = phyt_avg(ji,jj)  + ( &
+                  (zphn +zphd) * fthk * fq0 )
+               !!
+# endif
                !!======================================================================
                !! LOCAL GRID CELL TRENDS
@@ -5331 +5382 @@
       endif
       
+# if defined key_foam_medusa
+      !!----------------------------------------------------------------------
+      !! Dianostics 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( ln_diatrc ) THEN
          !!----------------------------------------------------------------------