TALK=f;RUN(1,1)
  DISPLAY
  This case differs from caseno=4 of case 204 in respect of the 
  boundary condition for U1. Here U1 is fixed at a single point,
  whereas in s204 U1 was fixed along the whole boundary.
  The solutions differ in  respect of U1, but not of V1, epsx, 
  epsy, etc.
  ENDDIS
  PHOTON USE
  p;;;;
  
  

  set prop off
  cl
  msg x-displacement fields
  gr ou z 1
  cont U1 z 1  fil;.0001
  pause
  cl
  msg y-displacement fields
  gr ou z 1
  cont V1 z 1 fil;.0001
  vec z 1 col 0
  pause
  cl
  msg drh1 fields
  gr ou z 1
  cont drh1 z 1 fil;.0001
  vec z 1  col 0
  pause
  cl
  msg enul fields
  gr ou z 1
  cont enul z 1 fil;.0001
  vec z 1 col 0
  pause

  ENDUSE

 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 TEXT(2D xy 2-material plate in tension; s208
LIBREF 
TITLE
integer(caseno)     
mesg(caseno = 1 : z free
mesg(caseno = 2 : z fixed
mesgm(Note: "Theoretical" values not provided
caseno=1
label ask
mesg(caseno=:caseno: Enter 1, 2 or blank
readvdu(caseno,int,1)
if(caseno.lt.1) then
 goto ask
endif  
if(caseno.gt.2) then
 goto ask
endif 
caseno
libref=caseno 
  Declarations and settings
REAL(FY,LX,LY,POISSON,YOUNG) 
FY= 40.0e6  ! H/m^2 = 40 N/mm^2
LX=90.e-3
LY=120.e-3 
  LX=LY/10
INTEGER(NXBODY,NYBODY)

 ************************************************************
  Group 2. Time dependence
 STEADY  =    T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES  =    T
 NXBODY = 10
 GRDPWR(X,NXBODY,LX,1)  

 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NYBODY = 10
 GRDPWR(Y,NYBODY,LY,1)  
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 NZ=1
 ZWLAST  = 0.001
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    T
 SOLVE(P1,V1,U1)
 SOLUTN(P1  ,Y,Y,Y,N,N,N)
 SOLUTN(U1  ,Y,Y,Y,N,N,Y)
 SOLUTN(V1  ,Y,Y,Y,N,N,Y)
 
 STORE(PRPS,DRH1,ENUL)
 STORE(STRX,STRY,STRZ,STXY)
 STORE(EPSY,EPSX,EPSZ)
 ************************************************************
  GROUP 8. ITERATION NUMBERS ETC
RESFAC=1.e-7
RESREF(V1)=0.0 
RESREF(U1)=0.0
LITER(V1) = 50
LITER(U1) = 50 
LITER(P1) = 50 

 ************************************************************
  GROUP 9. PROPERTIES
  
 CSG10='Q1'                  ! materials with various POISSON ratios
  MATFLG=T;NMAT=2         
  160    7800.0    0.3       473.0   43.0      1.0e-5   0.5e-11 
  161    7800.0    0.3       473.0   43.0      1.0e-5   0.5e-10
 
 ************************************************************
  GROUP 11. INITIAL VALUES
fiinit(p1)=0.0
fiinit(u1)=0.0
fiinit(v1)=0.0
 
 FIINIT(PRPS)=160
 PATCH(BODY2,INIVAL,1,NX,1,NY/2,1,1,1,1)
 INIT(BODY2,PRPS,FIXVAL,161)
 ************************************************************
  GROUP 13. BOUNDARY & SPECIAL SOURCES
 
PATCH(UP,NWALL,1,NX,NY,NY,1,1,1,1) ! top - fixed V1
COVAL(UP,V1,1,0.0)
PATCH(POINTU,CELL,1,1,NY,NY,1,1,1,1)  !  top-left - fixed U1
COVAL(POINTU,U1,FIXVAL,0)
   *** bottom - upward force
PATCH(FORC01,NORTH,1,NX,1,1,1,1,1,1)
COVAL(FORC01,V1,FIXFLU,FY)

IF(CASENO.EQ.2) THEN
SPEDAT(BOUNDARY,ZCONST,R,1.e20)
ELSE
SPEDAT(BOUNDARY,ZCONST,R,0.0)
ENDIF 
 ************************************************************
  GROUP 15. TERMINATE SWEEPS
 LSWEEP  =      1800
 ISG21=LSWEEP
 ************************************************************
  GROUP 17. RELAXATION
     #CONPROM
 RELAX(P1  ,LINRLX, 1.000000E+00)

      spedat(rlxfac,rlxu1d,r,0.3)  
      spedat(rlxfac,rlxv1d,r,0.3)  
 
 ************************************************************
  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
IXMON = NX-2
IYMON = 2
IZMON = 1
#maxmin
#endpause
  nprint=1;lsweep=20
STOP