for varn in var_map['2_data']: var=ds.variables[varn] try: tag=mesh.getTagHandle(pack_data_tag_name(varn, var.dtype.char, 'S2')) except Exception as ex: print "No tag found for variable '%s'" % varn continue var.set_auto_maskandscale(False) if len(var.shape) == 4: arr=var[ti,0,:,:].reshape(ntopo) else: arr=var[ti,:,:].reshape(ntopo) utils.set_packed_data(tag, tsvert, arr) # tag[set]=arr # Delete the 'default' tags that we don't need/want mesh.destroyTag(mesh.getTagHandle('DIRICHLET_SET'), True) mesh.destroyTag(mesh.getTagHandle('GEOM_DIMENSION'), True) mesh.destroyTag(mesh.getTagHandle('GLOBAL_ID'), True) mesh.destroyTag(mesh.getTagHandle('MATERIAL_SET'), True) mesh.destroyTag(mesh.getTagHandle('NEUMANN_SET'), True) mesh.save(out_path) print "Saved to %s" % out_path ds.close()
for ti in range(ntimes): set=mesh.createEntSet(ordered=False) set.add(quads) for varn in var_map['data']: var=ds.variables[varn] try: tag=mesh.getTagHandle(utils.pack_data_tag_name(varn, var.dtype.char)) except Exception as ex: print "No tag found for variable '%s'" % varn continue var.set_auto_maskandscale(False) if len(var.shape) == 4: arr=var[ti,0,:,:].reshape(nquads) else: arr=var[ti,:,:].reshape(nquads) # tag[set]=arr utils.set_packed_data(tag, set, arr) time_set.addChild(set) mesh.save(out_path) print "Saved to %s" % out_path ds.close()