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.