Changeset 8131 for branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcdms_medusa.F90

Timestamp:

2017-06-02T17:14:01+02:00 (7 years ago)

Author:

jpalmier

Message:

JPALM -- DMS add an updated version of ANDERSON sheme

File:

• branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcdms_medusa.F90 (modified) (5 diffs)

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

@@ -6 +6 @@
    !! History :
    !!  -   !  2014-08  (J. Palmieri - A. Yool)    added for UKESM1 project
+   !!  -   !  2017-05  (A. Yool)                  add extra Anderson scheme
    !!----------------------------------------------------------------------
 #if defined key_medusa && defined key_roam
@@ -38 +39 @@
 !=======================================================================
 !
-   SUBROUTINE trc_dms_medusa( chn, chd, mld, xqsr, xdin,   &  !! inputs
-     &  dms_andr, dms_simo, dms_aran, dms_hall )              !! outputs
+   SUBROUTINE trc_dms_medusa( chn, chd, mld, xqsr, xdin, xlim,  &  !! inputs
+     &  dms_andr, dms_simo, dms_aran, dms_hall dms_andm)           !! outputs
 !      
 !=======================================================================
@@ -69 +70 @@
       !!                    and is included because it is formally
       !!                    published (and different from the above)
+      !!
+      !! AXY (25/05/17): add extra "corrected" Anderson scheme
+      !!
+      !!                 5. As Anderson et al. (2001) but modified to
+      !!                    more accurately reflect nutrient limitation
+      !!                    status of phytoplankton community
       !!
       !! AXY (08/07/15): amend to remove Julien's original calculation
@@ -84 +91 @@
       REAL(wp), INTENT( in )    :: xqsr                 !! surface irradiance        (W/m2)
       REAL(wp), INTENT( in )    :: xdin                 !! surface DIN               (mmol N/m3)
+      REAL(wp), INTENT( in )    :: xlim                 !! surface DIN limitation    (mmol N/m3)
       REAL(wp), INTENT( inout ) :: dms_andr             !! DMS surface concentration (nmol/L) 
       REAL(wp), INTENT( inout ) :: dms_simo             !! DMS surface concentration (nmol/L) 
       REAL(wp), INTENT( inout ) :: dms_aran             !! DMS surface concentration (nmol/L) 
       REAL(wp), INTENT( inout ) :: dms_hall             !! DMS surface concentration (nmol/L) 
+      REAL(wp), INTENT( inout ) :: dms_andm             !! DMS surface concentration (nmol/L) 
 !
       REAL(wp) :: CHL, cmr, sw_dms
@@ -147 +156 @@
            dms_hall = fq3
         endif
+!
+! AXY (25/05/17): modified Anderson et al. (2001)
+        Jterm = xqsr + 1.0e-6
+        !! this version fixes the hard-coded assumption above
+        Qterm = xlim
+        fq1 = log10(CHL * Jterm * Qterm)
+        if (fq1 > 1.72) then
+           dms_andm = (8.24 * (fq1 - 1.72)) + 2.29
+        else
+           dms_andm = 2.29
+        endif
 
   END SUBROUTINE trc_dms_medusa