From 78f74e899177f4a0d125a68888d6a5b2d6a63208 Mon Sep 17 00:00:00 2001 From: Christoph Kords Date: Thu, 26 May 2011 09:23:13 +0000 Subject: [PATCH] forgot to use "mesh_ipCenterOfGravity" --- code/DAMASK_spectral.f90 | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/code/DAMASK_spectral.f90 b/code/DAMASK_spectral.f90 index a36e0b05d..ba8794822 100644 --- a/code/DAMASK_spectral.f90 +++ b/code/DAMASK_spectral.f90 @@ -47,6 +47,7 @@ program DAMASK_spectral use prec, only: pInt, pReal use IO use math + use mesh, only: mesh_ipCenterOfGravity use CPFEM, only: CPFEM_general, CPFEM_initAll use numerics, only: err_div_tol, err_stress_tol, err_stress_tolrel, err_defgrad_tol,& itmax, memory_efficient, mpieNumThreadsInt @@ -422,7 +423,7 @@ program DAMASK_spectral do k = 1, resolution(3); do j = 1, resolution(2); do i = 1, resolution(1) ielem = ielem + 1_pInt call CPFEM_general(CPFEM_mode,& ! first element in first iteration retains CPFEM_mode 1, - coordinates(1:3,i,j,k), + coordinates(1:3,i,j,k),& defgradold(i,j,k,:,:), defgrad(i,j,k,:,:),& ! others get 2 (saves winding forward effort) temperature,timeinc,ielem,1_pInt,& cstress,dsde, pstress, dPdF) @@ -481,7 +482,7 @@ program DAMASK_spectral do k = 1, resolution(3); do j = 1, resolution(2); do i = 1, resolution(1) ielem = ielem + 1_pInt call CPFEM_general(CPFEM_mode,& ! first element in first iteration retains CPFEM_mode 1, - coordinates(1:3,i,j,k), + coordinates(1:3,i,j,k),& defgradold(i,j,k,:,:), defgrad(i,j,k,:,:),& temperature,timeinc,ielem,1_pInt,& cstress,dsde, pstress, dPdF)