diff --git a/code/crystallite.f90 b/code/crystallite.f90 index 841d75e83..b64b97129 100644 --- a/code/crystallite.f90 +++ b/code/crystallite.f90 @@ -601,17 +601,27 @@ subroutine crystallite_stressAndItsTangent(updateJaco,rate_sensitivity) real(pReal) :: counter - if(iand(debug_level(debug_crystallite), debug_levelBasic) /= 0_pInt) then + if (iand(debug_level(debug_crystallite),debug_levelBasic) /= 0_pInt & + .and. ((e == debug_e .and. i == debug_i .and. g == debug_g) & + .or. .not. iand(debug_level(debug_crystallite), debug_levelSelective) /= 0_pInt)) then !$OMP CRITICAL (write2out) - write(6,'(/,a,i8,1x,i2,1x,i3)') '<< CRYST >> crystallite start at el ip g ', debug_e, debug_i, debug_g + write(6,'(/a,i8,1x,i8)') '<< CRYST >> element range ', FEsolving_execElem(1),FEsolving_execElem(2) + !$OMP END CRITICAL (write2out) + endif + + if (iand(debug_level(debug_crystallite),debug_levelExtensive) /= 0_pInt & + .and. ((e == debug_e .and. i == debug_i .and. g == debug_g) & + .or. .not. iand(debug_level(debug_crystallite), debug_levelSelective) /= 0_pInt)) then + !$OMP CRITICAL (write2out) + write(6,'(/,a,i8,1x,i2,1x,i3)') '<< CRYST >> values at el ip g ', e, i, g write(6,'(a,/,3(12x,3(f14.9,1x)/))') '<< CRYST >> F0 ', & - math_transpose33(crystallite_partionedF0(1:3,1:3,debug_g,debug_i,debug_e)) + math_transpose33(crystallite_partionedF0(1:3,1:3,g,i,e)) write(6,'(a,/,3(12x,3(f14.9,1x)/))') '<< CRYST >> Fp0', & - math_transpose33(crystallite_partionedFp0(1:3,1:3,debug_g,debug_i,debug_e)) + math_transpose33(crystallite_partionedFp0(1:3,1:3,g,i,e)) write(6,'(a,/,3(12x,3(f14.9,1x)/))') '<< CRYST >> Lp0', & - math_transpose33(crystallite_partionedLp0(1:3,1:3,debug_g,debug_i,debug_e)) + math_transpose33(crystallite_partionedLp0(1:3,1:3,g,i,e)) write(6,'(a,/,3(12x,3(f14.9,1x)/))') '<< CRYST >> F ', & - math_transpose33(crystallite_partionedF(1:3,1:3,debug_g,debug_i,debug_e)) + math_transpose33(crystallite_partionedF(1:3,1:3,g,i,e)) !$OMP END CRITICAL (write2out) endif @@ -2812,7 +2822,7 @@ subroutine crystallite_integrateStateFPI() - stateResiduum(1:mySizeDotState) ! need to copy to local variable, since we cant flush a pointer in openmp #ifndef _OPENMP - if (iand(debug_level(debug_crystallite), debug_levelBasic) /= 0_pInt & + if (iand(debug_level(debug_crystallite), debug_levelExtensive) /= 0_pInt & .and. ((e == debug_e .and. i == debug_i .and. g == debug_g) & .or. .not. iand(debug_level(debug_crystallite), debug_levelSelective) /= 0_pInt)) then write(6,'(a,i8,1x,i2,1x,i3,/)') '<< CRYST >> update state at el ip g ',e,i,g