TALK=T;RUN(1,1) DISPLAY A hydrocarbon fuel enters steadily a duct of uniform cross- section steadily, through the porous wall of which is supplied a steady stream of air, at a low rate for the first 3/4 of the length, and at an increased rate over the last 1/4 . The air and fuel are supposed to react in a mixed-is-burned manner; but smoke is presumed to be generated at a rate which is proportional to the concentration of unburned fuel times the absolute temperature raised to the power SMOEXP. The multi-fuel model of turbulence is used, but with arbitrarily chosen values of the micro-mixing rate, CONMIX * RATE. The situation is discussed in the report by DB Spalding, dated 1998 ENDDIS PHOTON USE AUTOPLOT file; phi 5 cl; da 1; prod; da 1; oxid; da 1; fuel; scale; cola 1; col3 2 colf 3 msg mass fractions of prod - yellow, oxid - blue, free fuel -red msg msg flow is from left to right pause cl; da 1; smos; da 1; smom; f11; scale; cola 1; col3 2; msg mass fractions of smoke according to: msg msg single-fluid - yellow; multi-fluid - blue pause cl; da 1; f1; da 1; f10; da 1; f20; scale; cola 1; col3 2; colf 3 msg mass fractions of f1 - yellow, f10 - blue, f20 - red pause cl; da 1; f1; da 1; f3; da 1; f5; da 1; f7; da 1; f9; da 1; f11; da 1; f13; da 1; f15; da 1; f17; da 1; f19; da 1; f20; scale; cola 1; col1 2; col1 3; col1 4; col1 5; col3 6; col1 7; col1 8; col1 9; col1 10; colf 11 msg mass fractions of some intermediate fluids also pause cl; da 1; mnsq; cola 1 msg root-mean-square fluctuations pause cl; da 1; tmp1; cola 1 ; da 1; tmpm; col3 2 msg temperatures: single fluid - yellow; multi-fluid - blue pause cl; da 1; fuel; cola 1; da 1; frfm; col3 2; msg free fuel from single- (yellow) and multi-fluid models (blue) pause enduse DISPLAY Notes: (1) four macros are provided in order reduce the repetiveness of MFM q1s. (2) they are named mfm, mfm1, mfm2 and mfm3, and are normally called in that order at the top of the q1 file, which then appears as: #mfm settings of nflr, nflf, etc #mfm1 settings of #mfm2 #mfm3 (3) the character variable mfm is declared and set in the always-loaded core-library macro 014.htm, which, in its turn, declares and sets the character variables mfm1, mfm2, mfm3. (4) the values are: mfm = $L004 mfm1 = $L003 mfm2 = $L002 mfm3 = $L001 ENDDIS 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 ************************************************************ Group 1. Run Title and Number ************************************************************ ************************************************************ TEXT(MFM; 1dsmoke; conmix=5. ) ************************************************************ ************************************************************ IRUNN = 1 ;LIBREF = 230 ************************************************************ Group 2. Time dependence STEADY = T ************************************************************ Group 3. X-Direction Grid Spacing CARTES = T NX = 50 XULAST =1. XFRAC(1)=0.02 ;XFRAC(2)=0.04 XFRAC(3)=0.06 ;XFRAC(4)=0.08 XFRAC(5)=0.1 ;XFRAC(6)=0.12 XFRAC(7)=0.14 ;XFRAC(8)=0.16 XFRAC(9)=0.18 ;XFRAC(10)=0.2 XFRAC(11)=0.22 ;XFRAC(12)=0.24 XFRAC(13)=0.26 ;XFRAC(14)=0.28 XFRAC(15)=0.3 ;XFRAC(16)=0.32 XFRAC(17)=0.34 ;XFRAC(18)=0.36 XFRAC(19)=0.38 ;XFRAC(20)=0.4 XFRAC(21)=0.42 ;XFRAC(22)=0.44 XFRAC(23)=0.46 ;XFRAC(24)=0.48 XFRAC(25)=0.5 ;XFRAC(26)=0.52 XFRAC(27)=0.54 ;XFRAC(28)=0.56 XFRAC(29)=0.58 ;XFRAC(30)=0.6 XFRAC(31)=0.62 ;XFRAC(32)=0.64 XFRAC(33)=0.66 ;XFRAC(34)=0.68 XFRAC(35)=0.7 ;XFRAC(36)=0.72 XFRAC(37)=0.74 ;XFRAC(38)=0.76 XFRAC(39)=0.78 ;XFRAC(40)=0.8 XFRAC(41)=0.82 ;XFRAC(42)=0.84 XFRAC(43)=0.86 ;XFRAC(44)=0.88 XFRAC(45)=0.9 ;XFRAC(46)=0.92 XFRAC(47)=0.94 ;XFRAC(48)=0.96 XFRAC(49)=0.98 ;XFRAC(50)=1. ************************************************************ Group 4. Y-Direction Grid Spacing NY = 1 YVLAST =0.1 YFRAC(1)=1. ************************************************************ Group 5. Z-Direction Grid Spacing PARAB = F NZ = 1 ZWLAST =1. ZFRAC(1)=1. ************************************************************ Group 6. Body-Fitted Coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T NAME(1)=P1 ;NAME(3)=U1 NAME(14)=H1 ;NAME(75)=T1C NAME(76)=T2C ;NAME(77)=T3C NAME(78)=T4C ;NAME(79)=T5C NAME(80)=T6C ;NAME(81)=T7C NAME(82)=T8C ;NAME(83)=T9C NAME(84)=T10C ;NAME(85)=T11C NAME(86)=T12C ;NAME(87)=T13C NAME(88)=T14C ;NAME(89)=T15C NAME(90)=T16C ;NAME(91)=T17C NAME(92)=T18C ;NAME(93)=T19C NAME(94)=T20C ;NAME(95)=S1C NAME(96)=S2C ;NAME(97)=S3C NAME(98)=S4C ;NAME(99)=S5C NAME(100)=S6C ;NAME(101)=S7C NAME(102)=S8C ;NAME(103)=S9C NAME(104)=S10C ;NAME(105)=S11C NAME(106)=S12C ;NAME(107)=S13C NAME(108)=S14C ;NAME(109)=S15C NAME(110)=S16C ;NAME(111)=S17C NAME(112)=S18C ;NAME(113)=S19C NAME(114)=S20C ;NAME(115)=PRDM NAME(116)=OXDM ;NAME(117)=FRFM NAME(118)=TMPM ;NAME(119)=SMOS NAME(120)=SMOM ;NAME(121)=F1 NAME(122)=F2 ;NAME(123)=F3 NAME(124)=F4 ;NAME(125)=F5 NAME(126)=F6 ;NAME(127)=F7 NAME(128)=F8 ;NAME(129)=F9 NAME(130)=F10 ;NAME(131)=F11 NAME(132)=F12 ;NAME(133)=F13 NAME(134)=F14 ;NAME(135)=F15 NAME(136)=F16 ;NAME(137)=F17 NAME(138)=F18 ;NAME(139)=F19 NAME(140)=F20 ;NAME(141)=MNSQ NAME(142)=AVEF ;NAME(143)=SPH1 NAME(144)=RHO1 ;NAME(145)=TMP1 NAME(146)=PROD ;NAME(147)=OXID NAME(148)=FUEL ;NAME(149)=RATE NAME(150)=MIXF * Y in SOLUTN argument list denotes: * 1-stored 2-solved 3-whole-field * 4-point-by-point 5-explicit 6-harmonic averaging SOLUTN(P1,Y,Y,N,N,N,Y) SOLUTN(U1,Y,Y,N,N,N,Y) SOLUTN(H1,Y,Y,N,N,N,Y) SOLUTN(T1C,Y,Y,N,N,N,Y) SOLUTN(T2C,Y,Y,N,N,N,Y) SOLUTN(T3C,Y,Y,N,N,N,Y) SOLUTN(T4C,Y,Y,N,N,N,Y) SOLUTN(T5C,Y,Y,N,N,N,Y) SOLUTN(T6C,Y,Y,N,N,N,Y) SOLUTN(T7C,Y,Y,N,N,N,Y) SOLUTN(T8C,Y,Y,N,N,N,Y) SOLUTN(T9C,Y,Y,N,N,N,Y) SOLUTN(T10C,Y,Y,N,N,N,Y) SOLUTN(T11C,Y,Y,N,N,N,Y) SOLUTN(T12C,Y,Y,N,N,N,Y) SOLUTN(T13C,Y,Y,N,N,N,Y) SOLUTN(T14C,Y,Y,N,N,N,Y) SOLUTN(T15C,Y,Y,N,N,N,Y) SOLUTN(T16C,Y,Y,N,N,N,Y) SOLUTN(T17C,Y,Y,N,N,N,Y) SOLUTN(T18C,Y,Y,N,N,N,Y) SOLUTN(T19C,Y,Y,N,N,N,Y) SOLUTN(T20C,Y,Y,N,N,N,Y) SOLUTN(S1C,Y,Y,N,N,N,Y) SOLUTN(S2C,Y,Y,N,N,N,Y) SOLUTN(S3C,Y,Y,N,N,N,Y) SOLUTN(S4C,Y,Y,N,N,N,Y) SOLUTN(S5C,Y,Y,N,N,N,Y) SOLUTN(S6C,Y,Y,N,N,N,Y) SOLUTN(S7C,Y,Y,N,N,N,Y) SOLUTN(S8C,Y,Y,N,N,N,Y) SOLUTN(S9C,Y,Y,N,N,N,Y) SOLUTN(S10C,Y,Y,N,N,N,Y) SOLUTN(S11C,Y,Y,N,N,N,Y) SOLUTN(S12C,Y,Y,N,N,N,Y) SOLUTN(S13C,Y,Y,N,N,N,Y) SOLUTN(S14C,Y,Y,N,N,N,Y) SOLUTN(S15C,Y,Y,N,N,N,Y) SOLUTN(S16C,Y,Y,N,N,N,Y) SOLUTN(S17C,Y,Y,N,N,N,Y) SOLUTN(S18C,Y,Y,N,N,N,Y) SOLUTN(S19C,Y,Y,N,N,N,Y) SOLUTN(S20C,Y,Y,N,N,N,Y) SOLUTN(PRDM,Y,N,N,N,N,Y) SOLUTN(OXDM,Y,N,N,N,N,Y) SOLUTN(FRFM,Y,N,N,N,N,Y) SOLUTN(TMPM,Y,N,N,N,N,Y) SOLUTN(SMOS,Y,Y,N,N,N,Y) SOLUTN(SMOM,Y,Y,N,N,N,Y) SOLUTN(F1,Y,Y,N,N,N,Y) SOLUTN(F2,Y,Y,N,N,N,Y) SOLUTN(F3,Y,Y,N,N,N,Y) SOLUTN(F4,Y,Y,N,N,N,Y) SOLUTN(F5,Y,Y,N,N,N,Y) SOLUTN(F6,Y,Y,N,N,N,Y) SOLUTN(F7,Y,Y,N,N,N,Y) SOLUTN(F8,Y,Y,N,N,N,Y) SOLUTN(F9,Y,Y,N,N,N,Y) SOLUTN(F10,Y,Y,N,N,N,Y) SOLUTN(F11,Y,Y,N,N,N,Y) SOLUTN(F12,Y,Y,N,N,N,Y) SOLUTN(F13,Y,Y,N,N,N,Y) SOLUTN(F14,Y,Y,N,N,N,Y) SOLUTN(F15,Y,Y,N,N,N,Y) SOLUTN(F16,Y,Y,N,N,N,Y) SOLUTN(F17,Y,Y,N,N,N,Y) SOLUTN(F18,Y,Y,N,N,N,Y) SOLUTN(F19,Y,Y,N,N,N,Y) SOLUTN(F20,Y,Y,N,N,N,Y) SOLUTN(MNSQ,Y,N,N,N,N,Y) SOLUTN(AVEF,Y,N,N,N,N,Y) SOLUTN(SPH1,Y,N,N,N,N,Y) SOLUTN(RHO1,Y,N,N,N,N,Y) SOLUTN(TMP1,Y,N,N,N,N,Y) SOLUTN(PROD,Y,N,N,N,N,Y) SOLUTN(OXID,Y,N,N,N,N,Y) SOLUTN(FUEL,Y,N,N,N,N,Y) SOLUTN(RATE,Y,N,N,N,N,Y) SOLUTN(MIXF,Y,Y,N,N,N,Y) DEN1 = 144 TEMP1 = 145 ************************************************************ Group 8. Terms & Devices * Y in TERMS argument list denotes: * 1-built-in source 2-convection 3-diffusion 4-transient * 5-first phase variable 6-interphase transport TERMS(P1,Y,Y,Y,N,Y,Y) TERMS(U1,Y,Y,Y,Y,Y,Y) TERMS(H1,N,Y,Y,N,Y,N) TERMS(T1C,N,Y,Y,Y,N,Y) TERMS(T2C,N,Y,Y,Y,Y,Y) TERMS(T3C,N,Y,Y,Y,N,Y) TERMS(T4C,N,Y,Y,Y,Y,Y) TERMS(T5C,N,Y,Y,Y,N,Y) TERMS(T6C,N,Y,Y,Y,Y,Y) TERMS(T7C,N,Y,Y,Y,N,Y) TERMS(T8C,N,Y,Y,Y,Y,Y) TERMS(T9C,N,Y,Y,Y,N,Y) TERMS(T10C,N,Y,Y,Y,Y,Y) TERMS(T11C,N,Y,Y,Y,N,Y) TERMS(T12C,N,Y,Y,Y,Y,Y) TERMS(T13C,N,Y,Y,Y,N,Y) TERMS(T14C,N,Y,Y,Y,Y,Y) TERMS(T15C,N,Y,Y,Y,N,Y) TERMS(T16C,N,Y,Y,Y,Y,Y) TERMS(T17C,N,Y,Y,Y,N,Y) TERMS(T18C,N,Y,Y,Y,Y,Y) TERMS(T19C,N,Y,Y,Y,N,Y) TERMS(T20C,N,Y,Y,Y,Y,Y) TERMS(S1C,N,Y,Y,Y,N,Y) TERMS(S2C,N,Y,Y,Y,Y,Y) TERMS(S3C,N,Y,Y,Y,N,Y) TERMS(S4C,N,Y,Y,Y,Y,Y) TERMS(S5C,N,Y,Y,Y,N,Y) TERMS(S6C,N,Y,Y,Y,Y,Y) TERMS(S7C,N,Y,Y,Y,N,Y) TERMS(S8C,N,Y,Y,Y,Y,Y) TERMS(S9C,N,Y,Y,Y,N,Y) TERMS(S10C,N,Y,Y,Y,Y,Y) TERMS(S11C,N,Y,Y,Y,N,Y) TERMS(S12C,N,Y,Y,Y,Y,Y) TERMS(S13C,N,Y,Y,Y,N,Y) TERMS(S14C,N,Y,Y,Y,Y,Y) TERMS(S15C,N,Y,Y,Y,N,Y) TERMS(S16C,N,Y,Y,Y,Y,Y) TERMS(S17C,N,Y,Y,Y,N,Y) TERMS(S18C,N,Y,Y,Y,Y,Y) TERMS(S19C,N,Y,Y,Y,N,Y) TERMS(S20C,N,Y,Y,Y,Y,Y) TERMS(SMOS,N,Y,Y,Y,N,Y) TERMS(SMOM,N,Y,Y,Y,Y,Y) TERMS(F1,N,Y,Y,Y,N,Y) TERMS(F2,N,Y,Y,Y,Y,Y) TERMS(F3,N,Y,Y,Y,N,Y) TERMS(F4,N,Y,Y,Y,Y,Y) TERMS(F5,N,Y,Y,Y,N,Y) TERMS(F6,N,Y,Y,Y,Y,Y) TERMS(F7,N,Y,Y,Y,N,Y) TERMS(F8,N,Y,Y,Y,Y,Y) TERMS(F9,N,Y,Y,Y,N,Y) TERMS(F10,N,Y,Y,Y,Y,Y) TERMS(F11,N,Y,Y,Y,N,Y) TERMS(F12,N,Y,Y,Y,Y,Y) TERMS(F13,N,Y,Y,Y,N,Y) TERMS(F14,N,Y,Y,Y,Y,Y) TERMS(F15,N,Y,Y,Y,N,Y) TERMS(F16,N,Y,Y,Y,Y,Y) TERMS(F17,N,Y,Y,Y,N,Y) TERMS(F18,N,Y,Y,Y,Y,Y) TERMS(F19,N,Y,Y,Y,N,Y) TERMS(F20,N,Y,Y,Y,Y,Y) TERMS(MIXF,N,Y,Y,Y,Y,Y) DIFCUT =0.5 ;ZDIFAC =1. GALA = F ;ADDDIF = F NEWRH1 = T ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1 DENPCO = T ************************************************************ Group 9. Properties used if PRPS is not stored, and where PRPS = -1.0 if it is! RHO1 = GRND6 ;TMP1 = GRND7 ;EL1 =0. TSURR =0. ;TEMP0 =0. ;PRESS0 =1.0E+05 DVO1DT =0. ;DRH1DP =0. RHO1A =16. ;RHO1B =29. ;RHO1C =28. TMP1A =1500. ;TMP1B =1500. ;TMP1C =1500. TMP2A =0.054825 ;TMP2B =4.9E+07 ;TMP2C =1. EMISS =0. ;SCATT =0. RADIA =0. ;RADIB =0. ENUL =1.0E-05 ;ENUT =0. PRNDTL(U1)=1. ;PRNDTL(H1)=1. PRNDTL(T1C)=1. ;PRNDTL(T2C)=1. PRNDTL(T3C)=1. ;PRNDTL(T4C)=1. PRNDTL(T5C)=1. ;PRNDTL(T6C)=1. PRNDTL(T7C)=1. ;PRNDTL(T8C)=1. PRNDTL(T9C)=1. ;PRNDTL(T10C)=1. PRNDTL(T11C)=1. ;PRNDTL(T12C)=1. PRNDTL(T13C)=1. ;PRNDTL(T14C)=1. PRNDTL(T15C)=1. ;PRNDTL(T16C)=1. PRNDTL(T17C)=1. ;PRNDTL(T18C)=1. PRNDTL(T19C)=1. ;PRNDTL(T20C)=1. PRNDTL(S1C)=1. ;PRNDTL(S2C)=1. PRNDTL(S3C)=1. ;PRNDTL(S4C)=1. PRNDTL(S5C)=1. ;PRNDTL(S6C)=1. PRNDTL(S7C)=1. ;PRNDTL(S8C)=1. PRNDTL(S9C)=1. ;PRNDTL(S10C)=1. PRNDTL(S11C)=1. ;PRNDTL(S12C)=1. PRNDTL(S13C)=1. ;PRNDTL(S14C)=1. PRNDTL(S15C)=1. ;PRNDTL(S16C)=1. PRNDTL(S17C)=1. ;PRNDTL(S18C)=1. PRNDTL(S19C)=1. ;PRNDTL(S20C)=1. PRNDTL(SMOS)=1. ;PRNDTL(SMOM)=1. PRNDTL(F1)=1. ;PRNDTL(F2)=1. PRNDTL(F3)=1. ;PRNDTL(F4)=1. PRNDTL(F5)=1. ;PRNDTL(F6)=1. PRNDTL(F7)=1. ;PRNDTL(F8)=1. PRNDTL(F9)=1. ;PRNDTL(F10)=1. PRNDTL(F11)=1. ;PRNDTL(F12)=1. PRNDTL(F13)=1. ;PRNDTL(F14)=1. PRNDTL(F15)=1. ;PRNDTL(F16)=1. PRNDTL(F17)=1. ;PRNDTL(F18)=1. PRNDTL(F19)=1. ;PRNDTL(F20)=1. PRNDTL(MIXF)=1. PRT(U1)=1. ;PRT(H1)=1. PRT(T1C)=1. ;PRT(T2C)=1. PRT(T3C)=1. ;PRT(T4C)=1. PRT(T5C)=1. ;PRT(T6C)=1. PRT(T7C)=1. ;PRT(T8C)=1. PRT(T9C)=1. ;PRT(T10C)=1. PRT(T11C)=1. ;PRT(T12C)=1. PRT(T13C)=1. ;PRT(T14C)=1. PRT(T15C)=1. ;PRT(T16C)=1. PRT(T17C)=1. ;PRT(T18C)=1. PRT(T19C)=1. ;PRT(T20C)=1. PRT(S1C)=1. ;PRT(S2C)=1. PRT(S3C)=1. ;PRT(S4C)=1. PRT(S5C)=1. ;PRT(S6C)=1. PRT(S7C)=1. ;PRT(S8C)=1. PRT(S9C)=1. ;PRT(S10C)=1. PRT(S11C)=1. ;PRT(S12C)=1. PRT(S13C)=1. ;PRT(S14C)=1. PRT(S15C)=1. ;PRT(S16C)=1. PRT(S17C)=1. ;PRT(S18C)=1. PRT(S19C)=1. ;PRT(S20C)=1. PRT(SMOS)=1. ;PRT(SMOM)=1. PRT(F1)=1. ;PRT(F2)=1. PRT(F3)=1. ;PRT(F4)=1. PRT(F5)=1. ;PRT(F6)=1. PRT(F7)=1. ;PRT(F8)=1. PRT(F9)=1. ;PRT(F10)=1. PRT(F11)=1. ;PRT(F12)=1. PRT(F13)=1. ;PRT(F14)=1. PRT(F15)=1. ;PRT(F16)=1. PRT(F17)=1. ;PRT(F18)=1. PRT(F19)=1. ;PRT(F20)=1. PRT(MIXF)=1. CP1 = GRND10 ;CP2 =1. CP1A =1500. ;CP1B =1500. ;CP1C =1500. ;CP1D =0. ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initial field variables (PHIs) FIINIT(P1)=1.0E-10 ;FIINIT(U1)=2. FIINIT(H1)=1.1595E+06 ;FIINIT(T1C)=1.0E-10 FIINIT(T2C)=1.0E-10 ;FIINIT(T3C)=1.0E-10 FIINIT(T4C)=1.0E-10 ;FIINIT(T5C)=1.0E-10 FIINIT(T6C)=1.0E-10 ;FIINIT(T7C)=1.0E-10 FIINIT(T8C)=1.0E-10 ;FIINIT(T9C)=1.0E-10 FIINIT(T10C)=1.0E-10 ;FIINIT(T11C)=1.0E-10 FIINIT(T12C)=1.0E-10 ;FIINIT(T13C)=1.0E-10 FIINIT(T14C)=1.0E-10 ;FIINIT(T15C)=1.0E-10 FIINIT(T16C)=1.0E-10 ;FIINIT(T17C)=1.0E-10 FIINIT(T18C)=1.0E-10 ;FIINIT(T19C)=1.0E-10 FIINIT(T20C)=1.0E-10 ;FIINIT(S1C)=1.0E-10 FIINIT(S2C)=1.0E-10 ;FIINIT(S3C)=1.0E-10 FIINIT(S4C)=1.0E-10 ;FIINIT(S5C)=1.0E-10 FIINIT(S6C)=1.0E-10 ;FIINIT(S7C)=1.0E-10 FIINIT(S8C)=1.0E-10 ;FIINIT(S9C)=1.0E-10 FIINIT(S10C)=1.0E-10 ;FIINIT(S11C)=1.0E-10 FIINIT(S12C)=1.0E-10 ;FIINIT(S13C)=1.0E-10 FIINIT(S14C)=1.0E-10 ;FIINIT(S15C)=1.0E-10 FIINIT(S16C)=1.0E-10 ;FIINIT(S17C)=1.0E-10 FIINIT(S18C)=1.0E-10 ;FIINIT(S19C)=1.0E-10 FIINIT(S20C)=1.0E-10 ;FIINIT(PRDM)=1.0E-10 FIINIT(OXDM)=1.0E-10 ;FIINIT(FRFM)=1.0E-10 FIINIT(TMPM)=1.0E-10 ;FIINIT(SMOS)=0. FIINIT(SMOM)=0. ;FIINIT(F1)=0.05 FIINIT(F2)=0.05 ;FIINIT(F3)=0.05 FIINIT(F4)=0.05 ;FIINIT(F5)=0.05 FIINIT(F6)=0.05 ;FIINIT(F7)=0.05 FIINIT(F8)=0.05 ;FIINIT(F9)=0.05 FIINIT(F10)=0.05 ;FIINIT(F11)=0.05 FIINIT(F12)=0.05 ;FIINIT(F13)=0.05 FIINIT(F14)=0.05 ;FIINIT(F15)=0.05 FIINIT(F16)=0.05 ;FIINIT(F17)=0.05 FIINIT(F18)=0.05 ;FIINIT(F19)=0.05 FIINIT(F20)=0.05 ;FIINIT(MNSQ)=1.0E-10 FIINIT(AVEF)=1.0E-10 ;FIINIT(SPH1)=1.0E-10 FIINIT(RHO1)=0.353169 ;FIINIT(TMP1)=1.0E-10 FIINIT(PROD)=1.0E-10 ;FIINIT(OXID)=1.0E-10 FIINIT(FUEL)=1. ;FIINIT(RATE)=100. FIINIT(MIXF)=1. No PATCHes yet used for this Group INIADD = F FSWEEP = 1 NAMFI =CHAM ************************************************************ Group 12. Patchwise adjustment of terms Patches for this group are printed with those for Group 13. Their names begin either with GP12 or & ************************************************************ Group 13. Boundary & Special Sources PATCH(FUIN ,WEST , 1, 1, 1, 1, 1, 1, 1, 1) COVAL(FUIN ,P1 , FIXFLU ,0.389704 ) COVAL(FUIN ,U1 ,0. ,2. ) COVAL(FUIN ,H1 ,0. ,5.01595E+07 ) COVAL(FUIN ,S1C ,0. ,0. ) COVAL(FUIN ,S2C ,0. ,0. ) COVAL(FUIN ,S3C ,0. ,0. ) COVAL(FUIN ,S4C ,0. ,0. ) COVAL(FUIN ,S5C ,0. ,0. ) COVAL(FUIN ,S6C ,0. ,0. ) COVAL(FUIN ,S7C ,0. ,0. ) COVAL(FUIN ,S8C ,0. ,0. ) COVAL(FUIN ,S9C ,0. ,0. ) COVAL(FUIN ,S10C,0. ,0. ) COVAL(FUIN ,S11C,0. ,0. ) COVAL(FUIN ,S12C,0. ,0. ) COVAL(FUIN ,S13C,0. ,0. ) COVAL(FUIN ,S14C,0. ,0. ) COVAL(FUIN ,S15C,0. ,0. ) COVAL(FUIN ,S16C,0. ,0. ) COVAL(FUIN ,S17C,0. ,0. ) COVAL(FUIN ,S18C,0. ,0. ) COVAL(FUIN ,S19C,0. ,0. ) COVAL(FUIN ,S20C,0. ,0. ) COVAL(FUIN ,F20 ,0. ,1. ) COVAL(FUIN ,MIXF,0. ,1. ) PATCH(OUTLET ,EAST , 50, 50, 1, 1, 1, 1, 1, 1) COVAL(OUTLET ,P1 ,1. ,0. ) COVAL(OUTLET ,SMOS,0. , SAME ) COVAL(OUTLET ,SMOM,0. , SAME ) COVAL(OUTLET ,F1 ,0. , SAME ) COVAL(OUTLET ,F2 ,0. , SAME ) COVAL(OUTLET ,F3 ,0. , SAME ) COVAL(OUTLET ,F4 ,0. , SAME ) COVAL(OUTLET ,F5 ,0. , SAME ) COVAL(OUTLET ,F6 ,0. , SAME ) COVAL(OUTLET ,F7 ,0. , SAME ) COVAL(OUTLET ,F8 ,0. , SAME ) COVAL(OUTLET ,F9 ,0. , SAME ) COVAL(OUTLET ,F10 ,0. , SAME ) COVAL(OUTLET ,F11 ,0. , SAME ) COVAL(OUTLET ,F12 ,0. , SAME ) COVAL(OUTLET ,F13 ,0. , SAME ) COVAL(OUTLET ,F14 ,0. , SAME ) COVAL(OUTLET ,F15 ,0. , SAME ) COVAL(OUTLET ,F16 ,0. , SAME ) COVAL(OUTLET ,F17 ,0. , SAME ) COVAL(OUTLET ,F18 ,0. , SAME ) COVAL(OUTLET ,F19 ,0. , SAME ) COVAL(OUTLET ,F20 ,0. , SAME ) PATCH(AIRIN ,NORTH , 1, 50, 1, 1, 1, 1, 1, 1) COVAL(AIRIN ,P1 , FIXFLU ,0.988874 ) COVAL(AIRIN ,H1 ,0. ,1.1595E+06 ) COVAL(AIRIN ,S1C ,0. ,0. ) COVAL(AIRIN ,S2C ,0. ,0. ) COVAL(AIRIN ,S3C ,0. ,0. ) COVAL(AIRIN ,S4C ,0. ,0. ) COVAL(AIRIN ,S5C ,0. ,0. ) COVAL(AIRIN ,S6C ,0. ,0. ) COVAL(AIRIN ,S7C ,0. ,0. ) COVAL(AIRIN ,S8C ,0. ,0. ) COVAL(AIRIN ,S9C ,0. ,0. ) COVAL(AIRIN ,S10C,0. ,0. ) COVAL(AIRIN ,S11C,0. ,0. ) COVAL(AIRIN ,S12C,0. ,0. ) COVAL(AIRIN ,S13C,0. ,0. ) COVAL(AIRIN ,S14C,0. ,0. ) COVAL(AIRIN ,S15C,0. ,0. ) COVAL(AIRIN ,S16C,0. ,0. ) COVAL(AIRIN ,S17C,0. ,0. ) COVAL(AIRIN ,S18C,0. ,0. ) COVAL(AIRIN ,S19C,0. ,0. ) COVAL(AIRIN ,S20C,0. ,0. ) COVAL(AIRIN ,F1 ,0. ,1. ) PATCH(AIRIN2 ,NORTH , 37, 50, 1, 1, 1, 1, 1, 1) COVAL(AIRIN2 ,P1 , FIXFLU ,28.253557 ) COVAL(AIRIN2 ,H1 ,0. ,1.1595E+06 ) COVAL(AIRIN2 ,S1C ,0. ,0. ) COVAL(AIRIN2 ,S2C ,0. ,0. ) COVAL(AIRIN2 ,S3C ,0. ,0. ) COVAL(AIRIN2 ,S4C ,0. ,0. ) COVAL(AIRIN2 ,S5C ,0. ,0. ) COVAL(AIRIN2 ,S6C ,0. ,0. ) COVAL(AIRIN2 ,S7C ,0. ,0. ) COVAL(AIRIN2 ,S8C ,0. ,0. ) COVAL(AIRIN2 ,S9C ,0. ,0. ) COVAL(AIRIN2 ,S10C,0. ,0. ) COVAL(AIRIN2 ,S11C,0. ,0. ) COVAL(AIRIN2 ,S12C,0. ,0. ) COVAL(AIRIN2 ,S13C,0. ,0. ) COVAL(AIRIN2 ,S14C,0. ,0. ) COVAL(AIRIN2 ,S15C,0. ,0. ) COVAL(AIRIN2 ,S16C,0. ,0. ) COVAL(AIRIN2 ,S17C,0. ,0. ) COVAL(AIRIN2 ,S18C,0. ,0. ) COVAL(AIRIN2 ,S19C,0. ,0. ) COVAL(AIRIN2 ,S20C,0. ,0. ) COVAL(AIRIN2 ,F1 ,0. ,1. ) PATCH(SMOKE ,PHASEM, 1, 50, 1, 1, 1, 1, 1, 1) COVAL(SMOKE ,T1C , GRND , GRND ) COVAL(SMOKE ,T2C , GRND , GRND ) COVAL(SMOKE ,T3C , GRND , GRND ) COVAL(SMOKE ,T4C , GRND , GRND ) COVAL(SMOKE ,T5C , GRND , GRND ) COVAL(SMOKE ,T6C , GRND , GRND ) COVAL(SMOKE ,T7C , GRND , GRND ) COVAL(SMOKE ,T8C , GRND , GRND ) COVAL(SMOKE ,T9C , GRND , GRND ) COVAL(SMOKE ,T10C, GRND , GRND ) COVAL(SMOKE ,T11C, GRND , GRND ) COVAL(SMOKE ,T12C, GRND , GRND ) COVAL(SMOKE ,T13C, GRND , GRND ) COVAL(SMOKE ,T14C, GRND , GRND ) COVAL(SMOKE ,T15C, GRND , GRND ) COVAL(SMOKE ,T16C, GRND , GRND ) COVAL(SMOKE ,T17C, GRND , GRND ) COVAL(SMOKE ,T18C, GRND , GRND ) COVAL(SMOKE ,T19C, GRND , GRND ) COVAL(SMOKE ,T20C, GRND , GRND ) COVAL(SMOKE ,S1C , GRND , GRND ) COVAL(SMOKE ,S2C , GRND , GRND ) COVAL(SMOKE ,S3C , GRND , GRND ) COVAL(SMOKE ,S4C , GRND , GRND ) COVAL(SMOKE ,S5C , GRND , GRND ) COVAL(SMOKE ,S6C , GRND , GRND ) COVAL(SMOKE ,S7C , GRND , GRND ) COVAL(SMOKE ,S8C , GRND , GRND ) COVAL(SMOKE ,S9C , GRND , GRND ) COVAL(SMOKE ,S10C, GRND , GRND ) COVAL(SMOKE ,S11C, GRND , GRND ) COVAL(SMOKE ,S12C, GRND , GRND ) COVAL(SMOKE ,S13C, GRND , GRND ) COVAL(SMOKE ,S14C, GRND , GRND ) COVAL(SMOKE ,S15C, GRND , GRND ) COVAL(SMOKE ,S16C, GRND , GRND ) COVAL(SMOKE ,S17C, GRND , GRND ) COVAL(SMOKE ,S18C, GRND , GRND ) COVAL(SMOKE ,S19C, GRND , GRND ) COVAL(SMOKE ,S20C, GRND , GRND ) COVAL(SMOKE ,SMOS, FIXFLU , GRND ) COVAL(SMOKE ,SMOM, FIXFLU , GRND ) XCYCLE = F EGWF = T WALLCO = GRND2 ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 500 ;ISWC1 = 1 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1 SELREF = T RESFAC =1.0E-03 ************************************************************ Group 16. Terminate Iterations LITER(P1)=20 ;LITER(U1)=20 LITER(H1)=20 ;LITER(T1C)=20 LITER(T2C)=20 ;LITER(T3C)=20 LITER(T4C)=20 ;LITER(T5C)=20 LITER(T6C)=20 ;LITER(T7C)=20 LITER(T8C)=20 ;LITER(T9C)=20 LITER(T10C)=20 ;LITER(T11C)=20 LITER(T12C)=20 ;LITER(T13C)=20 LITER(T14C)=20 ;LITER(T15C)=20 LITER(T16C)=20 ;LITER(T17C)=20 LITER(T18C)=20 ;LITER(T19C)=20 LITER(T20C)=20 ;LITER(S1C)=20 LITER(S2C)=20 ;LITER(S3C)=20 LITER(S4C)=20 ;LITER(S5C)=20 LITER(S6C)=20 ;LITER(S7C)=20 LITER(S8C)=20 ;LITER(S9C)=20 LITER(S10C)=20 ;LITER(S11C)=20 LITER(S12C)=20 ;LITER(S13C)=20 LITER(S14C)=20 ;LITER(S15C)=20 LITER(S16C)=20 ;LITER(S17C)=20 LITER(S18C)=20 ;LITER(S19C)=20 LITER(S20C)=20 ;LITER(SMOS)=20 LITER(SMOM)=20 ;LITER(F1)=20 LITER(F2)=20 ;LITER(F3)=20 LITER(F4)=20 ;LITER(F5)=20 LITER(F6)=20 ;LITER(F7)=20 LITER(F8)=20 ;LITER(F9)=20 LITER(F10)=20 ;LITER(F11)=20 LITER(F12)=20 ;LITER(F13)=20 LITER(F14)=20 ;LITER(F15)=20 LITER(F16)=20 ;LITER(F17)=20 LITER(F18)=20 ;LITER(F19)=20 LITER(F20)=20 ;LITER(MIXF)=20 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03 ENDIT(H1)=1.0E-03 ;ENDIT(T1C)=1.0E-03 ENDIT(T2C)=1.0E-03 ;ENDIT(T3C)=1.0E-03 ENDIT(T4C)=1.0E-03 ;ENDIT(T5C)=1.0E-03 ENDIT(T6C)=1.0E-03 ;ENDIT(T7C)=1.0E-03 ENDIT(T8C)=1.0E-03 ;ENDIT(T9C)=1.0E-03 ENDIT(T10C)=1.0E-03 ;ENDIT(T11C)=1.0E-03 ENDIT(T12C)=1.0E-03 ;ENDIT(T13C)=1.0E-03 ENDIT(T14C)=1.0E-03 ;ENDIT(T15C)=1.0E-03 ENDIT(T16C)=1.0E-03 ;ENDIT(T17C)=1.0E-03 ENDIT(T18C)=1.0E-03 ;ENDIT(T19C)=1.0E-03 ENDIT(T20C)=1.0E-03 ;ENDIT(S1C)=1.0E-03 ENDIT(S2C)=1.0E-03 ;ENDIT(S3C)=1.0E-03 ENDIT(S4C)=1.0E-03 ;ENDIT(S5C)=1.0E-03 ENDIT(S6C)=1.0E-03 ;ENDIT(S7C)=1.0E-03 ENDIT(S8C)=1.0E-03 ;ENDIT(S9C)=1.0E-03 ENDIT(S10C)=1.0E-03 ;ENDIT(S11C)=1.0E-03 ENDIT(S12C)=1.0E-03 ;ENDIT(S13C)=1.0E-03 ENDIT(S14C)=1.0E-03 ;ENDIT(S15C)=1.0E-03 ENDIT(S16C)=1.0E-03 ;ENDIT(S17C)=1.0E-03 ENDIT(S18C)=1.0E-03 ;ENDIT(S19C)=1.0E-03 ENDIT(S20C)=1.0E-03 ;ENDIT(SMOS)=1.0E-03 ENDIT(SMOM)=1.0E-03 ;ENDIT(F1)=1.0E-03 ENDIT(F2)=1.0E-03 ;ENDIT(F3)=1.0E-03 ENDIT(F4)=1.0E-03 ;ENDIT(F5)=1.0E-03 ENDIT(F6)=1.0E-03 ;ENDIT(F7)=1.0E-03 ENDIT(F8)=1.0E-03 ;ENDIT(F9)=1.0E-03 ENDIT(F10)=1.0E-03 ;ENDIT(F11)=1.0E-03 ENDIT(F12)=1.0E-03 ;ENDIT(F13)=1.0E-03 ENDIT(F14)=1.0E-03 ;ENDIT(F15)=1.0E-03 ENDIT(F16)=1.0E-03 ;ENDIT(F17)=1.0E-03 ENDIT(F18)=1.0E-03 ;ENDIT(F19)=1.0E-03 ENDIT(F20)=1.0E-03 ;ENDIT(MIXF)=1.0E-03 ************************************************************ Group 17. Relaxation RELAX(P1,LINRLX,0.5) RELAX(U1,FALSDT,1.) RELAX(H1,FALSDT,10.) RELAX(T1C,FALSDT,1.0E+09) RELAX(T2C,FALSDT,1.0E+09) RELAX(T3C,FALSDT,1.0E+09) RELAX(T4C,FALSDT,1.0E+09) RELAX(T5C,FALSDT,1.0E+09) RELAX(T6C,FALSDT,1.0E+09) RELAX(T7C,FALSDT,1.0E+09) RELAX(T8C,FALSDT,1.0E+09) RELAX(T9C,FALSDT,1.0E+09) RELAX(T10C,FALSDT,1.0E+09) RELAX(T11C,FALSDT,1.0E+09) RELAX(T12C,FALSDT,1.0E+09) RELAX(T13C,FALSDT,1.0E+09) RELAX(T14C,FALSDT,1.0E+09) RELAX(T15C,FALSDT,1.0E+09) RELAX(T16C,FALSDT,1.0E+09) RELAX(T17C,FALSDT,1.0E+09) RELAX(T18C,FALSDT,1.0E+09) RELAX(T19C,FALSDT,1.0E+09) RELAX(T20C,FALSDT,1.0E+09) RELAX(S1C,FALSDT,1.0E+09) RELAX(S2C,FALSDT,1.0E+09) RELAX(S3C,FALSDT,1.0E+09) RELAX(S4C,FALSDT,1.0E+09) RELAX(S5C,FALSDT,1.0E+09) RELAX(S6C,FALSDT,1.0E+09) RELAX(S7C,FALSDT,1.0E+09) RELAX(S8C,FALSDT,1.0E+09) RELAX(S9C,FALSDT,1.0E+09) RELAX(S10C,FALSDT,1.0E+09) RELAX(S11C,FALSDT,1.0E+09) RELAX(S12C,FALSDT,1.0E+09) RELAX(S13C,FALSDT,1.0E+09) RELAX(S14C,FALSDT,1.0E+09) RELAX(S15C,FALSDT,1.0E+09) RELAX(S16C,FALSDT,1.0E+09) RELAX(S17C,FALSDT,1.0E+09) RELAX(S18C,FALSDT,1.0E+09) RELAX(S19C,FALSDT,1.0E+09) RELAX(S20C,FALSDT,1.0E+09) RELAX(PRDM,LINRLX,1.) RELAX(OXDM,LINRLX,1.) RELAX(FRFM,LINRLX,1.) RELAX(TMPM,LINRLX,1.) RELAX(SMOS,FALSDT,1.0E+09) RELAX(SMOM,FALSDT,1.0E+09) RELAX(F1,LINRLX,0.25) RELAX(F2,LINRLX,0.25) RELAX(F3,LINRLX,0.25) RELAX(F4,LINRLX,0.25) RELAX(F5,LINRLX,0.25) RELAX(F6,LINRLX,0.25) RELAX(F7,LINRLX,0.25) RELAX(F8,LINRLX,0.25) RELAX(F9,LINRLX,0.25) RELAX(F10,LINRLX,0.25) RELAX(F11,LINRLX,0.25) RELAX(F12,LINRLX,0.25) RELAX(F13,LINRLX,0.25) RELAX(F14,LINRLX,0.25) RELAX(F15,LINRLX,0.25) RELAX(F16,LINRLX,0.25) RELAX(F17,LINRLX,0.25) RELAX(F18,LINRLX,0.25) RELAX(F19,LINRLX,0.25) RELAX(F20,LINRLX,0.25) RELAX(MNSQ,LINRLX,1.) RELAX(AVEF,LINRLX,1.) RELAX(SPH1,LINRLX,1.) RELAX(RHO1,LINRLX,0.1) RELAX(TMP1,LINRLX,1.) RELAX(PROD,LINRLX,1.) RELAX(OXID,LINRLX,1.) RELAX(FUEL,LINRLX,1.) RELAX(RATE,LINRLX,1.) RELAX(MIXF,FALSDT,1.0E+09) OVRRLX =0. EXPERT = F ;NNORSL = F ************************************************************ Group 18. Limits VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10 VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06 VARMAX(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10 VARMAX(T1C)=1.0E+10 ;VARMIN(T1C)=-1.0E+10 VARMAX(T2C)=1.0E+10 ;VARMIN(T2C)=-1.0E+10 VARMAX(T3C)=1.0E+10 ;VARMIN(T3C)=-1.0E+10 VARMAX(T4C)=1.0E+10 ;VARMIN(T4C)=-1.0E+10 VARMAX(T5C)=1.0E+10 ;VARMIN(T5C)=-1.0E+10 VARMAX(T6C)=1.0E+10 ;VARMIN(T6C)=-1.0E+10 VARMAX(T7C)=1.0E+10 ;VARMIN(T7C)=-1.0E+10 VARMAX(T8C)=1.0E+10 ;VARMIN(T8C)=-1.0E+10 VARMAX(T9C)=1.0E+10 ;VARMIN(T9C)=-1.0E+10 VARMAX(T10C)=1.0E+10 ;VARMIN(T10C)=-1.0E+10 VARMAX(T11C)=1.0E+10 ;VARMIN(T11C)=-1.0E+10 VARMAX(T12C)=1.0E+10 ;VARMIN(T12C)=-1.0E+10 VARMAX(T13C)=1.0E+10 ;VARMIN(T13C)=-1.0E+10 VARMAX(T14C)=1.0E+10 ;VARMIN(T14C)=-1.0E+10 VARMAX(T15C)=1.0E+10 ;VARMIN(T15C)=-1.0E+10 VARMAX(T16C)=1.0E+10 ;VARMIN(T16C)=-1.0E+10 VARMAX(T17C)=1.0E+10 ;VARMIN(T17C)=-1.0E+10 VARMAX(T18C)=1.0E+10 ;VARMIN(T18C)=-1.0E+10 VARMAX(T19C)=1.0E+10 ;VARMIN(T19C)=-1.0E+10 VARMAX(T20C)=1.0E+10 ;VARMIN(T20C)=-1.0E+10 VARMAX(S1C)=1.0E+10 ;VARMIN(S1C)=-1.0E+10 VARMAX(S2C)=1.0E+10 ;VARMIN(S2C)=-1.0E+10 VARMAX(S3C)=1.0E+10 ;VARMIN(S3C)=-1.0E+10 VARMAX(S4C)=1.0E+10 ;VARMIN(S4C)=-1.0E+10 VARMAX(S5C)=1.0E+10 ;VARMIN(S5C)=-1.0E+10 VARMAX(S6C)=1.0E+10 ;VARMIN(S6C)=-1.0E+10 VARMAX(S7C)=1.0E+10 ;VARMIN(S7C)=-1.0E+10 VARMAX(S8C)=1.0E+10 ;VARMIN(S8C)=-1.0E+10 VARMAX(S9C)=1.0E+10 ;VARMIN(S9C)=-1.0E+10 VARMAX(S10C)=1.0E+10 ;VARMIN(S10C)=-1.0E+10 VARMAX(S11C)=1.0E+10 ;VARMIN(S11C)=-1.0E+10 VARMAX(S12C)=1.0E+10 ;VARMIN(S12C)=-1.0E+10 VARMAX(S13C)=1.0E+10 ;VARMIN(S13C)=-1.0E+10 VARMAX(S14C)=1.0E+10 ;VARMIN(S14C)=-1.0E+10 VARMAX(S15C)=1.0E+10 ;VARMIN(S15C)=-1.0E+10 VARMAX(S16C)=1.0E+10 ;VARMIN(S16C)=-1.0E+10 VARMAX(S17C)=1.0E+10 ;VARMIN(S17C)=-1.0E+10 VARMAX(S18C)=1.0E+10 ;VARMIN(S18C)=-1.0E+10 VARMAX(S19C)=1.0E+10 ;VARMIN(S19C)=-1.0E+10 VARMAX(S20C)=1.0E+10 ;VARMIN(S20C)=-1.0E+10 VARMAX(PRDM)=1.0E+10 ;VARMIN(PRDM)=-1.0E+10 VARMAX(OXDM)=1.0E+10 ;VARMIN(OXDM)=-1.0E+10 VARMAX(FRFM)=1.0E+10 ;VARMIN(FRFM)=-1.0E+10 VARMAX(TMPM)=1.0E+10 ;VARMIN(TMPM)=-1.0E+10 VARMAX(SMOS)=1.0E+10 ;VARMIN(SMOS)=-1.0E+10 VARMAX(SMOM)=1.0E+10 ;VARMIN(SMOM)=-1.0E+10 VARMAX(F1)=1. ;VARMIN(F1)=0. VARMAX(F2)=1. ;VARMIN(F2)=0. VARMAX(F3)=1. ;VARMIN(F3)=0. VARMAX(F4)=1. ;VARMIN(F4)=0. VARMAX(F5)=1. ;VARMIN(F5)=0. VARMAX(F6)=1. ;VARMIN(F6)=0. VARMAX(F7)=1. ;VARMIN(F7)=0. VARMAX(F8)=1. ;VARMIN(F8)=0. VARMAX(F9)=1. ;VARMIN(F9)=0. VARMAX(F10)=1. ;VARMIN(F10)=0. VARMAX(F11)=1. ;VARMIN(F11)=0. VARMAX(F12)=1. ;VARMIN(F12)=0. VARMAX(F13)=1. ;VARMIN(F13)=0. VARMAX(F14)=1. ;VARMIN(F14)=0. VARMAX(F15)=1. ;VARMIN(F15)=0. VARMAX(F16)=1. ;VARMIN(F16)=0. VARMAX(F17)=1. ;VARMIN(F17)=0. VARMAX(F18)=1. ;VARMIN(F18)=0. VARMAX(F19)=1. ;VARMIN(F19)=0. VARMAX(F20)=1. ;VARMIN(F20)=0. VARMAX(MNSQ)=1. ;VARMIN(MNSQ)=0. VARMAX(AVEF)=1. ;VARMIN(AVEF)=0. VARMAX(SPH1)=1.0E+10 ;VARMIN(SPH1)=-1.0E+10 VARMAX(RHO1)=0.706339 ;VARMIN(RHO1)=0.035317 VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=773. VARMAX(PROD)=1.0E+10 ;VARMIN(PROD)=-1.0E+10 VARMAX(OXID)=1.0E+10 ;VARMIN(OXID)=-1.0E+10 VARMAX(FUEL)=1.0E+10 ;VARMIN(FUEL)=-1.0E+10 VARMAX(RATE)=1.0E+10 ;VARMIN(RATE)=-1.0E+10 VARMAX(MIXF)=1.0E+10 ;VARMIN(MIXF)=-1.0E+10 ************************************************************ Group 19. Data transmitted to GROUND PARSOL = F ISG62 = 1 SPEDAT(SET,MFM,FFUIN,R,1.) SPEDAT(SET,MFM,TAIRIN,R,773.) SPEDAT(SET,MFM,TFUIN,R,773.) SPEDAT(SET,MFM,TSTOIC,R,2606.314941) SPEDAT(SET,MFM,FACST,R,0.054825) SPEDAT(SET,MFM,GASCON,R,8314.) SPEDAT(SET,MFM,SPVLMX,R,100.) SPEDAT(SET,MFM,SPVLMN,R,1.0E-02) SPEDAT(SET,MFM,SMOCON,R,1.) SPEDAT(SET,MFM,SMOEXP,R,7.) SPEDAT(SET,MFM,MFMMOD,C,SCRS) SPEDAT(SET,MFM,NFLUIDS,I,20) SPEDAT(SET,MFM,NFLR,I,1) SPEDAT(SET,MFM,NFLF,I,20) SPEDAT(SET,MFM,CONREA,R,0.) SPEDAT(SET,MFM,CONMIX,R,5.) SPEDAT(SET,MFM,VISCON,R,0.) SPEDAT(SET,MFM,RADFACT,R,0.5) SPEDAT(SET,GXMONI,PLOTALL,L,T) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=3.786E+04 ;EX(U1)=77.07 EX(H1)=9.408E+06 ;EX(T1C)=695.700012 EX(T2C)=2280. ;EX(T3C)=2258. EX(T4C)=2166. ;EX(T5C)=2074. EX(T6C)=1982. ;EX(T7C)=1890. EX(T8C)=1798. ;EX(T9C)=1706. EX(T10C)=1614. ;EX(T11C)=1523. EX(T12C)=1431. ;EX(T13C)=1339. EX(T14C)=1247. ;EX(T15C)=1155. EX(T16C)=1063. ;EX(T17C)=971.299988 EX(T18C)=879.5 ;EX(T19C)=787.599976 EX(T20C)=695.700012 ;EX(S1C)=0. EX(S2C)=1.805E-05 ;EX(S3C)=1.53E-04 EX(S4C)=2.297E-04 ;EX(S5C)=2.581E-04 EX(S6C)=2.542E-04 ;EX(S7C)=2.309E-04 EX(S8C)=1.978E-04 ;EX(S9C)=1.616E-04 EX(S10C)=1.267E-04 ;EX(S11C)=9.567E-05 EX(S12C)=6.961E-05 ;EX(S13C)=4.881E-05 EX(S14C)=3.297E-05 ;EX(S15C)=2.142E-05 EX(S16C)=1.337E-05 ;EX(S17C)=8.002E-06 EX(S18C)=4.586E-06 ;EX(S19C)=2.451E-06 EX(S20C)=8.065E-07 ;EX(PRDM)=0.3831 EX(OXDM)=0. ;EX(FRFM)=0.1397 EX(TMPM)=1314. ;EX(SMOS)=1.696E-04 EX(SMOM)=7.553E-05 ;EX(F1)=0. EX(F2)=0.1142 ;EX(F3)=0.07434 EX(F4)=0.05543 ;EX(F5)=0.043 EX(F6)=0.03415 ;EX(F7)=0.02767 EX(F8)=0.02282 ;EX(F9)=0.01912 EX(F10)=0.01626 ;EX(F11)=0.01401 EX(F12)=0.01224 ;EX(F13)=0.01083 EX(F14)=9.705E-03 ;EX(F15)=8.832E-03 EX(F16)=8.191E-03 ;EX(F17)=7.813E-03 EX(F18)=7.901E-03 ;EX(F19)=0.01022 EX(F20)=0.03056 ;EX(MNSQ)=0.1773 EX(AVEF)=0.1683 ;EX(SPH1)=1500. EX(RHO1)=0.2311 ;EX(TMP1)=1997. EX(PROD)=0.6835 ;EX(OXID)=0.1856 EX(FUEL)=0.1309 ;EX(RATE)=100. EX(MIXF)=0.1683 ************************************************************ Group 21. Print-out of Variables INIFLD = F ;SUBWGR = F * Y in OUTPUT argument list denotes: * 1-field 2-correction-eq. monitor 3-selective dumping * 4-whole-field residual 5-spot-value table 6-residual table OUTPUT(P1,Y,N,Y,Y,Y,Y) OUTPUT(U1,Y,N,Y,Y,Y,Y) OUTPUT(H1,Y,N,Y,Y,Y,Y) OUTPUT(T1C,Y,N,Y,Y,Y,Y) OUTPUT(T2C,Y,N,Y,Y,Y,Y) OUTPUT(T3C,Y,N,Y,Y,Y,Y) OUTPUT(T4C,Y,N,Y,Y,Y,Y) OUTPUT(T5C,Y,N,Y,Y,Y,Y) OUTPUT(T6C,Y,N,Y,Y,Y,Y) OUTPUT(T7C,Y,N,Y,Y,Y,Y) OUTPUT(T8C,Y,N,Y,Y,Y,Y) OUTPUT(T9C,Y,N,Y,Y,Y,Y) OUTPUT(T10C,Y,N,Y,Y,Y,Y) OUTPUT(T11C,Y,N,Y,Y,Y,Y) OUTPUT(T12C,Y,N,Y,Y,Y,Y) OUTPUT(T13C,Y,N,Y,Y,Y,Y) OUTPUT(T14C,Y,N,Y,Y,Y,Y) OUTPUT(T15C,Y,N,Y,Y,Y,Y) OUTPUT(T16C,Y,N,Y,Y,Y,Y) OUTPUT(T17C,Y,N,Y,Y,Y,Y) OUTPUT(T18C,Y,N,Y,Y,Y,Y) OUTPUT(T19C,Y,N,Y,Y,Y,Y) OUTPUT(T20C,Y,N,Y,Y,Y,Y) OUTPUT(S1C,Y,N,Y,Y,Y,Y) OUTPUT(S2C,Y,N,Y,Y,Y,Y) OUTPUT(S3C,Y,N,Y,Y,Y,Y) OUTPUT(S4C,Y,N,Y,Y,Y,Y) OUTPUT(S5C,Y,N,Y,Y,Y,Y) OUTPUT(S6C,Y,N,Y,Y,Y,Y) OUTPUT(S7C,Y,N,Y,Y,Y,Y) OUTPUT(S8C,Y,N,Y,Y,Y,Y) OUTPUT(S9C,Y,N,Y,Y,Y,Y) OUTPUT(S10C,Y,N,Y,Y,Y,Y) OUTPUT(S11C,Y,N,Y,Y,Y,Y) OUTPUT(S12C,Y,N,Y,Y,Y,Y) OUTPUT(S13C,Y,N,Y,Y,Y,Y) OUTPUT(S14C,Y,N,Y,Y,Y,Y) OUTPUT(S15C,Y,N,Y,Y,Y,Y) OUTPUT(S16C,Y,N,Y,Y,Y,Y) OUTPUT(S17C,Y,N,Y,Y,Y,Y) OUTPUT(S18C,Y,N,Y,Y,Y,Y) OUTPUT(S19C,Y,N,Y,Y,Y,Y) OUTPUT(S20C,Y,N,Y,Y,Y,Y) OUTPUT(PRDM,Y,N,Y,N,N,N) OUTPUT(OXDM,Y,N,Y,N,N,N) OUTPUT(FRFM,Y,N,Y,N,N,N) OUTPUT(TMPM,Y,N,Y,N,N,N) OUTPUT(SMOS,Y,N,Y,Y,Y,Y) OUTPUT(SMOM,Y,N,Y,Y,Y,Y) OUTPUT(F1,Y,N,Y,Y,Y,Y) OUTPUT(F2,Y,N,Y,Y,Y,Y) OUTPUT(F3,Y,N,Y,Y,Y,Y) OUTPUT(F4,Y,N,Y,Y,Y,Y) OUTPUT(F5,Y,N,Y,Y,Y,Y) OUTPUT(F6,Y,N,Y,Y,Y,Y) OUTPUT(F7,Y,N,Y,Y,Y,Y) OUTPUT(F8,Y,N,Y,Y,Y,Y) OUTPUT(F9,Y,N,Y,Y,Y,Y) OUTPUT(F10,Y,N,Y,Y,Y,Y) OUTPUT(F11,Y,N,Y,Y,Y,Y) OUTPUT(F12,Y,N,Y,Y,Y,Y) OUTPUT(F13,Y,N,Y,Y,Y,Y) OUTPUT(F14,Y,N,Y,Y,Y,Y) OUTPUT(F15,Y,N,Y,Y,Y,Y) OUTPUT(F16,Y,N,Y,Y,Y,Y) OUTPUT(F17,Y,N,Y,Y,Y,Y) OUTPUT(F18,Y,N,Y,Y,Y,Y) OUTPUT(F19,Y,N,Y,Y,Y,Y) OUTPUT(F20,Y,N,Y,Y,Y,Y) OUTPUT(MNSQ,Y,N,Y,N,N,N) OUTPUT(AVEF,Y,N,Y,N,N,N) OUTPUT(SPH1,Y,N,Y,N,N,N) OUTPUT(RHO1,Y,N,Y,N,N,N) OUTPUT(TMP1,Y,N,Y,N,N,N) OUTPUT(PROD,Y,N,Y,N,N,N) OUTPUT(OXID,Y,N,Y,N,N,N) OUTPUT(FUEL,Y,N,Y,N,N,N) OUTPUT(RATE,Y,N,Y,N,N,N) OUTPUT(MIXF,Y,N,Y,Y,Y,Y) ************************************************************ Group 22. Monitor Print-Out IXMON = 37 ;IYMON = 1 ;IZMON = 1 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1 UWATCH = T ;USTEER = T HIGHLO = F ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 1000 ;NUMCLS = 5 NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000 IPLTF = 2 ;IPLTL = 1000 ;NPLT = 5 ISWPRF = 1 ;ISWPRL = 100000 ITABL = 3 ;IPROF = 1 ABSIZ =0.5 ;ORSIZ =0.2 NTZPRF = 1 ;NCOLPF = 50 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 40 PATCH(LONGPLOT,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(LONGPLOT,PROD,0. ,1. ) PLOT(LONGPLOT,OXID,0. ,1. ) PLOT(LONGPLOT,FUEL,0. ,1. ) PLOT(LONGPLOT,MIXF,0. ,1. ) PATCH(LONGPLT2,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(LONGPLT2,P1 ,0. ,0. ) PLOT(LONGPLT2,U1 ,0. ,0. ) PLOT(LONGPLT2,TMP1,0. ,0. ) PATCH(FPROF1 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(FPROF1 ,F1 ,0. ,0. ) PATCH(FPROF2 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(FPROF2 ,F2 ,0. ,0. ) PATCH(FPROF3 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(FPROF3 ,F3 ,0. ,0. ) PATCH(FPROF4 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(FPROF4 ,F4 ,0. ,0. ) PATCH(FPROF5 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(FPROF5 ,F5 ,0. ,0. ) PATCH(FPROF6 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(FPROF6 ,F6 ,0. ,0. ) PATCH(FPROF7 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(FPROF7 ,F7 ,0. ,0. ) PATCH(FPROF8 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(FPROF8 ,F8 ,0. ,0. ) PATCH(FPROF9 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(FPROF9 ,F9 ,0. ,0. ) PATCH(FPROF10 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(FPROF10 ,F10 ,0. ,0. ) PATCH(FPROF11 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(FPROF11 ,F11 ,0. ,0. ) PATCH(FPROF12 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(FPROF12 ,F12 ,0. ,0. ) PATCH(FPROF13 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(FPROF13 ,F13 ,0. ,0. ) PATCH(FPROF14 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(FPROF14 ,F14 ,0. ,0. ) PATCH(FPROF15 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(FPROF15 ,F15 ,0. ,0. ) PATCH(FPROF16 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(FPROF16 ,F16 ,0. ,0. ) PATCH(FPROF17 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(FPROF17 ,F17 ,0. ,0. ) PATCH(FPROF18 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(FPROF18 ,F18 ,0. ,0. ) PATCH(FPROF19 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(FPROF19 ,F19 ,0. ,0. ) PATCH(FPROF20 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(FPROF20 ,F20 ,0. ,0. ) PATCH(TPROF1 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(TPROF1 ,T1C ,0. ,0. ) PATCH(TPROF2 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(TPROF2 ,T2C ,0. ,0. ) PATCH(TPROF3 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(TPROF3 ,T3C ,0. ,0. ) PATCH(TPROF4 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(TPROF4 ,T4C ,0. ,0. ) PATCH(TPROF5 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(TPROF5 ,T5C ,0. ,0. ) PATCH(TPROF6 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(TPROF6 ,T6C ,0. ,0. ) PATCH(TPROF7 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(TPROF7 ,T7C ,0. ,0. ) PATCH(TPROF8 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(TPROF8 ,T8C ,0. ,0. ) PATCH(TPROF9 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(TPROF9 ,T9C ,0. ,0. ) PATCH(TPROF10 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(TPROF10 ,T10C,0. ,0. ) PATCH(TPROF11 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(TPROF11 ,T11C,0. ,0. ) PATCH(TPROF12 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(TPROF12 ,T12C,0. ,0. ) PATCH(TPROF13 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(TPROF13 ,T13C,0. ,0. ) PATCH(TPROF14 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(TPROF14 ,T14C,0. ,0. ) PATCH(TPROF15 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(TPROF15 ,T15C,0. ,0. ) PATCH(TPROF16 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(TPROF16 ,T16C,0. ,0. ) PATCH(TPROF17 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(TPROF17 ,T17C,0. ,0. ) PATCH(TPROF18 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(TPROF18 ,T18C,0. ,0. ) PATCH(TPROF19 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(TPROF19 ,T19C,0. ,0. ) PATCH(TPROF20 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(TPROF20 ,T20C,0. ,0. ) PATCH(SPROF1 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(SPROF1 ,S1C ,0. ,0. ) PATCH(SPROF2 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(SPROF2 ,S2C ,0. ,0. ) PATCH(SPROF3 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(SPROF3 ,S3C ,0. ,0. ) PATCH(SPROF4 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(SPROF4 ,S4C ,0. ,0. ) PATCH(SPROF5 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(SPROF5 ,S5C ,0. ,0. ) PATCH(SPROF6 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(SPROF6 ,S6C ,0. ,0. ) PATCH(SPROF7 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(SPROF7 ,S7C ,0. ,0. ) PATCH(SPROF8 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(SPROF8 ,S8C ,0. ,0. ) PATCH(SPROF9 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(SPROF9 ,S9C ,0. ,0. ) PATCH(SPROF10 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(SPROF10 ,S10C,0. ,0. ) PATCH(SPROF11 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(SPROF11 ,S11C,0. ,0. ) PATCH(SPROF12 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(SPROF12 ,S12C,0. ,0. ) PATCH(SPROF13 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(SPROF13 ,S13C,0. ,0. ) PATCH(SPROF14 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(SPROF14 ,S14C,0. ,0. ) PATCH(SPROF15 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(SPROF15 ,S15C,0. ,0. ) PATCH(SPROF16 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(SPROF16 ,S16C,0. ,0. ) PATCH(SPROF17 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(SPROF17 ,S17C,0. ,0. ) PATCH(SPROF18 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(SPROF18 ,S18C,0. ,0. ) PATCH(SPROF19 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(SPROF19 ,S19C,0. ,0. ) PATCH(SPROF20 ,PROFIL, 1, 50, 1, 1, 1, 1, 1, 1) PLOT(SPROF20 ,S20C,0. ,0. ) ************************************************************ Group 24. Dumps For Restarts SAVE = T ;NOWIPE = F NSAVE =CHAM STOP