TALK=F;RUN( 1, 1)
 CPVNAM=VDI;SPPNAM=Core

 TEXT(USP. Temperature in uniform flow
 STEADY  =    T
 RSET(M,32,32,1)
 RSET(D,DOM,1.,1.,1.E-2)
 
  Group 7. Variables: STOREd,SOLVEd,NAMEd
    * Non-default variable names
 NAME(150) =TEM1
    * Solved variables list
 SOLVE(TEM1)
    * Stored variables list
 STORE(U1,V1,W1)
    * Additional solver options
 SOLUTN(TEM1,Y,Y,Y,N,N,Y)

    GROUP 8. Terms (in differential equations) & devices
INTEGER(LEVEL)
LEVEL = 3
mesg(Level of coarsening is :LEVEL:
mesg(Do you want to change it? (y/n)
readvdu(ans,char,n)
if(:ans:.eq.y)then
 mesg(Enter level of coarsening (0,1,2):
 readvdu(LEVEL,int,LEVEL)
 mesg(New level of coarsening is :LEVEL:
endif

USP    = T
UAUTO  = F
USPDBG = F
UTCPLT = T
USPVTK = T
USPIMB = F
MXLEV  = LEVEL 
MYLEV  = LEVEL
MZLEV  = LEVEL
DOMAT  = -1
MINPRP = -1
MAXPRP = 100
CELLST = 10
FACEST = 10
PARSOL = F
 
  Group 9. Properties
REAL(FLAM,PE)
FLAM = 0.02
 ENUL=1.E-1
 RHO1=1.
 CP1 = 1.
 PRNDTL(TEM1) = -FLAM

REAL(Uin)
Uin = 0.1
mesg(Default value of uniform velocity is :Uin: m/s
mesg(Do you want to change it? (y/n)
readvdu(ans,char,n)
if(:ans:.eq.y)then
 mesg(Enter value of uniform velocity
 readvdu(Uin,real,Uin)
 mesg(New value of uniform velocity is :Uin: m/s
endif
PE = Uin/FLAM
mesg(Pecle's number is :PE:

FIINIT(U1) = Uin
FIINIT(V1) = 0.
FIINIT(W1) = 0.
 
PATCH(LEFT,WEST,1,1,1,NY,1,NZ,1,1)
COVAL(LEFT,TEM1,1.,0.)

PATCH(RIGHT,EAST,NX,NX,1,NY,1,NZ,1,1)
COVAL(RIGHT,TEM1,1.,100.)
 
  SAVE7BEGIN
(STORED of EXT1 is 100.*(exp(PE*XG)-1.)/(exp(PE)-1.))
(STORED of DTEM is EXT1-TEM1) 
  SAVE7END

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


mesg(Do you want to solve task without calculation of gradients ? (y/n)
readvdu(ans,char,n)
if(:ans:.eq.y)then
SPEDAT(SET,USP,CALCGRAD,L,F)
endif

 
  Group 15. Terminate Sweeps
  
 LSWEEP  =       100
 RESFAC  = 1.000000E-09
 
  Group 17. Relaxation
 RELAX(TEM1,LINRLX, 1.000000E+00)
 
  Group 20. Preliminary Printout
 ECHO    =    T
 
  Group 22. Monitor Print-Out
 IXMON   =         3 ;IYMON  =         3 ;IZMON  =         1
 NPRMON  =    100000
 NPRMNT  =         1
 TSTSWP  =        -1

  Group 23.Field Print-Out & Plot Control
 NPRINT  =    100000
 ISWPRF  =         1 ;ISWPRL =    100000

> DOM,    SIZE,        1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    MONIT,       5.000000E-01, 0.500000E+00, 0.500000E+00
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02

if(LEVEL.gt.0)then
> OBJ,    NAME,        B1
> OBJ,    POSITION,    0.500000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.500000E+00, 1.000000E+00, 1.000000E+00
> OBJ,    GEOMETRY,    fine
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        FINE_GRID_VOL
> OBJ,    LEVELS,       0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    AXES,        OXY
endif

if(LEVEL.gt.1)then
> OBJ,    NAME,        B2
> OBJ,    POSITION,    6.250000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        3.750000E-01, 1.000000E+00, 1.000000E+00
> OBJ,    GEOMETRY,    fine
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        FINE_GRID_VOL
> OBJ,    LEVELS,       1.000000E+00, 1.000000E+00, 1.000000E+00
> OBJ,    AXES,        OXY
endif

if(LEVEL.gt.2)then
> OBJ,    NAME,        B3
> OBJ,    POSITION,    7.500000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        2.500000E-01, 1.000000E+00, 1.000000E+00
> OBJ,    GEOMETRY,    fine
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        FINE_GRID_VOL
> OBJ,    LEVELS,       2.000000E+00, 2.000000E+00, 2.000000E+00
> OBJ,    AXES,        OXY
endif

STOP