Hi,
I've located an error in one of your FORTRAN 77 interface to the LAPACK Single Value Decomposition routines (DGESVD). This is located in https://software.intel.com/sites/products/documentation/doclib/mkl_sa/11...
The following two lines must be added. (1) in variable declaration and (2) just before calling DGESVD.
(1)
DOUBLE PRECISION Asvd( size(A,1), size(A,2) )
(2)
Asvd = A
Best wishes,
Mohamed Ali Al-Badri