662
662
call element_quality_pain_p0(output_positions, new_metric, quality)
664
664
call vtk_write_fields("adapted_quality", index = output_quality_index, &
665
665
& position = output_positions, model = output_positions%mesh, &
666
666
& sfields = (/quality/), tfields = (/new_metric/))