diff --git a/trunk/CPFEM.f90 b/trunk/CPFEM.f90
index 0a4c53d01..d5af213f5 100644
--- a/trunk/CPFEM.f90
+++ b/trunk/CPFEM.f90
@@ -106,6 +106,7 @@
 ! three dimensional stress state ?
     call math_init()
     call mesh_init()
+	call crystal_init()
     call constitutive_init()
     call CPFEM_init()
     CPFEM_first_call = .false.