PHOTON USE
  p;;;;;
  
  vi y; up z
  msg flow is from left to right
  msg contours of horizontal-direction velocity
  con u1 y 21 fi;0.001; pause
  msg contours of vertical-direction velocity
  con w1 y 21 fi;0.001; pause
  msg contours of pressure
  con p1 y 21 fi;0.001

  ENDUSE
 TEXT(F-1 Car; Balsa Wood Blank; Pot. Flow        
#cls
TITLE 
  DISPLAY
  
  The flow around a model Formula-1 racing-car is calculated by
  way of the velocity-potential method, which gives a solution
  rapidly, but which is realistic only for the forward part of 
  the vehicle.
  
  SOLVE(POT) and the setting of POT [ i.e. the velocity potential ]
  to constant values at upstream and downstream boundaries of the 
  domain activates the solution.
  
  A single sweep suffices; but it is necessary to use a large
  number of iterations in the linear-equation solver: 
  liter(pot) = - 999. [The minus, combined with UWATCH = T and
  TSTSWP = 1, cause the progress of the convergence to be shown
  on the VDU].
  
  The PHOTON USE commands show that the potential-flow solution
  predicts a pressure RISE at the rear of the vehicle which does
  not occur in practice.
  
  ENDDIS
#pause
   ************************************************************
   Q1 created by VDI menu, Version 3.4, Date 17/12/01
CPVNAM=VDI;SPPNAM=Core
 ************************************************************
 LIBREF =  128
 ************************************************************
  Group 1. Run Title
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,60,42,30,1.000000E-04)
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Solved variables list
         SOLVE(P1  ,U1  ,V1  ,W1  )
 STORE(P1,  U1  ,V1  ,W1  )
 SOLVE(POT)
  
 ************************************************************
  Group 13. Boundary & Special Sources
 PATCH(IN,EAST,NX,NX,1,NY,1,NZ,1,1)
 COVAL(IN,POT,1.E3,9)

 PATCH(OUT,WEST,1,1,1,NY,1,NZ,1,1)
 COVAL(OUT,POT,1.E3,-9)
 
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =    1
 ************************************************************
  Group 16. Terminate Iterations
 LITER(POT)=-500
 UWATCH=T
 ENDIT(POT)=1.E-9
 RESREF(POT)=1.E-20
 ISOLX=1;ISOLZ=1;ISOLY=1
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 POTVEL  =    T
 ************************************************************
 GVIEW(P,8.069060E-01,-5.902608E-01,2.224065E-02)
 GVIEW(UP,3.534978E-03,4.247757E-02,9.990911E-01)
 
> DOM,    SIZE,        1.000000E+00, 5.000000E-01, 5.000000E-01
> DOM,    MONIT,       3.993520E-01, 2.485000E-01, 5.072500E-02
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
> GRID,   RSET_X_1,     25,-1.070000E+00,G
> GRID,   RSET_X_6,     10, 1.100000E+00,G
> GRID,   RSET_Y_1,     11,-1.250000E+00,G
> GRID,   RSET_Y_9,     11, 1.250000E+00,G
> GRID,   RSET_Z_5,     10, 1.250000E+00,G
 
> OBJ,    NAME,        BODY
> OBJ,    POSITION,    5.000000E-01, 2.290000E-01, 5.000000E-03
> OBJ,    SIZE,        3.000000E-01, 4.200000E-02, 6.800000E-02
> OBJ,    CLIPART,     wednew
> OBJ,    ROTATION24,       19
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      198
 
> OBJ,    NAME,        R-FRONT
> OBJ,    POSITION,    7.325000E-01, 2.175000E-01, 0.000000E+00
> OBJ,    SIZE,        3.500000E-02, 5.000000E-03, 3.500000E-02
> OBJ,    CLIPART,     cylinder
> OBJ,    ROTATION24,        6
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      198
 
> OBJ,    NAME,        L-FRONT
> OBJ,    POSITION,    7.325000E-01, 2.775000E-01, 0.000000E+00
> OBJ,    SIZE,        3.500000E-02, 5.000000E-03, 3.500000E-02
> OBJ,    CLIPART,     cylinder
> OBJ,    ROTATION24,        5
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      198
 
> OBJ,    NAME,        R-BACK
> OBJ,    POSITION,    5.000000E-01, 2.025000E-01, 0.000000E+00
> OBJ,    SIZE,        4.000000E-02, 2.000000E-02, 4.000000E-02
> OBJ,    CLIPART,     cylinder
> OBJ,    ROTATION24,        5
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      198
 
> OBJ,    NAME,        L-BACK
> OBJ,    POSITION,    5.000000E-01, 2.775000E-01, 0.000000E+00
> OBJ,    SIZE,        4.000000E-02, 2.000000E-02, 4.000000E-02
> OBJ,    CLIPART,     cylinder
> OBJ,    ROTATION24,        5
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,      198
DISTIL=T
 EX(P1  )=   3.357E+02 
 EX(U1  )=   1.787E+01 
 EX(V1  )=   4.074E-01 
 EX(W1  )=   2.895E-01 
 store(prps);EX(PRPS)=   7.421E+00 
 EX(POT )=   3.718E+00