From e8628996ff1d88678cf095b27de0f32ebbfcdf03 Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Fri, 2 Sep 2011 13:50:05 +0000 Subject: [PATCH] added IO error in case mixed BC that allow rotation are given --- code/DAMASK_spectral.f90 | 5 ++++- code/IO.f90 | 2 ++ 2 files changed, 6 insertions(+), 1 deletion(-) diff --git a/code/DAMASK_spectral.f90 b/code/DAMASK_spectral.f90 index 2d7596419..fe4818c89 100644 --- a/code/DAMASK_spectral.f90 +++ b/code/DAMASK_spectral.f90 @@ -252,8 +252,11 @@ program DAMASK_spectral print '(a,i5)', 'Loadcase:', loadcase if (.not. followFormerTrajectory(loadcase)) & print '(a)', 'drop guessing along trajectory' - if (any(bc_mask(:,:,1,loadcase) .eqv. bc_mask(:,:,2,loadcase)))& ! exclusive or masking only + if (any(bc_mask(:,:,1,loadcase) .eqv. bc_mask(:,:,2,loadcase)))& ! exclusive or masking only call IO_error(31,loadcase) + if (any(bc_mask(1:3,1:3,2,loadcase).and.transpose(bc_mask(1:3,1:3,2,loadcase)).and.& + reshape((/.false.,.true.,.true.,.true.,.false.,.true.,.true.,.true.,.false./),(/3,3/))))& + call IO_error(38,loadcase) if (velGradApplied(loadcase)) then do j = 1, 3 if (any(bc_mask(j,:,1,loadcase) .eqv. .true.) .and.& diff --git a/code/IO.f90 b/code/IO.f90 index c432bd82c..eb4f2ce60 100644 --- a/code/IO.f90 +++ b/code/IO.f90 @@ -1198,6 +1198,8 @@ endfunction msg = 'non-positive result frequency in spectral loadcase' case (37) msg = 'incomplete loadcase' + case (38) + msg = 'mixed boundary conditions allow rotation' case (40) msg = 'path rectification error' case (41)