diff --git a/matlab/TCV_IMAS/tcv_get_ids_core_sources.m b/matlab/TCV_IMAS/tcv_get_ids_core_sources.m
index 445fdfe3b97a44730994778fb67f9e94f1fb8a02..96e8def3e50a3ccb07ac6afba13406f48bfb9375 100644
--- a/matlab/TCV_IMAS/tcv_get_ids_core_sources.m
+++ b/matlab/TCV_IMAS/tcv_get_ids_core_sources.m
@@ -173,7 +173,9 @@ if ~isempty(ec_gdat.ec.data) % if EC data available, fill sources
   end
 
   % load LIUQE data to convert
-  % \tilde{j}_// = <jdotB>/(R0<Bphi/R>) to j_//0 = <jdotB>/B0; Bphi=T(psi)/R
+  % \tilde{V,TORAY}_// = dI_\phi/dV = (1/2\pi)<j_\phi/R> to j_//0 = <jdotB>/B0;
+  % conversion using j_//,cd = <j_cd/B> B, thus
+  % j_// = 2\pi <B^2> / (R_0 <B_\phi/R>^2)
   [L,~,LY] = liuqe(shot,ec_tgrid_toray,'iterq',50,'ilim',3,'icsint',true);
   rho_pol_norm_liu = L.pQ; % normalized sqrt poloidal flux including axis
   mu0     = 4*pi*10^-7;