Changeset 9258


Ignore:
Timestamp:
2018-01-18T12:24:31+01:00 (3 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)

Location:
branches/UKMO/dev_r5518_GO6_package/NEMOGCM
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r5518_GO6_package/NEMOGCM/CONFIG/SHARED/namelist_medusa_ref

    r8280 r9258  
    140140!! 
    141141!! UKESM1 - new diagnostics  !! Jpalm 
     142!!       chl_out     : (IF key_oasis) use surf Chl (1) or mld Chl (2) 
     143!!       dmscut      : (if jdms_model=1) modify Anderson et al. 2001 param 
     144!!                     default = 1.72 
     145!!       dmsmin      : (if jdms_model=1) modify Anderson et al. 2001 dms minimum 
     146!!                     default = 2.29 
     147!!       dmsslp      : (if jdms_model=1) modify Anderson et al. 2001 dms slope 
     148!!                     default = 8.24 
    142149!!       jdms        :  include dms diagnostics 
    143150!!       jdms_input  :  use instant (0) or diel-avg (1) inputs 
    144151!!       jdms_model  :  which DMS model is passed to atmosphere 
    145 !! 
     152!!       scl_chl     : (IF key_oasis) scaling the exported to UM chlorophyll field 
    146153!!        
    147154&natbio 
     
    285292!! 
    286293!! UKESM1 - new diagnostics !! Jpalm 
     294      chl_out    = 1        !! (IF key_oasis) use surf Chl (1) or mld Chl (2) 
     295      dmscut     = 1.72     !! (if jdms_model=1) modify Anderson et al. 2001 param 
     296!!                             default = 1.72 
     297      dmsmin     = 2.29     !! (if jdms_model=1) modify Anderson et al. 2001 dms minimum 
     298!!                             default = 2.29 
     299      dmsslp     = 8.24     !! (if jdms_model=1) modify Anderson et al. 2001 dms slope 
     300!!                             default = 8.24 
    287301      jdms       = 1        !! include dms diagnostics 
    288302      jdms_input = 1        !! use instant (0) or diel-avg (1) inputs 
    289303      jdms_model = 3        !! choice of DMS model passed to atmosphere 
    290304!!                             1 = ANDR, 2 = SIMO, 3 = ARAN, 4 = HALL 
     305      scl_chl    = 1.0      !! (IF key_oasis) scaling the exported to UM chlorophyll field 
    291306!! 
    292307/ 
  • branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/TOP_SRC/MEDUSA/sms_medusa.F90

    r9114 r9258  
    177177   INTEGER  ::  jdms_model   !: choice of DMS model passed to atmosphere 
    178178!!                              1 = ANDR, 2 = SIMO, 3 = ARAN, 4 = HALL 
     179!! JPALM --19-12-2017 -- UM people need to tune the Anderson DMS scheme     
     180   REAL(wp) ::  dmsmin       !: Anderson DMS scheme - DMS minimum value 
     181   REAL(wp) ::  dmscut       !: Anderson DMS scheme - DMS cutoff value  
     182   REAL(wp) ::  dmsslp       !: Anderson DMS scheme - DMS slope  
    179183!! FOR UKESM    
    180184   REAL(wp) ::  scl_chl      !: scaling factor for tuned Chl passed to the UM  
  • branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcdms_medusa.F90

    r8132 r9258  
    113113! 
    114114! AXY (13/03/15): Anderson et al. (2001) 
     115!! JPALM --19-12-2017-- Tunable through the namelist 
     116!!                      within dmsmin - dmscut - dmsslp 
    115117        Jterm = xqsr + 1.0e-6 
    116118        !! this next line makes a hard-coded assumption about the  
     
    120122        Qterm = xdin / (xdin + 0.5) 
    121123        fq1 = log10(CHL * Jterm * Qterm) 
    122         if (fq1 > 1.72) then 
    123            dms_andr = (8.24 * (fq1 - 1.72)) + 2.29 
    124         else 
    125            dms_andr = 2.29 
     124        if (fq1 > dmscut) then 
     125           dms_andr = (dmsslp * (fq1 - dmscut)) + dmsmin 
     126        else 
     127           dms_andr = dmsmin 
    126128        endif 
    127129! 
  • branches/UKMO/dev_r5518_GO6_package/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcnam_medusa.F90

    r9114 r9258  
    8787      &  xsdiss,                                              & 
    8888      &  sedlam,sedlostpoc,jpkb,jdms,jdms_input,jdms_model,   & 
    89       &  scl_chl, chl_out 
     89      &  scl_chl, chl_out, dmsmin, dmscut, dmsslp 
    9090#if defined key_roam 
    9191      NAMELIST/natroam/ xthetaphy,xthetazoo,xthetanit,        & 
     
    127127      ! 1.4 namelist natbio : biological parameters 
    128128      ! ------------------------------------------- 
     129      !! Note: the default values below will all be overwritten by the 
     130      !!       input in the namelist natbio. 
    129131       
    130132      xxi         = 0. 
     
    250252      scl_chl     = 1. 
    251253      chl_out     = 1 
     254      dmsmin      = 2.29 !! Anderson DMS default 
     255      dmscut      = 1.72 !! Anderson DMS default 
     256      dmsslp      = 8.24 !! Anderson DMS default 
    252257             
    253258      !REWIND(numnatm) 
     
    403408!!       jdms_model  :  choice of DMS model passed to atmosphere 
    404409!!                      1 = ANDR, 2 = SIMO, 3 = ARAN, 4 = HALL, 5 = ANDM 
    405 !! 
     410!!       dmsmin      : DMS minimum value for DMS Anderson (ANDR) sheme ONLY 
     411!!       dmscut      : DMS cutoff value for DMS Anderson (ANDR) sheme ONLY 
     412!!       dmsslp      : DMS slope value for DMS Anderson (ANDR) sheme ONLY 
    406413!! UKESM1 - exported Chl to UM 
    407414!!       scl_chl     : scaling factor to tune the chl field sent to the UM 
     
    948955               &   ' choice of DMS model passed to atmosphere,                  jdms_model  = Anderson et al. (2001; modified)' 
    949956            endif 
    950           
     957            if (jdms_model .eq. 1) then 
     958               WRITE(numout,*)     & 
     959               &   ' Anderson DMS model tuned parameters:                       DMS minimum = ',dmsmin,'. -- Default = 2.29 ' 
     960               WRITE(numout,*)     & 
     961               &   ' Anderson DMS model tuned parameters:                       DMS cutoff  = ',dmscut,'. -- Default = 1.72 ' 
     962               WRITE(numout,*)     & 
     963               &   ' Anderson DMS model tuned parameters:                       DMS slope   = ',dmsslp,'. -- Default = 8.24 ' 
     964            endif 
     965 
    951966            WRITE(numout,*) '=== UKESM1 --  coupled Chl to the atmosphere' 
    952967            WRITE(numout,*)        & 
Note: See TracChangeset for help on using the changeset viewer.