ghdtheory_mod.f

Go to the documentation of this file.

  MODULE ghdtheory


      Use param
      Use param1

!
!     Zeroth order dissipation term
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: zeta0
!
!     cooling rate transport coefficient (1st order)
      DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: zetau
!
!     Thermal diffusivity DiT
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: dit
!
!     Mass mobility coefficient
      DOUBLE PRECISION, DIMENSION(:, :, :), ALLOCATABLE :: dijf
!
!     Thermal mobility
      DOUBLE PRECISION, DIMENSION(:, :, :), ALLOCATABLE :: lij
!
!     Ordinary diffusion
      DOUBLE PRECISION, DIMENSION(:, :, :), ALLOCATABLE :: dij
!
!     Dufour coefficient
      DOUBLE PRECISION, DIMENSION(:, :, :), ALLOCATABLE :: dijq
!
!     Species mass flux in X-direction
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: joix
!
!     Species mass flux in Y-direction
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: joiy
!
!     Species mass flux in Z-direction
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: joiz
!
!     external force in X-direction
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: fix
!
!     external force in Y-direction
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: fiy
!
!     external force in Z-direction
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: fiz

!     external force in X-direction--flux
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: fixvel
!
!     external force in Y-direction--flux
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: fiyvel
!
!     external force in Z-direction--flux
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: fizvel
!
!     external force Minus Drag term in Y-direction
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: fiminusdragx
!
!     Species mass flux Without Drag term in Y-direction
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: joiminusdragx
!
!     external force Minus Drag term in Y-direction
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: fiminusdragy
!
!     Species mass flux Without Drag term in Y-direction
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: joiminusdragy
!
!     external force Minus Drag term in Z-direction
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: fiminusdragz
!
!     Species mass flux Without Drag term in Z-direction
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: joiminusdragz


!     Everything that does not depend on the velocity of ith particle
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: deltau
!
!     Everything that does not depend on the velocity of the ith particle
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: deltav
!
!     Everything that does not depned on the velocity of the ith particle
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: deltaw

!     Everything that does not depend on the velocity of ith particle
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: del_dot_j

!     drag force in X-direction
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: dragfx
!
!     drag force in Y-direction
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: dragfy
!
!     drag force in Z-direction
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: dragfz
!     drag force in X-direction
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: dragfxflux
!
!     drag force in Y-direction
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: dragfyflux
!
!     drag force in Z-direction
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: dragfzflux
!     drag force in X-direction
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: beta_cell_x
!
!     drag force in Y-direction
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: beta_cell_y
!
!     drag force in Z-direction
      DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: beta_cell_z
!     drag force in X-direction
      DOUBLE PRECISION, DIMENSION(:, :, :), ALLOCATABLE :: beta_ij_cell_x
!
!     drag force in Y-direction
      DOUBLE PRECISION, DIMENSION(:, :, :), ALLOCATABLE :: beta_ij_cell_y
!
!     drag force in Z-direction
      DOUBLE PRECISION, DIMENSION(:, :, :), ALLOCATABLE :: beta_ij_cell_z

!       decide whether to do harmonic or arithmetic thermal diffusivity east
       LOGICAL, DIMENSION(:), ALLOCATABLE :: dit_harme
!       decide whether to do harmonic or arithmetic thermal diffusivity north
       LOGICAL, DIMENSION(:), ALLOCATABLE :: dit_harmn
!       decide whether to do harmonic or arithmetic thermal diffusivity top
       LOGICAL, DIMENSION(:), ALLOCATABLE :: dit_harmt

!       decide whether to do harmonic or arithmetic mass mobility east
       LOGICAL, DIMENSION(:,:), ALLOCATABLE :: dijf_harme
!       decide whether to do harmonic or arithmetic mass mobility north
       LOGICAL, DIMENSION(:,:), ALLOCATABLE :: dijf_harmn
!       decide whether to do harmonic or arithmetic mass mobility top
       LOGICAL, DIMENSION(:,:), ALLOCATABLE :: dijf_harmt

!       decide whether to do harmonic or arithmetic ordinary diffusivity east
       LOGICAL, DIMENSION(:,:), ALLOCATABLE :: dij_harme
!       decide whether to do harmonic or arithmetic ordinary diffusivity north
       LOGICAL, DIMENSION(:,:), ALLOCATABLE :: dij_harmn
!       decide whether to do harmonic or arithmetic ordinary diffusivity top
       LOGICAL, DIMENSION(:,:), ALLOCATABLE :: dij_harmt

      END MODULE ghdtheory