Exemplo n.º 1
0
    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()

Exemplo n.º 2
0
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()