From 45344cff6b6e0543820760ac0cb0ca4ceabae3f7 Mon Sep 17 00:00:00 2001 From: Franz Roters Date: Mon, 4 Feb 2013 14:34:01 +0000 Subject: [PATCH] added check if nGrains==1 when non-local plasticity is used --- code/IO.f90 | 2 ++ code/constitutive.f90 | 1 + 2 files changed, 3 insertions(+) diff --git a/code/IO.f90 b/code/IO.f90 index 478495ace..28a48f02c 100644 --- a/code/IO.f90 +++ b/code/IO.f90 @@ -1367,6 +1367,8 @@ subroutine IO_error(error_ID,e,i,g,ext_msg) case (212_pInt) msg = 'unknown plasticity output:' + case (252_pInt) + msg = 'nonlocal plasticity works only for direct CPFEM, i.e. one grain per integration point' case (253_pInt) msg = 'element type not supported for nonlocal plasticity' diff --git a/code/constitutive.f90 b/code/constitutive.f90 index 0330e8d9a..ac431879a 100644 --- a/code/constitutive.f90 +++ b/code/constitutive.f90 @@ -382,6 +382,7 @@ subroutine constitutive_init constitutive_sizePostResults(g,i,e) = constitutive_dislotwin_sizePostResults(myInstance) case (constitutive_nonlocal_label) + if(myNgrains/=1_pInt) call IO_error(252_pInt, e,i,g) select case(FE_geomtype(mesh_element(2,e))) case (7_pInt,8_pInt,9_pInt,10_pInt) ! all fine