223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

corrected indices in flux density output

This commit is contained in:
Christoph Kords 2011-01-28 12:49:58 +00:00
parent 0cf1404db0
commit b72d75ed05
1 changed files with 9 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
code/constitutive_nonlocal.f90
View File

@ -2375,39 +2375,39 @@ do o = 1,phase_Noutput(material_phase(g,ip,el))
cs = cs + ns cs = cs + ns
case ('fluxdensity_edge_neg_x') case ('fluxdensity_edge_neg_x')
constitutive_nonlocal_postResults(cs+1:cs+ns) = rhoSgl(:,2) * constitutive_nonlocal_v(:,2,g,ip,el) * m_currentconf(1,:,2) constitutive_nonlocal_postResults(cs+1:cs+ns) = - rhoSgl(:,2) * constitutive_nonlocal_v(:,2,g,ip,el) * m_currentconf(1,:,1)
cs = cs + ns cs = cs + ns
case ('fluxdensity_edge_neg_y') case ('fluxdensity_edge_neg_y')
constitutive_nonlocal_postResults(cs+1:cs+ns) = rhoSgl(:,2) * constitutive_nonlocal_v(:,2,g,ip,el) * m_currentconf(2,:,2) constitutive_nonlocal_postResults(cs+1:cs+ns) = - rhoSgl(:,2) * constitutive_nonlocal_v(:,2,g,ip,el) * m_currentconf(2,:,1)
cs = cs + ns cs = cs + ns
case ('fluxdensity_edge_neg_z') case ('fluxdensity_edge_neg_z')
constitutive_nonlocal_postResults(cs+1:cs+ns) = rhoSgl(:,2) * constitutive_nonlocal_v(:,2,g,ip,el) * m_currentconf(3,:,2) constitutive_nonlocal_postResults(cs+1:cs+ns) = - rhoSgl(:,2) * constitutive_nonlocal_v(:,2,g,ip,el) * m_currentconf(3,:,1)
cs = cs + ns cs = cs + ns
case ('fluxdensity_screw_pos_x') case ('fluxdensity_screw_pos_x')
constitutive_nonlocal_postResults(cs+1:cs+ns) = rhoSgl(:,3) * constitutive_nonlocal_v(:,3,g,ip,el) * m_currentconf(1,:,3) constitutive_nonlocal_postResults(cs+1:cs+ns) = rhoSgl(:,3) * constitutive_nonlocal_v(:,3,g,ip,el) * m_currentconf(1,:,2)
cs = cs + ns cs = cs + ns
case ('fluxdensity_screw_pos_y') case ('fluxdensity_screw_pos_y')
constitutive_nonlocal_postResults(cs+1:cs+ns) = rhoSgl(:,3) * constitutive_nonlocal_v(:,3,g,ip,el) * m_currentconf(2,:,3) constitutive_nonlocal_postResults(cs+1:cs+ns) = rhoSgl(:,3) * constitutive_nonlocal_v(:,3,g,ip,el) * m_currentconf(2,:,2)
cs = cs + ns cs = cs + ns
case ('fluxdensity_screw_pos_z') case ('fluxdensity_screw_pos_z')
constitutive_nonlocal_postResults(cs+1:cs+ns) = rhoSgl(:,3) * constitutive_nonlocal_v(:,3,g,ip,el) * m_currentconf(3,:,3) constitutive_nonlocal_postResults(cs+1:cs+ns) = rhoSgl(:,3) * constitutive_nonlocal_v(:,3,g,ip,el) * m_currentconf(3,:,2)
cs = cs + ns cs = cs + ns
case ('fluxdensity_screw_neg_x') case ('fluxdensity_screw_neg_x')
constitutive_nonlocal_postResults(cs+1:cs+ns) = rhoSgl(:,4) * constitutive_nonlocal_v(:,4,g,ip,el) * m_currentconf(1,:,4) constitutive_nonlocal_postResults(cs+1:cs+ns) = - rhoSgl(:,4) * constitutive_nonlocal_v(:,4,g,ip,el) * m_currentconf(1,:,2)
cs = cs + ns cs = cs + ns
case ('fluxdensity_screw_neg_y') case ('fluxdensity_screw_neg_y')
constitutive_nonlocal_postResults(cs+1:cs+ns) = rhoSgl(:,4) * constitutive_nonlocal_v(:,4,g,ip,el) * m_currentconf(2,:,4) constitutive_nonlocal_postResults(cs+1:cs+ns) = - rhoSgl(:,4) * constitutive_nonlocal_v(:,4,g,ip,el) * m_currentconf(2,:,2)
cs = cs + ns cs = cs + ns
case ('fluxdensity_screw_neg_z') case ('fluxdensity_screw_neg_z')
constitutive_nonlocal_postResults(cs+1:cs+ns) = rhoSgl(:,4) * constitutive_nonlocal_v(:,4,g,ip,el) * m_currentconf(3,:,4) constitutive_nonlocal_postResults(cs+1:cs+ns) = - rhoSgl(:,4) * constitutive_nonlocal_v(:,4,g,ip,el) * m_currentconf(3,:,2)
cs = cs + ns cs = cs + ns
case ('d_upper_edge') case ('d_upper_edge')