TALK=F;RUN( 1, 1)
 
 ************************************************************
   Q1 created by VDI menu, Version 2009, Date 17/02/09
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
  Group 1. Run Title
 TEXT(USP labirinth with conjugate heat transfer.                          )
 ************************************************************
title
  DISPLAY
  
  This case solves a two-dimensional steady hydrodynamics
  problem with conjugate heat transfer in labirinth channel.

  The analytical solution is absent. The solution can be compared 
  with solution of Structured PHOENICS.
  
  Cases changes inlet velocity as Uin=0.1*caseno
  
  Domain material is air.
  Upper material is epoxy. Lower material is gold.
              |------------------------------------|
              |                                    |
              |------|           EPOXY             |
   In Flow -> |----| |                             |
              |    | |                             |
              |    |  -------|                     |
              |    |-------| |                     |
              |            | |                     |
              |            | |                     |
              |            | ---------|            |
              |            |--------| |            |
              |                     | |            |
              |                     | |            |
              |     GOLD            | ---------|   |
              |                     |--------| |   |
              |                              | |   |
              |                              | |   |
              |                              | |   |
              |------------------------------------|                             
  
  ENDDIS
    
REAL(Uin)
INTEGER(caseno)
caseno=1
Uin=0.1*caseno
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
boolean(AGG)
AGG=F
mesg(Do you want to use AGG (y) or CGG (n)? (y/n)
readvdu(ans,char,n)
if(:ans:.eq.y)then
AGG=T
endif
if(AGG)then
 RSET(M,10,10,1)
else
 INTEGER(NCUBE,NCHAN)
 NCUBE = 80
 NCHAN = NCUBE/20
 RSET(M,NCUBE,NCUBE,1)
endif
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
 NAME(148) = PRPS;  NAME(147) = KOND
 NAME(149) = TEM1
    * Solved variables list
 SOLVE(P1  ,U1  ,V1, TEM1)
    * Stored variables list
 STORE(PRPS,KOND)
    * Additional solver options
 SOLUTN(P1  ,Y,Y,Y,N,N,Y)
 SOLUTN(TEM1  ,Y,Y,Y,N,N,Y)
 
 ************************************************************
  Group 8. Terms & Devices
 USP = T ; USPDBG= F
 UTCPLT = F ; USPIMB= F
 UAUTO = AGG ; UGTGRD= F
 USPVTK = T ; USPTER= F
 CELLST  =        10 ;FACEST =        10
if(AGG)then
 USCANX = T ; USCANY= T
 USCANZ = T
 MXLEV = 0 ;MYLEV = 0 ;MZLEV = 0
 DOMAT = -1
 MINPRP = -1 ; MAXPRP= 100
 UMXLEV = 4 ; UMXCLS= 10000
 UNOLAY = 2 ; UIOLEV= 0
 UWALEV = 0 ; UGRITR= 2
 UMXCRS = 16
else
 MXLEV   =         4 ;MYLEV  =         4 ;MZLEV  =         4
 DOMAT   =        0
 MINPRP  =        -1 ;MAXPRP =       250
endif
 ************************************************************
  Group 9. Properties
 PRESS0 =1.0E+05 ;TEMP0 =273.
    * Domain material index is   0 signifying:
    * Air at 20 deg C, 1 atm, treated as incompressible
 SETPRPS(1,  0)
 ENUT =0.
 DVO1DT =3.41E-03
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(PRPS)=-1.
   No PATCHes used for this Group
 if(AGG)then
else
 PATCH(B1,INIVAL, 4*NX/5+NCHAN+1,NX,1,NY,1,NZ,1,1)
 COVAL(B1, PRPS, 0.0, 104)

 PATCH(B2,INIVAL, 3*NX/5+NCHAN+1,4*NX/5+NCHAN,NY/5+NCHAN+1,NY,1,NZ,1,1)
 COVAL(B2, PRPS, 0.0, 104)
 PATCH(B3,INIVAL, 2*NX/5+NCHAN+1,3*NX/5+NCHAN,2*NY/5+NCHAN+1,NY,1,NZ,1,1)
 COVAL(B3, PRPS, 0.0, 104)
 PATCH(B4,INIVAL, NX/5+NCHAN+1,2*NX/5+NCHAN,3*NY/5+NCHAN+1,NY,1,NZ,1,1)
 COVAL(B4, PRPS, 0.0, 104)
 PATCH(B5,INIVAL, 1,NX/5+NCHAN,4*NY/5+NCHAN+1,NY,1,NZ,1,1)
 COVAL(B5, PRPS, 0.0, 104)
 PATCH(B6,INIVAL, 1,4*NX/5,1,NY/5,1,NZ,1,1)
 COVAL(B6, PRPS, 0.0, 107) 
 PATCH(B7,INIVAL, 1,3*NX/5,NY/5+1,2*NY/5,1,NZ,1,1)
 COVAL(B7, PRPS, 0.0, 107) 
 PATCH(B8,INIVAL, 1,2*NX/5,2*NY/5+1,3*NY/5,1,NZ,1,1)
 COVAL(B8, PRPS, 0.0, 107) 
 PATCH(B9,INIVAL, 1,NX/5,3*NY/5+1,4*NY/5,1,NZ,1,1)
 COVAL(B9, PRPS, 0.0, 107) 
endif 
 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
if(AGG)then
 PATCH(AA11,CELL, 1,8,1,1,1,1,1,1)
 COVAL(AA11,TEM1, FIXVAL, 100.) 
 PATCH(AA12,CELL, 1,NX,NY,NY,1,1,1,1)
 COVAL(AA12,TEM1, FIXVAL, 0.) 
else
 PATCH(AA10,NORTH, 1,NX,NY,NY,1,NZ,1,1)
 COVAL(AA10,TEM1, FIXVAL, 0.) 
 PATCH(AA11,SOUTH, 1,4*NX/5,1,1,1,NZ,1,1)
 COVAL(AA11,TEM1, FIXVAL, 100.) 
endif
 
 EGWF    =    T
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =       3000
 RESFAC  = 1.000000E-07
 ************************************************************
  Group 16. Terminate Iterations
 ************************************************************
  Group 17. Relaxation
 RELAX(P1 ,LINRLX,0.7 )
 RELAX(U1  ,FALSDT, 1.000000E-01)
 RELAX(V1  ,FALSDT, 1.000000E-01)
mesg(Do you want to use collocated arrangement (y) or staggered one (n)? (y/n)
readvdu(ans,char,n)
if(:ans:.eq.y)then
 SPEDAT(SET,USP,METHOD,I,1)
 RELAX(P1  ,LINRLX, 3.000000E-01)
 RELAX(U1  ,FALSDT, 5.000000E-02)
 RELAX(V1  ,FALSDT, 5.000000E-02)
mesg(Do you want to use SIMPLEST (y) or SIMPLE (n)? (y/n)
readvdu(ans,char,n)
if(:ans:.eq.y)then
 SPEDAT(SET,USP,SIMPLEST,L,T)
endif

endif


 ************************************************************
  Group 18. Limits
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 PARSOL = F
 ISG62 = 0
 ISG60 = 1
 SPEDAT(SET,USP,METHOD,I,1)
 SPEDAT(SET,USP,SIMPLEST,L,T)
 SPEDAT(SET,USPIO,VERTCENT,L,F)
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 SPEDAT(SET,SOLPHI,P1,I,24)
 SPEDAT(SET,SOLPHI,U1,I,0)
 SPEDAT(SET,PREPHI,U1,I,0)
 SPEDAT(SET,SOLPHI,V1,I,0)
 SPEDAT(SET,PREPHI,V1,I,0)
 SPEDAT(SET,SOLPHI,TEM1,I,24)
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  Group 21. Print-out of Variables
 ************************************************************
  Group 22. Monitor Print-Out
if(AGG)then
 IXMON = 5 ;IYMON = 4 ;IZMON = 1
else
 IXMON   =        50 ;IYMON  =        42 ;IZMON  =         1
endif
 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

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


 ************************************************************
  Group 24. Dumps For Restarts
if(AGG)then 
 GVIEW(P,-0.042817,0.054706,0.997584)
 GVIEW(UP,0.017703,0.998385,-0.05399)
 
> DOM,    SIZE,        1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    MONIT,       4.750000E-01, 3.750000E-01, 5.000000E-01
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
  > GRID,   RSET_X_1,     10, 1.000000E+00
  > GRID,   RSET_Y_1,     10, 1.000000E+00
  > GRID,   RSET_Z_1,      1, 1.000000E+00
 
> OBJ,    NAME,        IN
> OBJ,    POSITION,    0.000000E+00, 8.000000E-01, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 5.000000E-02, 1.000000E+00
> OBJ,    GEOMETRY,    cube3t
> OBJ,    TYPE,        INLET
> OBJ,    GRID,        NO
> OBJ,    PRESSURE,    0.
> OBJ,    VELOCITY,    0.1 ,0. ,0.
 
> OBJ,    NAME,        OUT
> OBJ,    POSITION,    8.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        5.000000E-02, 0.000000E+00, 1.000000E+00
> OBJ,    GEOMETRY,    cube12t
> OBJ,    TYPE,        OUTLET
> OBJ,    GRID,        NO
> OBJ,    PRESSURE,    0.
> OBJ,    COEFFICIENT, 1000.
 
> OBJ,    NAME,        B3
> OBJ,    POSITION,    0.000000E+00, 8.500000E-01, 0.000000E+00
> OBJ,    SIZE,        TO_END,       TO_END,       TO_END
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    cube14
> OBJ,    GRID,        NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    104,Epoxy
 
> OBJ,    NAME,        B4
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        2.000000E-01, 8.000000E-01, 1.000000E+00
> OBJ,    GEOMETRY,    cube14
> OBJ,    GRID,        NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    107, Gold
 
> OBJ,    NAME,        B5
> OBJ,    POSITION,    2.500000E-01, 6.000000E-01, 0.000000E+00
> OBJ,    SIZE,        TO_END,       2.500000E-01, 1.000000E+00
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    cube14
> OBJ,    GRID,        NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    104,Epoxy
 
> OBJ,    NAME,        B6
> OBJ,    POSITION,    2.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        2.500000E-01, 5.500000E-01, TO_END
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    cube14
> OBJ,    GRID,        NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    107, Gold
 
> OBJ,    NAME,        B7
> OBJ,    POSITION,    5.000000E-01, 4.000000E-01, 0.000000E+00
> OBJ,    SIZE,        TO_END,       2.000000E-01, 1.000000E+00
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    cube14
> OBJ,    GRID,        NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    104,Epoxy
 
> OBJ,    NAME,        B8
> OBJ,    POSITION,    4.500000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        2.000000E-01, 3.500000E-01, 1.000000E+00
> OBJ,    GEOMETRY,    cube14
> OBJ,    GRID,        NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    107, Gold
 
> OBJ,    NAME,        B9
> OBJ,    POSITION,    7.000000E-01, 2.000000E-01, 0.000000E+00
> OBJ,    SIZE,        TO_END,       2.000000E-01, 1.000000E+00
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    cube14
> OBJ,    GRID,        NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    104,Epoxy
 
> OBJ,    NAME,        B10
> OBJ,    POSITION,    6.500000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        1.500000E-01, 1.500000E-01, 1.000000E+00
> OBJ,    GEOMETRY,    cube14
> OBJ,    GRID,        NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    107, Gold
 
> OBJ,    NAME,        B11
> OBJ,    POSITION,    8.500000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        TO_END,       2.000000E-01, 1.000000E+00
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    cube14
> OBJ,    GRID,        NO
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    104,Epoxy
else
 GVIEW(P,0.000000E+00,0.000000E+00,-1.000000E+00)
 GVIEW(UP,0.000000E+00,1.000000E+00,0.000000E+00)
 
> DOM,    SIZE,        1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    MONIT,       4.950000E-01, 4.150000E-01, 5.000000E-02
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
 
> OBJ,    NAME,        IN
> OBJ,    POSITION,    0.000000E+00, 8.000000E-01, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 5.000000E-02, 1.000000E+00
> OBJ,    GEOMETRY,    cube3t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        INLET
> OBJ,    PRESSURE,     0.000000E+00
> OBJ,    VELOCITY,     1.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    TEMPERATURE,  0.000000E+00
 
> OBJ,    NAME,        OUT
> OBJ,    POSITION,    8.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        5.000000E-02, 0.000000E+00, 1.000000E+00
> OBJ,    GEOMETRY,    cube12t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        OUTLET
> OBJ,    PRESSURE,     0.000000E+00
> OBJ,    TEMPERATURE,  SAME
> OBJ,    COEFFICIENT,  1.000000E+03

endif
STOP