Changeset 7905 for branches/UKMO/r6232_HZG_WAVE-coupling/NEMOGCM/CONFIG/SHARED/namelist_ref

Timestamp:

2017-04-13T12:07:16+02:00 (7 years ago)

Author:

jcastill

Message:

Series of small bug fixes and stetic changes:

-Fix possible bug in the calculation of Stokes-Coriolis
-Move all the wave control variables to namelist namsbc_wave
-Use one namelist variable instead of two to set Stokes drift velocity coupling
-Cap the values of the Craig and Banner constant as calculated from wave input fields to take into account small values of the friction velocity
-Add new Phillips parametrization for Stokes drift vertical velocity, using the inverse depth scale as in Breivik 2015, instead of the peak wave number as calculated from wave input fields
-Better control of the wave fields that are read from file depending on the wave parameters

File:

• branches/UKMO/r6232_HZG_WAVE-coupling/NEMOGCM/CONFIG/SHARED/namelist_ref (modified) (7 diffs)

branches/UKMO/r6232_HZG_WAVE-coupling/NEMOGCM/CONFIG/SHARED/namelist_ref

@@ -266 +266 @@
                            !     =2 annual global mean of e-p-r set to zero
    ln_wave = .false.       !  Activate coupling with wave (T => fill namsbc_wave)  
-   ln_cdgw = .false.       !  Neutral drag coefficient read from wave model (T => ln_wave=.true. & fill namsbc_wave)  
-   ln_sdw  = .false.       !  Read 2D Surf Stokes Drift & Computation of 3D stokes drift (T => ln_wave=.true. & fill namsbc_wave)   
-   ln_tauoc= .false.       !  Activate ocean stress modified by external wave induced stress (T => ln_wave=.true. & fill namsbc_wave)  
-   ln_phioc= .false.       !  Activate wave to ocean energy (T => ln_wave=.true. & fill namsbc_wave)  
-   ln_stcor= .false.       !  Activate Stokes Coriolis term (T => ln_wave=.true. & ln_sdw=.true. & fill namsbc_wave)  
    nn_lsm  = 0             !  =0 land/sea mask for input fields is not applied (keep empty land/sea mask filename field) ,
                            !  =1:n number of iterations of land/sea mask application for input fields (fill land/sea mask filename field)
@@ -278 +273 @@
                            !  = 1  Average and redistribute per-category fluxes, forced mode only, not yet implemented coupled
                            !  = 2  Redistribute a single flux over categories (coupled mode only)
-   nn_drag   = 0           !  formula to calculate momentum from the wind components
-                           !  = 0 UKMO SHELF formulation
-                           !  = 1 standard formulation with forced of coupled drag coefficient
-                           !  = 2 standard formulation with constant drag coefficient
-                           !  = 3 momentum calculated from core forcing fields
-   nn_sdrift = 0           !  Parameterization for the calculation of 3D-Stokes drift from the surface Stokes drift
-                           !  = 0 Breivik 2015 parameterization: v_z=v_0*[exp(2*k*z)/(1-8*k*z)]
-                           !  = 1 Phillips:                      v_z=v_o*[exp(2*k*z)-beta*sqrt(-2*k*pi*z)*erfc(sqrt(-2*k*z))]
 /
 !-----------------------------------------------------------------------
@@ -394 +381 @@
    sn_rcv_mslp   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
    sn_rcv_phioc  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
-   sn_rcv_sdrfx  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
-   sn_rcv_sdrfy  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
+   sn_rcv_sdrft  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
    sn_rcv_wper   =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
    sn_rcv_wfreq  =       'none'                 ,    'no'    ,     ''      ,         ''          ,   '' 
@@ -948 +934 @@
    ln_zdfexp   = .false.   !  time-stepping: split-explicit (T) or implicit (F) time stepping
    nn_zdfexp   =    3            !  number of sub-timestep for ln_zdfexp=T
-   ln_zdfqiao  = .false.   !  Enhanced wave vertical mixing Qiao (2010) (T => ln_wave=.true. & ln_sdw=.true. & fill namsbc_wave)
 /
 !-----------------------------------------------------------------------
@@ -1022 +1007 @@
    nn_stab_func  =     2   !  stability function (0=Galp, 1= KC94, 2=CanutoA, 3=CanutoB)
    nn_clos       =     1   !  predefined closure type (0=MY82, 1=k-eps, 2=k-w, 3=Gen)
-   nn_wmix       =     0   !  type of wave breaking mixing
-           !                             ! = 0 Craig and Banner formulation (original NEMO formulation) 
-           !                             ! = 1 Janssen formulation (no assumption of direct energy conversion) 
 /
 !-----------------------------------------------------------------------
@@ -1302 +1284 @@
 !              !  file name  ! frequency (hours) ! variable     ! time interp. !  clim   ! 'yearly'/ ! weights  ! rotation ! land/sea mask !
 !              !             !  (if <0  months)  !   name       !   (logical)  !  (T/F)  ! 'monthly' ! filename ! pairing  ! filename      !
+   ln_cdgw     = .false.   !  Neutral drag coefficient read from wave model
+   ln_sdw      = .false.   !  Read 2D Surf Stokes Drift & Computation of 3D stokes drift
+   ln_stcor    = .false.   !  Activate Stokes Coriolis term
+   ln_tauoc    = .false.   !  Activate ocean stress modified by external wave induced stress
+   ln_phioc    = .false.   !  Activate wave to ocean energy
+   ln_rough    = .false.   !  Wave roughness equals the significant wave height
+   ln_zdfqiao  = .false.   !  Enhanced wave vertical mixing Qiao (2010)
+   nn_drag   = 0      !  formula to calculate momentum from the wind components
+                           ! = 0 UKMO SHELF formulation
+                           ! = 1 standard formulation with forced of coupled drag coefficient
+                           ! = 2 standard formulation with constant drag coefficient
+                           ! = 3 momentum calculated from core forcing fields
+   nn_sdrift   =  0   !  Parameterization for the calculation of 3D-Stokes drift from the surface Stokes drift
+                           ! = 0 Breivik 2015 parameterization: v_z=v_0*[exp(2*k*z)/(1-8*k*z)]
+                           ! = 1 Phillips:                      v_z=v_o*[exp(2*k*z)-beta*sqrt(-2*k*pi*z)*erfc(sqrt(-2*k*z))]
+   nn_wmix     =  0   !  type of wave breaking mixing
+                           ! = 0 Craig and Banner formulation (original NEMO formulation) 
+                           ! = 1 Janssen formulation (no assumption of direct energy conversion) 
    sn_cdg      =  'sdw_wave' ,        1          , 'drag_coeff' ,     .true.   , .false. , 'daily'   ,  ''      , ''       , ''
    sn_usd      =  'sdw_wave' ,        1          , 'u_sd2d'     ,     .true.   , .false. , 'daily'   ,  ''      , ''       , ''
@@ -1311 +1311 @@
    sn_tauoc    =  'sdw_wave' ,        1          , 'wave_stress',     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
    sn_phioc    =  'sdw_wave' ,        1          , 'wave_energy',     .true.   , .false. , 'daily'   ,  ''      , ''       , '' 
+   cn_dir      = './'  !  root directory for the location of wave forcing files
 ! 
-   cn_dir      = './'  !  root directory for the location of drag coefficient files
 /
 !-----------------------------------------------------------------------