diff --git a/matlab/TCV_IMAS/tcv_get_ids_equilibrium.m b/matlab/TCV_IMAS/tcv_get_ids_equilibrium.m
index 54f872b86726933316a3eec34036fc6a2bf403f8..10dcc0ccda6199e9380f4d28212375b4eff7e0f6 100644
--- a/matlab/TCV_IMAS/tcv_get_ids_equilibrium.m
+++ b/matlab/TCV_IMAS/tcv_get_ids_equilibrium.m
@@ -340,7 +340,7 @@ for it = 1:ntime
     end
   end
   constraints.pf_current = pf_current;
-  
+
   ids_equilibrium.time_slice{it}.constraints = constraints;
 end
 
@@ -518,9 +518,9 @@ params_eff.data_request = 'psi';
 profiles_2d.psi = gdat(params_equilibrium.shot,params_eff); % add psi_bound in a second step in special cases
 profiles_2d_desc.psi = [params_eff.data_request ' adding psi_bound in a 2nd step'];
 %profiles_2d.r = profiles_2d.psi;
-profiles_2d.r.data = repmat(repmat(profiles_2d.psi.dim{1},1,65),1,1,1299);
+profiles_2d.r.data = repmat(repmat(profiles_2d.psi.dim{1},1,numel(profiles_2d.psi.dim{2})),1,1,numel(profiles_2d.psi.dim{3}));
 profiles_2d_desc.r = 'from dim{1} of ''psi'' repeated';
-profiles_2d.z.data = repmat(repmat(profiles_2d.psi.dim{2}',28,1),1,1,1299);
+profiles_2d.z.data = repmat(repmat(profiles_2d.psi.dim{2}',numel(profiles_2d.psi.dim{1}),1),1,1,numel(profiles_2d.psi.dim{3}));
 profiles_2d_desc.z = 'from dim{2} of ''psi'' repeated';
 
 % theta = gdat(params_equilibrium.shot,'theta','machine',gdat_params.machine);