FDS Source: Iterate BI-DIR for dynamic calibration

Author: drjfloyd

Source/dump.f90

@@ -7926,10 +7926,16 @@ REAL(EB) RECURSIVE FUNCTION GAS_PHASE_OUTPUT(T,DT,NM,II,JJ,KK,IND,IND2,Y_INDEX,Z
          ! LJ AIR viscosity fit 100 K to 5000 K
          MU_G = 1.5205E-22_EB*PROBE_TMP**5 - 2.1417E-18_EB*PROBE_TMP**4 + 1.1402E-14_EB*PROBE_TMP**3 - &
                   2.9846E-11_EB*PROBE_TMP**2 + 5.9898E-8_EB*PROBE_TMP + 0.000002352_EB
-         RE_D = MIN(3800._EB,MAX(40._EB,PROBE_RHO*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
-         GAS_PHASE_OUTPUT_RES = SIGN(1._EB,COSTHETA)*1._EB/FAC*SQRT(2*PROBE_DELTA_P/PROBE_RHO)
+         JJJ = 1
+         BP_LOOP: DO
+            RE_D = MIN(3800._EB,MAX(40._EB,PROBE_RHO*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
+            GAS_PHASE_OUTPUT_RES = SIGN(1._EB,COSTHETA)*1._EB/FAC*SQRT(2._EB*PROBE_DELTA_P/PROBE_RHO)
+            IF (JJJ > 9 .OR. ABS(VEL-GAS_PHASE_OUTPUT_RES)/(GAS_PHASE_OUTPUT_RES+TWO_EPSILON_EB) < 0.001_EB) EXIT BP_LOOP
+            VEL = 0.2_EB*VEL+0.8_EB*GAS_PHASE_OUTPUT_RES
+            JJJ = JJJ + 1
+         ENDDO BP_LOOP
       ENDIF
 
    CASE(130) ! EXTINCTION