@ -139,7 +139,7 @@ contains
f_ => d_null_func_3d
end if
deltah = 1.d0/(idim+2)
deltah = 1.d0/(idim+1)
sqdeltah = deltah*deltah
deltah2 = 2.d0* deltah
f_ => s_null_func_3d
@ -231,7 +231,7 @@ contains
f_ => d_null_func_2d
deltah = done/(idim+2)
deltah = done/(idim+1)
deltah2 = (2*done)* deltah
@ -248,7 +248,7 @@ contains
f_ => s_null_func_2d
deltah = sone/(idim+2)
deltah = sone/(idim+1)
deltah2 = (2*sone)* deltah