```
TALK=T;RUN( 1, 1)

************************************************************
Q1 created by VDI menu, Version 3.6, Date 12/09/05
CPVNAM=VDI;SPPNAM=Core
Echo DISPLAY / USE settings
DISPLAY
tests symmetrical solution.

|-----------------|
|\\\\\\\\\\\\\ T=1|
|  \\\\\\ \\\\\\\ |
|    \\\\ k2 \\\\\|
|      \\\\ \\\\\\|
|        \\\\\\\\\|
|          \\\\\\\|
|            \\\\\|
|      k1      \\\|
Y |T=0             \|
| |-----------------|
+-- X

The solution should be lineary for k1=k2.
The 161 number property is the solid with incompressible
air properties.
ENDDIS
PHOTON USE
p;;;;

vi z;up y
msg Temperature contours
con TEM1 z 1 fil;.001
msg Press e to END
use patgeo
ENDUSE
************************************************************
IRUNN   =         1 ;LIBREF =       709
************************************************************
REAL(XLEN,K1,K2)
INTEGER(NCEL)
NCEL=3
mesg(NX number of cells is :ncel: OK? If not, insert uneven desired number,
ncel
XLEN=3.0
K1=1.0;K2=K1 ! heat conductivity
Group 1. Run Title
TEXT(2D PARSOL heat-conduction, xy plane
************************************************************
Group 2. Transience
************************************************************
Groups 3, 4, 5  Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,NCEL,NCEL,1)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS  =    T
* Non-default variable names
NAME(148) =TEM1 ; NAME(150) =PRPS
* Solved variables list
SOLVE(TEM1)
* Stored variables list
STORE(PRPS,KOND,DVAR,DRES,DCOR)
STORE(S_EX,D_EX) ! exact solution at sunny and dark
STORE(SE/S,DE/D) ! exact/newparsol_solution for sunny and dark
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
* Domain material is: 0   Air at 20 deg C, 1 atm, treated as i
SETPRPS(1,  0)
ENUT    = 0.000000E+00
DVO1DT  = 3.410000E-03
PRNDTL(TEM1) = -K1
CSG10='Q1'           ! signal use of the following properties line
! which correspond to steel
MATFLG=T;NMAT=1
<161>   artificial solid material having same properties as material 0
161      1.189      1.544E-05   1005.0   1.  3.41E-3   0.0
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(PRPS) = -1.000000E+00
FIINIT(DVAR) = -1.234E-11
FIINIT(DRES) = -1.234E-11
FIINIT(DCOR) = -1.234E-11
FIINIT(S_EX) = 0.; FIINIT(D_EX) = 0.
FIINIT(SE/S) = 0.; FIINIT(DE/D) = 0.
No PATCHes used for this Group
************************************************************
Group 12. Convection and diffusion adjustments
No PATCHes used for this Group
************************************************************
Group 13. Boundary & Special Sources

PATCH(HOT,VOLUME,NX,NX,NY,NY,NZ,NZ,1,1)
COVAL(HOT,TEM1,1.E6,1.0)
PATCH(COLD,VOLUME,1,1,1,1,1,1,1,1)
COVAL(COLD,TEM1,1.E6,0.0)

inform7begin
(stored var rst1 is resi(tem1))
(stored var crt1 is corr(tem1))
(stored s_ex at hot is 1.0)
patch(struct1,cell,2,2,1,1,1,nz,1,1)
(stored s_ex at struct1 is .325)
patch(struct2,cell,1,1,2,2,1,nz,1,1)
(stored s_ex at struct2 is .325)
patch(struct3,cell,nx-1,nx-1,ny,ny,1,nz,1,1)
(stored s_ex at struct3 is .675)
patch(struct4,cell,nx-1,nx-1,ny,ny,1,nz,1,1)
(stored s_ex at struct4 is .675)
patch(cut1,cell,nx,nx,1,1,1,nz,1,1)
(stored s_ex at cut1 is .4708)
(stored d_ex at cut1 is .5292)
patch(cut2,cell,1,1,ny,ny,1,nz,1,1)
(stored s_ex at cut2 is .4708)
(stored d_ex at cut2 is .5292)
patch(cut3,cell,nx/2+1,nx/2+1,ny/2+1,ny/2+1,1,nz,1,1)
(stored s_ex at cut3 is .45)
(stored d_ex at cut3 is .55)

(stored var SE/S is s_ex/tem1 with if(tem1!=0.))
(stored var DE/D is d_ex/dvar with if(dvar!=0.))
inform7end
EGWF    =    T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP  = 100
************************************************************
Group 16. Terminate Iterations
LITER (TEM1) = 100
RESREF(TEM1)=0.
ENDIT(TEM1)=0.

************************************************************
Group 17. Relaxation
relax(tem1,linrlx,1.)
************************************************************
Group 18. Limits
************************************************************
Group 19. EARTH Calls To GROUND Station
USEGRD  =    T  ;USEGRX =    T
ASAP    =    T
PARSOL  =    T
************************************************************
Group 20. Preliminary Printout
ECHO    =    T
DISTIL  =    T  ;NULLPR =    F
NDST    =         0
DSTTOL  = 1.000000E-02
************************************************************
Group 21. Print-out of Variables
************************************************************
Group 22. Monitor Print-Out
IXMON   =         1 ;IYMON  =         1 ;IZMON  =        -3
NPRMON  =    100000
NPRMNT  =         1
TSTSWP  =        -1
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT  =    100000
NXPRIN  =    1      ; NYPRIN  =    1
ISWPRF  =         1 ;ISWPRL =    100000
No PATCHes used for this Group
patch(plot,profil,1,nx,1,1,1,1,1,1)
coval(plot,tem1,0,0)
************************************************************
Group 24. Dumps For Restarts

GVIEW(P,0.000000E+00,0.000000E+00,1.000000E+00)
GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)

> DOM,    SIZE,        XLEN, XLEN, 1.000000E+00
> DOM,    MONIT,       5.000000E-01, 5.000000E-01, 5.000000E-01
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
> GRID,   RSET_X_1,      NCEL, 1.000000E+00
> GRID,   RSET_Y_1,      NCEL, 1.000000E+00
> GRID,   RSET_Z_1,      1, 1.000000E+00

> OBJ,    NAME,        BLOCK
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        XLEN, XLEN, 1.000000E+00
> OBJ,    GEOMETRY,    wedge1
> OBJ,    ROTATION24,        8
> OBJ,    GRID,        NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    161, air
> OBJ,    INI_TEMP,     0.000000E+00
STOP
isg60=1      ! set active NewParsol
```