TagIterateF.F

Go to the documentation of this file.

c TagIterate: Get a pointer to tag memory and interact directly with that 
c 
c This program shows how to get a pointer to tag memory and assign it to
c a fortran array, allowing an application to bypass the API and/or share
c memory with MOAB.
 
c Usage: TagIterate <mesh_filename>
 
 program tag_iterate
#include "iMesh_f.h"
 
c declarations
 imesh_instance mesh
 ibase_entityhandle root_set
 ibase_entityarriterator iter
 ibase_taghandle tagh
 ibase_handle_t rpents, ents, rpdata
 real*8 tag_data
 pointer(rpdata, tag_data(0:*))
 pointer(rpents, ents(0:*))
 pointer(ipoffsets, ioffsets(0,*))
 integer ierr, ents_alloc, ents_size, num_regions, atend, count
 
c create the Mesh instance
 call imesh_newmesh("MOAB", mesh, ierr)
 
c get the root set
 call imesh_getrootset(%VAL(mesh), root_set, ierr)
 
c load the mesh
 call imesh_load(%VAL(mesh), %VAL(root_set), "125hex.vtk", "",ierr)
 
c get the number of regions
 call imesh_getnumoftype(%VAL(mesh), %VAL(root_set), 
 1 %VAL(ibase_region), num_regions, ierr)
 
c get an iterator over regions
 call imesh_initentarriter(%VAL(mesh), %VAL(root_set), 
 1 %VAL(ibase_region), %VAL(imesh_all_topologies),%VAL(num_regions),
 1 %VAL(0), iter, ierr) 
 
c create a tag to put on the regions
 call imesh_createtagwithoptions(%VAL(mesh), "dumtag", 
 1 "moab:TAG_STORAGE_TYPE=DENSE moab:TAG_DEFAULT_VALUE=0.0", 
 1 %VAL(1), %VAL(ibase_double), tagh, ierr)
 
c iterate over tag memory
 10 call imesh_tagiterate(%VAL(mesh), %VAL(tagh), %VAL(iter), 
 1 rpdata, count, ierr)
 
c step the iterator over count entities
 call imesh_stepentarriter(%VAL(mesh), %VAL(iter), %VAL(count), 
 1 atend, ierr)
 
 if (atend .eq. 0) go to 10
 
 
 call imesh_dtor(%VAL(mesh), ierr)
 
 end