Merge pull request #15026 from drjfloyd/master

FDS Source: Edits to BI-DIRECTIONAL PROBE output.

Author: drjfloyd

Source/dump.f90

@@ -7897,16 +7897,15 @@ REAL(EB) RECURSIVE FUNCTION GAS_PHASE_OUTPUT(T,DT,NM,II,JJ,KK,IND,IND2,Y_INDEX,Z
       COSTHETA = (UU*ORIENTATION_VECTOR(1,DV%ORIENTATION_INDEX)+VV*ORIENTATION_VECTOR(2,DV%ORIENTATION_INDEX)+&
                   WW*ORIENTATION_VECTOR(3,DV%ORIENTATION_INDEX))/SQRT(VEL2)
       FAC = -2.308_EB*ABS(COSTHETA)**3 + 2.533_EB*ABS(COSTHETA)**2 + 0.7847_EB*ABS(COSTHETA)
-      DP = DP * FAC
-      VEL = SQRT(2._EB*DP/RHO(II,JJ,KK))
+      VEL = FAC*SQRT(VEL2)
       ZZ_GET(1:N_TRACKED_SPECIES) = ZZ(II,JJ,KK,1:N_TRACKED_SPECIES)
       CALL GET_VISCOSITY(ZZ_GET,MU_G,TMP(II,JJ,KK))
       RE_D = MIN(3800._EB,MAX(40._EB,RHO(II,JJ,KK)*VEL*PY%PROBE_DIAMETER/MU_G))
       FAC = 1.533_EB-0.001366_EB*RE_D+0.000001688_EB*RE_D**2-0.0000000009706_EB*RE_D**3+&
              0.0000000000002555_EB*RE_D**4-2.484E-17_EB*RE_D**5
       RHO_AIR = 350.9736_EB/PROBE_TMP !350 is 0.0288 101325/ 8.314472
-      GAS_PHASE_OUTPUT_RES = SIGN(1._EB,COSTHETA)*PY%CALIBRATION_CONSTANT*0.5_EB*RHO_AIR*(VEL*FAC)**2
-            
+      GAS_PHASE_OUTPUT_RES = SIGN(1._EB,COSTHETA)*VEL*PY%CALIBRATION_CONSTANT*FAC
+
    CASE(130) ! EXTINCTION
       ZZ_GET(1:N_TRACKED_SPECIES) = ZZ(II,JJ,KK,1:N_TRACKED_SPECIES)
       ZZ_FUEL = 0._EB