From 3eeb6bee2b565faa5193335505eca502188f2913 Mon Sep 17 00:00:00 2001
From: Laura Nervo <lnervo@esrf.fr>
Date: Mon, 10 Dec 2012 10:43:15 +0000
Subject: [PATCH] gtCrystCalculateReflections : fix small bug
Signed-off-by: Laura Nervo <laura.nervo@esrf.fr>
git-svn-id: https://svn.code.sf.net/p/dct/code/trunk@954 4c865b51-4357-4376-afb4-474e03ccb993
---
zUtil_Cryst/gtCrystCalculateReflections.m | 9 ++++++---
1 file changed, 6 insertions(+), 3 deletions(-)
diff --git a/zUtil_Cryst/gtCrystCalculateReflections.m b/zUtil_Cryst/gtCrystCalculateReflections.m
index 54addce0..8a315132 100644
--- a/zUtil_Cryst/gtCrystCalculateReflections.m
+++ b/zUtil_Cryst/gtCrystCalculateReflections.m
@@ -58,13 +58,13 @@ disp(msg)
[C,Cmat] = gtReadTextFile([filename typelist '.txt'],'%f %f %f %f','\n','#',[1 3]);
hkl = Cmat(:,1:3);
-sintheta = Cmat(:,4);
+sinth = Cmat(:,4);
if nargout > 1
load parameters;
energy = parameters.acq.energy;
- theta = asind(sintheta)*gtConvEnergyToWavelength(energy);
- if nargout == 3
+ theta = asind(sinth)*gtConvEnergyToWavelength(energy);
+ if nargout > 2
% for planes with Miller Indexes hkl -> hkil
hkil = hkl;
hkil(:,4) = hkl(:,3);
@@ -75,6 +75,9 @@ if nargout > 1
for ii=1:size(hkl,1)
fprintf('%3d %3d %3d %3d %7.3f\n',hkil(ii,:),theta(ii))
end
+ if nargout > 3
+ sintheta = sinth;
+ end
end
else
disp(' h k l sin(theta)/lambda')
--
GitLab