TALK=F;RUN(1,1)                                           
  DISPLAY
  2DXY Piezoelectric: Membrana Only  
  ENDDIS

  PHOTON USE 
  p;phi 
  1 1 1
  
  clear
  set prop off
  gr ou z 1
  con disy Z 1 fi;1.e-8
  pause
  
  clear
  gr ou z 1
  con disx Z 1 fi;1.e-8
  pause


  clear
  gr ou z 1
  SET VEC COMP - - -
  SET VEC COMP DISX DISY -
  VECT IZ 1 SH DASH 0
  msg Displacement vector
  pause

  autoplot 
  FILE 
  phi 5 

  clear
  data 1                                                                          
  DISY Y 1 X 1 m                                                                 
  colA 1          
  MU X 1.e6 1 1          
  MU Y 1.e6 1 1          
  scale
  plot 1                                                                      
  msg DISY(X) for Y=0
  pause

  ENDUSE 

 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
TEXT(Steady 2DXY Piezoelectric; S901
libref=901
TITLE
 
  Declarations and settings
REAL(LY,LX) 
LY=4.75e-6            ! 4.55 mkM
LX=30.e-6             !  half of membrane 
INTEGER(NYBODY, NXBODY)
 NYBODY = 16
 NXBODY = 50

real(StrX0, E21, V0, EY0 ,TimeMax)
E21 =  6.75      ! Cuol/m**2
V0  = -6.7    ! Volt
EY0 = V0/LY 
StrX0 = E21*EY0


 ************************************************************
  Group 2. Time dependence
 STEADY  =    T

 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES  =    T
 GRDPWR(Y,NYBODY,LY,1)  

 ************************************************************
  Group 4. Y-Direction Grid Spacing
 GRDPWR(X,NXBODY,LX,1)  
 
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 NZ=1
 ZWLAST  = LX
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
 SOLVE(DISY,DISX)
 TERMS (DISY,N,N,Y,Y,Y,N)
 TERMS (DISX,N,N,Y,Y,Y,N)
 STORE(STRY,STRX)
 
 STORE(PRPS,DEN1,ENUL,DVO1,DRH1)
 STORE(EPSY,EPSX,EPSZ)

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

 ************************************************************
  GROUP 9. PROPERTIES
  
 CSG10='Q1'                  ! materials with various POISSON ratios
  MATFLG=T;NMAT=1         
  160    7500.0    0.29  473.0   43.0    1.0e-5   1.65e-11 
   
 ************************************************************
  GROUP 11. INITIAL VALUES

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

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

 ************************************************************
  GROUP 13. BOUNDARY & SPECIAL SOURCES
 
PATCH(Left,WWALL,1,1,1,NY,1,1,1,LSTEP)    ! Left end is symmetrical 
COVAL(Left,DISX,1,0)
PATCH(BR,CELL,NX,NX,1,1,1,1,1,LSTEP)    ! bottom Right end fixed 
COVAL(BR,DISY,FIXVAL,0.0)
COVAL(BR,DISX,FIXVAL,0.0)

PATCH(FORCE,EAST,NX,NX,1,NY,1,1,1,LSTEP)
COVAL(FORCE,DISX,FIXFLU,StrX0)


 ************************************************************
  GROUP 15. TERMINATE SWEEPS
 LSWEEP  =    600  
 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 = NY/2
IYMON = 1
IZMON = 1
  inform7begin
 
  inform7end
  
  
STOP