visu.F90

Go to the documentation of this file.

 
 
PROGRAM main
 
  IMPLICIT NONE
 
  ! if the header file is in the "include" directory
  ! #include "libmmg2df.h"
 
  ! if the header file is in "include/mmg/mmg2d"
#include "mmg/mmg2d/libmmg2df.h"
 
  mmg5_data_ptr_t    :: mesh
  mmg5_data_ptr_t    :: sol,mmgmet
  INTEGER            :: ier,argc,i,j,opt
 
  CHARACTER(len=300) :: exec_name,filename,fileout
 
  print*,"  -- 2D MESH GENERATION FOR VISUALIZATION"
 
  argc =  command_argument_count();
  CALL get_command_argument(0, exec_name)
 
 
  IF ( argc /= 2 ) THEN
     print*," Usage: ",trim(adjustl(exec_name)),&
          " input_file_name output_file_name" &
     "  Generation of a triangular mesh for solution vizualisation (at VTK &
     & file format) from a Medit mesh file containing only points and the &
     & associated solution file."
     CALL exit(1);
  ENDIF
 
  ! Name and path of the mesh file
  CALL get_command_argument(1, filename)
  CALL get_command_argument(2, fileout)
 
  !!> ------------------------------ STEP   I --------------------------
  !! 1) Initialisation of mesh and sol structures */
  !! args of InitMesh:
  !! MMG5_ARG_start: we start to give the args of a variadic func
  !! MMG5_ARG_ppMesh: next arg will be a pointer over a MMG5_pMesh
  !! &mmgMesh: pointer toward your MMG5_pMesh (that store your mesh)
  !! MMG5_ARG_ppMet: next arg will be a pointer over a MMG5_pSol storing a metric
  !! &mmgSol: pointer toward your MMG5_pSol (that store your metric)
 
  mesh    = 0
  sol     = 0
  mmgmet  = 0
 
  CALL mmg2d_init_mesh(mmg5_arg_start, &
       mmg5_arg_ppmesh,mmgmesh,mmg5_arg_ppmet,mmgmet, &
       mmg5_arg_ppmet,sol, &
       mmg5_arg_end);
 
 
  !!> 2) Build initial mesh and solutions in MMG5 format
  !! Two solutions: just use the MMG2D_loadMesh function that will read a .mesh(b)
  !! file formatted or manually set your mesh using the MMG2D_Set* functions
 
  !!> Automatic loading of the mesh and multiple solutions
  CALL mmg2d_loadmesh(mesh,trim(adjustl(filename)),&
       len(trim(adjustl(filename))),ier)
  IF ( ier /= 1 )  CALL exit(102)
 
  CALL mmg2d_loadsol(mesh,sol,trim(adjustl(filename)),&
       len(trim(adjustl(filename))),ier)
  IF ( ier /= 1 )  CALL exit(103)
 
  !!> ------------------------------ STEP II ---------------------------
  !! Mesh generation
  CALL mmg2d_set_iparameter(mesh,mmgmet,mmg2d_iparam_noinsert,1,ier);
  IF ( ier /= 1 )  CALL exit(104)
 
  CALL mmg2d_set_iparameter(mesh,mmgmet,mmg2d_iparam_noswap,1,ier);
  IF ( ier /= 1 )  CALL exit(105)
 
 CALL mmg2d_set_iparameter(mesh,mmgmet,mmg2d_iparam_nomove,1,ier);
  IF ( ier /= 1 )  CALL exit(106)
 
  CALL mmg2d_mmg2dmesh(mesh,mmgmet,ier)
  IF ( ier /= mmg5_success )  CALL exit(200)
 
  !!> ------------------------------ STEP III --------------------------
  !! Save the new data
  !! Use the MMG2D_saveMesh/MMG2D_saveAllSols functions
  !! save the mesh
  CALL mmg2d_savevtkmesh(mesh,sol,trim(adjustl(fileout)),len(trim(adjustl(fileout))),ier)
  IF ( ier /= 1 ) CALL exit(300)
 
  !!> 3) Free the MMG2D structures
  CALL mmg2d_free_all(mmg5_arg_start, &
       mmg5_arg_ppmesh,mesh,mmg5_arg_ppmet,mmgmet, &
       mmg5_arg_ppmet,sol, &
       mmg5_arg_end)
 
END PROGRAM main