Changeset 6715 for branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90

Timestamp:

2016-06-16T18:58:14+02:00 (8 years ago)

Author:

jpalmier

Message:

JPALM -- 16-06-2016 -- MEDUSA branch update :

-- pass co2 flux and dms_surf through restart for atm coupling.
-- introduce CFC cycle for dynamic evolution comparison
-- add Tim Graham Age tracer
-- include MEDUSA Q10 modif
-- svn-key removed
-- still need debug stage

File:

• branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90 (modified) (11 diffs)

branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcbio_medusa.F90

@@ -175 +175 @@
       REAL(wp) ::    loc_T,fun_T,xvpnT,xvpdT
       INTEGER  ::    ieppley
+      !! AXY (16/05/11): per Katya's prompting, add in new T-dependence
+      !!                 for phytoplankton growth only (i.e. no change
+      !!                 for remineralisation)
+      REAL(wp) ::    fun_Q10
       !! AXY (01/03/10): add in mixed layer PP diagnostics
       REAL(wp), DIMENSION(jpi,jpj) ::  fprn_ml,fprd_ml
@@ -316 +320 @@
       !! AXY (13/03/15): add in other DMS calculations
       REAL(wp) ::    dms_andr, dms_simo, dms_aran, dms_hall
+#  if defined key_oasis3
+      REAL(wp), DIMENSION(jpi,jpj) :: pco2a_2d(ji,jj)        !! use 2D atm pCO2 from atm coupling
+#  endif
+
 # endif
       !! 
@@ -353 +361 @@
       REAL(wp), POINTER, DIMENSION(:,:  ) :: mdetc2d, gmidc2d, gmedc2d, f_pco2a2d, f_pco2w2d, f_co2flux2d
       REAL(wp), POINTER, DIMENSION(:,:  ) :: f_TDIC2d, f_TALK2d, f_kw6602d, f_pp02d, f_o2flux2d, f_o2sat2d
-      REAL(wp), POINTER, DIMENSION(:,:  ) :: dms_surf2d, dms_andr2d, dms_simo2d, dms_aran2d, dms_hall2d
+      REAL(wp), POINTER, DIMENSION(:,:  ) :: dms_andr2d, dms_simo2d, dms_aran2d, dms_hall2d, dms_surf2d
       REAL(wp), POINTER, DIMENSION(:,:  ) :: iben_n2d, iben_fe2d, iben_c2d, iben_si2d, iben_ca2d, oben_n2d
       REAL(wp), POINTER, DIMENSION(:,:  ) :: oben_fe2d, oben_c2d, oben_si2d, oben_ca2d, sfr_ocal2d
       REAL(wp), POINTER, DIMENSION(:,:  ) :: sfr_oarg2d, lyso_ca2d 
+
 # endif
       !! 2D var for diagnostics.
@@ -814 +823 @@
       ENDIF
       IF( med_diag%CO2FLUX%dgsave ) THEN
-         CALL wrk_alloc( jpi, jpj,    f_co2flux2d   )
+         CALL wrk_alloc( jpi, jpj,   f_co2flux2d    )
          f_co2flux2d(:,:)      = 0.0      !!
       ENDIF
@@ -897 +906 @@
       IF (jdms .eq. 1) THEN
          IF( med_diag%DMS_SURF%dgsave ) THEN
-            CALL wrk_alloc( jpi, jpj,  dms_surf2d     )
-         dms_surf2d(:,:)      = 0.0      !!
+            CALL wrk_alloc( jpi, jpj,   dms_surf2d    )
+            dms_surf2d(:,:)      = 0.0      !!
          ENDIF
          IF( med_diag%DMS_ANDR%dgsave ) THEN
             CALL wrk_alloc( jpi, jpj,   dms_andr2d    )
-         dms_andr2d(:,:)      = 0.0      !!
+            dms_andr2d(:,:)      = 0.0      !!
          ENDIF
          IF( med_diag%DMS_SIMO%dgsave ) THEN
             CALL wrk_alloc( jpi, jpj,  dms_simo2d     )
-         dms_simo2d(:,:)      = 0.0      !!
+            dms_simo2d(:,:)      = 0.0      !!
          ENDIF
          IF( med_diag%DMS_ARAN%dgsave ) THEN
             CALL wrk_alloc( jpi, jpj,   dms_aran2d    )
-         dms_aran2d(:,:)      = 0.0      !!
+            dms_aran2d(:,:)      = 0.0      !!
          ENDIF
          IF( med_diag%DMS_HALL%dgsave ) THEN
             CALL wrk_alloc( jpi, jpj,   dms_hall2d    )
-         dms_hall2d(:,:)      = 0.0      !!
+            dms_hall2d(:,:)      = 0.0      !!
          ENDIF
       ENDIF   
@@ -1042 +1051 @@
                   !! OPEN wet point IF..THEN loop
                   if (tmask(ji,jj,jk).eq.1) then
+#  if defined key_oasis3
+                     f_pco2a = pco2a_2d(ji,jj)        !! use 2D atm pCO2 from atm coupling
+#  endif
                      !! do carbonate chemistry
                      !!
@@ -1371 +1383 @@
                   !!
                   f_wind  = wndm(ji,jj)
+#  if defined key_oasis3
+                  f_pco2a = pco2a_2d(ji,jj)        !! use 2D atm pCO2 from atm coupling
+#  endif
+
                   !!
                   !! AXY (23/06/15): as part of an effort to update the carbonate chemistry
@@ -1706 +1722 @@
                loc_T   = tsn(ji,jj,jk,jp_tem)
                fun_T   = 1.066**(1.0 * loc_T)
+               !! AXY (16/05/11): add in new Q10 (1.5, not 2.0) for
+               !phytoplankton
+               !!                 growth; remin. unaffected
+               fun_Q10 = jq10**((loc_T - 0.0) / 10.0)
                if (jphy.eq.1) then
                   xvpnT = xvpn * fun_T
                   xvpdT = xvpd * fun_T
+               elseif (jphy.eq.2) then
+                  xvpnT = xvpn * fun_Q10
+                  xvpdT = xvpd * fun_Q10
                else
                   xvpnT = xvpn
@@ -2126 +2149 @@
                   fddc = xmdc * fun_T * zdtc
 # endif
+               elseif (jmd.eq.2) then
+                  !! AXY (16/05/13): add in Q10-based parameterisation (def in nmlst)
+                  !! temperature-dependent
+                  fdd  = xmd  * fun_Q10 * zdet
+#if defined key_roam
+                  fddc = xmdc * fun_Q10 * zdtc
+#endif
                else
                   !! temperature-independent
@@ -4327 +4357 @@
                      IF( med_diag%CO2FLUX%dgsave ) THEN
                          CALL iom_put( "CO2FLUX"  , f_co2flux2d )
-                          CALL wrk_dealloc( jpi, jpj,   f_co2flux2d   )
+                         zb_co2_flx = zn_co2_flx
+                         zn_co2_flx = f_co2flux2d
+                         CALL wrk_dealloc( jpi, jpj,   f_co2flux2d   )
                      ENDIF
                      IF( med_diag%OM_CAL%dgsave ) THEN
@@ -4368 +4400 @@
                        IF( med_diag%DMS_SURF%dgsave ) THEN
                          CALL iom_put( "DMS_SURF"  , dms_surf2d )
-                          CALL wrk_dealloc( jpi, jpj,   dms_surf2d   )
+                         zb_dms_srf = zn_dms_srf
+                         zn_dms_srf = dms_surf2d
+                         CALL wrk_dealloc( jpi, jpj,   dms_surf2d   ) 
                        ENDIF
                        IF( med_diag%DMS_ANDR%dgsave ) THEN