diff --git a/TCV/loadTCVdata.m b/TCV/loadTCVdata.m
index 02cce761f9d80821d7ae86e3c43719e8aa436874..770bbf4802d4a97054e980b2cdddee43cc5b1349 100644
--- a/TCV/loadTCVdata.m
+++ b/TCV/loadTCVdata.m
@@ -160,7 +160,7 @@ end
 liuqe23=[{'\results::i_p'} {'\results::z_axis'} {'\results::r_axis'} {'\results::q_psi'} ...
       {'\results::beta_tor'} {'\results::beta_pol'} {'\results::q_95'} {'\results::l_i'}  {'\results::delta_95'} ...
       {'\results::kappa_95'} {'\results::r_contour'} {'\results::z_contour'} {'\results::psi_axis'} ...
-      {'\results::thomson:psiscatvol'} {'\results::thomson:psi_max'} {'\results::rms_error'}];
+      {'\results::thomson:psiscatvol'} {'\results::thomson:psi_max'} {'\results::rms_error'} {'\results::total_energy'}];
 
 % all keywords and corresponding case to run below
 TCVkeywrdall=[{'Ip'} {'zmag'} {'rmag'} {'rcont'} {'zcont'} {'vol'} {'rhovol'} {'qrho'} {'q95'} {'kappa'} ...
@@ -448,9 +448,11 @@ switch TCVkeywrdcase{index}
         tracefirrat=tdi('\results::thomson:fir_thom_rat');
         tracefirrat.dim{1}=time;
       end
-      itim=iround(time,tracefirrat.dim{1});
       tracefirrat_data=NaN*ones(size(tracetdi.dim{1}));
-      tracefirrat_data(itim)=tracefirrat.data;
+      if ~isempty(tracefirrat.data)
+        itim=iround(time,tracefirrat.dim{1});
+        tracefirrat_data(itim)=tracefirrat.data;
+      end
     else
       nodenameeff='\results::thomson:te';
       tracetdi=tdi(nodenameeff);