diff --git a/crpptbx/TCV/gdat_tcv.m b/crpptbx/TCV/gdat_tcv.m
index 8c3d00030af63b01797c14073835ce93d1a2935a..5dec42a3ecf7ce3d117cff409233e47d647694ef 100644
--- a/crpptbx/TCV/gdat_tcv.m
+++ b/crpptbx/TCV/gdat_tcv.m
@@ -836,8 +836,8 @@ elseif strcmp(mapping_for_tcv.method,'switchcase')
       edge_str_ = '_edge';
       edge_str_dot = '.edge';
     end
-    psi_max=tdi(['\results::thomson' edge_str_dot ':psi_max' substr_liuqe]);
-    psiscatvol=tdi(['\results::thomson' edge_str_dot ':psiscatvol' substr_liuqe]);
+    psi_max=gdat([],['\results::thomson' edge_str_dot ':psi_max' substr_liuqe]);
+    psiscatvol=gdat([],['\results::thomson' edge_str_dot ':psiscatvol' substr_liuqe]);
     if abs(zshift)>1e-5
       % calculate new psiscatvol
       psitdi=tdi(['\results::psi' substr_liuqe]);