diff --git a/crpptbx/TCV/gdat_tcv.m b/crpptbx/TCV/gdat_tcv.m
index ecebc95f05cc64893563dd44e0f33f3e04e66ba5..32ee5b0ec47589824fcdb6cf342b9fc7ea7602a6 100644
--- a/crpptbx/TCV/gdat_tcv.m
+++ b/crpptbx/TCV/gdat_tcv.m
@@ -901,8 +901,8 @@ elseif strcmp(mapping_for_tcv.method,'switchcase')
       edge_str_ = '_edge';
       edge_str_dot = '.edge';
     end
-    psi_max=gdat([],['\results::thomson' edge_str_dot ':psi_max' substr_liuqe]);
-    psiscatvol=gdat([],['\results::thomson' edge_str_dot ':psiscatvol' substr_liuqe]);
+    psi_max=gdat([],['\results::thomson' edge_str_dot ':psi_max' substr_liuqe],'nverbose',gdat_params.nverbose);
+    psiscatvol=gdat([],['\results::thomson' edge_str_dot ':psiscatvol' substr_liuqe],'nverbose',gdat_params.nverbose);
     if abs(zshift)>1e-5
       % calculate new psiscatvol
       psitdi=tdi(['\results::psi' substr_liuqe]);
@@ -936,7 +936,9 @@ elseif strcmp(mapping_for_tcv.method,'switchcase')
         rho(ir,:)= sqrt(1.-psiscatvol.data(:,ir)./psi_max.data(:))';
       end
     else
-      if gdat_params.nverbose>=1; warning(['psiscatvol empty?, no rho calculated for data_request = ' data_request_eff]); end
+      if gdat_params.nverbose>=1 && gdat_data.gdat_params.edge==0
+        warning(['psiscatvol empty?, no rho calculated for data_request = ' data_request_eff]);
+      end
       rho=[];
     end
     gdat_data.dim{1}=rho;
@@ -1825,7 +1827,7 @@ function [gdat_data] = get_grids_1d(gdat_data,nbdim_x,nopt,nverbose);
 %
 
 gdat_data.grids_1d.rhopolnorm = gdat_data.x;
-if (nopt == 0)
+if (nopt == 0) || isempty(gdat_data.x) || isempty(gdat_data.t) || isempty(gdat_data.data) || ischar(gdat_data.data)
   gdat_data.grids_1d.rhotornorm = [];
   gdat_data.grids_1d.rhovolnorm = [];
   gdat_data.grids_1d.psi = [];
@@ -1839,6 +1841,7 @@ params_eff.data_request='rhotor';
 rhotor = gdat(gdat_data.shot,params_eff);
 params_eff.data_request='rhovol';
 rhovol = gdat(gdat_data.shot,params_eff);
+
 psi0 = interpos(rhotor.t',rhotor.psi_axis,gdat_data.t,-0.01);
 if (nbdim_x == 1)
   gdat_data.grids_1d.psi = gdat_data.grids_1d.rhopolnorm.^2*reshape(psi0,1,length(psi0));