@@ -99,7 +99,7 @@ END SUBROUTINE UPDATE_GLOBAL_OUTPUTS
9999SUBROUTINE DUMP_MESH_OUTPUTS (T ,DT ,NM ,FAKEWRITE )
100100
101101USE COMP_FUNCTIONS, ONLY : CURRENT_TIME
102- USE TURBULENCE, ONLY: SANDIA_OUT, SPECTRAL_OUTPUT
102+ USE TURBULENCE, ONLY: SANDIA_OUT
103103REAL (EB) :: TNOW
104104REAL (EB), INTENT (IN ) :: T,DT
105105INTEGER , INTENT (IN ) :: NM
@@ -649,9 +649,7 @@ END SUBROUTINE ASSIGN_FILE_NAMES
649649SUBROUTINE INITIALIZE_GLOBAL_DUMPS (T ,DT )
650650
651651USE COMP_FUNCTIONS, ONLY: CURRENT_TIME,GET_FILE_NUMBER,APPEND_FILE
652- USE PHYSICAL_FUNCTIONS, ONLY: GET_MASS_FRACTION_ALL,GET_MOLECULAR_WEIGHT
653652USE HVAC_ROUTINES, ONLY: N_DUCT_QUANTITY,N_NODE_QUANTITY
654- USE CONTROL_VARIABLES
655653REAL (EB) :: TNOW
656654REAL (EB), INTENT (IN ) :: T,DT
657655INTEGER :: NN,I,N,N_OUT,N_ZONE_TMP,LU,J, N_NODE_OUT, N_DUCT_OUT, NS
@@ -4213,10 +4211,10 @@ SUBROUTINE WRITE_DIAGNOSTICS(T,DT)
42134211
42144212USE COMP_FUNCTIONS, ONLY : CURRENT_TIME,GET_DATE,GET_DATE_ISO_8601
42154213REAL (EB), INTENT (IN ) :: T,DT
4216- INTEGER :: NM,II,JJ,KK,OUT_DIGITS,SOUT_DIGITS
4214+ INTEGER :: NM,II,JJ,KK,OUT_DIGITS,SOUT_DIGITS,MAX_VN_IJK( 3 ),MAX_CFL_IJK( 3 ),MAX_CFL_MESH,MAX_VN_MESH
42174215CHARACTER (120 ) :: SIMPLE_OUTPUT,SIMPLE_OUTPUT_ERR,OUT_FORMAT
42184216CHARACTER (LABEL_LENGTH) :: DATE
4219- REAL (EB) :: TNOW,CPUTIME,STIME,DTS
4217+ REAL (EB) :: TNOW,CPUTIME,STIME,DTS,MAX_CFL,MAX_VN
42204218
42214219TNOW = CURRENT_TIME()
42224220
@@ -4274,6 +4272,24 @@ SUBROUTINE WRITE_DIAGNOSTICS(T,DT)
42744272
42754273WRITE (LU_ERR,' (A)' ) TRIM (SIMPLE_OUTPUT_ERR)
42764274
4275+ ! Determine the mesh where the maximum CFL, VN, etc, occur
4276+
4277+ MAX_CFL = - 1._EB
4278+ MAX_VN = - 1._EB
4279+ DO NM= 1 ,NMESHES
4280+ M = > MESHES(NM)
4281+ IF (M% CFL> MAX_CFL) THEN
4282+ MAX_CFL = MAX (M% CFL,MAX_CFL)
4283+ MAX_CFL_MESH = NM
4284+ MAX_CFL_IJK = (/ M% ICFL,M% JCFL,M% KCFL/ )
4285+ ENDIF
4286+ IF (CHECK_VN .AND. M% VN> MAX_VN) THEN
4287+ MAX_VN = MAX (M% VN,MAX_VN)
4288+ MAX_VN_MESH = NM
4289+ MAX_VN_IJK = (/ M% I_VN,M% J_VN,M% K_VN/ )
4290+ ENDIF
4291+ ENDDO
4292+
42774293! Header for .out file
42784294
42794295IF (ICYC== 1 ) WRITE (LU_OUTPUT,100 )
@@ -4326,6 +4342,13 @@ SUBROUTINE WRITE_DIAGNOSTICS(T,DT)
43264342 ' on Mesh ' ,NM,' at (' ,II,' ,' ,JJ,' ,' ,KK,' )'
43274343ENDIF
43284344
4345+ WRITE (LU_OUTPUT,' (7X,A,E9.2,A,4(I0,A))' ) ' Maximum CFL Number : ' ,MAX_CFL,' on Mesh ' ,MAX_CFL_MESH,&
4346+ ' at (' ,MAX_CFL_IJK(1 ),' ,' ,MAX_CFL_IJK(2 ),' ,' ,MAX_CFL_IJK(3 ),' )'
4347+ IF (CHECK_VN) THEN
4348+ WRITE (LU_OUTPUT,' (7X,A,E9.2,A,4(I0,A))' ) ' Maximum VN Number : ' ,MAX_VN,' on Mesh ' ,MAX_VN_MESH,&
4349+ ' at (' ,MAX_VN_IJK(1 ),' ,' ,MAX_VN_IJK(2 ),' ,' ,MAX_VN_IJK(3 ),' )'
4350+ ENDIF
4351+
43294352WRITE (LU_OUTPUT,' (7X,A)' ) ' ---------------------------------------------------------------'
43304353
43314354DO NM= 1 ,NMESHES
@@ -4341,7 +4364,6 @@ SUBROUTINE WRITE_DIAGNOSTICS(T,DT)
43414364 WRITE (LU_OUTPUT,121 ) M% DT_RESTRICT_STORE
43424365 M% DT_RESTRICT_STORE= 0
43434366 ENDIF
4344-
43454367ENDDO
43464368
43474369WRITE (LU_OUTPUT,* )
@@ -8464,7 +8486,7 @@ REAL(EB) RECURSIVE FUNCTION GAS_PHASE_OUTPUT(T,DT,NM,II,JJ,KK,IND,IND2,Y_INDEX,Z
84648486 PROP_INDEX,REAC_INDEX,MATL_INDEX,ICC_IN,JCC_IN) RESULT(GAS_PHASE_OUTPUT_RES)
84658487
84668488USE MEMORY_FUNCTIONS, ONLY: REALLOCATE
8467- USE MATH_FUNCTIONS, ONLY: INTERPOLATE1D,INTERPOLATE1D_UNIFORM,EVALUATE_RAMP, UPDATE_HISTOGRAM
8489+ USE MATH_FUNCTIONS, ONLY: INTERPOLATE1D,INTERPOLATE1D_UNIFORM,UPDATE_HISTOGRAM
84688490USE PHYSICAL_FUNCTIONS, ONLY: GET_MASS_FRACTION,FED,FIC,GET_SPECIFIC_HEAT,RELATIVE_HUMIDITY, &
84698491 GET_CONDUCTIVITY,GET_MOLECULAR_WEIGHT,GET_MASS_FRACTION_ALL,GET_ENTHALPY,GET_SENSIBLE_ENTHALPY, &
84708492 GET_VISCOSITY,GET_POTENTIAL_TEMPERATURE,GET_SPECIFIC_GAS_CONSTANT,&
0 commit comments