add comment

matlab/TCV/gdat_tcv.m

@@ -2770,7 +2770,7 @@ elseif strcmp(mapping_for_tcv.method,'switchcase')
         nodenameeff_rho = strrep(nodenameeff,'q_psi','rho');
         rho_eff = mdsvalue(nodenameeff_rho);
         if numel(tracetdi.dim{1}) == numel(rho_eff)
-          tracetdi.dim{1}(:) = rho_eff;
+          tracetdi.dim{1}(:) = rho_eff; % it is not sqrt(linspace...) anymore, hence take rho itself
         end
       end
     end