************************************************************ ************************************************************ ---------------------------------------------------------- PHOENICS 2019 v1.0 - EARTH (C) Copyright 2019 Concentration Heat and Momentum Limited All rights reserved. Address: Bakery House, 40 High Street Wimbledon, London, SW19 5AU, England tel: +44 (0) 20-8947-7651 fax: +44 (0) 20-8879-3497 e-mail: phoenics@cham.co.uk web: www.cham.co.uk ---------------------------------------------------------- Code expiry date is the end : may 2019 ---------------------------------------------------------- Running with 64-bit Single Precision executable Working directory: P:\test\library\492 ************************************************************ Initial estimated storage requirement is 10000000 Information about material properties Total number of SPEDATs is 14 number of materials specified by SPEDATs is 1 solprp = 100 porprp = 198 vacprp = 199 !!!! The properties file is PROPS Properties being read from *props Properties have been read from PROPS Formula used for setting RHO1 Formula used for setting TMP1 PRPS is stored with initial value= =-1. Material properties used for phase 1 are density laminar viscosity turbulent viscosity mixing length specific heat >>> End of property-related data <<< ************************************************************ >>> Data read in from the Q1 file <<< The following text is provided as an example of what, by use of the readq1 command, can be transmitted to EARTH for writing near the top of the RESULT file. Such text should leave columns 1 and 2 blank; and it should not extend beyond the 68th column. **************************************************************** * This library case dates from the earliest days of PHOENICS, * * when Professor WU Chung-Hua ('turbomachinery Wu', who had * * returned to China from the USA) visited CHAM in 1982. * * * * The configuration of the combustion chamber, and its being * * supplied with premixed fuel vapour and air, was proposed by * * Professor Wu's accompanying assistant. * **************************************************************** >>> End of data read in from Q1 <<< ************************************************************ Number of F-array locations available is 10000000 Number used before BFC allowance is 68263 Number used after BFC allowance is 68263 ************************************************************ Standard k-e model constants AK =0.41 EWAL =8.6 CMU =0.5478 CD =0.1643 CMUCD =0.090004 C1E =1.44 C2E =1.92 The foregoing constants may be changed from the Q1 file by setting: SPEDAT(KECONST,name of constant,R,value) ************************************************************ biggest cell volume divided by average is 2.469999 at: ix = 1 iy = 10 iz = 7 xg =0.05236 yg =0.06175 zg =0.08575 smallest cell volume divided by average is 2.1224484E-02 at: ix = 1 iy = 1 iz = 3 xg =0.05236 yg =3.25E-03 zg =0.027 ratio of smallest to biggest is 8.5929111E-03 ************************************************************ Group 1. Run Title and Number ************************************************************ ************************************************************ TEXT(Idealised Gas-Turbine Combustion Chamber) ************************************************************ ************************************************************ IRUNN = 1 ;LIBREF = 492 ************************************************************ Group 2. Time dependence STEADY = T ************************************************************ Group 3. X-Direction Grid Spacing CARTES = F NX = 6 XULAST =0.628318 Method of pairs used for grid setting. XFRAC(1)=-6. ;XFRAC(2)=0.166666 ************************************************************ Group 4. Y-Direction Grid Spacing NY = 10 YVLAST =1. Method of pairs used for grid setting. YFRAC(1)=-5. ;YFRAC(3)=1. YFRAC(5)=1. ;YFRAC(7)=3. ************************************************************ Group 5. Z-Direction Grid Spacing PARAB = F NZ = 13 ZWLAST =1. Method of pairs used for grid setting. ZFRAC(1)=-2. ;ZFRAC(3)=1. ZFRAC(5)=1. ;ZFRAC(7)=1. ZFRAC(9)=1. ;ZFRAC(11)=7. ************************************************************ 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(14)=H1 ;NAME(136)=PRPS NAME(137)=DEGF ;NAME(138)=RRAT NAME(139)=EPKE ;NAME(140)=HPOR NAME(141)=NPOR ;NAME(142)=EPOR NAME(143)=VPOR ;NAME(144)=ENUT NAME(145)=RHO1 ;NAME(146)=TMP1 NAME(147)=PROD ;NAME(148)=OXID NAME(149)=FUEL ;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,Y,N,N,N) 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(H1,Y,Y,N,N,N,Y) SOLUTN(PRPS,Y,N,N,N,N,N) SOLUTN(DEGF,Y,N,N,N,N,N) SOLUTN(RRAT,Y,N,N,N,N,N) SOLUTN(EPKE,Y,N,N,N,N,N) SOLUTN(HPOR,Y,N,N,N,N,N) SOLUTN(NPOR,Y,N,N,N,N,N) SOLUTN(EPOR,Y,N,N,N,N,N) SOLUTN(VPOR,Y,N,N,N,N,N) SOLUTN(ENUT,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,Y,N,N,N,Y) SOLUTN(MIXF,Y,Y,N,N,N,Y) DEN1 = 145 VIST = 144 EPOR = 142 ;HPOR = 140 ;NPOR = 141 ;VPOR = 143 TEMP1 = 146 PRPS = 136 ************************************************************ 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,N) TERMS(U1,Y,Y,Y,N,Y,N) TERMS(V1,Y,Y,Y,N,Y,N) TERMS(W1,Y,Y,Y,N,Y,N) TERMS(KE,N,Y,Y,N,Y,N) TERMS(EP,N,Y,Y,N,Y,N) TERMS(H1,N,Y,Y,N,Y,N) TERMS(FUEL,N,Y,Y,N,Y,N) TERMS(MIXF,N,Y,Y,N,Y,N) DIFCUT =0.5 ;ZDIFAC =1. GALA = F ;ADDDIF = F NEWRH1 = T NEWENT = 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 = GRND8 ;EL1 = GRND4 TSURR =0. ;TEMP0 =0. ;PRESS0 =8.0E+05 DVO1DT =0. ;DRH1DP =0. RHO1A =16. ;RHO1B =29. ;RHO1C =28. TMP1A =0. ;TMP1B =0. ;TMP1C =0. TMP2A =0.054825 ;TMP2B =4.9E+07 ;TMP2C =0. EMISS =0. ;SCATT =0. RADIA =0. ;RADIB =0. EL1A =0. ;EL1B =0. ;EL1C =0. ENUL =9.983361E-06 ;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(H1)=1. PRNDTL(FUEL)=1. ;PRNDTL(MIXF)=1. PRT(U1)=1. ;PRT(V1)=1. PRT(W1)=1. ;PRT(KE)=1. PRT(EP)=1.314 ;PRT(H1)=1. PRT(FUEL)=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)=1.0E-10 FIINIT(V1)=1.0E-10 ;FIINIT(W1)=10. FIINIT(KE)=11.25 ;FIINIT(EP)=353.255707 FIINIT(H1)=1.1595E+06 ;FIINIT(PRPS)=-1. FIINIT(DEGF)=1.0E-10 ;FIINIT(RRAT)=1.0E-10 FIINIT(EPKE)=1.0E-10 ;FIINIT(HPOR)=1. FIINIT(NPOR)=1. ;FIINIT(EPOR)=1. FIINIT(VPOR)=1. ;FIINIT(ENUT)=1.0E-10 FIINIT(RHO1)=3.606 ;FIINIT(TMP1)=773. FIINIT(PROD)=1.0E-10 ;FIINIT(OXID)=1.0E-10 FIINIT(FUEL)=0.054825 ;FIINIT(MIXF)=0.054825 PATCH(CMP1 ,INIVAL, 1, 6, 1, 2, 1, 3, 1, 1) INIT(CMP1 ,PRPS,0. ,199. ) PATCH(CMP2 ,INIVAL, 1, 5, 7, 7, 1, 1, 1, 1) INIT(CMP2 ,EPOR,0. ,0.709 ) PATCH(CMP3 ,INIVAL, 1, 6, 7, 7, 1, 1, 1, 1) INIT(CMP3 ,VPOR,0. ,0.709 ) PATCH(CMP4 ,INIVAL, 1, 6, 7, 7, 1, 1, 1, 1) INIT(CMP4 ,NPOR,0. ,0.418 ) PATCH(CMP5 ,INIVAL, 1, 6, 7, 7, 1, 1, 1, 1) INIT(CMP5 ,HPOR,0. ,1. ) PATCH(CMP6 ,INIVAL, 1, 5, 8, 8, 1, 1, 1, 1) INIT(CMP6 ,EPOR,0. ,0.185 ) PATCH(CMP7 ,INIVAL, 1, 6, 8, 8, 1, 1, 1, 1) INIT(CMP7 ,VPOR,0. ,0.185 ) PATCH(CMP8 ,INIVAL, 1, 6, 8, 8, 1, 1, 1, 1) INIT(CMP8 ,V1 ,0. ,0. ) INIT(CMP8 ,NPOR,0. ,0. ) PATCH(CMP9 ,INIVAL, 1, 6, 8, 8, 1, 1, 1, 1) INIT(CMP9 ,HPOR,0. ,0.877 ) PATCH(CMP10 ,INIVAL, 1, 5, 9, 10, 1, 1, 1, 1) INIT(CMP10 ,U1 ,0. ,0. ) INIT(CMP10 ,EPOR,0. ,0. ) PATCH(CMP11 ,INIVAL, 1, 6, 9, 10, 1, 1, 1, 1) INIT(CMP11 ,PRPS,0. ,199. ) PATCH(CMP12 ,INIVAL, 1, 6, 9, 9, 1, 1, 1, 1) INIT(CMP12 ,V1 ,0. ,0. ) INIT(CMP12 ,NPOR,0. ,0. ) PATCH(CMP13 ,INIVAL, 1, 6, 9, 10, 1, 1, 1, 1) INIT(CMP13 ,W1 ,0. ,0. ) INIT(CMP13 ,HPOR,0. ,0. ) PATCH(CMP14 ,INIVAL, 1, 5, 8, 8, 2, 2, 1, 1) INIT(CMP14 ,EPOR,0. ,0.997 ) PATCH(CMP15 ,INIVAL, 1, 6, 8, 8, 2, 2, 1, 1) INIT(CMP15 ,VPOR,0. ,0.997 ) PATCH(CMP16 ,INIVAL, 1, 6, 8, 8, 2, 2, 1, 1) INIT(CMP16 ,NPOR,0. ,0.946 ) PATCH(CMP17 ,INIVAL, 1, 6, 8, 8, 2, 2, 1, 1) INIT(CMP17 ,HPOR,0. ,1. ) PATCH(CMP18 ,INIVAL, 1, 5, 9, 9, 2, 2, 1, 1) INIT(CMP18 ,EPOR,0. ,0.706 ) PATCH(CMP19 ,INIVAL, 1, 6, 9, 9, 2, 2, 1, 1) INIT(CMP19 ,VPOR,0. ,0.706 ) PATCH(CMP20 ,INIVAL, 1, 6, 9, 9, 2, 2, 1, 1) INIT(CMP20 ,NPOR,0. ,0.473 ) PATCH(CMP21 ,INIVAL, 1, 6, 9, 9, 2, 2, 1, 1) INIT(CMP21 ,HPOR,0. ,1. ) PATCH(CMP22 ,INIVAL, 1, 5, 10, 10, 2, 2, 1, 1) INIT(CMP22 ,EPOR,0. ,0.236 ) PATCH(CMP23 ,INIVAL, 1, 6, 10, 10, 2, 2, 1, 1) INIT(CMP23 ,VPOR,0. ,0.236 ) PATCH(CMP24 ,INIVAL, 1, 6, 10, 10, 2, 2, 1, 1) PATCH(CMP25 ,INIVAL, 1, 6, 10, 10, 2, 2, 1, 1) INIT(CMP25 ,HPOR,0. ,1. ) PATCH(WALL4 ,INIVAL, 1, 6, 2, 2, 1, 3, 1, 1) INIT(WALL4 ,V1 ,0. ,0. ) INIT(WALL4 ,NPOR,0. ,0. ) PATCH(WALL6 ,INIVAL, 1, 6, 1, 2, 3, 3, 1, 1) INIT(WALL6 ,W1 ,0. ,0. ) INIT(WALL6 ,HPOR,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(FOIN ,SOUTH , 1, 6, 3, 3, 3, 3, 1, 1) COVAL(FOIN ,P1 , FIXFLU ,360.600006 ) COVAL(FOIN ,U1 ,0. ,0. ) COVAL(FOIN ,V1 ,0. ,100. ) COVAL(FOIN ,W1 ,0. ,0. ) COVAL(FOIN ,KE ,0. ,0.25 ) COVAL(FOIN ,EP ,0. ,58.752609 ) COVAL(FOIN ,H1 ,0. ,3.845904E+06 ) COVAL(FOIN ,FUEL,0. ,0.054825 ) COVAL(FOIN ,MIXF,0. ,0.054825 ) PATCH(SOIN1 ,CELL , 1, 1, 10, 10, 7, 7, 1, 1) COVAL(SOIN1 ,P1 , FIXFLU ,8.157E-03 ) COVAL(SOIN1 ,U1 ,0. ,0. ) COVAL(SOIN1 ,V1 ,0. ,-40. ) COVAL(SOIN1 ,W1 ,0. ,0. ) COVAL(SOIN1 ,KE ,0. ,0.123 ) COVAL(SOIN1 ,EP ,0. ,8.429416 ) COVAL(SOIN1 ,H1 ,0. ,1.1595E+06 ) COVAL(SOIN1 ,FUEL,0. ,0. ) COVAL(SOIN1 ,MIXF,0. ,0. ) PATCH(SOIN2 ,CELL , 4, 4, 10, 10, 7, 7, 1, 1) COVAL(SOIN2 ,P1 , FIXFLU ,8.157E-03 ) COVAL(SOIN2 ,U1 ,0. ,0. ) COVAL(SOIN2 ,V1 ,0. ,-40. ) COVAL(SOIN2 ,W1 ,0. ,0. ) COVAL(SOIN2 ,KE ,0. ,0.123 ) COVAL(SOIN2 ,EP ,0. ,8.429416 ) COVAL(SOIN2 ,H1 ,0. ,1.1595E+06 ) COVAL(SOIN2 ,FUEL,0. ,0. ) COVAL(SOIN2 ,MIXF,0. ,0. ) PATCH(DILUIN ,CELL , 1, 1, 10, 10, 10, 10, 1, 1) COVAL(DILUIN ,P1 , FIXFLU ,0.0145 ) COVAL(DILUIN ,U1 ,0. ,0. ) COVAL(DILUIN ,V1 ,0. ,-40. ) COVAL(DILUIN ,W1 ,0. ,0. ) COVAL(DILUIN ,KE ,0. ,0.123 ) COVAL(DILUIN ,EP ,0. ,6.299674 ) COVAL(DILUIN ,H1 ,0. ,1.1595E+06 ) COVAL(DILUIN ,FUEL,0. ,0. ) COVAL(DILUIN ,MIXF,0. ,0. ) PATCH(OUTLET ,HIGH , 1, 6, 1, 10, 13, 13, 1, 1) COVAL(OUTLET ,P1 ,1000. ,0. ) COVAL(OUTLET ,U1 ,0. ,0. ) COVAL(OUTLET ,V1 ,0. ,0. ) COVAL(OUTLET ,W1 ,0. ,0. ) COVAL(OUTLET ,KE ,0. , SAME ) COVAL(OUTLET ,EP ,0. , SAME ) COVAL(OUTLET ,H1 ,0. , SAME ) COVAL(OUTLET ,FUEL,0. , SAME ) COVAL(OUTLET ,MIXF,0. , SAME ) PATCH(WALL1 ,NWALL , 1, 6, 10, 10, 3, 6, 1, 1) COVAL(WALL1 ,U1 , GRND2 ,0. ) COVAL(WALL1 ,W1 , GRND2 ,0. ) COVAL(WALL1 ,KE , GRND2 , GRND2 ) COVAL(WALL1 ,EP , GRND2 , GRND2 ) PATCH(WALL2 ,NWALL , 1, 6, 10, 10, 8, 9, 1, 1) COVAL(WALL2 ,U1 , GRND2 ,0. ) COVAL(WALL2 ,W1 , GRND2 ,0. ) COVAL(WALL2 ,KE , GRND2 , GRND2 ) COVAL(WALL2 ,EP , GRND2 , GRND2 ) PATCH(WALL3 ,NWALL , 1, 6, 10, 10, 11, 13, 1, 1) COVAL(WALL3 ,U1 , GRND2 ,0. ) COVAL(WALL3 ,W1 , GRND2 ,0. ) COVAL(WALL3 ,KE , GRND2 , GRND2 ) COVAL(WALL3 ,EP , GRND2 , GRND2 ) PATCH(WALL4-NW,NWALL , 1, 6, 2, 2, 1, 3, 1, 1) COVAL(WALL4-NW,U1 , GRND2 ,0. ) COVAL(WALL4-NW,W1 , GRND2 ,0. ) COVAL(WALL4-NW,KE , GRND2 , GRND2 ) COVAL(WALL4-NW,EP , GRND2 , GRND2 ) PATCH(WALL4-SW,SWALL , 1, 6, 3, 3, 1, 3, 1, 1) COVAL(WALL4-SW,U1 , GRND2 ,0. ) COVAL(WALL4-SW,W1 , GRND2 ,0. ) COVAL(WALL4-SW,KE , GRND2 , GRND2 ) COVAL(WALL4-SW,EP , GRND2 , GRND2 ) PATCH(WALL5 ,LWALL , 1, 6, 3, 6, 1, 1, 1, 1) COVAL(WALL5 ,U1 , GRND2 ,0. ) COVAL(WALL5 ,V1 , GRND2 ,0. ) COVAL(WALL5 ,KE , GRND2 , GRND2 ) COVAL(WALL5 ,EP , GRND2 , GRND2 ) PATCH(WALL6-HW,HWALL , 1, 6, 1, 2, 3, 3, 1, 1) COVAL(WALL6-HW,U1 , GRND2 ,0. ) COVAL(WALL6-HW,V1 , GRND2 ,0. ) COVAL(WALL6-HW,KE , GRND2 , GRND2 ) COVAL(WALL6-HW,EP , GRND2 , GRND2 ) PATCH(WALL6-LW,LWALL , 1, 6, 1, 2, 4, 4, 1, 1) COVAL(WALL6-LW,U1 , GRND2 ,0. ) COVAL(WALL6-LW,V1 , GRND2 ,0. ) COVAL(WALL6-LW,KE , GRND2 , GRND2 ) COVAL(WALL6-LW,EP , GRND2 , GRND2 ) PATCH(KESOURCE,PHASEM, 1, 6, 1, 10, 1, 13, 1, 1) COVAL(KESOURCE,KE , GRND4 , GRND4 ) COVAL(KESOURCE,EP , GRND4 , GRND4 ) PATCH(WHOLE ,PHASEM, 1, 6, 1, 10, 1, 13, 1, 1) COVAL(WHOLE ,FUEL,In-Form:COVAL(1.*EPKE&0.0)) 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. SELREF = T RESFAC =1.0E-05 ************************************************************ Group 16. Terminate Iterations LITER(P1)=20 ;LITER(U1)=20 LITER(V1)=20 ;LITER(W1)=20 LITER(KE)=20 ;LITER(EP)=20 LITER(H1)=20 ;LITER(FUEL)=20 LITER(MIXF)=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(H1)=1.0E-03 ;ENDIT(FUEL)=1.0E-03 ENDIT(MIXF)=1.0E-03 ************************************************************ Group 17. Relaxation RELAX(P1,LINRLX,1.) RELAX(U1,FALSDT,0.02) RELAX(V1,FALSDT,0.02) RELAX(W1,FALSDT,0.02) RELAX(KE,LINRLX,0.5) RELAX(EP,LINRLX,0.5) RELAX(H1,FALSDT,10.) RELAX(PRPS,LINRLX,1.) RELAX(DEGF,LINRLX,1.) RELAX(RRAT,LINRLX,1.) RELAX(EPKE,LINRLX,1.) RELAX(ENUT,LINRLX,1.) RELAX(RHO1,LINRLX,1.) RELAX(TMP1,LINRLX,1.) RELAX(PROD,LINRLX,1.) RELAX(OXID,LINRLX,1.) RELAX(FUEL,FALSDT,1.0E+09) RELAX(MIXF,FALSDT,1.0E+09) 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-10 VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10 VARMAX(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10 VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10 VARMAX(DEGF)=1.0E+10 ;VARMIN(DEGF)=-1.0E+10 VARMAX(RRAT)=1.0E+10 ;VARMIN(RRAT)=-1.0E+10 VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=0. VARMAX(HPOR)=1.0E+10 ;VARMIN(HPOR)=-1.0E+10 VARMAX(NPOR)=1.0E+10 ;VARMIN(NPOR)=-1.0E+10 VARMAX(EPOR)=1.0E+10 ;VARMIN(EPOR)=-1.0E+10 VARMAX(VPOR)=1.0E+10 ;VARMIN(VPOR)=-1.0E+10 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10 VARMAX(RHO1)=1.0E+10 ;VARMIN(RHO1)=-1.0E+10 VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=773. VARMAX(PROD)=1.0E+10 ;VARMIN(PROD)=0. VARMAX(OXID)=1.0E+10 ;VARMIN(OXID)=0. VARMAX(FUEL)=1.0E+10 ;VARMIN(FUEL)=0. VARMAX(MIXF)=1.0E+10 ;VARMIN(MIXF)=-1.0E+10 ************************************************************ Group 19. Data transmitted to GROUND GENK = T READQ1 = T PARSOL = F GEN1 = 681 ISG62 = 0 CHSOA =0.054825 ;CHSOB =1. SPEDAT(SET,PROPERTY,TMP1,C,=MAX(773.&(H1-4.9E+07*FUEL)/(1500.*FUEL+1500.*OXID+15$) SPEDAT(SET,PROPERTY,TMP1,C,00.*PROD))!IMAT>100) SPEDAT(SET,PROPERTY,RHO1,C,=(P1+8.0E+05)*1.0/(FUEL/16.+OXID/29.+PROD/28.)/(TMP1*$) SPEDAT(SET,PROPERTY,RHO1,C,8313.4)) SPEDAT(SET,STORED,PROD,C,=(MIXF-FUEL)*(1+17.24)!ZSLFIN) SPEDAT(SET,STORED,OXID,C,=1-FUEL-PROD!ZSLFIN) SPEDAT(SET,SOURCE,FUEL!WHOLE,C,=COVAL(1.*EPKE&0.0)) SPEDAT(SET,STORED,RRAT,C,=1.*EPKE*FUEL!ZSLFIN) SPEDAT(SET,LONGNAME,TMP1,C,absolute_temperature_of_the_gas_Kelvin) SPEDAT(SET,STORED,DEGF,C,=(9./5.)*(TMP1-273)!ZSLFIN) SPEDAT(SET,LONGNAME,DEGF,C,degrees_Fahrenheit) SPEDAT(SET,LONGNAME,RRAT,C,rate_of_consumption_of_fuel_kg/m^3sec) SPEDAT(SET,GXMONI,PLOTALL,L,T) SPEDAT(SET,MATERIAL,199,L,T) ************************************************************ Group 20. Preliminary Printout ************************************************************ 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(V1,Y,N,Y,Y,Y,Y) OUTPUT(W1,Y,N,Y,Y,Y,Y) OUTPUT(KE,Y,N,Y,Y,Y,Y) OUTPUT(EP,Y,N,Y,Y,Y,Y) OUTPUT(H1,Y,N,Y,Y,Y,Y) OUTPUT(PRPS,Y,N,Y,N,N,N) OUTPUT(DEGF,Y,N,Y,N,N,N) OUTPUT(RRAT,Y,N,Y,N,N,N) OUTPUT(EPKE,Y,N,Y,N,N,N) OUTPUT(HPOR,Y,N,Y,N,N,N) OUTPUT(NPOR,Y,N,Y,N,N,N) OUTPUT(EPOR,Y,N,Y,N,N,N) OUTPUT(VPOR,Y,N,Y,N,N,N) OUTPUT(ENUT,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,Y,Y,Y) OUTPUT(MIXF,Y,N,Y,Y,Y,Y) ************************************************************ Group 22. Monitor Print-Out IXMON = 3 ;IYMON = 5 ;IZMON = 5 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = 10005 UWATCH = F ;USTEER = F HIGHLO = F ************************************************************ Group 23.Field Print-Out & Plot Control NPRINT = 100 ;NUMCLS = 5 NXPRIN = 1 ;IXPRF = 3 ;IXPRL = 3 NYPRIN = 2 ;IYPRF = 1 ;IYPRL = 10 NZPRIN = 2 ;IZPRF = 1 ;IZPRL = 10000 XZPR = F ;YZPR = T IPLTF = 2 ;IPLTL = 100 ;NPLT = 5 ISWPRF = 1 ;ISWPRL = 100000 ITABL = 3 ;IPROF = 1 ABSIZ =0.5 ;ORSIZ =0.4 NTZPRF = 1 ;NCOLPF = 50 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 40 PATCH(IZ4 ,CONTUR, 1, 6, 1, 10, 4, 4, 1, 1) PLOT(IZ4 ,TMP1,1. ,10. ) PATCH(XSECIN1 ,CONTUR, 1, 1, 1, 10, 1, 13, 1, 1) PLOT(XSECIN1 ,W1 ,0. ,10. ) PLOT(XSECIN1 ,TMP1,0. ,10. ) PLOT(XSECIN1 ,FUEL,0. ,10. ) PLOT(XSECIN1 ,MIXF,0. ,10. ) PATCH(XSECIN2 ,CONTUR, 4, 4, 1, 10, 1, 13, 1, 1) PLOT(XSECIN2 ,W1 ,0. ,10. ) PLOT(XSECIN2 ,TMP1,0. ,10. ) PLOT(XSECIN2 ,FUEL,0. ,10. ) PLOT(XSECIN2 ,MIXF,0. ,10. ) ************************************************************ Group 24. Dumps For Restarts SAVE = T ;NOWIPE = F NSAVE =CHAM *** grid-geometry information *** X-coordinates of the cell centres 2.618E-01 Y-coordinates of the cell centres 3.250E-03 1.625E-02 2.925E-02 4.150E-02 5.525E-02 Z-coordinates of the cell centres 6.250E-03 2.700E-02 4.500E-02 8.575E-02 1.347E-01 1.837E-01 2.327E-01 X-coordinates of the (higher) cell faces 3.142E-01 Y-coordinates of the (higher) cell faces 6.500E-03 1.950E-02 3.250E-02 4.550E-02 5.850E-02 Z-coordinates of the (higher) cell faces 1.250E-02 2.900E-02 5.000E-02 9.800E-02 1.470E-01 1.960E-01 2.450E-01 Total number of F-array elements used is 71890 --- INTEGRATION OF EQUATIONS BEGINS --- ************************************************************ *** End of sweep 100 ************************************************************ Flow field at ITHYD= 1, ISWEEP= 100, ISTEP= 1 YZPR IX= 3 isweep = 100 istep = 1 Field Values of P1 IY= 10 none 9.796E+02 3.563E+02 -1.736E+00 2.127E+02 IY= 8 4.778E+02 8.468E+02 4.109E+01 3.230E+02 2.246E+02 IY= 6 4.387E+01 2.842E+02 -1.128E+02 4.074E+02 2.846E+02 IY= 4 4.883E+01 1.947E+02 1.127E+01 4.592E+02 3.650E+02 IY= 2 none none 1.904E+02 5.187E+02 4.912E+02 IZ= 1 3 5 7 9 IY= 10 -4.740E+01 3.571E-02 IY= 8 -8.850E+01 1.934E-02 IY= 6 -1.426E+01 1.974E-02 IY= 4 1.872E+02 4.030E-02 IY= 2 3.435E+02 3.719E-03 IZ= 11 13 Field Values of U1 IY= 10 none 2.915E-02 6.465E-02 2.576E+00 -2.157E+00 IY= 8 -1.370E-02 -1.915E-02 -1.383E-01 -2.574E+00 6.440E-01 IY= 6 -1.121E-02 -3.872E-02 -3.384E-01 1.045E+00 8.901E-01 IY= 4 4.080E-03 -2.520E-02 -2.514E-01 2.164E+00 -1.577E-02 IY= 2 none none -1.172E-02 3.564E-01 -5.558E-02 IZ= 1 3 5 7 9 IY= 10 -7.425E+00 -4.195E+00 IY= 8 -2.858E+00 -8.552E-01 IY= 6 2.923E+00 1.950E+00 IY= 4 3.442E+00 2.453E+00 IY= 2 5.230E-01 4.429E+00 IZ= 11 13 Field Values of V1 IY= 9 none 1.678E+01 2.217E+00 1.526E+00 5.191E-01 IY= 7 -1.105E+01 5.500E+01 1.579E+00 -8.835E+00 2.754E+00 IY= 5 -1.263E+01 8.063E+01 3.095E+00 -1.247E+01 4.359E+00 IY= 3 -5.432E+00 9.071E+01 3.671E+00 -9.781E+00 4.472E+00 IY= 1 none none 1.476E+00 -3.380E+00 2.003E+00 IZ= 1 3 5 7 9 IY= 9 6.906E-02 9.563E-01 IY= 7 -6.035E-01 4.644E-02 IY= 5 -2.390E+00 -9.799E-01 IY= 3 -1.691E+00 5.468E-01 IY= 1 -3.263E-01 7.081E-02 IZ= 11 13 Field Values of W1 IY= 10 none 1.358E+01 3.083E+01 2.373E+01 2.195E+01 IY= 8 -9.177E+00 7.890E+00 1.213E+01 1.214E+01 1.512E+01 IY= 6 -5.487E-01 9.266E-01 -7.051E+00 1.013E+01 1.153E+01 IY= 4 1.369E+01 -2.930E+00 -1.373E+01 9.001E+00 6.593E+00 IY= 2 none none -1.644E+01 9.174E+00 -1.834E+00 IZ= 1 3 5 7 9 IY= 10 1.861E+01 IY= 8 1.468E+01 IY= 6 1.224E+01 IY= 4 1.405E+01 IY= 2 1.836E+01 IZ= 11 Field Values of KE IY= 10 none 8.024E-01 6.826E+00 9.949E+01 3.719E+00 IY= 8 7.849E+02 4.477E+02 2.645E+02 1.206E+02 4.273E+01 IY= 6 1.171E+00 1.223E+02 1.120E+02 8.507E+01 4.908E+01 IY= 4 5.799E-01 4.279E+01 8.078E+01 8.835E+01 6.452E+01 IY= 2 none none 1.015E+02 1.692E+02 9.414E+01 IZ= 1 3 5 7 9 IY= 10 3.999E+00 3.110E+00 IY= 8 3.891E+01 2.581E+01 IY= 6 4.772E+01 2.787E+01 IY= 4 5.387E+01 4.265E+01 IY= 2 7.235E+01 1.041E+02 IZ= 11 13 Field Values of EP IY= 10 none 8.861E+01 2.199E+03 6.335E+04 8.841E+02 IY= 8 6.412E+05 3.917E+05 1.555E+05 6.443E+04 1.169E+04 IY= 6 8.121E+01 1.041E+05 3.458E+04 3.654E+04 1.118E+04 IY= 4 2.829E+01 4.201E+04 1.949E+04 3.354E+04 1.475E+04 IY= 2 none none 2.857E+04 6.470E+04 2.621E+04 IZ= 1 3 5 7 9 IY= 10 9.861E+02 6.760E+02 IY= 8 1.073E+04 4.995E+03 IY= 6 1.840E+04 6.428E+03 IY= 4 2.151E+04 1.013E+04 IY= 2 2.527E+04 3.302E+04 IZ= 11 13 Field Values of H1 IY= 10 none 3.471E+06 3.295E+06 2.990E+06 2.799E+06 IY= 8 3.496E+06 3.482E+06 2.976E+06 2.212E+06 2.150E+06 IY= 6 3.502E+06 3.568E+06 2.669E+06 2.070E+06 2.018E+06 IY= 4 3.502E+06 3.695E+06 2.286E+06 1.995E+06 1.974E+06 IY= 2 none none 2.086E+06 1.988E+06 1.916E+06 IZ= 1 3 5 7 9 IY= 10 2.437E+06 2.157E+06 IY= 8 1.966E+06 1.910E+06 IY= 6 1.771E+06 1.859E+06 IY= 4 1.646E+06 1.700E+06 IY= 2 1.555E+06 1.601E+06 IZ= 11 13 Field Values of PRPS IY= 10 empty pil prop pil prop pil prop pil prop IY= 8 pil prop pil prop pil prop pil prop pil prop IY= 6 pil prop pil prop pil prop pil prop pil prop IY= 4 pil prop pil prop pil prop pil prop pil prop IY= 2 empty empty pil prop pil prop pil prop IZ= 1 3 5 7 9 IY= 10 pil prop pil prop IY= 8 pil prop pil prop IY= 6 pil prop pil prop IY= 4 pil prop pil prop IY= 2 pil prop pil prop IZ= 11 13 Field Values of DEGF: degrees_Fahrenheit IY= 10 none 2.922E+03 2.949E+03 2.864E+03 2.749E+03 IY= 8 3.239E+03 2.424E+03 2.753E+03 2.078E+03 2.060E+03 IY= 6 3.313E+03 1.885E+03 2.614E+03 1.937E+03 1.911E+03 IY= 4 3.335E+03 1.478E+03 2.216E+03 1.865E+03 1.864E+03 IY= 2 none none 1.993E+03 1.870E+03 1.798E+03 IZ= 1 3 5 7 9 IY= 10 2.386E+03 2.079E+03 IY= 8 1.856E+03 1.794E+03 IY= 6 1.627E+03 1.735E+03 IY= 4 1.479E+03 1.546E+03 IY= 2 1.372E+03 1.429E+03 IZ= 11 13 Field Values of RRAT: rate_of_consumption_of_fuel_kg/m^3sec IY= 10 none 1.411E+00 2.815E+00 2.515E+00 4.802E-01 IY= 8 6.459E+00 1.878E+01 3.273E+00 7.714E-01 1.318E-01 IY= 6 4.695E-01 2.758E+01 5.087E-01 4.045E-01 7.327E-02 IY= 4 3.114E-01 4.116E+01 1.467E-01 2.442E-01 5.324E-02 IY= 2 none none 9.342E-02 1.556E-01 4.510E-02 IZ= 1 3 5 7 9 IY= 10 1.990E-01 6.576E-02 IY= 8 5.836E-02 2.027E-02 IY= 6 4.611E-02 1.759E-02 IY= 4 3.061E-02 1.051E-02 IY= 2 1.472E-02 7.229E-03 IZ= 11 13 Field Values of EPKE IY= 10 none 1.104E+02 3.220E+02 6.367E+02 2.377E+02 IY= 8 8.169E+02 8.748E+02 5.877E+02 5.342E+02 2.735E+02 IY= 6 6.935E+01 8.509E+02 3.087E+02 4.296E+02 2.279E+02 IY= 4 4.873E+01 9.818E+02 2.413E+02 3.796E+02 2.286E+02 IY= 2 none none 2.815E+02 3.824E+02 2.783E+02 IZ= 1 3 5 7 9 IY= 10 2.465E+02 2.174E+02 IY= 8 2.758E+02 1.935E+02 IY= 6 3.855E+02 2.306E+02 IY= 4 3.992E+02 2.375E+02 IY= 2 3.492E+02 3.166E+02 IZ= 11 13 Field Values of HPOR IY= 10 1.000E-11 1.0 1.0 1.0 1.0 IY= 8 8.770E-01 1.0 1.0 1.0 1.0 IY= 6 1.0 1.0 1.0 1.0 1.0 IY= 4 1.0 1.0 1.0 1.0 1.0 IY= 2 1.0 1.000E-11 1.0 1.0 1.0 IZ= 1 3 5 7 9 IY= 10 1.0 1.0 IY= 8 1.0 1.0 IY= 6 1.0 1.0 IY= 4 1.0 1.0 IY= 2 1.0 1.0 IZ= 11 13 Field Values of NPOR IY= 10 1.0 1.0 1.0 1.0 1.0 IY= 8 1.000E-11 1.0 1.0 1.0 1.0 IY= 6 1.0 1.0 1.0 1.0 1.0 IY= 4 1.0 1.0 1.0 1.0 1.0 IY= 2 1.000E-11 1.000E-11 1.0 1.0 1.0 IZ= 1 3 5 7 9 IY= 10 1.0 1.0 IY= 8 1.0 1.0 IY= 6 1.0 1.0 IY= 4 1.0 1.0 IY= 2 1.0 1.0 IZ= 11 13 Field Values of EPOR IY= 10 1.000E-11 1.0 1.0 1.0 1.0 IY= 8 1.850E-01 1.0 1.0 1.0 1.0 IY= 6 1.0 1.0 1.0 1.0 1.0 IY= 4 1.0 1.0 1.0 1.0 1.0 IY= 2 1.0 1.0 1.0 1.0 1.0 IZ= 1 3 5 7 9 IY= 10 1.0 1.0 IY= 8 1.0 1.0 IY= 6 1.0 1.0 IY= 4 1.0 1.0 IY= 2 1.0 1.0 IZ= 11 13 Field Values of VPOR IY= 10 1.000E-11 1.0 1.0 1.0 1.0 IY= 8 1.850E-01 1.0 1.0 1.0 1.0 IY= 6 1.0 1.0 1.0 1.0 1.0 IY= 4 1.0 1.0 1.0 1.0 1.0 IY= 2 1.000E-11 1.000E-11 1.0 1.0 1.0 IZ= 1 3 5 7 9 IY= 10 1.0 1.0 IY= 8 1.0 1.0 IY= 6 1.0 1.0 IY= 4 1.0 1.0 IY= 2 1.0 1.0 IZ= 11 13 Field Values of ENUT IY= 10 none 6.540E-04 1.907E-03 1.406E-02 1.408E-03 IY= 8 8.647E-02 4.606E-02 4.051E-02 2.032E-02 1.406E-02 IY= 6 1.519E-03 1.294E-02 3.267E-02 1.782E-02 1.938E-02 IY= 4 1.070E-03 3.922E-03 3.013E-02 2.095E-02 2.540E-02 IY= 2 none none 3.246E-02 3.981E-02 3.044E-02 IZ= 1 3 5 7 9 IY= 10 1.460E-03 1.287E-03 IY= 8 1.270E-02 1.200E-02 IY= 6 1.114E-02 1.088E-02 IY= 4 1.215E-02 1.616E-02 IY= 2 1.865E-02 2.953E-02 IZ= 11 13 Field Values of RHO1 IY= 10 none 1.427E+00 1.418E+00 1.461E+00 1.517E+00 IY= 8 1.305E+00 1.667E+00 1.511E+00 1.929E+00 1.944E+00 IY= 6 1.279E+00 2.039E+00 1.585E+00 2.044E+00 2.068E+00 IY= 4 1.272E+00 2.452E+00 1.828E+00 2.110E+00 2.111E+00 IY= 2 none none 1.998E+00 2.105E+00 2.173E+00 IZ= 1 3 5 7 9 IY= 10 1.716E+00 1.929E+00 IY= 8 2.117E+00 2.176E+00 IY= 6 2.352E+00 2.235E+00 IY= 4 2.534E+00 2.448E+00 IY= 2 2.683E+00 2.601E+00 IZ= 11 13 Field Values of TMP1: absolute_temperature_of_the_gas_Kelvin IY= 10 none 1.896E+03 1.911E+03 1.864E+03 1.800E+03 IY= 8 2.072E+03 1.620E+03 1.802E+03 1.427E+03 1.417E+03 IY= 6 2.113E+03 1.320E+03 1.725E+03 1.349E+03 1.335E+03 IY= 4 2.126E+03 1.094E+03 1.504E+03 1.309E+03 1.309E+03 IY= 2 none none 1.380E+03 1.312E+03 1.272E+03 IZ= 1 3 5 7 9 IY= 10 1.598E+03 1.428E+03 IY= 8 1.304E+03 1.270E+03 IY= 6 1.177E+03 1.237E+03 IY= 4 1.095E+03 1.132E+03 IY= 2 1.035E+03 1.067E+03 IZ= 11 13 Field Values of PROD IY= 10 none 6.272E-01 6.356E-01 6.093E-01 5.736E-01 IY= 8 7.255E-01 4.728E-01 5.748E-01 3.653E-01 3.598E-01 IY= 6 7.484E-01 3.055E-01 5.317E-01 3.217E-01 3.136E-01 IY= 4 7.553E-01 1.793E-01 4.083E-01 2.992E-01 2.990E-01 IY= 2 none none 3.390E-01 3.008E-01 2.787E-01 IZ= 1 3 5 7 9 IY= 10 4.609E-01 3.658E-01 IY= 8 2.965E-01 2.774E-01 IY= 6 2.255E-01 2.592E-01 IY= 4 1.795E-01 2.004E-01 IY= 2 1.464E-01 1.640E-01 IZ= 11 13 Field Values of OXID IY= 10 none 3.600E-01 3.556E-01 3.868E-01 4.244E-01 IY= 8 2.666E-01 5.057E-01 4.196E-01 6.332E-01 6.398E-01 IY= 6 2.448E-01 6.621E-01 4.666E-01 6.774E-01 6.861E-01 IY= 4 2.383E-01 7.788E-01 5.911E-01 7.002E-01 7.008E-01 IY= 2 none none 6.607E-01 6.987E-01 7.211E-01 IZ= 1 3 5 7 9 IY= 10 5.383E-01 6.339E-01 IY= 8 7.033E-01 7.225E-01 IY= 6 7.744E-01 7.408E-01 IY= 4 8.204E-01 7.995E-01 IY= 2 8.536E-01 8.359E-01 IZ= 11 13 Field Values of FUEL IY= 10 none 1.278E-02 8.741E-03 3.950E-03 2.020E-03 IY= 8 7.906E-03 2.147E-02 5.568E-03 1.444E-03 4.820E-04 IY= 6 6.770E-03 3.242E-02 1.648E-03 9.417E-04 3.215E-04 IY= 4 6.391E-03 4.192E-02 6.078E-04 6.433E-04 2.329E-04 IY= 2 none none 3.319E-04 4.068E-04 1.620E-04 IZ= 1 3 5 7 9 IY= 10 8.071E-04 3.026E-04 IY= 8 2.116E-04 1.047E-04 IY= 6 1.196E-04 7.628E-05 IY= 4 7.667E-05 4.424E-05 IY= 2 4.216E-05 2.283E-05 IZ= 11 13 Field Values of MIXF IY= 10 none 4.717E-02 4.359E-02 3.735E-02 3.347E-02 IY= 8 4.768E-02 4.739E-02 3.708E-02 2.147E-02 2.021E-02 IY= 6 4.780E-02 4.916E-02 3.080E-02 1.858E-02 1.751E-02 IY= 4 4.780E-02 5.175E-02 2.299E-02 1.705E-02 1.663E-02 IY= 2 none none 1.892E-02 1.690E-02 1.544E-02 IZ= 1 3 5 7 9 IY= 10 2.608E-02 2.036E-02 IY= 8 1.646E-02 1.532E-02 IY= 6 1.248E-02 1.429E-02 IY= 4 9.919E-03 1.103E-02 IY= 2 8.066E-03 9.016E-03 IZ= 11 13 ************************************************************ Whole-field residuals before solution at sweep 100 with resref values determined by EARTH & resfac=1.0E-05 variable resref (res sum)/resref (res sum) P1 1.837E-01 1.830E-03 3.360E-04 U1 1.273E+01 8.869E-04 1.129E-02 V1 1.273E+01 1.701E-04 2.166E-03 W1 1.273E+01 1.142E-04 1.454E-03 KE 9.955E+01 2.994E-04 2.981E-02 EP 1.240E+05 2.641E-04 3.274E+01 H1 3.852E+05 4.252E-05 1.638E+01 FUEL 1.146E-03 1.189E-04 1.362E-07 MIXF 3.515E-03 9.489E-05 3.336E-07 max abs corrections for solved-for variables for comparison with maximum & minimum values variable max-cor high low ( IX IY IZ) P1 1.959E-01 1.022E+03 -3.838E+02 ( 6, 1, 10) U1 7.009E-02 8.423E+00 -8.193E+00 ( 1, 1, 10) V1 1.196E-02 9.072E+01 -3.307E+01 ( 1, 9, 2) W1 8.753E-03 3.153E+01 -1.692E+01 ( 6, 1, 12) KE 8.340E-01 1.290E+03 1.596E-01 ( 1, 1, 13) EP 7.706E+02 1.469E+06 4.646E+00 ( 6, 1, 13) H1 1.327E+02 3.771E+06 1.287E+06 ( 1, 4, 4) FUEL 2.384E-06 4.833E-02 1.875E-05 ( 6, 5, 2) MIXF 2.708E-06 5.330E-02 2.609E-03 ( 1, 4, 4) SARAH = Self-Adjusting Relxtn. algoritHm DTFALS values set in EARTH by SARAH with relaxfac=1. U1 dtf=3.308664E-03 V1 dtf=3.308664E-03 W1 dtf=3.308664E-03 ************************************************************ Sources and sinks !! Zero nett sources are not printed !!! Nett Sources have units of mass_per_unit_time * variable Average values have units of the variable Typically the units of the sources are: U1,V1,W1 - Force - Newtons R1 - Mass - kg/s TEM1,H1 - Energy - Watts Nett source of U1 at patch named: OUTLET =-1.532227E-02 Nett source of U1 at patch named: WALL1 =-6.730359E-05 Nett source of U1 at patch named: WALL2 = 3.874044E-04 Nett source of U1 at patch named: WALL3 = 1.171143E-03 Nett source of U1 at patch named: WALL4-SW =-6.336369E-08 Nett source of U1 at patch named: WALL5 = 1.910003E-08 Nett source of U1 at patch named: WALL6-LW = 3.885600E-09 pos. sum=1.55857E-03 neg. sum=-0.01539 nett sum=-0.013831 Nett source of V1 at patch named: FOIN = 1.177630E+00 Nett source of V1 at patch named: SOIN1 =-3.262800E-01 Nett source of V1 at patch named: SOIN2 =-3.262800E-01 Nett source of V1 at patch named: DILUIN =-5.800000E-01 Nett source of V1 at patch named: OUTLET =-6.896576E-02 Nett source of V1 at patch named: WALL5 = 1.565584E-04 Nett source of V1 at patch named: WALL6-LW =-1.600011E-05 pos. sum=1.177786 neg. sum=-1.301542 nett sum=-0.123756 Nett source of W1 at patch named: OUTLET =-6.930737E-01 Nett source of W1 at patch named: WALL1 =-5.544975E-03 Nett source of W1 at patch named: WALL2 =-2.598631E-03 Nett source of W1 at patch named: WALL3 =-3.169282E-03 Nett source of W1 at patch named: WALL4-SW =-7.044840E-05 pos. sum=0. neg. sum=-0.704457 nett sum=-0.704457 Nett source of R1 at patch named: FOIN = 1.177630E-02 Nett source of R1 at patch named: SOIN1 = 8.157000E-03 Nett source of R1 at patch named: SOIN2 = 8.157000E-03 Nett source of R1 at patch named: DILUIN = 1.450000E-02 Nett source of R1 at patch named: OUTLET =-4.258938E-02 (Mass Out -4.803603E-02 In 5.446650E-03) pos. sum=0.04259 neg. sum=-0.042589 nett sum=9.164214E-07 Nett source of KE at patch named: FOIN = 2.944074E-03 (Average 2.500000E-01) Nett source of KE at patch named: SOIN1 = 1.003311E-03 (Average 1.230000E-01) Nett source of KE at patch named: SOIN2 = 1.003311E-03 (Average 1.230000E-01) Nett source of KE at patch named: DILUIN = 1.783500E-03 (Average 1.230000E-01) Nett source of KE at patch named: OUTLET =-1.140169E+00 (Ave Out 3.924890E+01 In 1.367566E+02) Nett source of KE at patch named: WALL1 = 2.719957E-09 Nett source of KE at patch named: WALL2 = 1.194936E-09 Nett source of KE at patch named: WALL3 = 1.838449E-09 Nett source of KE at patch named: WALL4-SW = 1.532919E-07 Nett source of KE at patch named: WALL5 =-5.163132E-08 Nett source of KE at patch named: WALL6-LW = 2.322345E-13 Nett source of KE at patch named: KESOURCE =-3.484712E+01 Non-linearised source for KE at: KESOURCE = 4.949973E+01 pos. sum=49.506458 neg. sum=-35.987286 nett sum=13.519173 Nett source of EP at patch named: FOIN = 6.918881E-01 (Average 5.875261E+01) Nett source of EP at patch named: SOIN1 = 6.875874E-02 (Average 8.429416E+00) Nett source of EP at patch named: SOIN2 = 6.875874E-02 (Average 8.429416E+00) Nett source of EP at patch named: DILUIN = 9.134528E-02 (Average 6.299675E+00) Nett source of EP at patch named: OUTLET =-2.680562E+02 (Ave Out 1.146409E+04 In 5.185394E+04) Nett source of EP at patch named: WALL1 = 1.403250E-06 Nett source of EP at patch named: WALL2 = 4.918402E-07 Nett source of EP at patch named: WALL3 = 7.374116E-07 Nett source of EP at patch named: WALL4-SW = 2.746960E-06 Nett source of EP at patch named: WALL5 =-9.032331E-07 Nett source of EP at patch named: WALL6-LW = 1.015301E-11 Nett source of EP at patch named: KESOURCE =-5.169827E+04 Non-linearised source for EP at: KESOURCE = 6.259446E+04 pos. sum=6.259538E+04 neg. sum=-5.196632E+04 nett sum=1.062905E+04 Nett source of H1 at patch named: FOIN = 4.529050E+04 (Average 3.845904E+06) Nett source of H1 at patch named: SOIN1 = 9.458041E+03 (Average 1.159500E+06) Nett source of H1 at patch named: SOIN2 = 9.458041E+03 (Average 1.159500E+06) Nett source of H1 at patch named: DILUIN = 1.681275E+04 (Average 1.159500E+06) Nett source of H1 at patch named: OUTLET =-8.102221E+04 (Ave Out 1.868434E+06 In 1.603186E+06) pos. sum=8.101934E+04 neg. sum=-8.102221E+04 nett sum=-2.875 Nett source of FUEL at patch named: FOIN = 6.456302E-04 (Average 5.482456E-02) Nett source of FUEL at patch named: OUTLET =-5.679269E-06 (Ave Out 1.206720E-04 In 2.161034E-05) Nett source of FUEL at patch named: WHOLE =-6.399860E-04 pos. sum=6.456302E-04 neg. sum=-6.456652E-04 nett sum=-3.504101E-08 Nett source of MIXF at patch named: FOIN = 6.456302E-04 (Average 5.482456E-02) Nett source of MIXF at patch named: OUTLET =-6.457104E-04 (Ave Out 1.446805E-02 In 9.054812E-03) pos. sum=6.456302E-04 neg. sum=-6.457104E-04 nett sum=-8.021016E-08 ************************************************************ Summary of sources from all Objects/Patches at sweep 100 Variable Inflow Outflow Nett P1 4.259029E-02 -4.258938E-02 9.164214E-07 U1 1.558570E-03 -1.538964E-02 -1.383107E-02 V1 1.177786E+00 -1.301542E+00 -1.237557E-01 W1 0.000000E+00 -7.044571E-01 -7.044571E-01 KE 4.950646E+01 -3.598729E+01 1.351917E+01 EP 6.259538E+04 -5.196632E+04 1.062905E+04 H1 8.101934E+04 -8.102221E+04 -2.875000E+00 FUEL 6.456302E-04 -6.456652E-04 -3.504101E-08 MIXF 6.456302E-04 -6.457104E-04 -8.021016E-08 ************************************************************ ************************************************************ spot values vs sweep or iteration number IXMON= 3 IYMON= 5 IZMON= 5 TIMESTEP= 1 Tabulation of abscissa and ordinates... ISWP P1 U1 V1 W1 KE 2 1.0000E-10 -6.0145E-04 9.1288E-01 7.8115E+00 6.8399E+01 7 -1.8484E+02 -1.8609E-02 3.3235E+00 -8.7431E-01 3.7902E+02 12 2.6677E+02 -1.4159E-03 2.3429E+00 -7.0781E-01 5.0303E+02 17 3.8250E+02 3.1274E-03 2.5609E+00 3.7452E-01 5.7661E+02 22 3.0377E+02 -2.4034E-02 2.6277E+00 -1.3393E-01 5.3634E+02 27 1.6231E+02 -6.9485E-02 3.3024E-01 -4.7835E+00 3.8281E+02 32 2.6402E+01 -1.3722E-01 -1.3004E-01 -9.5565E+00 2.4775E+02 37 -6.5102E+01 -2.1936E-01 1.2088E+00 -1.1088E+01 1.7520E+02 42 -1.0014E+02 -2.9860E-01 1.9217E+00 -1.1185E+01 1.3421E+02 47 -9.0141E+01 -3.2422E-01 2.3670E+00 -1.1094E+01 1.1838E+02 52 -7.0842E+01 -3.3762E-01 2.6731E+00 -1.1063E+01 1.0634E+02 57 -6.5606E+01 -3.4799E-01 2.8676E+00 -1.1075E+01 9.7709E+01 62 -6.6378E+01 -3.5335E-01 2.9672E+00 -1.1079E+01 9.3155E+01 67 -6.6303E+01 -3.5665E-01 2.9992E+00 -1.1087E+01 9.1054E+01 72 -6.5361E+01 -3.5834E-01 3.0244E+00 -1.1086E+01 8.9899E+01 77 -6.4357E+01 -3.5952E-01 3.0493E+00 -1.1083E+01 8.9100E+01 82 -6.3352E+01 -3.6047E-01 3.0692E+00 -1.1084E+01 8.8469E+01 87 -6.2544E+01 -3.6125E-01 3.0822E+00 -1.1086E+01 8.8001E+01 92 -6.2173E+01 -3.6173E-01 3.0896E+00 -1.1087E+01 8.7708E+01 97 -6.2050E+01 -3.6196E-01 3.0939E+00 -1.1087E+01 8.7547E+01 ISWP EP H1 FUEL MIXF 2 3.6344E+03 2.5235E+06 2.5391E-02 5.4825E-02 7 4.9449E+04 2.5187E+06 2.1840E-02 5.1603E-02 12 1.1186E+05 2.5161E+06 1.7484E-02 4.8345E-02 17 2.1753E+05 2.7194E+06 1.4144E-02 4.7402E-02 22 2.7305E+05 3.0084E+06 1.0318E-02 4.7271E-02 27 1.9197E+05 3.1353E+06 6.2367E-03 4.4806E-02 32 9.3780E+04 3.0522E+06 3.3393E-03 4.0122E-02 37 4.9503E+04 2.8547E+06 2.1287E-03 3.5014E-02 42 3.5945E+04 2.7014E+06 1.5676E-03 3.1588E-02 47 3.0400E+04 2.6217E+06 1.3575E-03 2.9876E-02 52 2.7040E+04 2.5718E+06 1.1935E-03 2.8833E-02 57 2.4681E+04 2.5409E+06 1.0878E-03 2.8195E-02 62 2.3323E+04 2.5223E+06 1.0354E-03 2.7814E-02 67 2.2688E+04 2.5130E+06 1.0136E-03 2.7622E-02 72 2.2347E+04 2.5070E+06 1.0022E-03 2.7500E-02 77 2.2108E+04 2.5029E+06 9.9475E-04 2.7416E-02 82 2.1918E+04 2.4998E+06 9.8874E-04 2.7353E-02 87 2.1778E+04 2.4977E+06 9.8446E-04 2.7311E-02 92 2.1689E+04 2.4966E+06 9.8215E-04 2.7289E-02 97 2.1638E+04 2.4960E+06 9.8109E-04 2.7276E-02 Variable 1 = P1 2 = U1 3 = V1 4 = W1 5 = KE Minval= -1.848E+02 -3.620E-01 -1.300E-01 -1.118E+01 6.840E+01 Maxval= 3.825E+02 3.127E-03 3.323E+00 7.812E+00 5.766E+02 Cellav= 2.628E+00 -2.325E-01 2.369E+00 -7.603E+00 2.020E+02 Variable 6 = EP 7 = H1 8 = FUEL 9 = MIXF Minval= 3.634E+03 2.496E+06 9.811E-04 2.728E-02 Maxval= 2.730E+05 3.135E+06 2.539E-02 5.482E-02 Cellav= 6.432E+04 2.640E+06 5.703E-03 3.537E-02 1.00 9..3.2..5.+6.7.+....+....+....+....+....+....+....+ . 2 2 3 3 3 3 0.90 + 9 5 3 3 3 3 + . 8 5 1 7 3 . 0.80 + 1 6 7 2 3 + . 9 9 9 . 0.70 + 8 6 3 + . 9 . 0.60 + 5 4 4 5 2 3 + . 4 4 8 7 . 0.50 + + . 9 . 0.40 + 6 8 3 + 1 7 4 6 . 0.30 3 9 7 + . . 0.20 + 8 5 7 1 1 1 1 1 1 1 1 1 1 . 6 3 6 9 1 . 0.10 + 8 6 9 7 6 + 7 7 7 8 9 9 7 7 6 6 6 6 6 6 0.00 6..1.+....+....+3.4.+8..8+8..8+.9.9+.9.9+.9..9.9..9 0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 the abscissa is ISWP. min= 2.00E+00 max= 9.70E+01 ************************************************************ ************************************************************ residuals vs sweep or iteration number Tabulation of abscissa and ordinates... ISWP P1 U1 V1 W1 KE 2 2.5134E+00 4.8309E-01 9.6465E-01 5.8502E-01 1.3826E+00 7 4.8977E-01 8.3765E-02 1.1859E-01 1.4450E-01 3.5410E-01 12 3.1541E-01 2.6530E-02 6.2704E-02 5.9395E-02 2.5169E-01 17 3.5625E-01 2.8698E-02 7.6136E-02 6.3623E-02 1.9848E-01 22 2.5227E-01 1.9514E-02 7.7383E-02 4.3322E-02 2.0167E-01 27 1.9140E-01 1.3922E-02 5.7320E-02 3.3826E-02 1.4969E-01 32 1.2314E-01 1.0534E-02 3.1752E-02 2.4000E-02 1.1002E-01 37 6.3463E-02 7.7605E-03 1.5687E-02 1.3153E-02 6.1520E-02 42 4.2425E-02 6.1118E-03 1.0825E-02 8.4759E-03 3.2670E-02 47 3.4363E-02 5.1858E-03 8.9368E-03 5.3174E-03 2.8353E-02 52 2.7723E-02 4.2326E-03 6.5319E-03 3.3286E-03 2.1691E-02 57 1.8717E-02 3.3639E-03 4.0875E-03 2.1325E-03 1.2491E-02 62 1.2948E-02 2.6555E-03 2.2894E-03 1.5258E-03 6.6611E-03 67 9.1903E-03 2.1402E-03 1.3248E-03 1.1213E-03 3.4568E-03 72 6.2290E-03 1.7992E-03 8.3996E-04 6.4523E-04 3.1925E-03 77 5.3209E-03 1.5409E-03 7.9553E-04 5.2368E-04 2.3287E-03 82 4.3347E-03 1.3474E-03 6.0217E-04 3.9356E-04 1.7158E-03 87 3.4089E-03 1.1851E-03 4.1921E-04 2.8916E-04 1.1376E-03 92 2.6771E-03 1.0437E-03 2.8537E-04 2.1354E-04 7.1128E-04 97 2.0303E-03 9.2707E-04 1.9593E-04 1.3662E-04 3.6481E-04 ISWP EP H1 FUEL MIXF 2 1.3855E+00 8.1325E-01 2.3431E+00 1.4889E+00 7 5.1516E-01 1.6145E-01 4.7991E-01 1.6556E-01 12 3.7591E-01 9.0658E-02 1.9346E-01 1.0252E-01 17 2.3377E-01 9.6627E-02 2.3286E-01 8.2960E-02 22 2.0712E-01 8.0901E-02 1.5480E-01 7.0680E-02 27 1.6261E-01 4.5421E-02 9.7026E-02 4.9694E-02 32 1.0875E-01 1.8611E-02 5.3382E-02 3.5859E-02 37 5.5057E-02 1.1424E-02 2.0333E-02 2.8749E-02 42 2.7625E-02 9.9958E-03 1.5624E-02 2.3510E-02 47 2.5247E-02 6.3496E-03 1.1697E-02 1.4443E-02 52 2.2808E-02 3.3579E-03 8.3200E-03 7.5386E-03 57 1.5055E-02 1.8103E-03 4.9456E-03 4.0424E-03 62 8.5829E-03 1.1212E-03 3.1001E-03 2.4996E-03 67 3.6298E-03 7.6357E-04 1.8215E-03 1.7011E-03 72 2.0316E-03 4.9389E-04 1.0276E-03 1.1008E-03 77 2.0409E-03 2.9075E-04 1.1891E-03 6.4830E-04 82 1.6750E-03 1.8585E-04 8.3388E-04 4.1475E-04 87 1.1043E-03 1.1354E-04 4.8384E-04 2.5356E-04 92 6.0735E-04 6.9930E-05 2.3402E-04 1.5612E-04 97 3.2055E-04 4.9089E-05 1.2511E-04 1.0958E-04 Variable 1 = P1 2 = U1 3 = V1 4 = W1 5 = KE Minval= -6.200E+00 -6.983E+00 -8.538E+00 -8.898E+00 -7.916E+00 Maxval= 9.216E-01 -7.276E-01 -3.599E-02 -5.361E-01 3.240E-01 Variable 6 = EP 7 = H1 8 = FUEL 9 = MIXF Minval= -8.045E+00 -9.922E+00 -8.986E+00 -9.119E+00 Maxval= 3.260E-01 -2.067E-01 8.515E-01 3.980E-01 1.00 9....+....+....+....+....+....+....+....+....+....+ . . 0.90 + 6 + . 8 6 . 0.80 + 5 7 + . 9 8 8 7 6 . 0.70 + 2 9 9 9 8 6 + . 9 . 0.60 + 9 9 + . 2 2 7 9 5 . 0.50 + 2 8 8 9 6 + . 2 3 8 9 6 . 0.40 + 2 1 4 9 6 + . 2 1 8 9 . 0.30 + 2 2 1 7 9 + . 2 1 8 9 8 . 0.20 + 2 3 8 9 8 + . 2 2 3 4 9 8 . 0.10 + 2 2 1 9 6 + . 2 2 9 . 0.00 +....+....+....+....+....+....+....+....+....+.2..9 0 .1 .2 .3 .4 .5 .6 .7 .8 .9 1.0 the abscissa is ISWP. min= 2.00E+00 max= 9.70E+01 ************************************************************ ************************************************************ PATCH(IZ4 ,CONTUR, 1, 6, 1, 10, 4, 4, 1, 1) PLOT(IZ4 ,TMP1, 1.000E+00, 1.000E+01) Contour plot of TMP1 Total=0.119741 Cell ave=1574.50415 Mass ave=1661.960693 Min valu=1365.696655 Max valu=1885.597046 PHOENICS polar plot: Angle X,radius Y Zero angle is at 3 o-clock 30 I::::::::::::::::::::::::::::::::::: .......I 31 I:::::::::::::::::::::::::::::::: 888 .....I 32 I::::::::::::::::::::::::::::666 888 ....I 33 I::::::::::::::::::::::::: 6666 888 ...I 34 I:::::::::::::::::::::44444 6666 888 ...I 35 I:::::::::::::::::: 4444 6666 888 ..I 36 I:::::::::::::: 2222 4444 666 88 .I 37 I::::::::::: 2222 4444 6666 888 .I 38 I:::::::00000 222 4444 666 888 .I 39 I.:::000000000 2222 4444 666 88 .I 40 I...::::::::::::::::::::::::::::::::::::::::::I 40 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 ************************************************************ PATCH(XSECIN1 ,CONTUR, 1, 1, 1, 10, 1, 13, 1, 1) PLOT(XSECIN1 ,W1 , 0.000E+00, 1.000E+01) Contour plot of W1 X-plane; Y horizontal, Z vertical Total=568.320007 Cell ave=4.371692 Min valu=-16.922176 Max valu=30.496031 N,S,E,W... denote: North, South, East .. 0 HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH 1 S N 2 S 4444444444444444444444444444444444444444N 3 S 44444444444444444444444444444444444444 N 4 S444444 66N 5 S4444 66666666 6666666666666N 6 S4 666666666666 6666666666666N 7 S 666666666666666 666666666N 8 S66666666666666666 666666N 9 S 66666666 666N 10 S 666666 N 11 S 66666666 N 12 S 66666666 4444444444444444N 13 S 666666666 4444444444444444444N 14 S66666666666 44444444444444444444 N 15 S 4444444444444444444N 16 S44444444 444444444444444444N 17 S 44444444 44444444444444444N 18 S22 444444 444444444444444N 19 S 444444 44444444444444N 20 S444444444 44444444444 N 21 S 444444444 N 22 S 4444444 N 23 S 444444444 N 24 S 4444444444N 25 S 4444444444N 26 S 444444444N 27 S4444444444444444444 444444N 28 S 4444444444 44N 29 S22222222222222 44444444 N 30 S000000 22222 4444444 666666N 31 S00000000000 2222 44444 6666666666N 32 S0000000000000 22222 444 6666 88N 33 S0000000000000000 2222 44 666 888 N 34 S00000000000000000 222 44 666 888 N 35 S0000000000000000 222 44 666 888 N 36 S222222222222222222222222 4444 66666 N 37 S 4444444 N 38 S 4444444444444 222222 444N 39 S 444 444 222222222222 N 40 LLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL 40 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 ************************************************************ PATCH(XSECIN1 ,CONTUR, 1, 1, 1, 10, 1, 13, 1, 1) PLOT(XSECIN1 ,TMP1, 0.000E+00, 1.000E+01) Contour plot of TMP1 X-plane; Y horizontal, Z vertical Total=1.821687E+05 Cell ave=1401.297485 Min valu=773. Max valu=2215.284912 N,S,E,W... denote: North, South, East .. 0 HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH 1 S222222222222222222 444444444444 N 2 S 2222222222222 444444444444 N 3 S 22222222222 444444444444 N 4 S 22222222222 44444444444 N 5 S 222222222 44444444444 N 6 S 2222222222 4444444444 N 7 S 22222222222 44444444 N 8 S 222222222222 4444444 N 9 S 22222222222222 444444 N 10 S 222222222222222222 44444N 11 S 2222222222222222222222 N 12 S 222222222N 13 S N 14 S N 15 S2222222222222222222 2222222222N 16 S222222222222222222222222222222222222 N 17 S 4444444 N 18 S 44444444 666N 19 S 44444444 66N 20 S 44444444 66N 21 S 44444444 6N 22 S 4444444 6N 23 S 444444 N 24 S 4N 25 S 22222222222222N 26 S 22222222222222222 N 27 S 222222N 28 S 444444444444444444444444444444444N 29 S 444444444444 6666666666 N 30 S 44444444444444 66666 N 31 S44444444444444 6666 888888N 32 S44444444444444 666666 8N 33 S44444444444444 66666666 N 34 S444444444444444 66666666666666 N 35 S44444444444444444 6666666666 N 36 S 222222 4444444 66666 N 37 S000000 2 444 666666666 8888N 38 S0000002468 888888888888888888888N 39 S0000002468 88888888 6 444444N 40 LLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL 40 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 ************************************************************ PATCH(XSECIN1 ,CONTUR, 1, 1, 1, 10, 1, 13, 1, 1) PLOT(XSECIN1 ,FUEL, 0.000E+00, 1.000E+01) Contour plot of FUEL X-plane; Y horizontal, Z vertical Total=0.454367 Cell ave=3.495128E-03 Min valu=0. Max valu=0.048325 N,S,E,W... denote: North, South, East .. 0 HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH 1 S000000000000000000000000000000000000000000000N 2 S000000000000000000000000000000000000000000000N 3 S000000000000000000000000000000000000000000000N 4 S000000000000000000000000000000000000000000000N 5 S000000000000000000000000000000000000000000000N 6 S000000000000000000000000000000000000000000000N 7 S000000000000000000000000000000000000000000000N 8 S000000000000000000000000000000000000000000000N 9 S000000000000000000000000000000000000000000000N 10 S000000000000000000000000000000000000000000000N 11 S000000000000000000000000000000000000000000000N 12 S000000000000000000000000000000000000000000000N 13 S000000000000000000000000000000000000000000000N 14 S000000000000000000000000000000000000000000000N 15 S000000000000000000000000000000000000000000000N 16 S000000000000000000000000000000000000000000000N 17 S000000000000000000000000000000000000000000000N 18 S000000000000000000000000000000000000000000000N 19 S000000000000000000000000000000000000000000000N 20 S000000000000000000000000000000000000000000000N 21 S000000000000000000000000000000000000000000000N 22 S000000000000000000000000000000000000000000000N 23 S000000000000000000000000000000000000000000000N 24 S000000000000000000000000000000000000000000000N 25 S000000000000000000000000000000000000000000000N 26 S000000000000000000000000000000000000000000000N 27 S000000000000000000000000000000000000000000000N 28 S000000000000000000000000000000000000000000000N 29 S00000000000000000000000000000000000000000000 N 30 S00000000000000000000000000000000000000000 N 31 S000000000000000000000000000000000000000 N 32 S00000000000000000000000000000000000 N 33 S000000000000000000000000000000000 N 34 S000000000000000000000000000000 N 35 S00000000000000000000000000 22N 36 S000000 24 6666 444444444 2222222N 37 S000000 24 66 4444444444 22222222N 38 S0000000000 00000 2222 N 39 S0000000000 000000N 40 LLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL 40 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 ************************************************************ PATCH(XSECIN1 ,CONTUR, 1, 1, 1, 10, 1, 13, 1, 1) PLOT(XSECIN1 ,MIXF, 0.000E+00, 1.000E+01) Contour plot of MIXF X-plane; Y horizontal, Z vertical Total=2.954735 Cell ave=0.022729 Min valu=0. Max valu=0.0533 N,S,E,W... denote: North, South, East .. 0 HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH 1 S 2222222222 4444N 2 S 2222222222 44444N 3 S 222222222 44444N 4 S 222222222 444444N 5 S 2222222222 444444N 6 S 2222222222 44444N 7 S 2222222222222 44444N 8 S 22222222222222 4444N 9 S 2222222222222222 4444N 10 S 22222222222222222 44N 11 S 2222222222222 N 12 S 222N 13 S N 14 S 00000000000000000N 15 S 222N 16 S22222222222222222222222222222222222222222 N 17 S22222222222222222222222222222222 4444N 18 S222222222 44444 N 19 S2222222 444444 N 20 S 44444 N 21 S 4444 N 22 S 44444 N 23 S 4444N 24 S 22222222222222222222 N 25 S 2222222222222222222222222222222N 26 S 2222222222222222222222222 N 27 S 22222222222N 28 S 44444444444444444444444N 29 S 4444444 N 30 S 444444 6666666666666666 N 31 S 444444 66666666 88N 32 S 444444 666666666 88N 33 S 4444444 6666666666 88N 34 S 4444444 66666666 88N 35 S 4444444 6666666 888N 36 S 24 6 8888888888888888888N 37 S0000002468 8888888888888888N 38 S0000002 888888888 88888888888888888N 39 S0000002 88888888888888 88888888888 6 444444N 40 LLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL 40 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 ************************************************************ PATCH(XSECIN2 ,CONTUR, 4, 4, 1, 10, 1, 13, 1, 1) PLOT(XSECIN2 ,W1 , 0.000E+00, 1.000E+01) Contour plot of W1 X-plane; Y horizontal, Z vertical Total=748.668701 Cell ave=5.75899 Min valu=-16.923037 Max valu=30.728359 N,S,E,W... denote: North, South, East .. 0 HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH 1 S N 2 S 44444444444444444444444444444444444444444N 3 S4444444444 4444444444444444 N 4 S4444 666N 5 S44 6666666666666666666666666666666666666 N 6 S 66666666666666666666666666666666666666 N 7 S 666666666666666666666666666666666666666 N 8 S66666666666666666666666666666666666666666 N 9 S 66666666666666666666666666666 N 10 S 66666666666666666666666666666 N 11 S 666666666666666666666666666666 N 12 S 666666666666666666666666666666 N 13 S 66666666666666666666666666666 N 14 S6666666666666666666666666666666666 N 15 S 6666666666666666666 N 16 S444444444 66666666666666666 N 17 S 444444444 6666666666666666 N 18 S222 444444 6666666666666666 N 19 S 44444444 666666666666666666 N 20 S44444444444 6666666666666666666N 21 S 6666N 22 S 6666666 66N 23 S N 24 S N 25 S 444N 26 S 44444N 27 S4444444444444444444 4444N 28 S 444444444 N 29 S2222222222222222 444444 N 30 S0000000 22222 444444 666666 N 31 S0000000000000 2222 4444 6666666666N 32 S000000000000000 2222 444 6666 8N 33 S000000000000000000 2222 44 666 88 N 34 S000000000000000000 222 44 66 888 N 35 S00000000000000000 222 44 666 8888 N 36 S222222222222222222222222 4444 66666 N 37 S 4444444 N 38 S 4444444444444 2222222 44N 39 S 444 444 222222222222 N 40 LLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL 40 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 ************************************************************ PATCH(XSECIN2 ,CONTUR, 4, 4, 1, 10, 1, 13, 1, 1) PLOT(XSECIN2 ,TMP1, 0.000E+00, 1.000E+01) Contour plot of TMP1 X-plane; Y horizontal, Z vertical Total=1.838872E+05 Cell ave=1414.516968 Min valu=773. Max valu=2214.753418 N,S,E,W... denote: North, South, East .. 0 HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH 1 S222222222222222222222222222 4N 2 S 2222222222222222222222 44N 3 S 22222222222222222222 444N 4 S 22222222222222222222 444N 5 S 222222222222222222 444N 6 S 222222222222222222 4444N 7 S 222222222222222222 4444N 8 S 222222222222222222 4444 N 9 S 222222222222222222 4444 N 10 S 2222222222222222 4444 N 11 S 2222222222222222 4444 N 12 S 22222222222222 44444 N 13 S 22222222222222 44444 N 14 S 222222222222222 444444 6N 15 S2222222222222222222222 444444 6N 16 S222222222222222 44444444 66N 17 S 444444444 66N 18 S 4444444444444 666N 19 S 4444444444444 666N 20 S 44444444444444 666N 21 S 4444444444444444 666N 22 S 4444444444444444 666N 23 S 444444444444444444 N 24 S 44444N 25 S N 26 S 22222222222N 27 S 44444444444444444 2N 28 S 4444444444444 44444444444N 29 S 44444444444 666666666666666666666N 30 S 44444444444444 66666 N 31 S44444444444444 66666 888888N 32 S444444444444444 66666 8N 33 S444444444444444 666666666 N 34 S4444444444444444 6666666666666 N 35 S44444444444444444 6666666666 N 36 S 222222 4444444 66666 N 37 S000000 2 444 666666666 8888N 38 S0000002468 888888888888888888888N 39 S0000002468 88888888 6 444444N 40 LLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL 40 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 ************************************************************ PATCH(XSECIN2 ,CONTUR, 4, 4, 1, 10, 1, 13, 1, 1) PLOT(XSECIN2 ,FUEL, 0.000E+00, 1.000E+01) Contour plot of FUEL X-plane; Y horizontal, Z vertical Total=0.455718 Cell ave=3.505524E-03 Min valu=0. Max valu=0.048324 N,S,E,W... denote: North, South, East .. 0 HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH 1 S000000000000000000000000000000000000000000000N 2 S000000000000000000000000000000000000000000000N 3 S000000000000000000000000000000000000000000000N 4 S000000000000000000000000000000000000000000000N 5 S000000000000000000000000000000000000000000000N 6 S000000000000000000000000000000000000000000000N 7 S000000000000000000000000000000000000000000000N 8 S000000000000000000000000000000000000000000000N 9 S000000000000000000000000000000000000000000000N 10 S000000000000000000000000000000000000000000000N 11 S000000000000000000000000000000000000000000000N 12 S000000000000000000000000000000000000000000000N 13 S000000000000000000000000000000000000000000000N 14 S000000000000000000000000000000000000000000000N 15 S000000000000000000000000000000000000000000000N 16 S000000000000000000000000000000000000000000000N 17 S000000000000000000000000000000000000000000000N 18 S000000000000000000000000000000000000000000000N 19 S000000000000000000000000000000000000000000000N 20 S000000000000000000000000000000000000000000000N 21 S000000000000000000000000000000000000000000000N 22 S000000000000000000000000000000000000000000000N 23 S000000000000000000000000000000000000000000000N 24 S000000000000000000000000000000000000000000000N 25 S000000000000000000000000000000000000000000000N 26 S000000000000000000000000000000000000000000000N 27 S000000000000000000000000000000000000000000000N 28 S000000000000000000000000000000000000000000000N 29 S00000000000000000000000000000000000000000000 N 30 S00000000000000000000000000000000000000000 N 31 S000000000000000000000000000000000000000 N 32 S000000000000000000000000000000000000 N 33 S000000000000000000000000000000000 N 34 S000000000000000000000000000000 N 35 S00000000000000000000000000 22N 36 S000000 24 6666 444444444 2222222N 37 S000000 24 66 4444444444 22222222N 38 S0000000000 00000 2222 N 39 S0000000000 000000N 40 LLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL 40 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 ************************************************************ PATCH(XSECIN2 ,CONTUR, 4, 4, 1, 10, 1, 13, 1, 1) PLOT(XSECIN2 ,MIXF, 0.000E+00, 1.000E+01) Contour plot of MIXF X-plane; Y horizontal, Z vertical Total=3.008694 Cell ave=0.023144 Min valu=0. Max valu=0.053297 N,S,E,W... denote: North, South, East .. 0 HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH 1 S 22222222222222222222222222 N 2 S 2222222222222222222222222 N 3 S 22222222222222222222222 N 4 S 22222222222222222222222 N 5 S 2222222222222222222222 N 6 S 222222222222222222222 4N 7 S 22222222222222222222 4N 8 S 222222222222222222 44N 9 S 22222222222222222 44N 10 S 2222222222222222 444N 11 S 222222222222222 4444N 12 S 22222222222222 4444N 13 S 2222222222222 4444 N 14 S 22222222222222 4444 N 15 S 2222222222222222222 4444 N 16 S22222222222222222222222222 4444 N 17 S222222222222222222222 44444 N 18 S222222222 444444 N 19 S222222 44444 N 20 S 444444 N 21 S 444444 N 22 S 444444 N 23 S 444444 N 24 S 4N 25 S 2222 N 26 S 22222222222222N 27 S 2222N 28 S 444444444444444444444444444N 29 S 4444444 N 30 S 444444 6666666666666666 N 31 S 4444444 66666666 88N 32 S 4444444 666666666 88N 33 S 4444444 666666666 88N 34 S 4444444 66666666 88N 35 S 4444444 6666666 888N 36 S 24 6 8888888888888888888N 37 S0000002468 8888888888888888N 38 S0000002 888888888 88888888888888888N 39 S0000002 88888888888888888888888888 6 444444N 40 LLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL 40 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 ************************************************************ ************************************************************ SATLIT RUN NUMBER = 1 ; LIBRARY REF.= 492 Run started at 12:04:40 on Wednesday, 17 April 2019 Run completed at 12:04:41 on Wednesday, 17 April 2019 CPU time of run 1 s This includes 1 seconds of user time and 0 seconds of system time. TIME/(VARIABLES*CELLS*TSTEPS*SWEEPS*ITS) = 1.425E-06 ************************************************************