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

lgw: first pass of usr1 with renorm only

parent 6bf3b1fe

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

+22 −21
Original line number Diff line number Diff line
@@ -50,8 +50,7 @@ 
!  Define local variables here

!

      INTEGER :: ijk, ijmk

      DOUBLE PRECISION :: Voltot, Ncells, renorm_test

      DOUBLE PRECISION :: sum_momy_g, sum_momy_s, sum_momy_tot

      DOUBLE PRECISION :: Vol_tot, flux_tot_y

!

!  Include files defining statement functions here

!
@@ -64,31 +63,33 @@ 
!     Normalize vertical velocity variables to Vk' = Vk - jy

!     ie ' variables are relative to the zero net flux plane



      !number of computational cells for avg

      Ncells = ZERO

      Voltot = ZERO

      tot_momy_sum = ZERO

      renorm_test  = ZERO

      sum_momy_

      flux_tot_y = ZERO

      Vol_tot = ZERO 

      DO IJK = IJKSTART3, IJKEND3

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

        IF (FLUID_AT(IJK)) THEN

        IF (.NOT. FLUID_AT(IJK)) CYCLE

          IJMK = JM_OF(IJK)

!!!          Ncells = Ncells + ONE

!!!          Voltot = Voltot + VOL(IJK)

          tot_momy_sum = tot_momy_sum &

                       + EP_S(IJK,1)*VOL(IJK)*(V_s(IJK,1) +

V_s(IJMK,1))/2.0d0 &

                       + EP_G(IJK)*VOL(IJK)*(V_g(IJK) + V_g(IJMK))/2.0d0

          flux_tot_y = flux_tot_y + 0.5d0 * VOL(IJK) * &

                     ( (1.0d0 - EP_G(IJK))*(V_s(IJK,1) + V_s(IJMK,1)) + &

                              + EP_G(IJK) *(V_g(IJK)   + V_g(IJMK)))

!         Vol_tot = Vol_tot + VOL(IJK) 

        ENDIF

      ENDDO



!!!      CALL GLOBAL_ALL_SUM(Ncells)

!!!      CALL GLOBAL_ALL_SUM(Voltot)

      Ncells = DBLE(IMAX*JMAX*KMAX)

      Voltot = DBLE(XLENGTH*YLENGTH*ZLENGTH)

      CALL GLOBAL_ALL_SUM(tot_momy_sum)

      CALL GLOBAL_ALL_SUM(flux_tot_y)

!     this only works for regular rect domains

!     CALL GLOBAL_ALL_SUM(Voltot) 

      Vol_tot = DBLE(XLENGTH*YLENGTH*ZLENGTH)

      flux_tot_y = flux_tot_y / Vol_tot

      

      IF (flux_tot_y .GT. 1.0d-15)

        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 

          V_g(IJK)   = V_g(IJK)   - flux_tot_y

          V_s(IJK,1) = V_s(IJK,1) - flux_tot_y 

        ENDDO

      ENDIF 





!     check for usr IO time
@@ -102,7 +103,7 @@ V_s(IJMK,1))/2.0d0 & 
!     print

      IF(myPE .NE. PE_IO) RETURN

      OPEN(UNIT=888,FILE='uio_stats.dat',STATUS='OLD',ACCESS='APPEND')

      WRITE(888,8882) TIME, sum_momy_tot

      WRITE(888,8882) TIME, flux_tot_y

      CLOSE(888)



      RETURN