module mod_parameters
    implicit none
    !!Constants
    real(8)   , parameter  :: Pi             = 4.0d0*atan(1.0d0)
    real(8)   , parameter  :: Grav_Const     = 6.67408 * (1.0d-11)
    !! Conversions
    real(8)   , parameter  :: gcm3_to_kgm3   = 1000.0d0
    real(8)   , parameter  :: kgm3_to_gcm3   = 1.d-3
    real(8)   , parameter  :: m_to_km        = 1.d-3
    real(8)   , parameter  :: km_to_m        = 1.d3
    real(8)   , parameter  :: GPa_to_Pa      = 1.0d9
    real(8)   , parameter  :: Pa_to_GPa      = 1.0d-9
    real(8)   , parameter  ::  W_to_GW       = 1.0d-9
    real(8)   , parameter  ::  W_to_TW       = 1.0d-12
    real(8)   , parameter  :: rad_to_deg     = 180.0d0/Pi
    real(8)   , parameter  :: pjmstr         = 1.0d0/dsqrt(3.0d0)
    real(8)   ,  parameter :: day2sec        = 86400d0
    !!Zeros
    complex(8), parameter  :: Im             = dcmplx(0.0d0, 1.0d0)
    complex(8), parameter  :: dzero          = dcmplx(0.0d0, 0.0d0)
    real(8)   , parameter  :: zero           = 0.0d0
    !! Numerical parameters
    integer   , parameter  :: Max_num_period = 2000
    integer   , parameter  :: Max_jmax       = 400
    integer   , parameter  :: sph_lo_diag    = 12
    integer   , parameter  :: sph_hi_diag    = 9
    integer   , parameter  :: tor_lo_diag    = 3
    integer   , parameter  :: tor_hi_diag    = 4
    integer   , parameter  :: sph_band       = sph_lo_diag + sph_hi_diag + 1
    integer   , parameter  :: tor_band       = tor_lo_diag + tor_hi_diag + 1
    integer   , parameter  :: sph_diag       = sph_lo_diag + 1
    integer   , parameter  :: tor_diag       = tor_lo_diag + 1
    real(8)   , parameter  :: CN_theta       = 1.0d0
    real(8)   , parameter  :: CFL            = 1.0d0

    !! EXPERIMENTAL: Bulk viscosity of the ocean for compressible models
    real(8)   , parameter  :: bulk_vis = 1d-3
    real(8)   , parameter  ::  semimajor_a = 1070400.d3
end module mod_parameters