Merge pull request #641 from yeganer/fix_j_blue_estimator_hack

Fix detailed radiative rates

tardis/montecarlo/base.py

@@ -137,8 +137,10 @@ def run(self, model, no_of_packets, no_of_virtual_packets=0, nthreads=1,last_run
         # Workaround so that j_blue_estimator is in the right ordering
         # They are written as an array of dimension (no_of_shells, no_of_lines)
         # but python expects (no_of_lines, no_of_shells)
-        self.j_blue_estimator = self.j_blue_estimator.flatten().reshape(
+        self.j_blue_estimator = np.ascontiguousarray(
+                self.j_blue_estimator.flatten().reshape(
                 self.j_blue_estimator.shape, order='F')
+                )
         self.Edotlu_estimator = self.Edotlu_estimator.flatten().reshape(
                 self.Edotlu_estimator.shape, order='F')