Skip to content

Commit 7471836

Browse files
drjfloyddrjfloyd
authored
Merge pull request #15024 from drjfloyd/master
FDS Source: Add BI-DIRECTIONAL PROBE output
2 parents 43eafbd + 93c3442 commit 7471836

File tree

5 files changed

+63
-12
lines changed

5 files changed

+63
-12
lines changed

‎Manuals/FDS_User_Guide/FDS_User_Guide.tex‎

Lines changed: 13 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -10558,7 +10558,15 @@ \subsection{Thermocouples}
1055810558
\label{fig:thermocouple_time_constant}
1055910559
\end{figure}
1056010560

10561+
\subsection{Bi-Directional Probe}
10562+
\label{info:bidir_probe}
1056110563

10564+
The output quantity \ct{BI-DIRECTIONAL PROBE} is the velocity of a modeled bi-directional probe. A bi-directional probe uses the following equation:
10565+
\be
10566+
C \sqrt{\frac{2 \Delta P}{\rho}}
10567+
\label{BDP}
10568+
\ee
10569+
where $C$ is a calibration constant (default value is 0.93), $\Delta P$ is the pressure difference across the probe, and $\rho$ is the gas density at the probe. In a typical experiment, the gas density is computed assuming standard pressure (101325 Pa), the molecular weight of air (28.8 g/mol), and the temperature as measured by a thermocouple near the probe. Bi-directional probes have biases due to both the Reynolds number (based on the probe diameter) of the flow and the angle of the flow with respect to the probe axis. This model accounts for those sensitivities and the impact of density differences from varied molecular weight at the probe. The orientation of the probe can be specified with either \ct{IOR} or \ct{ORIENTATION} on \ct{DEVC}. A probe with \ct{IOR}=-1 would have a positive velocity output when the flow is in the negative x direction. Parameters for the probe can be specified with a \ct{PROP_ID} on the \ct {DEVC}. The calibration constant and the probe diameter (default of 0.0254 m) can be set respectively with \ct{CALIBRATION_CONSTANT} and \ct{PROBE_DIAMETER} on \ct{PROP}. If the probe temperature is an aspirated thermocouple or other measurement not sensitive to the radiative environment, then set \ct{TC=F} on \ct{PROP}.
1056210570

1056310571
\subsection{Volume Flow}
1056410572
\label{info:volume_flow}
@@ -11726,6 +11734,7 @@ \section{Device, Control, and Other Miscellaneous Output Quantities}
1172611734
\endhead
1172711735
\ct{ACTUATED SPRINKLERS} & Section~\ref{info:TIMING} & & D \\ \hline
1172811736
\ct{ASPIRATION} & Section~\ref{info:aspiration_detector} & \%/m & D \\ \hline
11737+
\ct{BI-DIRECTIONAL PROBE} & Section~\ref{info:bidir_probe} & m/s & D \\ \hline
1172911738
\ct{CHAMBER OBSCURATION} & Section~\ref{info:smoke_detector} & \%/m & D \\ \hline
1173011739
\ct{CELL INDEX I} & Mesh cell index in x & & D,S \\ \hline
1173111740
\ct{CELL INDEX J} & Mesh cell index in y & & D,S \\ \hline
@@ -13171,6 +13180,7 @@ \section{\texorpdfstring{{\tt PROP}}{PROP} (Device Properties)}
1317113180
\ct{ALPHA_E} & Real & Section~\ref{info:smoke_detector} & & 0. \\ \hline
1317213181
\ct{BETA_C} & Real & Section~\ref{info:smoke_detector} & & 1. \\ \hline
1317313182
\ct{BETA_E} & Real & Section~\ref{info:smoke_detector} & & 1. \\ \hline
13183+
\ct{CALIBRATION_CONSTANT} & Real & Section~\ref{info:bidir_probe} & & 0.93 \\ \hline
1317413184
\ct{CHARACTERISTIC_VELOCITY} & Real & Section~\ref{info:pressure_coefficient} & m/s & 1. \\ \hline
1317513185
\ct{C_FACTOR} & Real & Section~\ref{info:sprinklers} & (m/s)$^{1/2}$ & 0. \\ \hline
1317613186
\ct{DENSITY} & Real & Section~\ref{info:THERMOCOUPLE} & kg/m$^3$ & 8908. \\ \hline
@@ -13207,6 +13217,7 @@ \section{\texorpdfstring{{\tt PROP}}{PROP} (Device Properties)}
1320713217
\ct{PDPA_RADIUS} & Real & Section~\ref{PDPA} & m & 0. \\ \hline
1320813218
\ct{PDPA_START} & Real & Section~\ref{PDPA} & s & 0. \\ \hline
1320913219
\ct{PRESSURE_RAMP} & Character & Section~\ref{info:sprinklers} & & \\ \hline
13220+
\ct{PROBE_DIAMETER} & Real & Section~\ref{info:bidir_probe} & m & 0.0254 \\ \hline
1321013221
\ct{P0} & Real & Section~\ref{info:velocity_patch} & m/s & 0. \\ \hline
1321113222
\ct{PX(3)} & Real & Section~\ref{info:velocity_patch} & m/s & 0. \\ \hline
1321213223
\ct{PXX(3,3)} & Real & Section~\ref{info:velocity_patch} & m/s & 0. \\ \hline
@@ -13222,6 +13233,7 @@ \section{\texorpdfstring{{\tt PROP}}{PROP} (Device Properties)}
1322213233
\ct{SPRAY_PATTERN_MU} & Real & Section~\ref{info:sprinklers} & degrees & 0. \\ \hline
1322313234
\ct{SPRAY_PATTERN_SHAPE} & Character & Section~\ref{info:sprinklers} & & \ct{'GAUSSIAN'} \\ \hline
1322413235
\ct{SPRAY_PATTERN_TABLE} & Character & Section~\ref{info:sprinklers} & & \\ \hline
13236+
\ct{TC} & Logical & Section~\ref{info:bidir_probe} & & \ct{T} \\ \hline
1322513237
\ct{TIME_CONSTANT} & Real & Section~\ref{info:THERMOCOUPLE} & s & \\ \hline
1322613238
\ct{VELOCITY_COMPONENT} & Integer & Section~\ref{info:velocity_patch} & & \\ \hline
1322713239
\ct{VIEW_ANGLE} & Real & Section~\ref{info:heat_flux} & degrees & 180. \\ \hline
@@ -14444,6 +14456,7 @@ \chapter{Error Codes}
1444414456
885 \> \ct{DEVC ... must have coordinates, even if not a point quantity.} \> Section~\ref{info:DEVC} \\
1444514457
886 \> \ct{DEVC ... must have either an output QUANTITY or PROP_ID.} \> Section~\ref{info:DEVC} \\
1444614458
887 \> \ct{DEVC ... must have an ORIENTATION.} \> Section~\ref{info:heat_flux} \\
14459+
\> \> or Section~\ref{info:bidir_probe} \\
1444714460
888 \> \ct{DEVC ... TEMPORAL_STATISTIC is not recognized.} \> Section~\ref{info:statistics} \\
1444814461
889 \> \ct{DEVC ... SPATIAL_STATISTIC is not recognized.} \> Section~\ref{info:statistics} \\
1444914462
890 \> \ct{DEVC ... STATISTICS is not recognized.} \> Section~\ref{info:statistics} \\

‎Source/data.f90‎

Lines changed: 6 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -608,6 +608,12 @@ SUBROUTINE DEFINE_OUTPUT_QUANTITIES
608608
OUTPUT_QUANTITY(113)%UNITS = 'kW/m2'
609609
OUTPUT_QUANTITY(113)%SHORT_NAME = 'hflux'
610610

611+
! Model of a Bi-Directional Probe
612+
613+
OUTPUT_QUANTITY(114)%NAME = 'BI-DIRECTIONAL PROBE'
614+
OUTPUT_QUANTITY(114)%UNITS = 'm/s'
615+
OUTPUT_QUANTITY(114)%SHORT_NAME = 'bdp'
616+
611617
! Special Outputs for Partially Stirred Batch Reactor Model
612618

613619
OUTPUT_QUANTITY(130)%NAME = 'EXTINCTION'

‎Source/devc.f90‎

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -11,13 +11,13 @@ MODULE DEVICE_VARIABLES
1111

1212
TYPE PROPERTY_TYPE
1313
REAL(EB) :: DENSITY,DIAMETER,EMISSIVITY,HEAT_TRANSFER_COEFFICIENT,SPECIFIC_HEAT,RTI,TIME_CONSTANT, &
14-
ACTIVATION_TEMPERATURE,ACTIVATION_OBSCURATION, &
14+
ACTIVATION_TEMPERATURE,ACTIVATION_OBSCURATION, CALIBRATION_CONSTANT,&
1515
ALPHA_E,ALPHA_C,BETA_E,BETA_C,CHARACTERISTIC_VELOCITY,PARTICLE_VELOCITY,MASS_FLOW_RATE,FLOW_RATE,FLOW_TAU, &
1616
GAUGE_EMISSIVITY,GAUGE_TEMPERATURE,INITIAL_TEMPERATURE,K_FACTOR,C_FACTOR,OPERATING_PRESSURE,OFFSET,&
17-
SPRAY_ANGLE(2,2),P0=0._EB,PX(3)=0._EB,PXX(3,3)=0._EB,VIEW_ANGLE
17+
SPRAY_ANGLE(2,2),P0=0._EB,PX(3)=0._EB,PXX(3,3)=0._EB,VIEW_ANGLE,PROBE_DIAMETER
1818
INTEGER :: PDPA_M=0,PDPA_N=0,N_SMOKEVIEW_PARAMETERS=0,N_SMOKEVIEW_IDS=0,N_INSERT,I_VEL=0,PARTICLES_PER_SECOND
1919
LOGICAL :: PDPA_INTEGRATE=.TRUE.,PDPA_NORMALIZE=.TRUE.,HISTOGRAM_NORMALIZE=.TRUE.,HISTOGRAM=.FALSE., &
20-
HISTOGRAM_CUMULATIVE=.FALSE.,SPARK=.FALSE.
20+
HISTOGRAM_CUMULATIVE=.FALSE.,SPARK=.FALSE.,TC=.TRUE.
2121
REAL(EB) :: PDPA_START=0._EB,PDPA_END=1.E6_EB,PDPA_RADIUS=0.1_EB
2222
REAL(EB), ALLOCATABLE, DIMENSION(:) :: TABLE_ROW, V_FACTOR
2323
INTEGER :: PART_INDEX=-1,FLOW_RAMP_INDEX,SPRAY_PATTERN_INDEX,Z_INDEX=-999,Y_INDEX=-999,PRESSURE_RAMP_INDEX

‎Source/dump.f90‎

Lines changed: 27 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -7328,13 +7328,13 @@ REAL(EB) RECURSIVE FUNCTION GAS_PHASE_OUTPUT(T,DT,NM,II,JJ,KK,IND,IND2,Y_INDEX,Z
73287328
INTEGER, INTENT(IN) :: II,JJ,KK,IND,IND2,NM,Y_INDEX,Z_INDEX,ELEM_INDX,PART_INDEX,VELO_INDEX,PIPE_INDEX,PROP_INDEX,REAC_INDEX, &
73297329
MATL_INDEX
73307330
INTEGER, INTENT(IN), OPTIONAL :: ICC_IN,JCC_IN
7331-
REAL(EB) :: H_TC,TMP_TC,RE_D,NUSSELT,VEL,K_G,MU_G,&
7331+
REAL(EB) :: H_TC,TMP_TC,RE_D,NUSSELT,VEL,K_G,MU_G,COSTHETA,FAC,DP,&
73327332
Q_SUM,TMP_G,UU,VV,WW,VEL2,Y_MF_INT,PATHLENGTH,EXT_COEF,MASS_EXT_COEF,ZZ_FUEL,ZZ_OX,&
73337333
VELSR,WATER_VOL_FRAC,RHS,DT_C,DT_E,T_RATIO,Y_E_LAG, H_G,H_G_SUM,CPBAR,CP,ZZ_GET(1:N_TRACKED_SPECIES),RCON,&
73347334
EXPON,Y_SPECIES,MEC,Y_SPECIES2,Y_H2O,R_Y_H2O,R_DN,SGN,Y_ALL(N_SPECIES),H_S,D_Z_N(0:I_MAX_TEMP),&
73357335
DISSIPATION_RATE,S11,S22,S33,S12,S13,S23,DUDX,DUDY,DUDZ,DVDX,DVDY,DVDZ,DWDX,DWDY,DWDZ,ONTHDIV,SS,ETA,DELTA,R_DX2,&
73367336
UVW,UODX,VODY,WODZ,XHAT,ZHAT,BBF,GAMMA_LOC,VC,VOL,PHI,GAS_PHASE_OUTPUT_CC,&
7337-
GAS_PHASE_OUTPUT_CFA,CFACE_AREA,VELOCITY_COMPONENT(1:3),ATOTV(1:3),TMP_F,R_D,MW
7337+
GAS_PHASE_OUTPUT_CFA,CFACE_AREA,VELOCITY_COMPONENT(1:3),ATOTV(1:3),TMP_F,R_D,MW,RHO_AIR,PROBE_TMP
73387338
INTEGER :: N,I,J,K,NN,IL,III,JJJ,KKK,IP,JP,KP,FED_ACTIVITY,IP1,JP1,KP1,IM1,JM1,KM1,IIM1,JJM1,KKM1,NR,NS,RAM,&
73397339
ICC,JCC,NCELL,AXIS,ICF,NFACE,JCF,JCC_LO,JCC_HI,PDPA_FORMULA,IC
73407340
REAL(FB) :: RN
@@ -7882,6 +7882,31 @@ REAL(EB) RECURSIVE FUNCTION GAS_PHASE_OUTPUT(T,DT,NM,II,JJ,KK,IND,IND2,Y_INDEX,Z
78827882
ENDIF
78837883
GAS_PHASE_OUTPUT_RES = (GAS_PHASE_OUTPUT_RES - K_G*(TMP(IP,JP,KP)-TMP(II,JJ,KK))*R_DN)*0.001
78847884

7885+
CASE(114) ! BI-DIRECTIONAL PROBE
7886+
IF (PY%TC) THEN
7887+
PROBE_TMP = GAS_PHASE_OUTPUT(T,DT,NM,II,JJ,KK,110,IND2,Y_INDEX,Z_INDEX,ELEM_INDX,PART_INDEX,VELO_INDEX,PIPE_INDEX,&
7888+
PROP_INDEX,REAC_INDEX,MATL_INDEX,ICC_IN,JCC_IN) + TMPM
7889+
ELSE
7890+
PROBE_TMP = TMP(II,JJ,KK)
7891+
ENDIF
7892+
UU = U(II,JJ,KK)
7893+
VV = V(II,JJ,KK)
7894+
WW = W(II,JJ,KK)
7895+
VEL2 = UU**2+VV**2+WW**2
7896+
DP = 0.5_EB*VEL2*RHO(II,JJ,KK)
7897+
COSTHETA = (UU*ORIENTATION_VECTOR(1,DV%ORIENTATION_INDEX)+VV*ORIENTATION_VECTOR(2,DV%ORIENTATION_INDEX)+&
7898+
WW*ORIENTATION_VECTOR(3,DV%ORIENTATION_INDEX))/SQRT(VEL2)
7899+
FAC = -2.308_EB*ABS(COSTHETA)**3 + 2.533_EB*ABS(COSTHETA)**2 + 0.7847_EB*ABS(COSTHETA)
7900+
DP = DP * FAC
7901+
VEL = SQRT(2._EB*DP/RHO(II,JJ,KK))
7902+
ZZ_GET(1:N_TRACKED_SPECIES) = ZZ(II,JJ,KK,1:N_TRACKED_SPECIES)
7903+
CALL GET_VISCOSITY(ZZ_GET,MU_G,TMP(II,JJ,KK))
7904+
RE_D = MIN(3800._EB,MAX(40._EB,RHO(II,JJ,KK)*VEL*PY%PROBE_DIAMETER/MU_G))
7905+
FAC = 1.533_EB-0.001366_EB*RE_D+0.000001688_EB*RE_D**2-0.0000000009706_EB*RE_D**3+&
7906+
0.0000000000002555_EB*RE_D**4-2.484E-17_EB*RE_D**5
7907+
RHO_AIR = 350.9736_EB/PROBE_TMP !350 is 0.0288 101325/ 8.314472
7908+
GAS_PHASE_OUTPUT_RES = SIGN(1._EB,COSTHETA)*PY%CALIBRATION_CONSTANT*0.5_EB*RHO_AIR*(VEL*FAC)**2
7909+
78857910
CASE(130) ! EXTINCTION
78867911
ZZ_GET(1:N_TRACKED_SPECIES) = ZZ(II,JJ,KK,1:N_TRACKED_SPECIES)
78877912
ZZ_FUEL = 0._EB

‎Source/read.f90‎

Lines changed: 14 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -6481,30 +6481,30 @@ SUBROUTINE READ_PROP
64816481
USE DEVICE_VARIABLES
64826482
USE PHYSICAL_FUNCTIONS, ONLY : SPRAY_ANGLE_DISTRIBUTION
64836483
USE MATH_FUNCTIONS, ONLY: GET_TABLE_INDEX
6484-
REAL(EB) :: ACTIVATION_OBSCURATION,ACTIVATION_TEMPERATURE,ALPHA_C,ALPHA_E,BETA_C,BETA_E, &
6484+
REAL(EB) :: ACTIVATION_OBSCURATION,ACTIVATION_TEMPERATURE,ALPHA_C,ALPHA_E,BETA_C,BETA_E,CALIBRATION_CONSTANT, &
64856485
HEAT_TRANSFER_COEFFICIENT,DIAMETER,DENSITY,SPECIFIC_HEAT, &
64866486
C_FACTOR,CHARACTERISTIC_VELOCITY,ORIFICE_DIAMETER,EMISSIVITY, &
64876487
PARTICLE_VELOCITY,FLOW_RATE,FLOW_TAU,GAUGE_EMISSIVITY,GAUGE_TEMPERATURE,INITIAL_TEMPERATURE,K_FACTOR,&
64886488
LENGTH,SPRAY_ANGLE(2,2),OFFSET,OPERATING_PRESSURE,RTI,PDPA_START,PDPA_END,PDPA_RADIUS,MASS_FLOW_RATE,&
6489-
SPRAY_PATTERN_MU,SPRAY_PATTERN_BETA,HISTOGRAM_LIMITS(2),P0,PX(3),PXX(3,3),TIME_CONSTANT,VIEW_ANGLE
6489+
SPRAY_PATTERN_MU,SPRAY_PATTERN_BETA,HISTOGRAM_LIMITS(2),P0,PX(3),PXX(3,3),TIME_CONSTANT,VIEW_ANGLE,PROBE_DIAMETER
64906490
INTEGER ::I,N,NN,PDPA_M,PDPA_N,PARTICLES_PER_SECOND,VELOCITY_COMPONENT,HISTOGRAM_NBINS,FED_ACTIVITY
6491-
LOGICAL :: PDPA_INTEGRATE,PDPA_NORMALIZE,HISTOGRAM_NORMALIZE,HISTOGRAM,HISTOGRAM_CUMULATIVE,SPARK
6491+
LOGICAL :: PDPA_INTEGRATE,PDPA_NORMALIZE,HISTOGRAM_NORMALIZE,HISTOGRAM,HISTOGRAM_CUMULATIVE,SPARK,TC
64926492
CHARACTER(LABEL_LENGTH) :: SMOKEVIEW_ID(SMOKEVIEW_OBJECTS_DIMENSION),QUANTITY='null',PART_ID='null',FLOW_RAMP='null', &
64936493
SPRAY_PATTERN_TABLE='null',SPEC_ID='null',&
64946494
PRESSURE_RAMP='null',SMOKEVIEW_PARAMETERS(SMOKEVIEW_OBJECTS_DIMENSION), &
64956495
SPRAY_PATTERN_SHAPE='GAUSSIAN'
64966496
TYPE (PROPERTY_TYPE), POINTER :: PY
64976497

6498-
NAMELIST /PROP/ ACTIVATION_OBSCURATION,ACTIVATION_TEMPERATURE,ALPHA_C,ALPHA_E,BETA_C,BETA_E, &
6498+
NAMELIST /PROP/ ACTIVATION_OBSCURATION,ACTIVATION_TEMPERATURE,ALPHA_C,ALPHA_E,BETA_C,BETA_E,CALIBRATION_CONSTANT,&
64996499
CHARACTERISTIC_VELOCITY,C_FACTOR,DENSITY,DIAMETER,EMISSIVITY,FED_ACTIVITY,FLOW_RAMP,FLOW_RATE,FLOW_TAU, &
65006500
GAUGE_EMISSIVITY,GAUGE_TEMPERATURE,HEAT_TRANSFER_COEFFICIENT,HISTOGRAM,HISTOGRAM_CUMULATIVE, &
65016501
HISTOGRAM_LIMITS,HISTOGRAM_NBINS,HISTOGRAM_NORMALIZE,ID, &
65026502
INITIAL_TEMPERATURE,K_FACTOR,LENGTH,MASS_FLOW_RATE,OFFSET,OPERATING_PRESSURE,ORIFICE_DIAMETER,P0,&
65036503
PARTICLES_PER_SECOND,PARTICLE_VELOCITY,PART_ID,PDPA_END,&
65046504
PDPA_INTEGRATE,PDPA_M,PDPA_N,PDPA_NORMALIZE,PDPA_RADIUS,&
6505-
PDPA_START,PRESSURE_RAMP,PX,PXX,QUANTITY,RTI,SMOKEVIEW_ID,SMOKEVIEW_PARAMETERS,SPARK,&
6505+
PDPA_START,PRESSURE_RAMP,PROBE_DIAMETER,PX,PXX,QUANTITY,RTI,SMOKEVIEW_ID,SMOKEVIEW_PARAMETERS,SPARK,&
65066506
SPEC_ID,SPECIFIC_HEAT,SPRAY_ANGLE,&
6507-
SPRAY_PATTERN_BETA,SPRAY_PATTERN_MU,SPRAY_PATTERN_SHAPE,SPRAY_PATTERN_TABLE,TIME_CONSTANT,VELOCITY_COMPONENT,&
6507+
SPRAY_PATTERN_BETA,SPRAY_PATTERN_MU,SPRAY_PATTERN_SHAPE,SPRAY_PATTERN_TABLE,TC,TIME_CONSTANT,VELOCITY_COMPONENT,&
65086508
VIEW_ANGLE
65096509

65106510
! Count the PROP lines in the input file. Note how many of these are cables.
@@ -6547,6 +6547,7 @@ SUBROUTINE READ_PROP
65476547
PY%ALPHA_E = ALPHA_E
65486548
PY%BETA_C = BETA_C
65496549
PY%BETA_E = BETA_E
6550+
PY%CALIBRATION_CONSTANT = CALIBRATION_CONSTANT
65506551
PY%DENSITY = DENSITY
65516552
PY%DIAMETER = DIAMETER
65526553
PY%EMISSIVITY = EMISSIVITY
@@ -6562,6 +6563,7 @@ SUBROUTINE READ_PROP
65626563
PY%OFFSET = OFFSET
65636564
PY%OPERATING_PRESSURE = OPERATING_PRESSURE
65646565
PY%PART_ID = PART_ID
6566+
PY%PROBE_DIAMETER = PROBE_DIAMETER
65656567
PY%QUANTITY = QUANTITY
65666568
IF (PY%PART_ID/='null' .AND. PY%QUANTITY == 'null' ) PY%QUANTITY = 'NOZZLE FLOW RATE'
65676569
PY%RTI = RTI
@@ -6620,6 +6622,7 @@ SUBROUTINE READ_PROP
66206622
CALL SPRAY_ANGLE_DISTRIBUTION(PY%SPRAY_LON,PY%SPRAY_LAT,PY%SPRAY_LON_CDF,PY%SPRAY_LAT_CDF, &
66216623
SPRAY_PATTERN_BETA,SPRAY_PATTERN_MU,PY%SPRAY_ANGLE,SPRAY_PATTERN_SHAPE,NDC2)
66226624
ENDIF
6625+
PY%TC = TC
66236626
PY%TIME_CONSTANT = TIME_CONSTANT
66246627

66256628
! PDPA model
@@ -6789,6 +6792,7 @@ SUBROUTINE SET_PROP_DEFAULTS
67896792
ALPHA_E = 0.0_EB
67906793
BETA_C = -1.0_EB
67916794
BETA_E = -1.0_EB
6795+
CALIBRATION_CONSTANT = 0.93_EB
67926796
DENSITY = 8908._EB ! kg/m3 (Nickel)
67936797
DIAMETER = 0.001 ! m
67946798
EMISSIVITY = 0.85_EB
@@ -6828,6 +6832,7 @@ SUBROUTINE SET_PROP_DEFAULTS
68286832
HISTOGRAM_NORMALIZE = .TRUE.
68296833
PRESSURE_RAMP = 'null'
68306834
P0 = -1.E10_EB
6835+
PROBE_DIAMETER = 0.0254_EB
68316836
PX = 0._EB
68326837
PXX = 0._EB
68336838
QUANTITY = 'null'
@@ -6844,6 +6849,7 @@ SUBROUTINE SET_PROP_DEFAULTS
68446849
SPRAY_PATTERN_SHAPE = 'GAUSSIAN'
68456850
SPRAY_PATTERN_MU = -1._EB
68466851
SPRAY_PATTERN_BETA = 5.0_EB
6852+
TC = .TRUE.
68476853
TIME_CONSTANT = -1._EB
68486854
FED_ACTIVITY = 2 ! light work
68496855
VELOCITY_COMPONENT = 0
@@ -14096,7 +14102,8 @@ SUBROUTINE READ_DEVC
1409614102
DV%QUANTITY(1)=='GAUGE HEAT FLUX GAS' .OR. &
1409714103
DV%QUANTITY(1)=='RADIANCE' .OR. &
1409814104
DV%QUANTITY(1)=='ADIABATIC SURFACE TEMPERATURE GAS' .OR. &
14099-
DV%QUANTITY(1)=='RADIOMETER GAS') THEN
14105+
DV%QUANTITY(1)=='RADIOMETER GAS' .OR. &
14106+
DV%QUANTITY(1)=='BI-DIRECTIONAL PROBE') THEN
1410014107
IF (DV%ORIENTATION_INDEX==0) THEN
1410114108
WRITE(MESSAGE,'(3A)') 'ERROR(887): DEVC ',TRIM(ID),' must have an ORIENTATION.'
1410214109
CALL SHUTDOWN(MESSAGE) ; RETURN

0 commit comments

Comments
 (0)