Changeset 7894 for branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/mocsy_wrapper.F90

Timestamp:

2017-04-11T13:58:38+02:00 (7 years ago)

Author:

jpalmier

Message:

JPALM -- 11-04-2017 -- MEDUSA spring tidy-up refreshning session

File:

• branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/mocsy_wrapper.F90 (modified) (4 diffs)

branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/MEDUSA/mocsy_wrapper.F90

@@ -6 +6 @@
    !! History :
    !!  -   !  2015-06  (A. Yool)             added for UKESM project
-   !!----------------------------------------------------------------------
-#if defined key_medusa && defined key_roam && defined key_mocsy
+   !!  -   !  2017-04  (A. Yool)             alter optK1K2 to 'w14' option
+   !!----------------------------------------------------------------------
+#if defined key_medusa && defined key_roam
    !!----------------------------------------------------------------------
    !!                                        MEDUSA carbonate chemistry
    !!----------------------------------------------------------------------
-   !!   trc_co2_medusa        :  
+   !!   mocsy_wrapper
    !!----------------------------------------------------------------------
 
@@ -218 +219 @@
 ! present)
 !
+!=======================================================================
+!
+! AXY (05/04/17): alter options to include optK1K2 = 'w14'
+!
+! In conversation with Jim Orr, Waters (2014) is now the 
+! preferred option for optK1K2 in the case of global scale 
+! simulations since this formulation works over broader ranges
+! of temperature (0 < T < 50) and salinity (1 < S < 50).
+!
+! NOTE: *contrary* to the notice above, MEDUSA has now been
+! revised to remove the PML carbonate chemistry routine as part
+! of a wider code review. This routine is not anticipated to be 
+! used again, and has several out of date parameterisations 
+! relative to MOCSY. As a result, key_mocsy is no longer required.
+!
         USE mocsy_singledouble
         IMPLICIT NONE
@@ -249 +265 @@
 !  CONFIGURE OPTIONS
 !> ======================================================================
-!> OPTIONS: see complete documentation in 'vars' subroutine
+!> OPTIONS: see complete documentation in 'vars' subroutine (in mocsy_mainmod.F90)
+!> AXY (05/04/17): optK1K2 switched from 'm10' to 'w14'
 !> Typical options for MODELS
    optCON  = 'mol/m3' ! input concentrations are in MOL/M3
@@ -255 +272 @@
    optP    = 'm'      ! input variable 'depth' is in METERS
    optB    = 'l10'    ! Lee et al. (2010) formulation for total boron
-   optK1K2 = 'l'      ! Lueker et al. (2000) formulations for K1 & K2 (best practices)
+   optK1K2 = 'w14'    ! Lueker et al. (2000) formulations for K1 & K2 (best practices)
    optKf   = 'dg'     ! Dickson & Riley (1979) formulation for Kf (recommended by Dickson & Goyet, 1994)
+!> optK1K2 = 'l'      ! Lueker et al. (2000) formulations for K1 & K2 (best practices)
+!> optKf   = 'dg'     ! Dickson & Riley (1979) formulation for Kf (recommended by Dickson & Goyet, 1994)
 
 !> ======================================================================