#*******************************************************************
#**
#**        v e m b l d e x m 0 1
#**
#**  two coupled Poisson equation on the 2-dimensional unit cube.
#**  The mesh is read from a vecfem input file (e.g. generated by
#**  vemgen2d).
#**
#**   by L. Grosz                           Karlsruhe, Jan. 1995
#**
#*******************************************************************
#**
#**  The data set of this examples has two parts (search for
#**  'cut here'). The first part specifies the problem
#**  (please copy it to 'vembldexm01.equation') and the second part
#**  defines the control parameters  (please copy it to
#**  'vembldexm01.resource'). The FORTRAN code for the solution
#**  of the problem is generated by entering
#**  'vembuild vembldexm01' into your shell.
#**
#*******************************************************************
####### cut here to get vembldexm01.equation <<<<<<<<<<<<<<<<<<<<<<<<<
#*******************************************************************
#**
#**  The problem is a system of two coupled 2-dimensional Poisson
#**  equations with Dirichlet and Neuman boundary conditions.
#**
#**  The domain is [0,1] x [0,1], where the mesh is givem by
#**  rec.vem or is generated by vemgen2d. The mesh uses
#**  triangular elements of order 2 and line elements of order 2.
#**  Dirichlet conditions are set on all boundary points for the
#**  first component, but there is no Dirichlet condition for
#**  the second component.
#**
#*******************************************************************
#**
#**  The functions b1, b2, f1, f2, g1 and g2 are selected, so that
#**
#**     u1 = 1
#**     u2 = 2
#**
#**  gets the exact solution of the problem.
#**
     b1=1
     b2=2
     f1=10*b1+b2
     f2=b1+10*b2
     g1=0
     g2=0
#**
#*******************************************************************
#**
#**    The Dirichlet conditions:
#**
     u1=b1
     u2=b2
#**
#**  'u2=b2' is not necessary, since no Dirichlet conditions
#**  for the second component are specified in the input mesh.
#**
#*******************************************************************
#**
#**  the functional equation :
#**
    area{v1x1*u1x1+v1x2*u1x2+(10.*u1+    u2)*v1 - f1*v1
       + v2x1*u2x1+v2x2*u2x2+(    u1+10.*u2)*v2 - f2*v2}
       + line{v1*g1+v2*g2} =0
#**
#*******************************************************************
>>>>>>> cut here to get vembldexm01.resource <<<<<<<<<<<<<<<<<<<<<<<<<
#*******************************************************************
#**
#**  The problem has a two dimensional domain and two solution
#**  component:
#**
     DIM=2
     NK=2
#**
#*******************************************************************
#**
#**  One processor with maximal 20 Mbytes are used. Maximal 1000
#**  nodes and 1000 elements are allowed:
#**
    PROCESS_STORAGE=20
    PROCESS_MAXNN=1000
    PROCESS_MAXNE=1000
#**
#*******************************************************************
#**
#**  the is read from the file rec.vem. other meshes can be
#**  created vemgen2d.
#**
    MESH_PREP=print
    MESH_FILEIN=rec.vem
#**
#*******************************************************************
#**
#**  The solution components are printed into the file solution.out
#**  and the indicator, which considers the error of all components
#**  is printed into file error.out.
#**
    MESH_POSTP=print
    OUTPUT_INDEX=11
    OUTPUT_FILE=solution.out
    OUTPUT_ERRINDEX=11
    OUTPUT_ERRFILE=error.out
#**
#*******************************************************************