TALK=F;RUN( 1, 1)
 
 ************************************************************
   Q1 created by VDI menu, Version 2007, Date 26/09/07
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
 IRUNN   =         1 ;LIBREF =         0
 ************************************************************
  Group 1. Run Title
 TEXT(Body with large Conductivity in uniform field
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 INTEGER(LEVEL)
 LEVEL = 4
 mesg(Refinement level = :LEVEL:. Enter other if you want
 readvdu(LEVEL, int, :LEVEL:)

 RSET(M,2,1,2)
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Non-default variable names
 NAME(147) =TEM1
    * Solved variables list
 SOLVE(TEM1)
    * Stored variables list
    * Additional solver options
 SOLUTN(TEM1,Y,Y,Y,N,N,Y)
 
 ************************************************************
  Group 8. Terms & Devices
 USP     =    T  ; USPDBG=    F
 UTCPLT  =    F  ; USPIMB=    F
 UAUTO   =    T  ; UGTGRD=    F
 USPVTK  =    T  ; USPTER=    F
 USCANX  =    T  ; USCANY=    F
 USCANZ  =    T
 MXLEV   =         0 ;MYLEV  =         0 ;MZLEV  =         0
 DOMAT   =        -1
 CELLST  =        10 ; FACEST=        10
 MINPRP  =        -1 ; MAXPRP=       250
 UMXLEV = LEVEL ; UMXCLS= 10000
 UNOLAY  =         2 ; UIOLEV=         1
 UWALEV = 0 ; UGRITR= 2
 UMXCRS  =        16
 USPREL  = 7.000000E-01 ; UMNCEL= 1.000000E-04
 UDVREL  = 1.000000E+00
PARSOL = T

mesg(PARSOL = T. Do you change it? (y/n)
readvdu(ans,char,n)
if(:ans:.eq.y)then
 PARSOL = F
 mesg(PARSOl = F
endif

mesg(Do you want to view results in the centres of cells? (y/n)
readvdu(ans,char,n)
if(:ans:.eq.y)then
SPEDAT(SET,USPIO,VERTCENT,L,F)
endif


INTEGER(caseno)
caseno = 0
mesg(Body is cube. Do you wand to change shape of body? (y/n)
readvdu(ans,char,n)
if(:ans:.eq.y)then
 mesg(caseno = 0: cube
 mesg(caseno = 1: cylinder
 mesg(caseno = 2: tallwedge
 mesg(caseno = 3: wedge
 mesg(caseno = ? (0/1/2/3/4)
 readvdu(caseno,int,:caseno:)
endif


 ************************************************************
  Group 9. Properties
 PRESS0  = 1.000000E+05 ;TEMP0  = 2.730000E+02
    * Domain material index is   0 signifying:
    * Air at 20 deg C, 1 atm, treated as incompressible
 SETPRPS(1,  0)
 ENUT    = 0.000000E+00
 DVO1DT  = 3.410000E-03
 PRNDTL(TEM1) = -0.1
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
   No PATCHes used for this Group
 
 
 INIADD  =    F
 ************************************************************
  Group 12. Convection and diffusion adjustments
   No PATCHes used for this Group
 ************************************************************
  Group 13. Boundary & Special Sources
   No PATCHes used for this Group
 
 EGWF    =    T
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =       300
 RESFAC  = 1.000000E-07
 ************************************************************
  Group 16. Terminate Iterations
 LITER(TEM1) = 100
 ENDIT(TEM1) = 1.E-7
 ************************************************************
  Group 17. Relaxation
            RELAX(TEM1,FALSDT,1000.)
 ************************************************************
  Group 18. Limits
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  Group 21. Print-out of Variables
 ************************************************************
  Group 22. Monitor Print-Out
 NPRMON  =    100000
 NPRMNT  =         1
 TSTSWP  =        -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT  =    100000
 ISWPRF  =         1 ;ISWPRL =    100000
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 
 GVIEW(P,9.983340E-02,-9.949917E-01,-4.974976E-03)
 GVIEW(UP,0.000000E+00,-4.999955E-03,9.999875E-01)
 
> DOM,    SIZE,        1.000000E+00, 1.000000E+00, 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
> GRID,   AUTO,         F T F
  > GRID,   RSET_X_1,      2, 1.000000E+00
  > GRID,   RSET_Y_1,      1, 1.000000E+00
  > GRID,   RSET_Z_1,      2, 1.000000E+00
 
> OBJ,    NAME,        BODY
> OBJ,    POSITION,    1.000000E-01, 0.000000E+00, 3.000000E-01
> OBJ,    SIZE,        8.000000E-01, 1.000000E+00, 4.000000E-01
if(caseno.eq.0)then
> OBJ,    GEOMETRY,    cube14
> OBJ,    ROTATION24,        1
endif
if(caseno.eq.1)then
> OBJ,    GEOMETRY,    cylinder
> OBJ,    ROTATION24,        5
endif
if(caseno.eq.2)then
> OBJ,    GEOMETRY,    tallwedge
> OBJ,    ROTATION24,        5
endif
if(caseno.eq.3)then
> OBJ,    GEOMETRY,    wedge
> OBJ,    ROTATION24,        5
endif
> OBJ,    GRID,        NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    110,SILVER at 27 deg C
 
> OBJ,    NAME,        PL1
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 1.000000E+00, 1.000000E+00
> OBJ,    GEOMETRY,    cube13
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE
> OBJ,    SURF_TEMP,    0.000000E+00, 0.000000E+00
 
> OBJ,    NAME,        PL2
> OBJ,    POSITION,    1.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 1.000000E+00, 1.000000E+00
> OBJ,    GEOMETRY,    cube13
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE
> OBJ,    SURF_TEMP,    0.000000E+00, 1.000000E+02
STOP