particles.f90

Go to the documentation of this file.

!********************************************************************************
!********************************************************************************
MODULE particles_module
  !
  IMPLICIT NONE

  REAL*8 :: min_phi

  REAL*8 :: max_phi

  REAL*8 :: diam2

  REAL*8 :: max_rho

  REAL*8 :: diam1

  REAL*8 :: min_rho

  INTEGER :: no_bins

  REAL*8 :: sigma_phi

  REAL*8 :: mean_phi

  REAL*8, ALLOCATABLE, DIMENSION(:) :: fracsolid 
  
  REAL*8, ALLOCATABLE, DIMENSION(:) :: init_fracsolid 

  REAL*8, ALLOCATABLE, DIMENSION(:) :: rhosol

  REAL*8, ALLOCATABLE, DIMENSION(:) :: diam

  REAL*8, ALLOCATABLE, DIMENSION(:) :: C_d

  REAL*8 :: C_d_input

  REAL*8, ALLOCATABLE, DIMENSION(:) :: sepsolid

  REAL*8, ALLOCATABLE, DIMENSION(:) :: solidmassflux_fract

  REAL*8, ALLOCATABLE, DIMENSION(:) :: v_s

  REAL*8, ALLOCATABLE, DIMENSION(:) :: S 
  
  REAL*8 :: Acc_rate
  
  REAL*8 :: Total_Acc_rate
  
  REAL*8 :: Pn

  REAL*8 :: k = 0.4D0

  REAL*8 :: ustar !> shear velocity of flow
  
  REAL*8 :: h_bl

  REAL*8 :: ks = 1.D0
  
  REAL*8 :: eta0
  
  REAL*8 :: wavelengthfreen 
  
  REAL*8 :: wavelengthshearn
  
  REAL*8 :: wavelengthfreex 
  
  REAL*8 :: wavelengthshearx
  
  REAL*8 :: wavelength_height 
  
  REAL*8, ALLOCATABLE, DIMENSION(:) :: eta
  
  REAL*8, ALLOCATABLE, DIMENSION(:) :: BV
  
  REAL*8, ALLOCATABLE, DIMENSION(:) :: SS0
  
  REAL*8, ALLOCATABLE, DIMENSION(:) :: norm_flux
  
  REAL*8, ALLOCATABLE, DIMENSION(:) :: ipf
  
  REAL*8, ALLOCATABLE, DIMENSION(:) :: fpf

  INTEGER :: iclass
  INTEGER :: final = 9.D0
   
  REAL*8 :: sumsed

  REAL*8 :: delta_phi
  REAL*8 :: phi_left
  REAL*8 :: phi_right
  REAL*8, ALLOCATABLE, DIMENSION(:) :: phi_dist
  REAL*8 :: left_tail_fract 
  REAL*8 :: right_tail_fract
  REAL*8 :: rho_one

  !
  SAVE
  !----------------------------------------------------------------------
CONTAINS
  !----------------------------------------------------------------------

  !******************************************************************************
  !******************************************************************************
  SUBROUTINE initialize_particles
    INTEGER :: i

    no_bins = iclass

    ! Define a normal distribution

    delta_phi= (max_phi - min_phi) / no_bins

    phi_left = min_phi 

    DO i=1,no_bins

       phi_dist(i) = phi_left + 0.5d0 * delta_phi

       phi_right = min_phi + (i) * delta_phi

       ! Convert grainsize to m   
       diam(i) = 1.d-3 * 2.0**(-phi_dist(i))
    
       IF ( diam(i) .LE. diam1 ) THEN

          rhosol(i) = min_rho

       ELSEIF ( diam(i) .LE. diam2 ) THEN

          rhosol(i) = min_rho + ( diam(i) - diam1 ) / ( diam2 - diam1 ) *       &
               ( max_rho - min_rho)

       ELSE
          
          rhosol(i) = max_rho
   
       END IF
       
       c_d(i) = c_d_input

       fracsolid(i) = cdf(phi_right) - cdf(phi_left)
       
       phi_left = phi_right

    END DO

    left_tail_fract= cdf(min_phi)
    right_tail_fract= 1 - cdf(max_phi)

    If (left_tail_fract .GT. fracsolid(1)) THEN

       fracsolid(1) = fracsolid(1) + left_tail_fract

    END IF

    If (right_tail_fract .GT. fracsolid(no_bins)) THEN

       fracsolid(no_bins) = fracsolid(no_bins) + right_tail_fract

    END IF

    fracsolid(1:no_bins) = fracsolid(1:no_bins)/sum(fracsolid(1:no_bins))
    
    RETURN
    
  END SUBROUTINE initialize_particles

  !******************************************************************************
  !
  !******************************************************************************

  FUNCTION cdf(phi_in)
    !
    IMPLICIT NONE

    REAL*8 :: cdf
    REAL*8, INTENT(IN) :: phi_in

    cdf = 0.5d0 * (1.d0 + erf((phi_in - mean_phi)/(sigma_phi * sqrt(2.d0))))

  END FUNCTION cdf

  !******************************************************************************
  !******************************************************************************
  SUBROUTINE allocate_particles
    !
    IMPLICIT NONE
    !
    ALLOCATE(solidmassflux_fract(iclass))
    ALLOCATE(v_s(iclass))
    ALLOCATE(s(iclass))
    ALLOCATE(diam(iclass))
    ALLOCATE(rhosol(iclass))
    ALLOCATE(fracsolid(iclass))
    ALLOCATE(phi_dist(iclass))
    ALLOCATE(c_d(iclass))
    ALLOCATE(norm_flux(iclass))
    ALLOCATE(ipf(iclass))
    ALLOCATE(fpf(iclass))
    ALLOCATE(eta(final))
    ALLOCATE(ss0(final))
    ALLOCATE(bv(final))    
    !
    RETURN
  END SUBROUTINE allocate_particles
  !----------------------------------------------------------------------
END MODULE particles_module
!----------------------------------------------------------------------