From 4a8bc166349f13a6d277f67fce2ff5bf4aa53609 Mon Sep 17 00:00:00 2001
From: GudniRos <gkrosen@gmail.com>
Date: Wed, 25 Nov 2015 16:22:06 -0800
Subject: [PATCH] Fix a bug in reshapeing in adjoint projFieldsDeriv

---
 simpegMT/SurveyMT.py | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/simpegMT/SurveyMT.py b/simpegMT/SurveyMT.py
index 6c08264c..7547ca5d 100644
--- a/simpegMT/SurveyMT.py
+++ b/simpegMT/SurveyMT.py
@@ -299,7 +299,8 @@ class RxMT(Survey.BaseRx):
                 # Calculate the complex derivative
                 PDeriv_real = ZijN_uV(aHd*v) - aHd_uV(Zij.T*aHd*v)#
                 # NOTE: Need to reshape the output to go from 2*nU array to a (nU,2) matrix for each polarization
-                PDeriv_real = PDeriv_real.reshape((mesh.nE,2))
+                # PDeriv_real = np.hstack((mkvc(PDeriv_real[:len(PDeriv_real)/2],2),mkvc(PDeriv_real[len(PDeriv_real)/2::],2)))
+                PDeriv_real = PDeriv_real.reshape((2,mesh.nE)).T
             # Extract the data
             if real_or_imag == 'imag':
                 Pv = 1j*PDeriv_real