Work around to avoid using the RANK intrinsic.
Work around to avoid using the RANK intrinsic.

Oops.  It turns out that the rank intrinsic is in TS 29113 (a technical
specification beyond Fortran 2008) and is not implemented in many
compilers.  I was misled by the fact that gfortran has it!

Fortunately rank(x) and size(shape(x)) give the same result...

hande-qmc/hande,ruthfranklin/hande,hande-qmc/hande,hande-qmc/hande,hande-qmc/hande,hande-qmc/hande