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