def save_region_field_meshes(self, filename_trunk): output("saving regions of fields...") for field in self.fields: fregion = self.domain.regions[field.region_name] output("field %s: saving regions..." % field.name) for region in self.domain.regions: if not fregion.contains(region): continue output(region.name) aux = Mesh.from_region_and_field(region, field) aux.write("%s_%s_%s.mesh" % (filename_trunk, region.name, field.name), io="auto") output("...done") output("...done")