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