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")
