Changeset 9187 for branches/NERC/dev_r5518_GO6_MEDUSA_DMS_TUNING/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcdms_medusa.F90

Timestamp:

2018-01-05T17:33:15+01:00 (6 years ago)

Author:

jpalmier

Message:

JPALM -- GMED #365 -- enable Anderson DMS tuning

File:

• branches/NERC/dev_r5518_GO6_MEDUSA_DMS_TUNING/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcdms_medusa.F90 (modified) (2 diffs)

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

@@ -113 +113 @@
 !
 ! AXY (13/03/15): Anderson et al. (2001)
+!! JPALM --19-12-2017-- Tunable through the namelist
+!!                      within dmsmin - dmscut - dmsslp
         Jterm = xqsr + 1.0e-6
         !! this next line makes a hard-coded assumption about the 
@@ -120 +122 @@
         Qterm = xdin / (xdin + 0.5)
         fq1 = log10(CHL * Jterm * Qterm)
-        if (fq1 > 1.72) then
-           dms_andr = (8.24 * (fq1 - 1.72)) + 2.29
-        else
-           dms_andr = 2.29
+        if (fq1 > dmscut) then
+           dms_andr = (dmsslp * (fq1 - dmscut)) + dmsmin
+        else
+           dms_andr = dmsmin
         endif
 !