talk=f;run(1,1) DISPLAY Steady state YZ test ENDDIS ************************************************************ Group 1. Run Title TEXT(YZ TEST: I129 ) ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,1,12,10) * Set overall domain extent: * xulast yvlast zwlast name XSI= 1.000E-01;YSI= 1.200E-01;ZSI= 1.000E+00;RSET(D,CHAM ) * Set objects: x0 y0 z0 * dx dy dz name XPO= 0.000E+00;YPO= 1.000E-02;ZPO= 0.000E+00 XSI= 1.000E-01;YSI= 1.100E-01;ZSI= 0.000E+00;RSET(B,B1 ) XPO= 0.000E+00;YPO= 0.000E+00;ZPO= 1.000E+00 XSI= 1.000E-01;YSI= 1.100E-01;ZSI= 0.000E+00;RSET(B,B2 ) XPO= 0.000E+00;YPO= 0.000E+00;ZPO= 0.000E+00 XSI= 1.000E-01;YSI= 0.000E+00;ZSI= 1.000E+00;RSET(B,B3 ) XPO= 0.000E+00;YPO= 1.200E-01;ZPO= 0.000E+00 XSI= 1.000E-01;YSI= 0.000E+00;ZSI= 1.000E+00;RSET(B,B4 ) XPO= 0.000E+00;YPO= 0.000E+00;ZPO= 0.000E+00 XSI= 1.000E-01;YSI= 1.000E-02;ZSI= 5.000E-01;RSET(B,B5 ) XPO= 0.000E+00;YPO= 1.100E-01;ZPO= 5.000E-01 XSI= 1.000E-01;YSI= 1.000E-02;ZSI= 5.000E-01;RSET(B,B6 ) ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd ONEPHS = T * Non-default variable names NAME( 44) =TRAD ; NAME( 45) =TRES NAME( 46) =ABSV ; NAME( 47) =BLOK NAME( 48) =TEM1 ; NAME( 49) =PRPS NAME( 50) =DEN1 * Solved variables list SOLVE(P1 ,V1 ,W1 ,KE ,EP ,TEM1) * Stored variables list STORE(DEN1,PRPS,BLOK,ABSV,TRES,TRAD) * Additional solver options SOLUTN(P1 ,Y,Y,Y,N,N,N) SOLUTN(KE ,Y,Y,N,N,N,N) SOLUTN(EP ,Y,Y,N,N,N,N) SOLUTN(TEM1,Y,Y,Y,N,N,Y) IVARBK = -1 ;ISOLBK = 1 ************************************************************ Group 8. Terms & Devices TERMS (KE ,N,Y,Y,Y,Y,N) TERMS (EP ,N,Y,Y,Y,Y,N) NEWRH1 = T NEWENL = T NEWENT = T ISOLX = 0 ;ISOLY = 0 ;ISOLZ = 0 ************************************************************ Group 9. Properties PRESS0 = 1.000E+05 TEMP0 = 2.730E+02 EL1 = KE15DEP ENUT = PRKOLM setprps(1,2) PRT (EP ) = 1.314E+00 ************************************************************ Group 11.Initialise Var/Porosity Fields FIINIT(KE ) = 2.500E-01 ;FIINIT(EP ) = 2.050E+00 FIINIT(BLOK) = 1.000E+00 ;FIINIT(TEM1) = 2.000E+01 CONPOR(B5 , -1.00,CELL ,-#1,-#1,-#1,-#1,-#1,-#1) INIT(B5 ,BLOK, 0.000E+00, 2.000E+00) INIT(B5 ,PRPS, 0.000E+00, 1.000E+02) CONPOR(B6 , -1.00,CELL ,-#1,-#1,-#3,-#3,-#2,-#2) INIT(B6 ,BLOK, 0.000E+00, 3.000E+00) INIT(B6 ,PRPS, 0.000E+00, 1.000E+02) INIADD = F ************************************************************ Group 13. Boundary & Special Sources PATCH (B2 ,HIGH ,#1,#1,#1,#2,#2,#2,1,1) COVAL (B2 ,P1 , 1.000E+00, 0.000E+00) COVAL (B2 ,KE , 0.000E+00, SAME ) COVAL (B2 ,EP , 0.000E+00, SAME ) COVAL (B2 ,TEM1, 0.000E+00, 2.000E+01) INLET (B1 ,LOW ,#1,#1,#2,#3,#1,#1,1,1) VALUE (B1 ,P1 , 1.189E+01) VALUE (B1 ,W1 , 1.000E+01) VALUE (B1 ,KE , 2.500E-01) VALUE (B1 ,EP , 2.050E+00) VALUE (B1 ,TEM1, 2.000E+01) PATCH (HOTB3 ,SWALL ,#1,#1,#1,#1,#1,#2,1,1) COVAL (HOTB3 ,TEM1, LOGLAW , 3.000E+01) PATCH (HOTB4 ,NWALL ,#1,#1,#3,#3,#1,#2,1,1) COVAL (HOTB4 ,TEM1, LOGLAW , 1.000E+02) PATCH (KESOURCE,PHASEM,#1,#NREGX,#1,#NREGY,#1,#NREGZ,#1,#NREGT) COVAL (KESOURCE,KE , GRND4 , GRND4 ) COVAL (KESOURCE,EP , GRND4 , GRND4 ) PATCH (BUOYANCY,PHASEM,#1,#NREGX,#1,#NREGY,#1,#NREGZ,#1,#NREGT) COVAL (BUOYANCY,V1 , FIXFLU , DENSDIFF ) BUOYA = 0.000E+00 ; BUOYB =-9.810E+00 ; BUOYC = 0.000E+00 BUOYD = 1.189E+00 EGWF = T ************************************************************ Group 15. Terminate Sweeps LSWEEP = 100; SELREF = T; RESFAC = 1.000E-02 ************************************************************ Group 16. Terminate Iterations LITER (TRAD) = 50 ;LITER (TEM1) = 15 ENDIT (TRAD) = GRND1 ;ENDIT (TEM1) = GRND1 ************************************************************ Group 17. Relaxation RELAX(P1 ,LINRLX, 1.000E+00) RELAX(V1 ,FALSDT, 1.000E-03) RELAX(W1 ,FALSDT, 1.000E-03) RELAX(KE ,FALSDT, 1.000E-03) RELAX(EP ,FALSDT, 1.000E-03) KELIN = 1 ************************************************************ Group 19. EARTH Calls To GROUND Station NAMGRD =FLAR GENK = T ************************************************************ Group 22. Monitor Print-Out IXMON = 1 ;IYMON = 5 ;IZMON = 5 TSTSWP = -1 ************************************************************ Group 23.Field Print-Out & Plot Control NYPRIN = 2; NPLT = 1 ************************************************************ Group 24. Dumps For Restarts spedat(set,flair,trad,r,25) STOP