Changeset 8908 for branches/2017/dev_METO_2017/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_oce.F90

Timestamp:

2017-12-06T10:36:02+01:00 (6 years ago)

Author:

timgraham

Message:

Merged branches/UKMO/r8727_WAVE-2_Clementi_add_coupling into branch

File:

• branches/2017/dev_METO_2017/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_oce.F90 (modified) (2 diffs)

branches/2017/dev_METO_2017/NEMOGCM/NEMO/OPA_SRC/SBC/sbc_oce.F90

@@ -65 +65 @@
    LOGICAL , PUBLIC ::   ln_sdw         !: true if 3d stokes drift from wave model
    LOGICAL , PUBLIC ::   ln_tauoc       !: true if normalized stress from wave is used
+   LOGICAL , PUBLIC ::   ln_tauw        !: true if ocean stress components from wave is used
    LOGICAL , PUBLIC ::   ln_stcor       !: true if Stokes-Coriolis term is used
+   !
+   INTEGER , PUBLIC ::   nn_sdrift      ! type of parameterization to calculate vertical Stokes drift
    !
    LOGICAL , PUBLIC ::   ln_icebergs    !: Icebergs
@@ -79 +82 @@
    INTEGER , PUBLIC, PARAMETER ::   jp_none    = 5        !: for OPA when doing coupling via SAS module
    
+   !!----------------------------------------------------------------------
+   !!           Stokes drift parametrization definition 
+   !!----------------------------------------------------------------------
+   INTEGER , PUBLIC, PARAMETER ::   jp_breivik  = 0     ! Breivik 2015: v_z=v_0*[exp(2*k*z)/(1-8*k*z)]
+   INTEGER , PUBLIC, PARAMETER ::   jp_phillips = 1     ! Phillips: v_z=v_o*[exp(2*k*z)-beta*sqrt(-2*k*pi*z)*erfc(sqrt(-2*k*z))]
+   INTEGER , PUBLIC, PARAMETER ::   jp_peakfr   = 2     ! Phillips using the peak wave number read from wave model instead of the inverse depth scale
+
    !!----------------------------------------------------------------------
    !!           component definition