TALK=F;RUN(1,1)                                           
  DISPLAY
  A FIXED rectangular beam of dimensions Lx 
  makes the longitudal fluctuations.
  ENDDIS

  PHOTON USE 
  p;parphi 
  1 1 4000
  
  clear
  set prop off
  view y
  gr ou y 1
  con disx y 1 fi;0.0001
  msg DISX(X,Z=Time)
  pause
  
  
  autoplot 
  FILE 
  parphi 5 

  clear
  data 1                                                                          
  DISX Z 10 X 1 m                                                                  
  col1 1          
  data 1                                                                          
  DISX Z 50 X 1 m                                                                  
  colA 2          
  data 1                                                                          
  DISX Z 100 X 1 m                                                                  
  colF 3          
  scale
  plot 1 2 3                                                                      
  msg DISX(X) for Time=10(W),50(Y),100(R)
  pause

  clear
  data 1                                                                          
  DISX X 50 Z 1 m                                                                  
  col1 1          
  scale
  plot 1                                                                      
  msg DISX(Time) for X=LX/2
  pause

  ENDUSE

 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
TEXT(Wave in fixed 1D x beam; S671
libref=677
TITLE
integer(caseno)
mesg(caseno 1 : Initial velocity of displacement 
mesg(caseno 2 : Initial displacement
caseno=1
mesgm(caseno = :caseno: Enter another if not OK
readvdu(caseno, int, 1)
caseno 

 
  Declarations and settings
REAL(FX,LYZ,LX,POISSON,YOUNG) 
FX= 40.0e6  ! H/m^2 = 40 N/mm^2
LYZ=40.e-3
LX=120.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(NXBODY,NXIB,NXIE)

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

 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES  =    T
 NY = 1
 YVLAST  = LYZ

 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NXBODY = 100
 NXIB = 42
 NXIE = 58
 
 GRDPWR(X,NXBODY,LX,1)  
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 NZ=1
 ZWLAST  = LYZ
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
 SOLVE(DISX)
 SOLUTN(DISX  ,Y,Y,Y,N,N,Y)

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

 STORE(VDSX)
 
 STORE(PRPS,DEN1,ENUL,DVO1,DRH1)
 STORE(STRX)
 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

FIINIT(DISX)=0.0
FIINIT(VDSX)=0.0
FIINIT(PRPS)=160

RELAX(DISX,LINRLX, 1.000000E+00)

 ************************************************************
  GROUP 13. BOUNDARY & SPECIAL SOURCES
 
PATCH(Left,WWALL,1,1,1,1,1,1,1,LSTEP)    ! left end fixed 
COVAL(Left,DISX,1,0.0)
PATCH(Right,EWALL,NX,NX,1,1,1,1,1,LSTEP)    ! right end fixed 
COVAL(Right,DISX,1,0.0)

PATCH(InBlock,INIVAL,NXIB,NXIE,1,1,1,1,1,1)

if(caseno.eq.1) then
 INIT (InBlock,VDSX,FIXVAL,1)
endif
if(caseno.eq.2) then
 INIT (InBlock,DISX,FIXVAL,1.e-6)
endif
 ************************************************************
  GROUP 15. TERMINATE SWEEPS
 LSWEEP  =    10  
 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
IXMON = NX/2
IYMON = 1
IZMON = 1

IDISPA = 1

  inform7begin
 
  inform7end
  
  
STOP