diff --git a/src/phase.f90 b/src/phase.f90 index 6bea38161..281f47de7 100644 --- a/src/phase.f90 +++ b/src/phase.f90 @@ -19,7 +19,7 @@ module phase implicit none private - type(Rotation), dimension(:,:,:), allocatable, protected :: & + type(Rotation), dimension(:,:,:), allocatable :: & material_orientation0 !< initial orientation of each grain,IP,element type(rotation), dimension(:,:,:), allocatable :: & diff --git a/src/phase_mechanical.f90 b/src/phase_mechanical.f90 index 383bf07a8..61782bab2 100644 --- a/src/phase_mechanical.f90 +++ b/src/phase_mechanical.f90 @@ -226,7 +226,7 @@ module subroutine mechanical_init(materials,phases) allocate(phase_mechanical_P(phases%length)) allocate(phase_mechanical_S0(phases%length)) - allocate(material_orientation0(homogenization_maxNconstituents,phases%length,maxval(material_phaseMemberAt))) + allocate(material_orientation0(homogenization_maxNconstituents,phases%length,maxVal(material_phaseMemberAt))) do ph = 1, phases%length Nconstituents = count(material_phaseAt == ph) * discretization_nIPs @@ -279,7 +279,7 @@ module subroutine mechanical_init(materials,phases) endif - !$OMP PARALLEL DO PRIVATE(ph,me) + !$OMP PARALLEL DO PRIVATE(ph,me,material,constituents,constituent) do el = 1, size(material_phaseMemberAt,3); do ip = 1, size(material_phaseMemberAt,2) do co = 1, homogenization_Nconstituents(material_homogenizationAt(el)) material => materials%get(discretization_materialAt(el))