Changeset 9258 for branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcnam_medusa.F90

Timestamp:

2018-01-18T12:24:31+01:00 (6 years ago)

Author:

marc

Message:

Adding Julien's GMED ticket 365 for DMS constants for Anderson scheme (incl fixes to a couple of other tickets)

File:

• branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcnam_medusa.F90 (modified) (5 diffs)

branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcnam_medusa.F90

@@ -87 +87 @@
       &  xsdiss,                                              &
       &  sedlam,sedlostpoc,jpkb,jdms,jdms_input,jdms_model,   &
-      &  scl_chl, chl_out
+      &  scl_chl, chl_out, dmsmin, dmscut, dmsslp
 #if defined key_roam
       NAMELIST/natroam/ xthetaphy,xthetazoo,xthetanit,        &
@@ -127 +127 @@
       ! 1.4 namelist natbio : biological parameters
       ! -------------------------------------------
+      !! Note: the default values below will all be overwritten by the
+      !!       input in the namelist natbio.
       
       xxi         = 0.
@@ -250 +252 @@
       scl_chl     = 1.
       chl_out     = 1
+      dmsmin      = 2.29 !! Anderson DMS default
+      dmscut      = 1.72 !! Anderson DMS default
+      dmsslp      = 8.24 !! Anderson DMS default
             
       !REWIND(numnatm)
@@ -403 +408 @@
 !!       jdms_model  :  choice of DMS model passed to atmosphere
 !!                      1 = ANDR, 2 = SIMO, 3 = ARAN, 4 = HALL, 5 = ANDM
-!!
+!!       dmsmin      : DMS minimum value for DMS Anderson (ANDR) sheme ONLY
+!!       dmscut      : DMS cutoff value for DMS Anderson (ANDR) sheme ONLY
+!!       dmsslp      : DMS slope value for DMS Anderson (ANDR) sheme ONLY
 !! UKESM1 - exported Chl to UM
 !!       scl_chl     : scaling factor to tune the chl field sent to the UM
@@ -948 +955 @@
                &   ' choice of DMS model passed to atmosphere,                  jdms_model  = Anderson et al. (2001; modified)'
             endif
-         
+            if (jdms_model .eq. 1) then
+               WRITE(numout,*)     &
+               &   ' Anderson DMS model tuned parameters:                       DMS minimum = ',dmsmin,'. -- Default = 2.29 '
+               WRITE(numout,*)     &
+               &   ' Anderson DMS model tuned parameters:                       DMS cutoff  = ',dmscut,'. -- Default = 1.72 '
+               WRITE(numout,*)     &
+               &   ' Anderson DMS model tuned parameters:                       DMS slope   = ',dmsslp,'. -- Default = 8.24 '
+            endif
+
             WRITE(numout,*) '=== UKESM1 --  coupled Chl to the atmosphere'
             WRITE(numout,*)        &