diff --git a/crpptbx/TCV_IMAS/tcv_get_ids_equilibrium.m b/crpptbx/TCV_IMAS/tcv_get_ids_equilibrium.m
index 3e59d78d330f285d14bda8c1dbd1388f898b7b59..af0e8a6fc8b85807d5bc955dfa3599e687a7290b 100644
--- a/crpptbx/TCV_IMAS/tcv_get_ids_equilibrium.m
+++ b/crpptbx/TCV_IMAS/tcv_get_ids_equilibrium.m
@@ -266,7 +266,7 @@ profiles_1d.j_tor.data = - 2.*pi.* j_tor; % 2pi sigma_bp * jtor above (Eq. (30)
 j_par = - profiles_1d.f_df_dpsi.data./profiles_1d.f.data./mu0.*gradpsi_over_R_sq./2./pi ...
         - profiles_1d.f.data .*2*pi .*profiles_1d.dpressure_dpsi.data ...
         - profiles_1d.f.data .*2*pi /mu0.*profiles_1d.f_df_dpsi.data.*profiles_1d.gm1.data;
-profiles_1d.j_parallel.data = j_par./repmat(ids_equilibrium.vacuum_toroidal_field.b0.data',size(profiles_1d.f.data,1),1);
+profiles_1d.j_parallel.data = j_par./repmat(ids_equilibrium.vacuum_toroidal_field.b0',size(profiles_1d.f.data,1),1);
 
 profiles_1d_fieldnames = fieldnames(profiles_1d);
 special_fields = {'geometric_axis', 'rho_tor_norm', 'psi'}; % fields needing non-automatic treatments