TALK=F;RUN( 1, 1)
 
 ************************************************************
   Q1 created by VDI menu, Version 2007, Date 26/09/07
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
  Echo DISPLAY / USE settings
  DISPLAY
    2D (x,y) Y-shape beam is free for Time <= 0.
 
   :----- fixed (DisX=DisY=0)
   :
   :    **************************
   V    *            :  Top      :
   /*****            :           :
   /*      ***********************
   /*      *
   /*      ***********************
   /*****            :           :
        *            :  Bottom   :
        **************************
 
    Initial condition Time = 0:
    caseno = 1:  in Top  VDSY = 0.1
    caseno = 2:  in Top and Bottom  VDSY = 0.1
    caseno = 3:  in Top VDSY = 0.1, in Bottom  VDSY = -0.1
  ENDDIS
  PHOTON USE
  p;parphi
  1 1 200
 
  clear
  set prop off
  view y
  gr ou y 30
  con disy y 30 fi;0.0001
  msg DISY (X,Z=Time) on axis of top beam
  pause
 
  clear
  view 1 1 1
  gr ou z 10 x 1 12 y 13 24
  gr ou z 10 x 9 m y 1 12
  gr ou z 10 x 9 m y 25 m
  SET VEC COMP - - -
  SET VEC COMP DISX DISY -
  VECT IZ 10 SH DASH 0
 
  gr ou z 100 x 1 12 y 13 24
  gr ou z 100 x 9 m y 1 12
  gr ou z 100 x 9 m y 25 m
  SET VEC COMP - - -
  SET VEC COMP DISX DISY -
  VECT IZ 100 SH DASH 0
 
  gr ou z 190 x 1 12 y 13 24
  gr ou z 190 x 9 m y 1 12
  gr ou z 190 x 9 m y 25 m
  SET VEC COMP - - -
  SET VEC COMP DISX DISY -
  VECT IZ 190 SH DASH 0
  msg Displacement vector (Z=Time=10, 100, 190)
 
  pause
 
  autoplot
  FILE
  parphi 5
 
  clear
  data 1
  DISY X m Y m Z 1 m
  col1 1
  data 1
  DISY X m Y 1 Z 1 m
  colA 2
  scale
  plot 1 2
  msg DISY(Time) on end of top(W) and bottom(Y) beam
  pause
 
 
  ENDUSE
 ************************************************************
 IRUNN   =         1 ;LIBREF =       673
 ************************************************************
  Group 1. Run Title
 TEXT(Free wave in 2D "tuning-fork" model; S67)
 ************************************************************
  Group 2. Transience
 STEADY=F
    * Set overall time and no. of steps
 RSET(U,0.000000E+00,3.000000E-03,200)
    * Modify regions
 RSET(T,1,200,1.000000E+00)
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,32,36,1,1.000000E-05)
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
    * Non-default variable names
 NAME(137) =EPSZ ; NAME(138) =EPSX
 NAME(139) =EPSY ; NAME(140) =STRY
 NAME(141) =STRX ; NAME(142) =DRH1
 NAME(143) =DVO1 ; NAME(144) =ENUL
 NAME(145) =DEN1 ; NAME(146) =PRPS
 NAME(147) =VDSX ; NAME(148) =VDSY
 NAME(149) =DISY ; NAME(150) =DISX
    * Solved variables list
 SOLVE(DISY,DISX)
    * Stored variables list
 STORE(VDSY,VDSX,PRPS,DEN1,ENUL,DVO1,DRH1,STRX)
 STORE(STRY,EPSY,EPSX,EPSZ)
    * Additional solver options
 SOLUTN(DISY,Y,Y,Y,N,N,Y)
 SOLUTN(DISX,Y,Y,Y,N,N,Y)
 
 ************************************************************
  Echo InForm settings for Group  7
  INFORM7BEGIN
integer(caseno)
    caseno 1 : Positive initial Velocity in TOP beam
    caseno 2 : Positive initial Velocity in BOTH beam
    caseno 3 : Symmetrical initial conditions
   ==================
caseno=3
   ==================
  INFORM7END
 ************************************************************
  Group 8. Terms & Devices
 TERMS (DISY,N,N,Y,Y,Y,N)
 TERMS (DISX,N,N,Y,Y,Y,N)
 ************************************************************
  Group 9. Properties
    * Domain material index is 160 signifying:
    * NOTSET
 SETPRPS(1,160,Q1)
 ENUT    = 0.000000E+00
 DRH1DP  = 5.000000E-12
 DVO1DT  = 1.000000E-05
   *  List of user-defined materials to be read by EARTH
    MATFLG=T;IMAT=1
   *  Name
   *Ind. Dens.  Viscos.  Spec.heat  Conduct.  Expans.  Compr.
   *  
    160 7800.0 0.3 473.0 43.0 1.0E-5 0.5E-11
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(EPSZ) =  1.001000E-10 ;FIINIT(EPSX) =  1.001000E-10
 FIINIT(EPSY) =  1.001000E-10 ;FIINIT(STRY) =  1.001000E-10
 FIINIT(STRX) =  1.001000E-10 ;FIINIT(DRH1) =  1.001000E-10
 FIINIT(DVO1) =  1.001000E-10 ;FIINIT(ENUL) =  1.001000E-10
 FIINIT(DEN1) =  1.001000E-10 ;FIINIT(PRPS) =  1.600000E+02
 FIINIT(VDSX) =  0.000000E+00 ;FIINIT(VDSY) =  0.000000E+00
 FIINIT(DISY) =  0.000000E+00 ;FIINIT(DISX) =  0.000000E+00
 
 PATCH (INITVEL ,INIVAL,4,0,0,0,0,0,1,1)
 
 PATCH (INITVEL1,INIVAL,5,0,0,0,0,0,1,1)
 
 
 INIADD  =    F
 ************************************************************
  Echo InForm settings for Group 11
  INFORM11BEGIN
(INITIAL of VDSY at INITVEL is 0.1)
if(caseno.eq.2) then
(INITIAL of VDSY at INITVEL1 is 0.1)
endif
if(caseno.eq.3) then
(INITIAL of VDSY at INITVEL1 is -0.1)
endif
  INFORM11END
 ************************************************************
  Group 12. Convection and diffusion adjustments
   No PATCHes used for this Group
 ************************************************************
  Group 13. Boundary & Special Sources
 
 PATCH (LEFT    ,WWALL ,6,0,0,0,0,0,1,200)
 COVAL (LEFT    ,DISY, 1.000000E+00, 0.000000E+00)
 COVAL (LEFT    ,DISX, 1.000000E+00, 0.000000E+00)
 
 EGWF    =    T
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =       120
 RESFAC  = 1.000000E-03
 ************************************************************
  Group 16. Terminate Iterations
 ************************************************************
  Group 17. Relaxation
 RELAX(DISY,LINRLX, 1.000000E+00)
 RELAX(DISX,LINRLX, 1.000000E+00)
 SPEDAT(SET,GXMONI,TRANSIENT,L,F)
 ************************************************************
  Group 18. Limits
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 USEGRD  =    T  ;USEGRX =    T
 STRA    =    T
 ASAP    =    T
 PARSOL  =    F
 ISG21   =       120
 ISG52   =         3
 CSG10   ='q1'
 ************************************************************
  Group 20. Preliminary Printout
 ECHO    =    T
 ************************************************************
  Group 21. Print-out of Variables
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =        32 ;IYMON  =        36 ;IZMON  =         1
 NPRMON  =    100000
 NPRMNT  =         1
 TSTSWP  =        -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT  =    100000
 NTPRIN  =    100000 ;ISTPRF =         1 ;ISTPRL =    100000
 NXPRIN  =         1
 NYPRIN  =         1
 ISWPRF  =         1 ;ISWPRL =    100000
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 IDISPA  =         1 ;IDISPB =         0 ;IDISPC =         0
 
 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.200000E-01, 3.600000E-02, 1.000000E-02
> DOM,    MONIT,       1.181250E-01, 3.550000E-02, 5.000000E-03
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    INCREMENT,   1.000000E-02, 1.000000E-02, 1.000000E-02
> GRID,   RSET_X_1,      8, 1.000000E+00
> GRID,   RSET_X_2,      4, 1.000000E+00
> GRID,   RSET_X_3,     13, 1.000000E+00
> GRID,   RSET_X_4,      7, 1.000000E+00
> GRID,   RSET_Y_1,     12, 1.000000E+00
> GRID,   RSET_Y_2,     12, 1.000000E+00
> GRID,   RSET_Y_3,     12, 1.000000E+00
  > GRID,   RSET_Z_1,      1, 1.000000E+00
 
> OBJ,    NAME,        FLUID1
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        3.000000E-02, 1.200000E-02, 1.000000E-02
> OBJ,    GEOMETRY,    cubet
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    0,  Air at 20 deg C, 1 atm, treated as i
> OBJ,    TIME_LIMITS,   ALWAYS_ACTIVE
> OBJ,    SCAL_FIXF,   Per-unit-volume
 
> OBJ,    NAME,        FLUID2
> OBJ,    POSITION,    0.000000E+00, 2.400000E-02, 0.000000E+00
> OBJ,    SIZE,        3.000000E-02, 1.200000E-02, 1.000000E-02
> OBJ,    GEOMETRY,    cubet
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    0,  Air at 20 deg C, 1 atm, treated as i
> OBJ,    TIME_LIMITS,   ALWAYS_ACTIVE
> OBJ,    SCAL_FIXF,   Per-unit-volume
 
> OBJ,    NAME,        FLUID3
> OBJ,    POSITION,    4.500000E-02, 1.200000E-02, 0.000000E+00
> OBJ,    SIZE,        7.500000E-02, 1.200000E-02, 1.000000E-02
> OBJ,    GEOMETRY,    cubet
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    0,  Air at 20 deg C, 1 atm, treated as i
> OBJ,    TIME_LIMITS,   ALWAYS_ACTIVE
> OBJ,    SCAL_FIXF,   Per-unit-volume
 
> OBJ,    NAME,        INITVEL
> OBJ,    POSITION,    9.375001E-02, 2.400000E-02, 0.000000E+00
> OBJ,    SIZE,        2.624998E-02, 1.200000E-02, 1.000000E-02
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 
> OBJ,    NAME,        INITVEL1
> OBJ,    POSITION,    9.375002E-02, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        2.624997E-02, 1.200000E-02, 1.000000E-02
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 
> OBJ,    NAME,        LEFT
> OBJ,    POSITION,    0.000000E+00, 1.200000E-02, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 1.200000E-02, 1.000000E-02
> OBJ,    GEOMETRY,    default
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        USER_DEFINED
 ************************************************************
STOP