plume.f90

Go to the documentation of this file.

!********************************************************************************
!
!********************************************************************************
MODULE plume_module
  !

  USE variables, ONLY : wp

  IMPLICIT NONE
  !
  REAL(wp) :: s
  REAL(wp) :: x
  REAL(wp) :: y
  REAL(wp) :: z
  REAL(wp) :: r
  REAL(wp) :: u
  REAL(wp) :: v
  REAL(wp) :: w
  REAL(wp) :: mag_u
  REAL(wp) :: phi
  REAL(wp) :: rp
  REAL(wp) :: alpha_inp
  REAL(wp) :: beta_inp
  REAL(wp) :: prob_factor
  LOGICAL :: particles_loss

  !
  REAL(wp) :: vent_height
  REAL(wp) :: w0
  REAL(wp) :: r0
  REAL(wp) :: log10_mfr
  !
  SAVE

CONTAINS

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

  SUBROUTINE initialize_plume

    IMPLICIT NONE

    x = 0.0_wp
    y = 0.0_wp
    z = vent_height
    s = 0.0_wp
    r = r0
    u = 1.0e-5_wp
    v = 1.0e-5_wp
    w = w0

    mag_u = sqrt(u*u+v*v+w*w)
    phi = atan(w/u)

    

    RETURN

  END SUBROUTINE initialize_plume

END MODULE plume_module