C conversion: Steven G. Johnson (stevenj@alum.mit.edu)
Original Fortran code: Joerg.M.Gablonsky (jmgablon@mailandnews.com)
