AIMAG() is compliant with Fortran77 standard - comments updated to reflect that