diff --git a/mc_slab.c b/mc_slab.c
index a127c0c21d804c33c316702fbe60a2513c779e83..99291417d3a22b459c7a07609c102e167ad0e306 100644
--- a/mc_slab.c
+++ b/mc_slab.c
@@ -301,6 +301,7 @@ int slab_mc(int nphot, int seed)
         Magfield_lab[0] = Magfield_lab_new[0];
         Magfield_lab[1] = Magfield_lab_new[1];
         Magfield_lab[2] = Magfield_lab_new[2];
+        
       }
       else
 
@@ -515,7 +516,8 @@ void comptonization(double E_lab,
   //================================================================== */
 
   status = depolarization_probability(E_erf, E_new_erf, costheta_erf, phik, PolarFlag);
-
+  
+ 
   if (status == 1)
   {
     exit(1);
@@ -537,12 +539,6 @@ void comptonization(double E_lab,
 
   elmag_lorentztransform(ElectField_new_erf, MagField_Erf, el_vel_vect, Efield_lab_new, Magfield_lab_new, 2);
 
-  // =========================================================================== */
-  // Back convert from (K1)' in ERF to lab*/
-  // =========================================================================== */
-
-  versor_boost(k_new_erf, el_vel_vect, k_new_lab, 2);
-
   // ================================================================== */
   // Back convert electron energy from ERF to lab */
   // ================================================================== */
diff --git a/solve_rte.c b/solve_rte.c
index 44139b8009e452f0c6bec189f3194dda16f05f92..9d7ee65aec981b62b8edfe33168fb844997ac5b1 100644
--- a/solve_rte.c
+++ b/solve_rte.c
@@ -32,7 +32,7 @@ int solve_rte(int k_iter, int seed_distribution, int Nstep_tau)
   printf("\n\nReceveing Nstep_tau %d\n", Nstep_tau);
   int ii, jj, kk, status;
 
-  Nstep_mu = 20;
+  Nstep_mu = 30;
   /* Nstep_tau = 100;*/
 
   tau0 = disktau / 2.;
@@ -191,15 +191,6 @@ int solve_rte(int k_iter, int seed_distribution, int Nstep_tau)
     ptr_data->Spline_I2d_r_downstream = Spline_I2d_r_downstream;
 
 
-
-
-   /* printf("CHECK HERE !\n");
-    double zio=0;
-        int status = gsl_spline2d_eval_extrap_e(Spline_I2d_l_downstream, 2*tau0, 0.5, xacc_2d, yacc_2d, &zio);
-
-       printf("================> status %d val %lf\n", status, zio);
-*/
-
     /*=======================================================================*/
     /*Solve the RTE at the quadrature points*/
     /*=======================================================================*/