TALK=F;RUN(1,1)                                           
  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 
 

 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
TEXT(Free wave in 2D "tuning-fork" model; S673
libref=673
TITLE

integer(caseno)
mesg(caseno 1 : Positive initial Velocity in TOP beam
mesg(caseno 2 : Positive initial Velocity in BOTH beam 
mesg(caseno 3 : Symmetrical initial conditions 
caseno=1
mesgm(caseno = :caseno: Enter another if not OK
readvdu(caseno, int, 1)
caseno 
 
  Declarations and settings
REAL(LZ,LX,LY,POISSON,YOUNG) 
LX = 120.e-3
LY = 36.e-3
LZ = 10.e-3 
YOUNG   = 1/0.5E-11   ! Young's modulus used for exact solutions
POISSON=0.3           ! Poisson's ratio used for exact solutions
                      ! The above must correspond to the data
                      ! below MATFLG
INTEGER(NYBODY,NXBODY)

 ************************************************************
  Group 2. Time dependence
 STEADY  =    F
 RSET(U,0,30.e-4,200)

 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES  =    T
 NXBODY = 32
 GRDPWR(X,NXBODY,LX,1)  
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NYBODY = 36
 GRDPWR(Y,NYBODY,LY,1)  
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 NZ=1
 ZWLAST  = LZ
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
 SOLVE(DISX,DISY)
 SOLUTN(DISY  ,Y,Y,Y,N,N,Y)
 SOLUTN(DISX  ,Y,Y,Y,N,N,Y)

 TERMS (DISY,N,N,Y,Y,Y,N)
 TERMS (DISX,N,N,Y,Y,Y,N)

 STORE(VDSY,VDSX)
 
 STORE(PRPS,DEN1,ENUL,DVO1,DRH1)
 STORE(STRX,STRY)
 STORE(EPSY,EPSX,EPSZ)

 ************************************************************
  GROUP 8. ITERATION NUMBERS ETC

 ************************************************************
  GROUP 9. PROPERTIES
  
 CSG10='Q1'                  ! materials with various POISSON ratios
  MATFLG=T;NMAT=1         
  160    7800.0    0.3       473.0   43.0      1.0e-5   0.5e-11 
   
 ************************************************************
  GROUP 11. INITIAL VALUES
 
   **** CALCULATE initial conditions ***

FIINIT(disx)=0
FIINIT(disy)=0
FIINIT(VDSX)=0.0
FIINIT(VDSY)=0.0
FIINIT(PRPS)=160

PATCH(FLUID1,INIVAL,1,8,1,12,1,1,1,1)
INIT(FLUID1,PRPS,FIXVAL,0)

PATCH(FLUID2,INIVAL,1,8,ny-11,ny,1,1,1,1)
INIT(FLUID2,PRPS,FIXVAL,0)

PATCH(FLUID3,INIVAL,13,nx,13,ny-12,1,1,1,1)
INIT(FLUID3,PRPS,FIXVAL,0)

  **** Initial Velocity ---
PATCH(INITVEL,INIVAL,26,nx,NY-11,NY,1,1,1,1)
INIT(INITVEL,VDSY,FIXVAL,0.1)
if(caseno.eq.2) then
 PATCH(INITVEL1,INIVAL,26,nx,1,12,1,1,1,1)
 INIT(INITVEL1,VDSY,FIXVAL, 0.1)
endif
if(caseno.eq.3) then
 PATCH(INITVEL1,INIVAL,26,nx,1,12,1,1,1,1)
 INIT(INITVEL1,VDSY,FIXVAL, -0.1)
endif


RELAX(DISY,LINRLX, 1.000000E+00)
RELAX(DISX,LINRLX, 1.000000E+00)

 ************************************************************
  GROUP 13. BOUNDARY & SPECIAL SOURCES
 
PATCH(Left,WWALL,1,1,13,24,1,1,1,LSTEP)    ! west edge fixed 
COVAL(Left,DISY,1,0.0)
COVAL(Left,DISX,1,0.0)


 ************************************************************
  GROUP 15. TERMINATE SWEEPS
 LSWEEP  =      120
 ISG21=LSWEEP

 ************************************************************
  GROUP 17. RELAXATION
 RELAX(P1  ,LINRLX, 1.000000E+00)
 ************************************************************
  GROUP 19. DATA TRANSMITTED TO GROUND
 STRA    =    T
 PARSOL  =    F
 ISG52   =    3   !   probe & res

 ************************************************************
  GROUP 23.FIELD PRINT-OUT & PLOT CONTROL
TSTSWP = - 1   ! graphic-mode
NYPRIN = 1
NXPRIN = 1
IXMON = NX
IYMON = NY
IZMON = 1

IDISPA = 1

  inform7begin
 
  inform7end
  
 ************************************************************
  
STOP