diff --git a/code/config/Phase_Dislokmc_Tungsten.config b/code/config/Phase_Dislokmc_Tungsten.config
new file mode 100644
index 000000000..e27ac3588
--- /dev/null
+++ b/code/config/Phase_Dislokmc_Tungsten.config
@@ -0,0 +1,43 @@
+### $Id$ ###
+[Tungsten]
+elasticity          hooke
+plasticity          dislokmc
+
+
+### Material parameters ###
+lattice_structure     bcc
+C11                   523.0e9       # From Marinica et al. Journal of Physics: Condensed Matter(2013)
+C12                   202.0e9
+C44                   161.0e9
+
+grainsize             2.0e-5        # Average grain size [m] 2.0e-5
+SolidSolutionStrength 0.0           # Strength due to elements in solid solution
+
+### Dislocation glide parameters ###
+#per family
+Nslip               12 0 
+slipburgers         2.72e-10        # Burgers vector of slip system [m]
+rhoedge0            1.0e12          # Initial edge dislocation density [m/m**3] 1.0e12
+rhoedgedip0         1.0             # Initial edged dipole dislocation density [m/m**3]
+Qedge               2.725e-19       # Activation energy for dislocation glide [J]
+v0                  3560.3          # Initial glide velocity [m/s] 1.0e-4 (kmC)
+p_slip              0.16            # p-exponent in glide velocity
+q_slip              1.00            # q-exponent in glide velocity
+u_slip              2.47            # u-exponent of stress pre-factor (kmC)                                                                     
+s_slip              0.97            # self hardening in glide velocity (kmC)
+tau_peierls         2.03e9          # peierls stress [Pa]
+#v_slip              1.58            # v-exponent in glide velocity (kmC)
+
+
+#hardening
+dipoleformationfactor 0             # to have hardening due to dipole formation off
+CLambdaSlip         10.0            # Adj. parameter controlling dislocation mean free path
+D0                  4.0e-5          # Vacancy diffusion prefactor [m**2/s]
+Qsd                 4.5e-19         # Activation energy for climb [J]
+Catomicvolume       1.0             # Adj. parameter controlling the atomic volume [in b]
+Cedgedipmindistance 1.0             # Adj. parameter controlling the minimum dipole distance [in b]
+#interaction_slipslip    0 0 0 0 0 0
+#interaction_slipslip    0.1 0.1 0.14 0.14 0.14 0.14 
+#interaction_slipslip 	 1 0.1 0.9 0.5 0.1 0.7
+interaction_slipslip	 0.2 0.11 0.19 0.15 0.11 0.17
+#interaction_slipslip    1 1 1.4 1.4 1.4 1.4
diff --git a/code/constitutive_titanmod.f90 b/code/constitutive_titanmod.f90
index c1fd23982..961de4323 100644
--- a/code/constitutive_titanmod.f90
+++ b/code/constitutive_titanmod.f90
@@ -669,7 +669,7 @@ subroutine constitutive_titanmod_init(fileUnit)
  sanityChecks: do phase = 1_pInt, size(phase_plasticity)
    myPhase: if (phase_plasticity(phase) == PLASTICITY_TITANMOD_ID) then
      instance = phase_plasticityInstance(phase) 
-     if (sum(constitutive_titanmod_Nslip(:,instance)) <= 0_pInt) &
+     if (sum(constitutive_titanmod_Nslip(:,instance)) < 0_pInt) &
        call IO_error(211_pInt,el=instance,ext_msg='nslip ('//PLASTICITY_TITANMOD_label//')')
      if (sum(constitutive_titanmod_Ntwin(:,instance)) < 0_pInt) &
        call IO_error(211_pInt,el=instance,ext_msg='ntwin ('//PLASTICITY_TITANMOD_label//')')
@@ -718,7 +718,7 @@ subroutine constitutive_titanmod_init(fileUnit)
      if (constitutive_titanmod_twinhpconstant(instance) <= 0.0_pReal) &
        call IO_error(211_pInt,el=instance,ext_msg='twinhpconstant ('//PLASTICITY_TITANMOD_label//')')
      if (constitutive_titanmod_aTolRho(instance) <= 0.0_pReal) &
-       call IO_error(211_pInt,el=instance,ext_msg='aTolRho ('//PLASTICITY_TITANMOD_label//')')
+       call IO_error(211_pInt,el=instance,ext_msg='aTol_rho ('//PLASTICITY_TITANMOD_label//')')
      
 !--------------------------------------------------------------------------------------------------
 ! determine total number of active slip or twin systems