diff --git a/TCV/loadTCVdata.m b/TCV/loadTCVdata.m
index ce4b2c870030d7fd49991507978d42ebbbcb66d7..e1a81a3f2252287299e485ad2e088a139c495ab4 100644
--- a/TCV/loadTCVdata.m
+++ b/TCV/loadTCVdata.m
@@ -20,9 +20,9 @@
 % 'terho'= Te profile on (rho=sqrt(psi),time) mesh.Note rho is a 2D array as depends on time. ADD error bars in .std
 % 'profnerho' =  ne smoothed or fitted , vs (rho,t) (from Thomson fit)
 % 'profterho' =  te smoothed or fitted , vs (rho,t) (from Thomson fit)
-% 'neft' =  ne fitted from data on rho mesh (from proffit.local_time:neft)
+% 'neft' =  ne fitted from data on rho mesh (from proffit.local_time:neft_abs)
 % 'teft' =  te fitted from data on rho mesh (from proffit.local_time:teft)
-% 'neftav' =  ne fitted from averaged over time data on rho mesh (from proffit.avg_time:neft)
+% 'neftav' =  ne fitted from averaged over time data on rho mesh (from proffit.avg_time:neft_abs)
 % 'teftav' =  te fitted from averaged over time data on rho mesh (from proffit.avg_time:teft)
 % 'ece'  =  electron cyclotron emission
 % 'sxr'  =  soft x-ray emission
@@ -357,6 +357,7 @@ switch TCVkeywrdcase{index}
       nodenameeff='\results::thomson:te';
       tracetdi=tdi(nodenameeff);
       tracestd=tdi('\results::thomson:te:error_bar');
+      trace_abs=tdi('\results::thomson:fir_thom_rat');
     end
     trace.data=tracetdi.data'; % Thomson data as (t,z)
     trace.std=tracestd.data';