VECFEM3 Reference Manual: vembldexm01.html

    #*******************************************************************
    #**
    #**        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
    #**
    #*******************************************************************