diff --git a/crpptbx/TCV_IMAS/tcv_get_ids_equilibrium.m b/crpptbx/TCV_IMAS/tcv_get_ids_equilibrium.m
index 6d30400d278e9460cba418ed1dcd2b54dc0a54b5..c9bdb3dafadda2dcb416810cc51454f84ba1f282 100644
--- a/crpptbx/TCV_IMAS/tcv_get_ids_equilibrium.m
+++ b/crpptbx/TCV_IMAS/tcv_get_ids_equilibrium.m
@@ -26,14 +26,14 @@ ids_equilibrium_description.vacuum_toroidal_field = vacuum_toroidal_field_desc;
 ids_equilibrium.time = vacuum_toroidal_field.b0.t;
 ids_equilibrium_description.time = '.t subfield from: [''b0'',''source'',''liuqe'']';
 
-ids_equilibrium.time_slice(1:length(ids_equilibrium.time)) = ids_equilibrium.time_slice(1);
+ids_equilibrium.time_slice(1:numel(ids_equilibrium.time)) = ids_equilibrium.time_slice(1);
 
 % load time array data to copy to time_slices
 
 % global_quantities data into local global_quantities.* structure with correct end names and global_quantities_desc.* with description. Use temp.* and temp_desc.* structures for temporary data
 
 % brute force solution load all eqdsks
-% $$$ for it=1:length(ids_equilibrium.time)
+% $$$ for it=1:numel(ids_equilibrium.time)
 % $$$   ids_equilibrium.time(it)
 % $$$   temp.eqdsks{it}=gdat(params_equilibrium.shot,'eqdsk','time',ids_equilibrium.time(it),'write',0,'machine',machine);
 % $$$ end
@@ -76,8 +76,8 @@ global_quantities_desc.w_mhd = 'w_mhd';
 
 global_quantities_fieldnames = fieldnames(global_quantities);
 special_fields = {'magnetic_axis', 'psi_axis', 'q_min'}; % fields needing non-automatic treatments
-for it=1:length(ids_equilibrium.time)
-  for i=1:length(global_quantities_fieldnames)
+for it=1:numel(ids_equilibrium.time)
+  for i=1:numel(global_quantities_fieldnames)
     if ~any(strcmp(global_quantities_fieldnames{i},special_fields))
       if ~isstruct(ids_equilibrium.time_slice{it}.global_quantities.(global_quantities_fieldnames{i}))
         ids_equilibrium.time_slice{it}.global_quantities.(global_quantities_fieldnames{i}) = ...
@@ -90,7 +90,7 @@ for it=1:length(ids_equilibrium.time)
 end
 
 % special case
-for it=1:length(ids_equilibrium.time)
+for it=1:numel(ids_equilibrium.time)
   ids_equilibrium.time_slice{it}.global_quantities.magnetic_axis.r = temp.r_magnetic_axis.data(it);
   ids_equilibrium.time_slice{it}.global_quantities.magnetic_axis.z = temp.z_magnetic_axis.data(it);
   ids_equilibrium.time_slice{it}.global_quantities.psi_axis = temp.psi_axis.data(it) + ...
@@ -136,8 +136,8 @@ temp_desc.z_lcfs = 'z_contour_edge';
 
 boundary_fieldnames = fieldnames(boundary);
 special_fields = {'lcfs', 'geometric_axis', 'x_point'}; % fields needing non-automatic treatments
-for it=1:length(ids_equilibrium.time)
-  for i=1:length(boundary_fieldnames)
+for it=1:numel(ids_equilibrium.time)
+  for i=1:numel(boundary_fieldnames)
     if ~any(strcmp(boundary_fieldnames{i},special_fields))
       if ~isstruct(ids_equilibrium.time_slice{it}.boundary.(boundary_fieldnames{i}))
         ids_equilibrium.time_slice{it}.boundary.(boundary_fieldnames{i}) = ...
@@ -150,7 +150,7 @@ for it=1:length(ids_equilibrium.time)
 end
 
 % special cases
-for it=1:length(ids_equilibrium.time)
+for it=1:numel(ids_equilibrium.time)
   ids_equilibrium.time_slice{it}.boundary.outline.r = temp.r_lcfs.data(:,it);
   ids_equilibrium.time_slice{it}.boundary.outline.z = temp.z_lcfs.data(:,it);
   ids_equilibrium.time_slice{it}.boundary.lcfs.r = ids_equilibrium.time_slice{it}.boundary.outline.r;
@@ -159,10 +159,12 @@ for it=1:length(ids_equilibrium.time)
   ids_equilibrium.time_slice{it}.boundary.geometric_axis.z = temp.zgeom.data(it);
   if temp.n_x_point.data(it) > 0
     ids_equilibrium.time_slice{it}.boundary.x_point(1:temp.n_x_point.data(it)) = ids_equilibrium.time_slice{it}.boundary.x_point(1);
-    for i=1:length(temp.n_x_point.data(it))
+    for i=1:temp.n_x_point.data(it)
       ids_equilibrium.time_slice{it}.boundary.x_point{i}.r = temp.r_x_point.data(i,it);
       ids_equilibrium.time_slice{it}.boundary.x_point{i}.z = temp.z_x_point.data(i,it);
     end
+  else
+    ids_equilibrium.time_slice{it}.boundary.x_point = {};
   end
 end
 
@@ -221,8 +223,8 @@ profiles_1d.volume = gdat(params_equilibrium.shot,'volume_rho','machine',machine
 
 profiles_1d_fieldnames = fieldnames(profiles_1d);
 special_fields = {'geometric_axis', 'rho_tor_norm', 'psi'}; % fields needing non-automatic treatments
-for it=1:length(ids_equilibrium.time)
-  for i=1:length(profiles_1d_fieldnames)
+for it=1:numel(ids_equilibrium.time)
+  for i=1:numel(profiles_1d_fieldnames)
     if ~any(strcmp(profiles_1d_fieldnames{i},special_fields))
       if ~isstruct(ids_equilibrium.time_slice{it}.profiles_1d.(profiles_1d_fieldnames{i}))
         if ~ischar(profiles_1d.(profiles_1d_fieldnames{i}).data) && ~isempty(profiles_1d.(profiles_1d_fieldnames{i}).data) ...
@@ -238,9 +240,9 @@ for it=1:length(ids_equilibrium.time)
 end
 
 % special cases
-nrho = length(profiles_1d.rho_tor.x);
-ntime = length(temp.psi_axis.data);
-for it=1:length(ids_equilibrium.time)
+nrho = numel(profiles_1d.rho_tor.x);
+ntime = numel(temp.psi_axis.data);
+for it=1:numel(ids_equilibrium.time)
   ids_equilibrium.time_slice{it}.profiles_1d.rho_tor_norm = ids_equilibrium.time_slice{it}.profiles_1d.rho_tor./ ...
       ids_equilibrium.time_slice{it}.profiles_1d.rho_tor(end);
   ids_equilibrium.time_slice{it}.profiles_1d.psi = temp.psi_axis.data(it) + ...
@@ -270,8 +272,8 @@ profiles_2d.psi = gdat(params_equilibrium.shot,'psi','machine',machine); % add p
 
 profiles_2d_fieldnames = fieldnames(profiles_2d);
 special_fields = {'grid', 'grid_type'}; % fields needing non-automatic treatments
-for it=1:length(ids_equilibrium.time)
-  for i=1:length(profiles_2d_fieldnames)
+for it=1:numel(ids_equilibrium.time)
+  for i=1:numel(profiles_2d_fieldnames)
     if ~any(strcmp(profiles_2d_fieldnames{i},special_fields))
       if ~isstruct(ids_equilibrium.time_slice{it}.profiles_2d{1}.(profiles_2d_fieldnames{i}))
         if ~ischar(profiles_2d.(profiles_2d_fieldnames{i}).data) && ~isempty(profiles_2d.(profiles_2d_fieldnames{i}).data) ...
@@ -287,7 +289,7 @@ for it=1:length(ids_equilibrium.time)
 end
 
 % special cases
-for it=1:length(ids_equilibrium.time)
+for it=1:numel(ids_equilibrium.time)
   ids_equilibrium.time_slice{it}.profiles_2d{1}.grid_type.name = profiles_2d.grid_type.name;
   ids_equilibrium.time_slice{it}.profiles_2d{1}.grid_type.index = profiles_2d.grid_type.index;
   ids_equilibrium.time_slice{it}.profiles_2d{1}.grid_type.description = profiles_2d.grid_type.description;
@@ -300,6 +302,6 @@ end
 
 % make arrays not filled in empty:
 ids_equilibrium.grids_ggd = {};
-for it=1:length(ids_equilibrium.time_slice)
+for it=1:numel(ids_equilibrium.time_slice)
   ids_equilibrium.time_slice{it}.ggd = {};
 end