TALK=T;RUN(1,1) PHOTON USE p;;;;; up z msg velocity vectors * set vec ref;0.5 vec x 1 sh; gr ou x 1 msg Pressto continue pause; cl msg contour of pressure con p1 x 1 fi;0.002; gr ou x 1 msg Press to continue pause; cl msg contour of SLTM, deduced from MFM con SLTM x 1 fi;0.002; gr ou x 1 pause; cl msg contour of SLTS, dedued from single-fluid model con SLTS x 1 fi;0.002; gr ou x 1 msg Press to continue pause; cl msg contour of c1 con c1 x 1 fi;0.002; gr ou x 1 msg Press to continue pause; cl msg contour of enut con enut x 1 fi;0.002; gr ou x 1 msg Press to continue pause; cl msg contour of ep/ke con RATE x 1 fi;0.002; gr ou x 1 msg Press to continue pause; cl msg contour of ke con ke x 1 fi;0.002; gr ou x 1 pause; cl msg contour of mnsq con mnsq x 1 fi;0.002; gr ou x 1 pause; cl msg contour of avef con avef x 1 fi;0.002; gr ou x 1 pause; cl msg contour of fluid 1 volume fraction con f1 x 1 fi;0.002; gr ou x 1 msg Press to continue pause; cl msg contour of fluid 2 volume fraction con f2 x 1 fi;0.002; gr ou x 1 msg Press to continue pause; cl msg contour of fluid 3 volume fraction con f3 x 1 fi;0.002; gr ou x 1 msg Press to continue pause; cl msg contour of fluid 4 volume fraction con f4 x 1 fi;0.002; gr ou x 1 msg Press to continue pause; cl msg contour of fluid 5 volume fraction con f5 x 1 fi;0.002; gr ou x 1 msg Press to continue pause; cl msg contour of fluid 6 volume fraction con f6 x 1 fi;0.002; gr ou x 1 msg Press to continue pause; cl msg contour of fluid 7 volume fraction con f7 x 1 fi;0.002; gr ou x 1 msg Press to continue pause; cl msg contour of fluid 8 volume fraction con f8 x 1 fi;0.002; gr ou x 1 msg Press to continue pause; cl msg contour of fluid 9 volume fraction con f9 x 1 fi;0.002; gr ou x 1 msg Press to continue pause; cl msg contour of fluid 10 volume fraction con f10 x 1 fi;0.002; gr ou x 1 pause; cl msg contour of fluid 11 volume fraction con f11 x 1 fi;0.002; gr ou x 1 msg Press to continue enduse DISPLAY The geometry of this case is similar to that of case W918; but, unlike that case, ONEPHS=T. There are therefore only 3 velocity components. On the other hand, 11 fluids are distinguished, formed by mixing an acidic fluid with an alkaline one in different proportions. __________|.|________ The sketch illustrates the apparatus and | ||| | the initial state of the two liquids. | upper |.| | | liquid ||| | They are both at rest, and are in contact | |.| acid | at a horizontal interface |..........|||........| | lower |.| alkali | The paddle is supposed to be suddenly set | liquid ||| | in motion. | --------- | |paddle ////\\\\\ | The computational task is to predict both | . | the macro-mixing, represented by the ------------|--------- subsequent distributions of velocity, .< axis of pressure and time-average concentration, symmetry but also the extent to which the two The stirred mixing tank liquids are mixed together at any point. ENDDIS 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 ************************************************************ Group 1. Run Title and Number ************************************************************ ************************************************************ TEXT(MFM; Stirred Tank; Nfluids = 21 ) ************************************************************ ************************************************************ IRUNN = 1 ;LIBREF = 501 ************************************************************ Group 2. Time dependence STEADY = F * Set overall time and no. of steps TFIRST =0. ;TLAST =10. FSTEP = 1 ;LSTEP = 10 TFRAC(1)=0.1 ;TFRAC(2)=0.2 TFRAC(3)=0.3 ;TFRAC(4)=0.4 TFRAC(5)=0.5 ;TFRAC(6)=0.6 TFRAC(7)=0.7 ;TFRAC(8)=0.8 TFRAC(9)=0.9 ;TFRAC(10)=1. ************************************************************ Group 3. X-Direction Grid Spacing CARTES = F NX = 1 XULAST =0.1 XFRAC(1)=1. ************************************************************ Group 4. Y-Direction Grid Spacing NY = 10 YVLAST =0.5 YFRAC(1)=0.1 ;YFRAC(2)=0.2 YFRAC(3)=0.3 ;YFRAC(4)=0.4 YFRAC(5)=0.5 ;YFRAC(6)=0.6 YFRAC(7)=0.7 ;YFRAC(8)=0.8 YFRAC(9)=0.9 ;YFRAC(10)=1. ************************************************************ Group 5. Z-Direction Grid Spacing PARAB = F NZ = 10 ZWLAST =1. ZFRAC(1)=0.1 ;ZFRAC(2)=0.2 ZFRAC(3)=0.3 ;ZFRAC(4)=0.4 ZFRAC(5)=0.5 ;ZFRAC(6)=0.6 ZFRAC(7)=0.7 ;ZFRAC(8)=0.8 ZFRAC(9)=0.9 ;ZFRAC(10)=1. ************************************************************ Group 6. Body-Fitted Coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T NAME(1)=P1 ;NAME(3)=U1 NAME(5)=V1 ;NAME(7)=W1 NAME(12)=KE ;NAME(13)=EP NAME(16)=C1 ;NAME(131)=F1 NAME(132)=F2 ;NAME(133)=F3 NAME(134)=F4 ;NAME(135)=F5 NAME(136)=F6 ;NAME(137)=F7 NAME(138)=F8 ;NAME(139)=F9 NAME(140)=F10 ;NAME(141)=F11 NAME(142)=MNSQ ;NAME(143)=AVEF NAME(144)=SLTS ;NAME(145)=SLTM NAME(146)=MIXL ;NAME(147)=EPKE NAME(148)=RATE ;NAME(149)=LEN1 NAME(150)=ENUT * 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,Y,N,N,Y) SOLUTN(U1,Y,Y,N,N,N,Y) SOLUTN(V1,Y,Y,N,N,N,Y) SOLUTN(W1,Y,Y,N,N,N,Y) SOLUTN(KE,Y,Y,N,N,N,N) SOLUTN(EP,Y,Y,N,N,N,N) SOLUTN(C1,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(MNSQ,Y,N,N,N,N,Y) SOLUTN(AVEF,Y,N,N,N,N,Y) SOLUTN(SLTS,Y,Y,N,N,N,Y) SOLUTN(SLTM,Y,Y,N,N,N,Y) SOLUTN(MIXL,Y,N,N,N,N,N) SOLUTN(EPKE,Y,N,N,N,N,N) SOLUTN(RATE,Y,N,N,N,N,Y) SOLUTN(LEN1,Y,N,N,N,N,Y) SOLUTN(ENUT,Y,N,N,N,N,Y) VIST = 150 LEN1 = 149 ************************************************************ 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(V1,Y,Y,Y,Y,Y,Y) TERMS(W1,Y,Y,Y,Y,Y,Y) TERMS(KE,N,Y,Y,Y,Y,N) TERMS(EP,N,Y,Y,Y,Y,N) TERMS(C1,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(SLTS,N,Y,Y,Y,Y,Y) TERMS(SLTM,N,Y,Y,Y,N,Y) DIFCUT =0.5 ;ZDIFAC =1. GALA = F ;ADDDIF = F NEWENT = T ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1 ************************************************************ Group 9. Properties used if PRPS is not stored, and where PRPS = -1.0 if it is! RHO1 =1000. ;TMP1 =0. ;EL1 = GRND4 TSURR =0. ;TEMP0 =0. ;PRESS0 =0. DVO1DT =0. ;DRH1DP =0. EMISS =0. ;SCATT =0. RADIA =0. ;RADIB =0. EL1A =0. ;EL1B =0. ;EL1C =0. ENUL =1.0E-02 ;ENUT = GRND3 ENUTA =0. ;ENUTB =0. ;ENUTC =0. IENUTA = 0 PRNDTL(U1)=1. ;PRNDTL(V1)=1. PRNDTL(W1)=1. ;PRNDTL(KE)=1. PRNDTL(EP)=1. ;PRNDTL(C1)=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(SLTS)=1. PRNDTL(SLTM)=1. PRT(U1)=1. ;PRT(V1)=1. PRT(W1)=1. ;PRT(KE)=1. PRT(EP)=1.314 ;PRT(C1)=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(SLTS)=1. PRT(SLTM)=1. CP1 =1. ;CP2 =1. ************************************************************ Group 10.Inter-Phase Transfer Processes ************************************************************ Group 11.Initial field variables (PHIs) FIINIT(P1)=1.0E-10 ;FIINIT(U1)=1.0E-10 FIINIT(V1)=1.0E-10 ;FIINIT(W1)=1.0E-10 FIINIT(KE)=1. ;FIINIT(EP)=0.1 FIINIT(C1)=0. ;FIINIT(F1)=0. FIINIT(F2)=0. ;FIINIT(F3)=0. FIINIT(F4)=0. ;FIINIT(F5)=0. FIINIT(F6)=0. ;FIINIT(F7)=0. FIINIT(F8)=0. ;FIINIT(F9)=0. FIINIT(F10)=0. ;FIINIT(F11)=1. FIINIT(MNSQ)=1.0E-10 ;FIINIT(AVEF)=1.0E-10 FIINIT(SLTS)=1.0E-10 ;FIINIT(SLTM)=0. FIINIT(MIXL)=1.0E-10 ;FIINIT(EPKE)=1.0E-10 FIINIT(RATE)=0.1 ;FIINIT(LEN1)=0.2 FIINIT(ENUT)=0.1 PATCH(START ,INIVAL, 1, 1, 1, 10, 1, 5, 1, 10) INIT(START ,C1 ,0. ,1. ) INIT(START ,F1 ,0. ,1. ) INIT(START ,F2 ,0. ,0. ) INIT(START ,F3 ,0. ,0. ) INIT(START ,F4 ,0. ,0. ) INIT(START ,F5 ,0. ,0. ) INIT(START ,F6 ,0. ,0. ) INIT(START ,F7 ,0. ,0. ) INIT(START ,F8 ,0. ,0. ) INIT(START ,F9 ,0. ,0. ) INIT(START ,F10 ,0. ,0. ) INIT(START ,F11 ,0. ,0. ) 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(KESOURCE,PHASEM, 0, 0, 0, 0, 0, 0, 1, 10) COVAL(KESOURCE,KE , GRND4 , GRND4 ) COVAL(KESOURCE,EP , GRND4 , GRND4 ) PATCH(CYLWALL ,NWALL , 1, 1, 10, 10, 1, 10, 1, 10) COVAL(CYLWALL ,U1 , GRND2 ,0. ) COVAL(CYLWALL ,W1 , GRND2 ,0. ) COVAL(CYLWALL ,KE , GRND2 , GRND2 ) COVAL(CYLWALL ,EP , GRND2 , GRND2 ) PATCH(TOP ,HWALL , 1, 1, 1, 10, 10, 10, 1, 10) COVAL(TOP ,U1 , GRND2 ,0. ) COVAL(TOP ,V1 , GRND2 ,0. ) COVAL(TOP ,KE , GRND2 , GRND2 ) COVAL(TOP ,EP , GRND2 , GRND2 ) PATCH(BOTTOM ,LWALL , 1, 1, 1, 10, 1, 1, 1, 10) COVAL(BOTTOM ,U1 , GRND2 ,0. ) COVAL(BOTTOM ,V1 , GRND2 ,0. ) COVAL(BOTTOM ,KE , GRND2 , GRND2 ) COVAL(BOTTOM ,EP , GRND2 , GRND2 ) PATCH(PRELIEF ,HIGH , 1, 1, 10, 10, 10, 10, 1, 100) COVAL(PRELIEF ,P1 ,1. ,0. ) COVAL(PRELIEF ,C1 ,0. ,0. ) COVAL(PRELIEF ,F1 ,0. ,0. ) COVAL(PRELIEF ,F2 ,0. ,0. ) COVAL(PRELIEF ,F3 ,0. ,0. ) COVAL(PRELIEF ,F4 ,0. ,0. ) COVAL(PRELIEF ,F5 ,0. ,0. ) COVAL(PRELIEF ,F6 ,0. ,0. ) COVAL(PRELIEF ,F7 ,0. ,0. ) COVAL(PRELIEF ,F8 ,0. ,0. ) COVAL(PRELIEF ,F9 ,0. ,0. ) COVAL(PRELIEF ,F10 ,0. ,0. ) COVAL(PRELIEF ,F11 ,0. ,0. ) PATCH(PADDLE ,CELL , 1, 1, 5, 5, 2, 2, 1, 10) COVAL(PADDLE ,U1 , FIXVAL ,10. ) PATCH(SALT1 ,PHASEM, 1, 1, 1, 10, 1, 10, 1, 10) COVAL(SALT1 ,SLTM, FIXFLU , GRND ) PATCH(SALT2 ,PHASEM, 1, 1, 1, 10, 1, 10, 1, 10) COVAL(SALT2 ,SLTS, FIXFLU , GRND ) IURVAL = -1 XCYCLE = F EGWF = T WALLCO = GRND2 ************************************************************ Group 14. Downstream Pressure For PARAB ************************************************************ Group 15. Terminate Sweeps LSWEEP = 100 ;ISWC1 = 1 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1 SARAH =1.0E-03 SELREF = T RESFAC =1.0E-03 ************************************************************ Group 16. Terminate Iterations LITER(P1)=20 ;LITER(U1)=10 LITER(V1)=10 ;LITER(W1)=10 LITER(KE)=20 ;LITER(EP)=20 LITER(C1)=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(SLTS)=20 ;LITER(SLTM)=20 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03 ENDIT(V1)=1.0E-03 ;ENDIT(W1)=1.0E-03 ENDIT(KE)=1.0E-03 ;ENDIT(EP)=1.0E-03 ENDIT(C1)=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(SLTS)=1.0E-03 ;ENDIT(SLTM)=1.0E-03 ************************************************************ Group 17. Relaxation RELAX(P1,LINRLX,1.) RELAX(U1,LINRLX,0.1) RELAX(V1,FALSDT,9.999999E-04) RELAX(W1,FALSDT,9.999999E-04) RELAX(KE,LINRLX,0.5) RELAX(EP,LINRLX,0.5) RELAX(C1,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(MNSQ,LINRLX,1.) RELAX(AVEF,LINRLX,1.) RELAX(SLTS,FALSDT,1.0E+09) RELAX(SLTM,FALSDT,1.0E+09) RELAX(MIXL,LINRLX,1.) RELAX(EPKE,LINRLX,1.) RELAX(RATE,LINRLX,1.) RELAX(LEN1,LINRLX,1.) RELAX(ENUT,LINRLX,0.1) KELIN = 3 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(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06 VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06 VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-04 VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10 VARMAX(C1)=1.0E+10 ;VARMIN(C1)=-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(MNSQ)=1. ;VARMIN(MNSQ)=0. VARMAX(AVEF)=1. ;VARMIN(AVEF)=0. VARMAX(SLTS)=1.0E+10 ;VARMIN(SLTS)=-1.0E+10 VARMAX(SLTM)=1.0E+10 ;VARMIN(SLTM)=-1.0E+10 VARMAX(MIXL)=1.0E+10 ;VARMIN(MIXL)=-1.0E+10 VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10 VARMAX(RATE)=1.0E+04 ;VARMIN(RATE)=-1.0E+10 VARMAX(LEN1)=1.0E+10 ;VARMIN(LEN1)=-1.0E+10 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10 ************************************************************ Group 19. Data transmitted to GROUND GENK = T PARSOL = F ISG62 = 1 SPEDAT(SET,STORED,MIXL,C,=LEN1!ZSLFIN) SPEDAT(SET,STORED,RATE,C,=EPKE!ZSLFIN) SPEDAT(SET,MFM,MFMMOD,C,MFM) SPEDAT(SET,MFM,NFLUIDS,I,11) SPEDAT(SET,MFM,NFLR,I,1) SPEDAT(SET,MFM,NFLF,I,11) SPEDAT(SET,MFM,CONREA,R,0.) SPEDAT(SET,MFM,CONMIX,R,5.) SPEDAT(SET,MFM,VISCON,R,0.) SPEDAT(SET,GXMONI,PLOTALL,L,T) ************************************************************ Group 20. Preliminary Printout DISTIL = T ;NULLPR = F NDST = 0 DSTTOL =1.0E-02 EX(P1)=74.32 ;EX(U1)=0. EX(V1)=0.09456 ;EX(W1)=0.1466 EX(KE)=0.05478 ;EX(EP)=0.2388 EX(C1)=0.4929 ;EX(F1)=1.383E-05 EX(F2)=2.812E-05 ;EX(F3)=8.503E-05 EX(F4)=7.441E-03 ;EX(F5)=0.3341 EX(F6)=0.3102 ;EX(F7)=0.1923 EX(F8)=0.119 ;EX(F9)=0.03107 EX(F10)=4.958E-03 ;EX(F11)=7.499E-04 EX(MNSQ)=0.07536 ;EX(AVEF)=0.5198 EX(SLTS)=2.013E-06 ;EX(SLTM)=1.702E-06 EX(MIXL)=0.01619 ;EX(EPKE)=2.002 EX(RATE)=2.002 ;EX(LEN1)=0.01619 EX(ENUT)=1.641E-03 ************************************************************ Group 21. Print-out of Variables IURPRN = -1 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,Y,Y,Y,Y,Y) OUTPUT(U1,N,N,Y,Y,N,N) OUTPUT(V1,Y,Y,Y,Y,Y,Y) OUTPUT(W1,Y,Y,Y,Y,Y,Y) OUTPUT(KE,Y,N,Y,Y,Y,Y) OUTPUT(EP,Y,N,Y,Y,Y,Y) OUTPUT(C1,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(MNSQ,Y,N,Y,N,N,N) OUTPUT(AVEF,Y,N,Y,N,N,N) OUTPUT(SLTS,Y,N,Y,Y,Y,Y) OUTPUT(SLTM,Y,N,Y,Y,Y,Y) OUTPUT(MIXL,Y,N,Y,N,N,N) OUTPUT(EPKE,Y,N,Y,N,N,N) OUTPUT(RATE,Y,N,Y,N,N,N) OUTPUT(LEN1,Y,N,Y,N,N,N) OUTPUT(ENUT,Y,N,Y,N,N,N) ************************************************************ Group 22. Monitor Print-Out IXMON = 1 ;IYMON = 5 ;IZMON = 5 NPRMON = 20 ;NPRMNT = 1 ;TSTSWP = -1 UWATCH = T ;USTEER = T HIGHLO = F ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100000 ;NUMCLS = 5 NTPRIN = 100000 ;ISTPRF = 1 ;ISTPRL = 100000 NYPRIN = 1 ;IYPRF = 1 ;IYPRL = 10000 NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000 XZPR = F ;YZPR = F IPLTF = 1 ;IPLTL = -1 ;NPLT = -1 ISWPRF = 1 ;ISWPRL = 100000 ITABL = 3 ;IPROF = 1 ABSIZ =0.5 ;ORSIZ =0.2 NTZPRF = 1 ;NCOLPF = 50 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20 PATCH(MAP ,CONTUR, 1, 1, 1, 10, 1, 10, 1, 10) PLOT(MAP ,U1 ,0. ,10. ) PLOT(MAP ,V1 ,0. ,10. ) PLOT(MAP ,W1 ,0. ,10. ) PLOT(MAP ,C1 ,0. ,10. ) PATCH(IZEQNZM1,PROFIL, 1, 1, 1, 10, 9, 9, 1, 10) PLOT(IZEQNZM1,W1 ,0. ,0. ) PATCH(IZEQNZD2,PROFIL, 1, 1, 1, 10, 1, 5, 1, 10) PLOT(IZEQNZD2,W1 ,-1. ,-1. ) PATCH(TIMEPLOT,PROFIL, 1, 1, 5, 5, 5, 5, 1, 10) PLOT(TIMEPLOT,U1 ,0. ,0. ) PLOT(TIMEPLOT,V1 ,0. ,0. ) PLOT(TIMEPLOT,W1 ,0. ,0. ) PLOT(TIMEPLOT,C1 ,0. ,0. ) PATCH(PROF1 ,PROFIL, 1, 1, 1, 10, 1, 10, 1, 1) PLOT(PROF1 ,F1 ,0. ,0. ) PATCH(PROF2 ,PROFIL, 1, 1, 1, 10, 1, 10, 1, 1) PLOT(PROF2 ,F2 ,0. ,0. ) PATCH(PROF3 ,PROFIL, 1, 1, 1, 10, 1, 10, 1, 1) PLOT(PROF3 ,F3 ,0. ,0. ) PATCH(PROF4 ,PROFIL, 1, 1, 1, 10, 1, 10, 1, 1) PLOT(PROF4 ,F4 ,0. ,0. ) PATCH(PROF5 ,PROFIL, 1, 1, 1, 10, 1, 10, 1, 1) PLOT(PROF5 ,F5 ,0. ,0. ) PATCH(PROF6 ,PROFIL, 1, 1, 1, 10, 1, 10, 1, 1) PLOT(PROF6 ,F6 ,0. ,0. ) PATCH(PROF7 ,PROFIL, 1, 1, 1, 10, 1, 10, 1, 1) PLOT(PROF7 ,F7 ,0. ,0. ) PATCH(PROF8 ,PROFIL, 1, 1, 1, 10, 1, 10, 1, 1) PLOT(PROF8 ,F8 ,0. ,0. ) PATCH(PROF9 ,PROFIL, 1, 1, 1, 10, 1, 10, 1, 1) PLOT(PROF9 ,F9 ,0. ,0. ) PATCH(PROF10 ,PROFIL, 1, 1, 1, 10, 1, 10, 1, 1) PLOT(PROF10 ,F10 ,0. ,0. ) PATCH(PROF11 ,PROFIL, 1, 1, 1, 10, 1, 10, 1, 1) PLOT(PROF11 ,F11 ,0. ,0. ) PATCH(MIDDLE ,PROFIL, 1, 1, 5, 5, 1, 10, 1, 10) PLOT(MIDDLE ,MNSQ,0. ,0. ) PLOT(MIDDLE ,AVEF,0. ,0. ) PLOT(MIDDLE ,ENUT,0. ,0. ) ************************************************************ Group 24. Dumps For Restarts SAVE = T ;NOWIPE = F NSAVE =CHAM STOP