all but the first loop: calculated.f

For the second and further loops the field d(*) containing the length of the edges and $d_{min}$ are recalculated. The h-field is obtained by interpolation within the unrefined mesh.