6 #define ERROR(rval) if (0 .ne. rval) call exit(1)
23 integer comm1, mysize,myproc,ier
26 ibase_entitysethandle handle
27 ibase_entityhandle root_set
29 ibase_entityarriterator iter
30 integer local_dims(6),global_dims(6)
31 integer geom_dim,num_regions, num_verts, num_quads, count
58 call imesh_newmesh(
'MOAB', mesh, ier)
62 call imesh_createstructuredmesh(%VAL(mesh), local_dims,
63 1 global_dims, %VAL(0),%VAL(0),%VAL(0), %VAL(1), %VAL(-1),
64 1 %VAL(-1), %VAL(-1), %VAL(0), %VAL(1), %VAL(1), handle, ier)
67 call imesh_getrootset(%VAL(mesh), root_set, ier)
70 call imesh_getgeometricdimension(%VAL(mesh), geom_dim, ier)
73 call imesh_getnumoftype(%VAL(mesh), %VAL(root_set),
74 1 %VAL(ibase_face), num_quads, ier)
77 call imesh_getnumoftype(%VAL(mesh), %VAL(root_set),
78 1 %VAL(ibase_vertex), num_verts, ier)
82 call imesh_initentarriter(%VAL(mesh), %VAL(root_set),
83 1 %VAL(ibase_face), %VAL(imesh_quadrilateral),%VAL(num_quads),
86 call imesh_createtagwithoptions(%VAL(mesh),
"XM1",
87 1
"moab:TAG_STORAGE_TYPE=DENSE; moab:TAG_DEFAULT_VALUE=0.0",
88 1 %VAL(5), %VAL(ibase_double), tagh, ier)
90 call imesh_tagiterate(%VAL(mesh),%VAL(tagh),%VAL(iter),rpxm1,
93 call imesh_endentarriter(%VAL(mesh), %VAL(iter), ier)
102 call imesh_dtor(%VAL(mesh), ier)
105 if (rsum .ne. 5*64*64)
call exit(1)