[patch, libfortran] Speed up cshift for dim > 1

[Dominique d'Humières]

Hi Thomas,

Your patch works as advertised! For the record, the following CSHIFT is still 4 times slower than the DO loop

    td(:,:) = cshift(array=t(:,:), shift=vect(:), dim=1)

Thanks for working on this issue.

Dominique