From 36d4bfe56646366714bea7299ee7a2ef105b7a8a Mon Sep 17 00:00:00 2001
From: Olivier Sauter <Olivier.Sauter@epfl.ch>
Date: Mon, 17 Jan 2022 17:19:15 +0100
Subject: [PATCH] try to finalize tcv_get_ids_equilibrium changes

to be tested
---
 matlab/TCV_IMAS/tcv_get_ids_equilibrium.m | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/matlab/TCV_IMAS/tcv_get_ids_equilibrium.m b/matlab/TCV_IMAS/tcv_get_ids_equilibrium.m
index 54f872b8..10dcc0cc 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);
-- 
GitLab