From effec313ee1c233fa0bc9f93cc0b745ce039ea77 Mon Sep 17 00:00:00 2001 From: Pratheek Shanthraj Date: Tue, 21 Apr 2015 12:19:41 +0000 Subject: [PATCH] fixed bug in cleavage system set up for ortho lattice structures --- code/lattice.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/code/lattice.f90 b/code/lattice.f90 index b02852536..bb3c9a3df 100644 --- a/code/lattice.f90 +++ b/code/lattice.f90 @@ -1515,8 +1515,8 @@ subroutine lattice_initializeStructure(myPhase,CoverA,aA,aM,cM) myNtrans = 0_pInt myNcleavage = lattice_ortho_Ncleavage do i = 1_pInt, myNcleavage ! assign cleavage system vectors - cd(1:3,i) = lattice_iso_systemCleavage(1:3,i)/math_norm3(lattice_iso_systemCleavage(1:3,i)) - cn(1:3,i) = lattice_iso_systemCleavage(4:6,i)/math_norm3(lattice_iso_systemCleavage(4:6,i)) + cd(1:3,i) = lattice_iso_systemCleavage(1:3,i)/math_norm3(LATTICE_ortho_systemCleavage(1:3,i)) + cn(1:3,i) = lattice_iso_systemCleavage(4:6,i)/math_norm3(LATTICE_ortho_systemCleavage(4:6,i)) ct(1:3,i) = math_vectorproduct(cd(1:3,i),cn(1:3,i)) enddo lattice_NcleavageSystem(1:lattice_maxNcleavageFamily,myPhase) = lattice_iso_NcleavageSystem