mesh_visu.f90

Go to the documentation of this file.

program mesh_visu

  use choral
  use choral_constants

  implicit none

  !! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  !! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  !!
  !!      VARAIBLE DEFINITION : BEGIN
  !!
  !!  

  !! verb    = verbosity level
  !!
  integer, parameter :: verb = 2

  !! msh_file  = mesh file name
  !! msh       = mesh 
  !! mesh_desc = a variable for the mesh analysis process
  !!
  character(len=150)      :: msh_file
  type(mesh)              :: msh
  integer, dimension(4,6) :: mesh_desc

  !!
  !!   VARAIBLE DEFINITION : END
  !!
  !! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  !! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

  call choral_init(verb=verb)

  write(*,*)
  write(*,*)'mesh_visu: start'
  write(*,*)""
  write(*,*)"    ASSEMBLE A MESH "
  write(*,*)""
  write(*,*)"    and visualise it with gmsh"
  write(*,*)""

  !!
  !!  settings: choose a mesh
  !!
  ! msh_file = trim(GMSH_DIR)//"disk/disk1_1.msh"
  msh_file = trim(gmsh_dir)//"disk/disk2_1.msh"
  ! msh_file = trim(GMSH_DIR)//"sphere/sphere1_1.msh"
  ! msh_file = trim(GMSH_DIR)//"sphere/sphere2_1.msh"


  write(*,*) ""
  write(*,*)"================================== ASSEMBLING THE MESH "
  msh = mesh(trim(msh_file), "gmsh")
  if (verb>1) call print(msh)

  write(*,*) ""
  write(*,*)"================================== ANALYSE THE MESH "
  mesh_desc = mesh_analyse(msh, verb-1)

  write(*,*) ""
  write(*,*)"================================== MESH DISPLAY "
  call system("gmsh -option "//trim(gmsh_dir)//"gmsh-options-mesh "//trim(msh_file))

  
end program mesh_visu