Changeset 9187 for branches/NERC/dev_r5518_GO6_MEDUSA_DMS_TUNING/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcnam_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/trcnam_medusa.F90 (modified) (4 diffs)

branches/NERC/dev_r5518_GO6_MEDUSA_DMS_TUNING/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,        &
@@ -250 +250 @@
       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 +406 @@
 !!       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 +953 @@
                &   ' 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,*)        &