%%%%%%%%%
Number of subdomains
%%%%%%%%%
*NEL
    1
    
%%%%%%%%%
Number of nodes of each subdomain
In the first column there is the subdomain number, in the second, its number of nodes
%%%%%%%%%
*NNOD
       1       40
      

%%%%%%%%%
The nodal cordinates of each subdomain are given below.
The coordinates of each vertex, its nodal number, and the number of elements between two adjacent vertices have to be given 
for all the subdomains
    *NODE
     1  --->subdomain number
     1  -0.100E+01  -0.100E+01  --->node number of the first vertex,  its x1 coordinate,   its x2 coordinate (node 1 x1=-0.100E+01 x2=-0.100E+01)
     2                          --->number of elements with equal length to be placed between vertex1 (node 1) and vertex2 (node 3)
     3  -0.600E+00  -0.100E+01  --->node number of the 2nd vertex, its x1 coordinate,   its x2 coordiante
     2
     5  -0.600E+00  -0.600E+00
     2
     7  -0.100E+01  -0.600E+00
     2
     9  -0.100E+01  -0.100E+01
     0

%%%%%%%%%
*NODE
     1
     1  -0.100E+01  -0.100E+01
     10
     11  0.100E+01  -0.100E+01
     10
     21  0.100E+01   0.100E+01
     10
     31  -0.100E+01  0.100E+01
     10
     41  -0.100E+01  -0.100E+01
     0


%%%%%%%%%
Boundary condition for each subdomain. 
Firstly, the prescribed displacements are specified, then, after a row of zeros, the tractions. In the first column the first node (element in case of presribed traction) 
for which the prescrption is valid, 
in the second column, the last node (or element). In the third and fourth column there is a zero if the component in x1 direction is not prescribed, 1 otherwise. 
In the fourth column the same, but for the component in x2 direction. Finally, in the 5th and 6th column, the prescribed values for the displacements or tractions in x1 and
x2 direction respectively.
%%%%%%%%%%%%%%%%%%%%%
    *BOUNDARY
     1   --->subdomain number
    1       3       1       1   0.000E+00   0.000E+00  --->zero diplacements in x1 and x2 direction prescribed from the first to the third node.
    0    0    0    0    0    0
    0    0    0    0    0    0
    0    0    0    0    0    0

%%%%%%%%%
*BOUNDARY
    1
    1      11       1       1   0.000E+00   0.000E+00
    13     21       1       1   0.000E+00   0.000E+00
    0    0    0    0    0    0
    29    30    1    1    0.d00    -1.d00
    0    0    0    0    0    0



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
In the section MATERIAL the shear modulus elastic parameter mu for the material is defined.
In the first two columns the subdomains having the same shear modulus are specified, in the third column the value of the shear modulus is defined.
*MATERIAL
       1       1       0.100E+01
       0       0          0


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%