Commit b051de0b authored by Gregor Michalicek's avatar Gregor Michalicek

Add error message in io/io_matrix.F90/read_matrix_DA if the matrix dimensions are wrong

parent 0f321e18
......@@ -87,7 +87,13 @@ CONTAINS
l_real=mat%l_real
READ(id,rec=rec,iostat=err) matsize1,matsize2
if (matsize1<1) call judft_error("Data not found in file")
IF (mat%matsize1.NE.matsize1.OR.mat%matsize2.NE.matsize2) CALL mat%alloc(l_real,matsize1,matsize2)
IF (mat%matsize1.NE.matsize1.OR.mat%matsize2.NE.matsize2) THEN
WRITE(*,*) 'mat%matsize1 = ', mat%matsize1
WRITE(*,*) 'mat%matsize2 = ', mat%matsize2
WRITE(*,*) 'stored matsize1 = ', matsize1
WRITE(*,*) 'stored matsize2 = ', matsize2
CALL juDFT_error("matrix size wrong", calledby="read_matrix_DA")
END IF
IF (mat%l_real) THEN
READ(id,rec=rec,iostat=err) matsize1,matsize2,mat%data_r
ELSE
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment