c====================== block data file "kmixbdta.h" =======================
c
c     data statements needed for "kmix" vertical mixing subroutines
c     in "kmixs.F"
c
c
c-----------------------------------------------------------------------
c     parameters for several subroutines
c-----------------------------------------------------------------------
c
c
      data epsilon_def                 /   0.1    /
      data vonk_def                    /   0.4    /
      data conc1_def                   /   5.0    /
      data conam_def,concm_def,conc2_def,zetam_def / 1.257,8.380,16.0,- 0.2 /
      data conas_def,concs_def,conc3_def,zetas_def /-28.86,98.96,16.0,- 1.0 /
c
c
c-----------------------------------------------------------------------
c     parameters for subroutine "bldepth"
c-----------------------------------------------------------------------
c
c
c     Ricr    = critical bulk Richardson Number            = 0.3
c     hbf     = fraction of bounadry layer depth to 
c               which absorbed solar radiation 
c               contributes to surface buoyancy forcing    = 1.0
c
      data Ricr_def      /  0.3  /
      data cekman_def    /  0.7  /
      data cmonob_def    /  1.0  /
      data concv_def     /  1.8  /
      data hbf_def       /  1.0  /
c
c
c-----------------------------------------------------------------------
c     limits for lookup table of wm and ws in subroutine "kmixinit"
c     (zmin and zmax are in m3/s3, umin and umax are in m/s)
c-----------------------------------------------------------------------
c
c
      data zmin_def,zmax_def / -4.e-7 ,  0.0  /
      data umin_def,umax_def /  0.    ,   .04 /
c
c
c-----------------------------------------------------------------------
c     parameters for subroutine "ri_iwmix"
c     to compute vertical mixing coefficients below boundary layer:
c-----------------------------------------------------------------------
c
c
c-----------------------------------------------------------------------
c     Riinfty         = local Richardson Number limit for shear instability
c     num_v_smooth_Ri = number of times Ri is vertically smoothed
c     BVSQcon         = Brunt-Vaisala squared    -- units: 1/s**2
c-----------------------------------------------------------------------
c
      data Riinfty_def               /  0.7   /
      data BVSQcon_def               / -0.2e-4/
      data num_v_smooth_Ri_def       /  0     /
c
c-----------------------------------------------------------------------
c   *  The vertical mixing coefficients are defined
c      in (m2/s) units in "kmixinit". they are initialized
c      and read via namelist "eddy" in (cm2/s)
c      units using their c-g-s names)
c      (m2/s)                                                 (cm2/s)
c      difm0   = viscosity max due to shear instability     = vvcric
c      difs0   = diffusivity ..                             = vdcric
c      difmiw  = viscosity background due to internal waves = fkpm
c      difsiw  = diffusivity ..                             = fkph
c      difmcon = viscosity due to convective instability    = vvclim
c      difscon = diffusivity ..                             = vdclim
c-----------------------------------------------------------------------
c
c    the following two parameters are resolution dependent: (blkdta.h)
      data fkpm_def      /   16.7  /
      data fkph_def      /    0.5  /

      data vvcric_def    /   50.0  /
      data vdcric_def    /   50.0  /
      data vvclim_def    / 1000.0  /
      data vdclim_def    / 1000.0  /


c
c-----------------------------------------------------------------------
c     parameters for subroutine "blmix"
c-----------------------------------------------------------------------
c
c
c-----------------------------------------------------------------------
c     to compute mixing within boundary layer:
c
c     cstar   = proportionality coefficient
c               for nonlocal transport
c-----------------------------------------------------------------------
      data cstar_def     /   10.0   /

c
c
c     jwatertype = jerlov water type: I  IA IB II III
c                                     1  2  3  4  5
c-------------------------------------------------------------------------
c
      data jwatertype_def   /   1   /