New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 10302 for branches/UKMO/dev_r5518_GO6_package_FOAMv14/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcnam_medusa.F90 – NEMO

Ignore:
Timestamp:
2018-11-13T18:21:16+01:00 (5 years ago)
Author:
dford
Message:

Merge in revisions 8447:10159 of dev_r5518_GO6_package.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/UKMO/dev_r5518_GO6_package_FOAMv14/NEMOGCM/NEMO/TOP_SRC/MEDUSA/trcnam_medusa.F90

    r8442 r10302  
    2424   USE sms_medusa      ! sms trends 
    2525   USE iom             ! I/O manager 
     26   USE sbc_oce, ONLY: lk_oasis 
    2627   !!USE trc_nam_dia     ! JPALM 13-11-2015 -- if iom_use for diag 
    2728 
    … …  
    8586      &  jriver_n,jriver_si,jriver_c,jriver_alk,jriver_dep,   & 
    8687      &  xsdiss,                                              & 
    87       &  sedlam,sedlostpoc,jpkb,jdms,jdms_input,jdms_model 
     88      &  sedlam,sedlostpoc,jpkb,jdms,jdms_input,jdms_model,   & 
     89      &  scl_chl, chl_out, dmsmin, dmscut, dmsslp,            & 
     90      &  ln_foam_medusa 
    8891#if defined key_roam 
    8992      NAMELIST/natroam/ xthetaphy,xthetazoo,xthetanit,        & 
     93      &    xobs_xco2a,                                        & 
    9094      &    xthetarem,xo2min  
    9195#endif 
    … …  
    125129      ! 1.4 namelist natbio : biological parameters 
    126130      ! ------------------------------------------- 
     131      !! Note: the default values below will all be overwritten by the 
     132      !!       input in the namelist natbio. 
    127133       
    128134      xxi         = 0. 
    … …  
    246252      jdms_input  = 0 
    247253      jdms_model  = 0 
     254      scl_chl     = 1. 
     255      chl_out     = 1 
     256      dmsmin      = 2.29 !! Anderson DMS default 
     257      dmscut      = 1.72 !! Anderson DMS default 
     258      dmsslp      = 8.24 !! Anderson DMS default 
     259!! 
     260      ln_foam_medusa = .FALSE. 
    248261             
    249262      !REWIND(numnatm) 
    … …  
    399412!!       jdms_model  :  choice of DMS model passed to atmosphere 
    400413!!                      1 = ANDR, 2 = SIMO, 3 = ARAN, 4 = HALL, 5 = ANDM 
    401 !! 
     414!!       dmsmin      : DMS minimum value for DMS Anderson (ANDR) sheme ONLY 
     415!!       dmscut      : DMS cutoff value for DMS Anderson (ANDR) sheme ONLY 
     416!!       dmsslp      : DMS slope value for DMS Anderson (ANDR) sheme ONLY 
     417!! UKESM1 - exported Chl to UM 
     418!!       scl_chl     : scaling factor to tune the chl field sent to the UM 
     419!!       chl_out     : select the chl field to send at the UM: 
     420!!                     1- Surf Chl ; 2- MLD Chl  
     421!! 
     422!! FOAM - observation operator and data assimilation 
     423!!       ln_foam_medusa : calculate required diagnostics 
     424 
    402425      IF(lwp) THEN 
    403426!! 
    … …  
    909932!! 
    910933!! UKESM1 - new diagnostics  !! Jpalm; AXY (08/07/15) 
    911          WRITE(numout,*) '=== UKESM1-related parameters' 
    912          WRITE(numout,*)     & 
    913          &   ' include DMS diagnostic?,                                   jdms        = ', jdms 
    914          if (jdms_input .eq. 0) then 
    915             WRITE(numout,*)     & 
    916             &   ' use instant (0) or diel-avg (1) inputs,                    jdms_input  = instantaneous' 
    917          else 
    918             WRITE(numout,*)     & 
    919             &   ' use instant (0) or diel-avg (1) inputs,                    jdms_input  = diel-average' 
    920          endif 
    921     if (jdms_model .eq. 1) then 
    922             WRITE(numout,*)     & 
    923             &   ' choice of DMS model passed to atmosphere,                  jdms_model  = Anderson et al. (2001)' 
    924     elseif (jdms_model .eq. 2) then 
    925             WRITE(numout,*)     & 
    926             &   ' choice of DMS model passed to atmosphere,                  jdms_model  = Simo & Dachs (2002)' 
    927     elseif (jdms_model .eq. 3) then 
    928             WRITE(numout,*)     & 
    929             &   ' choice of DMS model passed to atmosphere,                  jdms_model  = Aranami & Tsunogai (2004)' 
    930     elseif (jdms_model .eq. 4) then 
    931             WRITE(numout,*)     & 
    932             &   ' choice of DMS model passed to atmosphere,                  jdms_model  = Halloran et al. (2010)' 
    933     elseif (jdms_model .eq. 5) then 
    934             WRITE(numout,*)     & 
    935             &   ' choice of DMS model passed to atmosphere,                  jdms_model  = Anderson et al. (2001; modified)' 
    936          endif 
     934         WRITE(numout,*) '=== UKESM1-related parameters ===' 
     935         WRITE(numout,*) ' ---- --- ---' 
     936 
     937         IF (lk_oasis) THEN 
     938            WRITE(numout,*) '=== UKESM1 --  coupled DMS to the atmosphere' 
     939            WRITE(numout,*)     & 
     940            &   ' include DMS diagnostic?,                                   jdms        = ', jdms 
     941            if (jdms_input .eq. 0) then 
     942               WRITE(numout,*)     & 
     943               &   ' use instant (0) or diel-avg (1) inputs,                    jdms_input  = instantaneous' 
     944            else 
     945               WRITE(numout,*)     & 
     946               &   ' use instant (0) or diel-avg (1) inputs,                    jdms_input  = diel-average' 
     947            endif 
     948          if (jdms_model .eq. 1) then 
     949               WRITE(numout,*)     & 
     950               &   ' choice of DMS model passed to atmosphere,                  jdms_model  = Anderson et al. (2001)' 
     951       elseif (jdms_model .eq. 2) then 
     952               WRITE(numout,*)     & 
     953               &   ' choice of DMS model passed to atmosphere,                  jdms_model  = Simo & Dachs (2002)' 
     954       elseif (jdms_model .eq. 3) then 
     955               WRITE(numout,*)     & 
     956               &   ' choice of DMS model passed to atmosphere,                  jdms_model  = Aranami & Tsunogai (2004)' 
     957       elseif (jdms_model .eq. 4) then 
     958               WRITE(numout,*)     & 
     959               &   ' choice of DMS model passed to atmosphere,                  jdms_model  = Halloran et al. (2010)' 
     960       elseif (jdms_model .eq. 5) then 
     961               WRITE(numout,*)     & 
     962               &   ' choice of DMS model passed to atmosphere,                  jdms_model  = Anderson et al. (2001; modified)' 
     963            endif 
     964            if (jdms_model .eq. 1) then 
     965               WRITE(numout,*)     & 
     966               &   ' Anderson DMS model tuned parameters:                       DMS minimum = ',dmsmin,'. -- Default = 2.29 ' 
     967               WRITE(numout,*)     & 
     968               &   ' Anderson DMS model tuned parameters:                       DMS cutoff  = ',dmscut,'. -- Default = 1.72 ' 
     969               WRITE(numout,*)     & 
     970               &   ' Anderson DMS model tuned parameters:                       DMS slope   = ',dmsslp,'. -- Default = 8.24 ' 
     971            endif 
     972 
     973            WRITE(numout,*) '=== UKESM1 --  coupled Chl to the atmosphere' 
     974            WRITE(numout,*)        & 
     975               &   ' Scaling factor to export tuned Chl to the atmosphere       scl_chl  = ', scl_chl 
     976            IF (chl_out .eq. 1) THEN 
     977               WRITE(numout,*)        & 
     978               &   ' Chl field to be scaled and sent to the atmosphere:         chl_out  = Surface Chl field ' 
     979            ELSEIF (chl_out .eq. 2) THEN 
     980               WRITE(numout,*)        & 
     981               &   ' Chl field to be scaled and sent to the atmosphere:         chl_out  = MLD Chl field ' 
     982            ENDIF 
     983         ENDIF ! IF lk_oasis=true 
     984!! FOAM 
     985         WRITE(numout,*) '=== FOAM-related parameters' 
     986         WRITE(numout,*)     & 
     987         &   ' calculate diagnostics for data assimilation,            ln_foam_medusa = ', ln_foam_medusa 
    937988!! 
    938989      ENDIF 
    … …  
    9901041      xthetarem = 0. 
    9911042      xo2min    = 0. 
     1043      xobs_xco2a = 0. 
    9921044 
    9931045      !READ(numnatm,natroam) 
    … …  
    10091061!!       xthetarem :  oxygen consumption by carbon remineralisation 
    10101062!!       xo2min    :  oxygen minimum concentration 
     1063!!       xobs_xco2a : observed atmospheric xCO2 (not used if <= 0) 
    10111064 
    10121065      IF(lwp) THEN 
    … …  
    10261079          WRITE(numout,*)     & 
    10271080          &   ' oxygen minimum concentration                               xo2min      = ', xo2min 
     1081          WRITE(numout,*)     & 
     1082          &   ' observed atmospheric xCO2 (not used if <= 0)               xobs_xco2a  = ', xobs_xco2a 
    10281083       ENDIF 
    10291084 
    … …  
    20532108          med_diag%OCN_DPCO2%dgsave = .FALSE. 
    20542109      ENDIF 
    2055       !! 
     2110      !! UKESM additional 
     2111      IF  (iom_use("CHL_MLD")) THEN  
     2112          med_diag%CHL_MLD%dgsave = .TRUE. 
     2113      ELSE  
     2114          med_diag%CHL_MLD%dgsave = .FALSE. 
     2115      ENDIF 
     2116      IF  (iom_use("CHL_CPL")) THEN  
     2117          med_diag%CHL_CPL%dgsave = .TRUE. 
     2118      ELSE  
     2119          med_diag%CHL_CPL%dgsave = .FALSE. 
     2120      ENDIF 
     2121      !! 3D 
    20562122      IF  (iom_use("TPP3")) THEN  
    20572123          med_diag%TPP3%dgsave = .TRUE. 
Note: See TracChangeset for help on using the changeset viewer.