- checked with gfortran73, gfortran83 and NAG - yet TODO: still getting an error with ifort (dunno why)