Skip to main content

GmshModelGeo example 2

In this example we test following methods:

In this example we create a circle geometry and mesh.

Usage​

PROGRAM main
USE easifemBase
USE easifemClasses
IMPLICIT NONE
TYPE( gmsh_ ) :: gmsh
  • Instance of [[Gmsh_]].
INTEGER( I4B ) :: ierr, tag
REAL( DFP ), PARAMETER :: R=1.0, Xc=0.0, Yc=0.0, Zc=0.0
REAL( DFP ) :: x, y, z, meshSize
CHARACTER( LEN = GMSH_API_MAX_STR_LEN ), ALLOCATABLE :: names( : )
CHARACTER( LEN = GMSH_API_MAX_STR_LEN ) :: name
CHARACTER(LEN=*), PARAMETER :: PREFIX="t1"
  • [[Gmsh_#Initialize]]
  • [[GmshModel_#Add]]
  • [[GmshModel_#SetFileName]]
ierr = gmsh%initialize()
ierr = gmsh%model%add( PREFIX )
ierr = gmsh%model%setFileName( PREFIX//".geo" )

Points​

  • [[GmshModelGeo_#AddPoint]]

Point-1: Specifying tag explicitly

x = Xc-R; y = Yc; z = Zc; meshSize = 0.1; tag=1
ierr = gmsh%model%geo%addPoint(x,y,z,meshSize,tag)
x = Xc+R; y = Yc; z = Zc; meshSize = 0.1; tag=2
ierr = gmsh%model%geo%addPoint(x,y,z,meshSize,tag)
x = Xc; y = Yc; z = Zc; meshSize = 0.1; tag=3
ierr = gmsh%model%geo%addPoint(x,y,z,meshSize,tag)

Line and curves​

  • [[GmshModelGeo_#AddCircleArc]]

  • explicit tag

  • automatic tag

  • if tag=-1 it means automatic tag.

  • nx, ny, nz are optional, default value is zero.

ierr = gmsh%model%geo%AddCircleArc( &
  & startTag=1, &
  & centerTag=3, &
  & endTag=2, &
  & tag=1, &
  & nx=0.0, &
  & ny=0.0, &
  & nz=0.0 )
call ok(ierr .NE. 0_I4B, "AddCircleArc(1)")

Adding another circle arc.

ierr = gmsh%model%geo%AddCircleArc( &
  & startTag=2, &
  & centerTag=3, &
  & endTag=1, &
  & tag=2)
call ok(ierr .NE. 0_I4B, "AddCircleArc(2)")

CurveLoop​

  • [[GmshModelGeo_#AddCurveLoop]]

  • reorient is optional, default value is false.

  • tag is optional

ierr = gmsh%model%geo%addCurveLoop( &
  & curveTags=[1,2], tag=1, reorient=.TRUE.)

AddPlaneSurface​

  • [[GmshModelGeo_#AddPlaneSurface]]

  • tag is explicitly given.

ierr = gmsh%model%geo%addPlaneSurface(wireTags=[1], tag=1)

Synchronize​

  • [[GmshModelGeo_#Synchronize]]
ierr = gmsh%model%geo%Synchronize()

Write​

  • [[Gmsh_#Write]]
ierr = gmsh%Write(PREFIX//".msh")

Cleanup​

  • [[Gmsh_#Finalize]]
  ierr = gmsh%finalize()

END PROGRAM main