TALK=T;RUN(1,1)
L231
#L230
Case 231 starts
DISPLAY
This case extends case
230
by treating the temperatures and smoke concentrations as
continuously-varying attributes.
Their names are T1C, T2C etc and S1C, S2C etc
The scalar variable RADFACT is introduced in order that the
effects of radiative heat losses associated with smoke content
may be qualitatively simulated.
ENDDIS
SPEDAT( MFM, RADFACT, R, 0.5)
DO II=1,NFLUIDS
SOLVE(S:NFLUIDS+1-II:C) ! smoke concentrations
COVAL(SMOKE,S:NFLUIDS+1-II:C,GRND,GRND)
COVAL(FUIN,S:NFLUIDS+1-II:C,ONLYMS,0.0)
COVAL(AIRIN,S:NFLUIDS+1-II:C,ONLYMS,0.0)
COVAL(AIRIN2,S:NFLUIDS+1-II:C,ONLYMS,0.0)
ENDDO
DO II=1,NFLUIDS
SOLVE(T:NFLUIDS+1-II:C) ! temperatures
COVAL(SMOKE,T:NFLUIDS+1-II:C,GRND,GRND)
ENDDO
DO II=1,NFLUIDS ! line-printer plots
NAMPROF=TPROF:II:
PATCH(:NAMPROF:,PROFIL,1,NX,1,1,1,1,1,1)
COVAL(:NAMPROF:,T:II:C,0.0,0.0)
ENDDO
DO II=1,NFLUIDS
NAMPROF=SPROF:II:
PATCH(:NAMPROF:,PROFIL,1,NX,1,1,1,1,1,1)
COVAL(:NAMPROF:,S:II:C,0.0,0.0)
ENDDO
EX(T1C )= 6.957E+02; EX(S1C )= 0.0
EX(T2C )= 2.280E+03; EX(S2C )= 1.805E-05
EX(T3C )= 2.258E+03; EX(S3C )= 1.530E-04
EX(T4C )= 2.166E+03; EX(S4C )= 2.297E-04
EX(T5C )= 2.074E+03; EX(S5C )= 2.581E-04
EX(T6C )= 1.982E+03; EX(S6C )= 2.542E-04
EX(T7C )= 1.890E+03; EX(S7C )= 2.309E-04
EX(T8C )= 1.798E+03; EX(S8C )= 1.978E-04
EX(T9C )= 1.706E+03; EX(S9C )= 1.616E-04
EX(T10C)= 1.614E+03; EX(S10C)= 1.267E-04
EX(T11C)= 1.523E+03; EX(S11C)= 9.567E-05
EX(T12C)= 1.431E+03; EX(S12C)= 6.961E-05
EX(T13C)= 1.339E+03; EX(S13C)= 4.881E-05
EX(T14C)= 1.247E+03; EX(S14C)= 3.297E-05
EX(T15C)= 1.155E+03; EX(S15C)= 2.142E-05
EX(T16C)= 1.063E+03; EX(S16C)= 1.337E-05
EX(T17C)= 9.713E+02; EX(S17C)= 8.002E-06
EX(T18C)= 8.795E+02; EX(S18C)= 4.586E-06
EX(T19C)= 7.876E+02; EX(S19C)= 2.451E-06
EX(T20C)= 6.957E+02; EX(S20C)= 8.065E-07
EX(TMPM)= 1.314E+03
EX(SMOS)= 1.696E-04
EX(SMOM)= 7.553E-05
STOP