all but the first loop: calculated.f

For the second and further loops only the field d(*) containing the length of the edges and d_{min} has to be determined. The h-field is obtained by interpolation within the unrefined mesh.