! mod_wim_prams.F
      ! Author: Timothy Williams
      ! Date:   20141024, 14:03:11 CEST
      ! parameters needed for WIM

      module mod_wim_prams
      implicit none

      !!ice/water parameters
      real,    parameter   :: gravity      = 9.81   ! g (m/s^2)  
      real,    parameter   :: rhowtr       = 1025   ! density of water (kg/m^3)
      real,    parameter   :: rhoice_waves = 922.5  ! density of ice (kg/m^3)

      real,    parameter   :: brine_vol    = 0.1    !(no units) volume fraction of brine (Williams et al, 2013a) 
      real,    parameter   :: poisson      = 0.3    !(no units)
      real,    parameter   :: visc_RP      = 13     !(Pa.s/m) Robinson-Palmer viscosity (Williams et al, 2013a)
      real,    parameter   :: sigc0        = 1.76e6 ! Fresh water flexural strength (Pa)
      real,    parameter   :: young0       = 10.0e9 ! Fresh water Young's modulus (Pa)

      !!FSD parameters
      real,    parameter   :: mu           = 1.0    ! Fatigue coefficient
      integer, parameter   :: xi           = 2      ! fragmentation fraction
      real,    parameter   :: fragility    = 0.9    ! Toyota et al. (in press)
      
      !!options for determining breaking strain, Young's modulus
      integer, parameter   :: BRK_OPT     = 1   !1: Williams et al (2013);            2: Marchenko's COMSOL est
      integer, parameter   :: YOUNG_OPT   = 1   !1: Vernon's est for Young's modulus; 2: Marchenko's est

      !!need setting (waves_init, mod_waves_init.F)
      real                 :: epsc      ! breaking strain
      real                 :: sigc      ! flexural strength (Pa)
      real                 :: young     ! Young's modulus (Pa)
      !real, dimension(3)   :: prams_RP  ! for Robinson-Palmer dispersion relation

      end module mod_wim_prams