Commit 193b76fa authored by William D. Fullmer's avatar William D. Fullmer
Browse files

lgw bug fix

parent 645d4e32

liu.ge.wang/build.19.3.gnu/usr1.f

+1 −1
Original line number Diff line number Diff line
@@ -82,7 +82,7 @@ 
      Vol_tot = DBLE(XLENGTH*YLENGTH*ZLENGTH)

      flux_tot_y = flux_tot_y / Vol_tot

      

      IF (flux_tot_y .GT. 1.0d-15) THEN

      IF (DABS(flux_tot_y) .GT. 1.0d-15) THEN

        DO IJK = IJKSTART3, IJKEND3

          IF(.NOT.IS_ON_myPE_wobnd(I_OF(IJK), J_OF(IJK), K_OF(IJK))) CYCLE

          IF (.NOT.FLUID_AT(IJK)) CYCLE

tools/calc_dp_ufs.f

0 → 100755
+49 −0
Original line number Diff line number Diff line 
      PROGRAM calc_DP

!---------------------------------------------------------------------------------

!

! Purpose:

!       Calculate the body force and pressure drop for the gas-solid

!       sedimentation problems

!

!---------------------------------------------------------------------------------

      IMPLICIT NONE

!---------------------------------------------------------------------------------

!     Local variables

!---------------------------------------------------------------------------------

      INTEGER :: i, j

      DOUBLE PRECISION :: rho_g, rho_s, Drho, rrho, mu_g

      DOUBLE PRECISION :: phi, diap, Ar, grav, Ly, DP

!---------------------------------------------------------------------------------

!     set variables

!---------------------------------------------------------------------------------

      rho_s = 1.5d0

      rho_g = 1.3d-3

      diap  = 7.5d-3

      mu_g  = 1.8d-4

      phi   = 0.10d0

      Ly    = 2.304d0

      grav  = 9.8d+3

!---------------------------------------------------------------------------------

!     read in

!---------------------------------------------------------------------------------

!---------------------------------------------------------------------------------

!     Calc other variables

!       Ar = rho_g*Drho*grav*diap**3 / mu_g**2

!---------------------------------------------------------------------------------

      rrho  = rho_s/rho_g

      Drho  = rho_s - rho_g

      Ar    = rho_g*Drho*diap**3/mu_g**2

      DP    = (phi*rho_s + (1.0d0 - phi)*rho_g)*grav*Ly

!--------------------------------------------------------------------------------

!     Write out

!---------------------------------------------------------------------------------

    6 FORMAT (A9,F12.6)

   16 FORMAT (A9,F24.16)

      WRITE(*,*)   '  '

      WRITE(*,16)  ' DP    = ', DP

      WRITE(*,*)   '  '

!

!

!

      END PROGRAM calc_DP