FDS Source: use local VEL_G array for particle drag interpolation

rmcdermo · Sep 18, 2024 · 1933bec · 1933bec
1 parent 3f8a838
commit 1933bec
Showing 1 changed file with 6 additions and 6 deletions.
diff --git a/Source/part.f90 b/Source/part.f90
@@ -2522,8 +2522,8 @@ SUBROUTINE MOVE_IN_GAS
             EMBER_DENSITY,EMBER_VOLUME=0._EB,ACCEL_X,ACCEL_Y,ACCEL_Z,&
             LP_FORCE,FACE_VOLS(2,2,2),VEL_G_INT(3),VOL_WGT(2,2,2),&
             EMBER_PACKING_RATIO,LOCAL_PACKING_RATIO,LPC_GEOM_FACTOR
-REAL(EB) :: WGT(2,2,2,3)
-REAL(EB), POINTER, DIMENSION(:,:,:) :: FV_D=>NULL(),VEL_G=>NULL()
+REAL(EB) :: WGT(2,2,2,3),VEL_G(2,2,2)
+REAL(EB), POINTER, DIMENSION(:,:,:) :: FV_D=>NULL()
 REAL(EB), SAVE :: BETA
 INTEGER :: IIX,JJY,KKZ,IL,JL,KL,AXIS,N_LPC2
 LOGICAL :: STUCK=.FALSE.
@@ -2559,21 +2559,21 @@ SUBROUTINE MOVE_IN_GAS
 DO AXIS=IAXIS,KAXIS
    IL = IIX; JL = JJY; KL = KKZ
    IF (AXIS==IAXIS) THEN
-      VEL_G => U
       IL = FLOOR(XI)
+      VEL_G = U(IL:IL+1,JL:JL+1,KL:KL+1)
    ELSEIF (AXIS==JAXIS) THEN
-      VEL_G => V
       JL = FLOOR(YJ)
+      VEL_G = V(IL:IL+1,JL:JL+1,KL:KL+1)
    ELSEIF (AXIS==KAXIS) THEN
-      VEL_G => W
       KL = FLOOR(ZK)
+      VEL_G = W(IL:IL+1,JL:JL+1,KL:KL+1)
    ENDIF
    IF (ICC>0) THEN
       CALL GET_FACE_IDW(AXIS,IL,JL,KL,BC%X,BC%Y,BC%Z,WGT(:,:,:,AXIS))
    ELSE
       CALL GET_FACE_TLW(AXIS,IL,JL,KL,BC%X,BC%Y,BC%Z,WGT(:,:,:,AXIS))
    ENDIF
-   VEL_G_INT(AXIS) = SUM(VEL_G(IL:IL+1,JL:JL+1,KL:KL+1)*WGT(:,:,:,AXIS))
+   VEL_G_INT(AXIS) = SUM(VEL_G*WGT(:,:,:,AXIS))
 ENDDO
 UBAR = VEL_G_INT(IAXIS)
 VBAR = VEL_G_INT(JAXIS)