Changeset 9292 for branches/UKMO/dev_r5518_GO6_package_asm_3d_bgc/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90

Timestamp:

2018-01-30T19:41:58+01:00 (6 years ago)

Author:

dford

Message:

Merge in changes from dev_r5518_GO6_package_asm_surf_bgc_v2.

File:

• branches/UKMO/dev_r5518_GO6_package_asm_3d_bgc/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90 (modified) (3 diffs)

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

@@ -102 +102 @@
       !! JPALM (27-06-2016): add lk_oasis for CO2 and DMS coupling with atm
       USE sbc_oce,                    ONLY: lk_oasis
+# if defined key_foam_medusa
+      USE sms_medusa,                 ONLY: hist_pco2, xobs_xco2a,          &
+                                            pgrow_avg, ploss_avg,           &
+                                            phyt_avg, mld_max
+# else
       USE sms_medusa,                 ONLY: hist_pco2
+# endif
       USE trc,                        ONLY: ln_rsttr, nittrc000, trn
       USE bio_medusa_init_mod,        ONLY: bio_medusa_init
@@ -319 +325 @@
          f_xco2a(:,:) = fq4
       endif
+#  if defined key_foam_medusa
+      IF ( xobs_xco2a > 0.0 ) THEN
+         IF(lwp) WRITE(numout,*) ' using observed atm pCO2 = ', xobs_xco2a
+         f_xco2a(:,:) = xobs_xco2a
+      ELSE
+         IF(lwp) WRITE(numout,*) ' xobs_xco2a <= 0 so using default atm pCO2'
+      ENDIF
+#  endif
 #  if defined key_axy_pi_co2
       !! OCMIP pre-industrial pCO2
@@ -358 +372 @@
       !!          x * 30d + 1*rdt(i.e: mod = rdt)   
       !!          ++ need to pass carb-chem output var through restarts
+#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
          !! simulation timestep and every month subsequently; the resulting 3D