diff --git a/code/constitutive_nonlocal.f90 b/code/constitutive_nonlocal.f90 index 9c117eafc..9104eb6ea 100644 --- a/code/constitutive_nonlocal.f90 +++ b/code/constitutive_nonlocal.f90 @@ -1313,11 +1313,14 @@ if (verboseDebugger .and. selectiveDebugger) then !$OMPEND CRITICAL (write2out) endif -if (.not. ( mesh_element(2,el) == 1_pInt & - .or. mesh_element(2,el) == 6_pInt & - .or. mesh_element(2,el) == 7_pInt & - .or. mesh_element(2,el) == 8_pInt )) & - call IO_error(-1,el,ip,g,'element type not supported for nonlocal constitution') +select case(mesh_element(2,el)) + case (1,6,7,8,9) + ! all fine + case default + call IO_error(-1,el,ip,g,'element type not supported for nonlocal constitution') + + +end select myInstance = phase_constitutionInstance(material_phase(g,ip,el)) myStructure = constitutive_nonlocal_structure(myInstance)