FDS Source: allow locking time step during FREEZE_TIME

Author: Hodges

Source/main.f90

@@ -544,6 +544,7 @@ PROGRAM FDS
 IF (UNFREEZE_TIME > 0._EB) THEN 
    FREEZE_VELOCITY=.TRUE. 
    SOLID_PHASE_ONLY=.TRUE.
+   LOCK_TIME_STEP=.TRUE.
 ENDIF
 
 IF (MY_RANK==0 .AND. VERBOSE) CALL VERBOSE_PRINTOUT('Starting the time-stepping')
@@ -586,6 +587,16 @@ PROGRAM FDS
    IF ((UNFREEZE_TIME > 0._EB).AND.(T>UNFREEZE_TIME)) THEN 
       FREEZE_VELOCITY=.FALSE.
       SOLID_PHASE_ONLY=.FALSE.
+      LOCK_TIME_STEP=.FALSE.
+      
+      UNFREEZE_DT_LOOP: DO NM=LOWER_MESH_INDEX,UPPER_MESH_INDEX
+         CALL CHECK_STABILITY(DT,DT_NEW,T,NM)
+      ENDDO UNFREEZE_DT_LOOP
+      IF (N_MPI_PROCESSES>1) THEN
+         CALL MPI_ALLGATHERV(MPI_IN_PLACE,0,MPI_DATATYPE_NULL,DT_NEW(1:NMESHES),&
+                             COUNTS,DISPLS,MPI_DOUBLE_PRECISION,MPI_COMM_WORLD,IERR)
+      ENDIF
+      DT = MINVAL(DT_NEW)
    ENDIF
 
    !================================================================================================================================

Source/velo.f90

@@ -12,7 +12,8 @@ MODULE VELO
 PRIVATE
 
 PUBLIC VELOCITY_PREDICTOR,VELOCITY_CORRECTOR,NO_FLUX,BAROCLINIC_CORRECTION,MATCH_VELOCITY,MATCH_VELOCITY_FLUX,&
-       VELOCITY_BC,COMPUTE_VISCOSITY,VISCOSITY_BC,VELOCITY_FLUX,VELOCITY_FLUX_CYLINDRICAL
+       VELOCITY_BC,COMPUTE_VISCOSITY,VISCOSITY_BC,VELOCITY_FLUX,VELOCITY_FLUX_CYLINDRICAL,&
+       CHECK_STABILITY
 
 
 CONTAINS