From 882c44a61d0b9e23c789ef92c2fed1eb751ea9c7 Mon Sep 17 00:00:00 2001 From: Philip Eisenlohr Date: Tue, 24 Jul 2012 15:29:23 +0000 Subject: [PATCH] changed the majority of debugging output to level "extensive" --- code/CPFEM.f90 | 20 +++++++++++--------- code/config/debug.config | 2 +- 2 files changed, 12 insertions(+), 10 deletions(-) diff --git a/code/CPFEM.f90 b/code/CPFEM.f90 index 787095df1..43d57bf7d 100644 --- a/code/CPFEM.f90 +++ b/code/CPFEM.f90 @@ -115,7 +115,8 @@ subroutine CPFEM_init use prec, only: pInt use debug, only: debug_level, & debug_CPFEM, & - debug_levelBasic + debug_levelBasic, & + debug_levelExtensive use IO, only: IO_read_jobBinaryFile use FEsolving, only: parallelExecution, & symmetricSolver, & @@ -145,7 +146,7 @@ subroutine CPFEM_init ! *** restore the last converged values of each essential variable from the binary file if (restartRead) then - if (iand(debug_level(debug_CPFEM), debug_levelBasic) /= 0_pInt) then + if (iand(debug_level(debug_CPFEM), debug_levelExtensive) /= 0_pInt) then !$OMP CRITICAL (write2out) write(6,'(a)') '<< CPFEM >> Restored state variables of last converged step from binary files' !$OMP END CRITICAL (write2out) @@ -236,6 +237,7 @@ subroutine CPFEM_general(mode, coords, ffn, ffn1, Temperature, dt, element, IP, use debug, only: debug_level, & debug_CPFEM, & debug_levelBasic, & + debug_levelExtensive, & debug_levelSelective, & debug_e, & debug_i, & @@ -394,7 +396,7 @@ subroutine CPFEM_general(mode, coords, ffn, ffn1, Temperature, dt, element, IP, j = 1:mesh_maxNips, & k = 1:mesh_NcpElems ) & constitutive_state0(i,j,k)%p = constitutive_state(i,j,k)%p ! microstructure of crystallites - if (iand(debug_level(debug_CPFEM), debug_levelBasic) /= 0_pInt) then + if (iand(debug_level(debug_CPFEM), debug_levelExtensive) /= 0_pInt) then !$OMP CRITICAL (write2out) write(6,'(a)') '<< CPFEM >> Aging states' if (debug_e == cp_en .and. debug_i == IP) then @@ -416,7 +418,7 @@ subroutine CPFEM_general(mode, coords, ffn, ffn1, Temperature, dt, element, IP, ! * dump the last converged values of each essential variable to a binary file if (restartWrite) then - if (iand(debug_level(debug_CPFEM), debug_levelBasic) /= 0_pInt) then + if (iand(debug_level(debug_CPFEM), debug_levelExtensive) /= 0_pInt) then !$OMP CRITICAL (write2out) write(6,'(a)') '<< CPFEM >> Writing state variables of last converged step to binary files' !$OMP END CRITICAL (write2out) @@ -485,7 +487,7 @@ subroutine CPFEM_general(mode, coords, ffn, ffn1, Temperature, dt, element, IP, if (terminallyIll .or. outdatedFFN1 .or. any(abs(ffn1 - materialpoint_F(1:3,1:3,IP,cp_en)) > defgradTolerance)) then if (.not. terminallyIll .and. .not. outdatedFFN1) then - if (iand(debug_level(debug_CPFEM), debug_levelBasic) /= 0_pInt) then + if (iand(debug_level(debug_CPFEM), debug_levelExtensive) /= 0_pInt) then !$OMP CRITICAL (write2out) write(6,'(a,1x,i8,1x,i2)') '<< CPFEM >> OUTDATED at element ip',cp_en,IP write(6,'(a,/,3(12x,3(f10.6,1x),/))') '<< CPFEM >> FFN1 old:',math_transpose33(materialpoint_F(1:3,1:3,IP,cp_en)) @@ -512,7 +514,7 @@ subroutine CPFEM_general(mode, coords, ffn, ffn1, Temperature, dt, element, IP, FEsolving_execElem(2) = cp_en FEsolving_execIP(1,cp_en) = IP FEsolving_execIP(2,cp_en) = IP - if (iand(debug_level(debug_CPFEM), debug_levelBasic) /= 0_pInt) then + if (iand(debug_level(debug_CPFEM), debug_levelExtensive) /= 0_pInt) then !$OMP CRITICAL (write2out) write(6,'(a,i8,1x,i2)') '<< CPFEM >> Calculation for element ip ',cp_en,IP !$OMP END CRITICAL (write2out) @@ -523,7 +525,7 @@ subroutine CPFEM_general(mode, coords, ffn, ffn1, Temperature, dt, element, IP, !* parallel computation and calulation not yet done elseif (.not. CPFEM_calc_done) then - if (iand(debug_level(debug_CPFEM), debug_levelBasic) /= 0_pInt) then + if (iand(debug_level(debug_CPFEM), debug_levelExtensive) /= 0_pInt) then !$OMP CRITICAL (write2out) write(6,'(a,i8,a,i8)') '<< CPFEM >> Calculation for elements ',FEsolving_execElem(1),' to ',FEsolving_execElem(2) !$OMP END CRITICAL (write2out) @@ -533,7 +535,7 @@ subroutine CPFEM_general(mode, coords, ffn, ffn1, Temperature, dt, element, IP, call mesh_build_subNodeCoords() ! update subnodal coordinates call mesh_build_ipCoordinates() ! update ip coordinates #endif - if (iand(debug_level(debug_CPFEM), debug_levelBasic) /= 0_pInt) then + if (iand(debug_level(debug_CPFEM), debug_levelExtensive) /= 0_pInt) then !$OMP CRITICAL (write2out) write(6,'(a,i8,a,i8)') '<< CPFEM >> Start stress and tangent ',FEsolving_execElem(1),' to ',FEsolving_execElem(2) !$OMP END CRITICAL (write2out) @@ -637,7 +639,7 @@ subroutine CPFEM_general(mode, coords, ffn, ffn1, Temperature, dt, element, IP, Temperature = materialpoint_Temperature(IP,cp_en) ! homogenized result except for potentially non-isothermal starting condition. endif - if (mode < 3 .and. iand(debug_level(debug_CPFEM), debug_levelBasic) /= 0_pInt & + if (mode < 3 .and. iand(debug_level(debug_CPFEM), debug_levelExtensive) /= 0_pInt & .and. ((debug_e == cp_en .and. debug_i == IP) & .or. .not. iand(debug_level(debug_CPFEM), debug_levelSelective) /= 0_pInt)) then !$OMP CRITICAL (write2out) diff --git a/code/config/debug.config b/code/config/debug.config index f279a219d..a17cc35cc 100644 --- a/code/config/debug.config +++ b/code/config/debug.config @@ -15,7 +15,7 @@ lattice # lattice.f90, possible value: ba constitutive # constitutive_*.f90 possible values: basic, extensive, selective crystallite # crystallite.f90 possible values: basic, extensive, selective homogenization # homogenization_*.f90 possible values: basic, extensive, selective -CPFEM # CPFEM.f90 possible values: basic, selective +CPFEM # CPFEM.f90 possible values: basic, extensive, selective spectral # DAMASK_spectral.f90 possible values: basic, fft, restart, divergence abaqus # ABAQUS FEM solver possible values: basic #