From 7605aab46fee3ad97b9496f01a034be32614ea1b Mon Sep 17 00:00:00 2001
From: Andreas Suter <andreas.suter@psi.ch>
Date: Tue, 15 May 2012 07:57:06 +0000
Subject: [PATCH] shift the calculated field by the dead layer

---
 src/tests/nonlocal/PPippard.cpp | 10 ++++++++--
 1 file changed, 8 insertions(+), 2 deletions(-)

diff --git a/src/tests/nonlocal/PPippard.cpp b/src/tests/nonlocal/PPippard.cpp
index def47843..e278477c 100644
--- a/src/tests/nonlocal/PPippard.cpp
+++ b/src/tests/nonlocal/PPippard.cpp
@@ -420,12 +420,18 @@ void PPippard::SaveField()
   fprintf(fp, "%% Data --------------------------------------\n");
   fprintf(fp, "%% z (nm), B/B_0 \n");
   if (fParams.specular) {
+    for (Int_t i=0; i<(Int_t)(fParams.deadLayer/f_dz); i++) {
+      fprintf(fp, "%lf, %lf\n", f_dz*(Double_t)i, 1.0);
+    }
     for (Int_t i=0; i<PippardFourierPoints/2; i++) {
-      fprintf(fp, "%lf, %lf\n", f_dz*(Double_t)i, fFieldB[i+fShift][1]);
+      fprintf(fp, "%lf, %lf\n", fParams.deadLayer + f_dz*(Double_t)i, fFieldB[i+fShift][1]);
     }
   } else {
+    for (Int_t i=0; i<(Int_t)(fParams.deadLayer/(f_dz * XiP_T(fParams.t))); i++) {
+      fprintf(fp, "%lf, %lf\n", f_dz * XiP_T(fParams.t) * (Double_t)i, 1.0);
+    }
     for (Int_t i=0; i<PippardDiffusePoints; i++) {
-      fprintf(fp, "%lf, %lf\n", f_dz * XiP_T(fParams.t) * (Double_t)i, (*fFieldDiffuse)(i));
+      fprintf(fp, "%lf, %lf\n", fParams.deadLayer + f_dz * XiP_T(fParams.t) * (Double_t)i, (*fFieldDiffuse)(i));
     }
   }