make empty cells as {}

git-svn-id: https://spcsvn.epfl.ch/repos/TCV/gdat/trunk@11661 d63d8f72-b253-0410-a779-e742ad2e26cf

crpptbx/TCV_IMAS/tcv2ids2database.m

@@ -13,7 +13,7 @@ ids2get = {'equilibrium', 'magnetics', 'pf_active','wall'}; % default will load
 %% can plot with: [plotids_H] = plotids([],idsok.ids{1},idsok.ids{2});
 if isfield(ids_from_tcv,'tf')
   % avoid present problem of not saving data if ids_from_tcv.tf.field_map{1}.time is not set
-  if abs(ids_from_tcv.tf.field_map{1}.time+9.e40)<1e-5
+  if length(ids_from_tcv.tf.field_map)>0 && abs(ids_from_tcv.tf.field_map{1}.time+9.e40)<1e-5
     ids_from_tcv.tf.field_map{1}.time = -8.9000e+40;
     disp(['WARNING: ids_from_tcv.tf.field_map{1}.time changed to ' num2str(ids_from_tcv.tf.field_map{1}.time) ' to avoid writing error']);
   end
@@ -21,3 +21,11 @@ end
 %%
 [ids_put_status] = ids2database(shot,run_out,occurence,ids_from_tcv);
 
+
+% [plotids_H] = plotids([],idsok.ids);
+
+% tic;[TTTe,idse]=evalc('ids_get(expIdx, ''equilibrium'')');toc
+% tic;[TTTm,idsm]=evalc('ids_get(expIdx, ''magnetics'')');toc
+% tic;[TTTp,idsp]=evalc('ids_get(expIdx, ''pf_active'')');toc
+% tic;[TTTt,idst]=evalc('ids_get(expIdx, ''tf'')');toc
+% tic;[TTTw,idsw]=evalc('ids_get(expIdx, ''wall'')');toc

crpptbx/TCV_IMAS/tcv_get_ids_equilibrium.m

@@ -294,3 +294,10 @@ for it=1:length(ids_equilibrium.time)
   ids_equilibrium.time_slice{it}.profiles_2d{1}.psi(:,:) = ids_equilibrium.time_slice{it}.profiles_2d{1}.psi(:,:) + ...
       global_quantities.psi_boundary.data(it);
 end
+
+
+% make arrays not filled in empty:
+ids_equilibrium.grids_ggd = {};
+for it=1:length(ids_equilibrium.time_slice)
+  ids_equilibrium.time_slice{it}.ggd = {};
+end

crpptbx/TCV_IMAS/tcv_get_ids_tf.m

@@ -22,3 +22,7 @@ else
   % leave defaults
   ids_tf_out_description = struct([]);
 end
+
+% make arrays not filled in empty: not the case for magnetics
+ids_tf_out.coil = {};
+ids_tf_out.field_map = {};

crpptbx/TCV_IMAS/tcv_get_ids_wall.m

 function [ids_wall,ids_description,varargout] = tcv_get_ids_wall(shot, ids_wall_empty, varargin)
 
 % Input pharser
-[ids_wall, params] = tcv_ids_headpart(shot, ids_wall_empty,'wall',varargin{:});
+[ids_wall, params] = tcv_ids_headpart(shot, ids_wall_empty,'wall','homogeneous_time',0,varargin{:});
 
 % Get data of outline
 [ids_wall.description_2d]= tcv_get_ids_descprition_2d(params.shot, ids_wall.description_2d(1));
 
 % make arrays not filled in empty: not the case for magnetics
-% put does not work: ids_wall.global_quantities = struct([]);
-% put does not work: ids_wall.global_quantities.electrons = struct([]);
 ids_wall.global_quantities.neutral = {}; % wall does not work with Antoine's addpath
 ids_wall.description_ggd = {};
-% put does not work: ids_wall.description_2d{1}.type = struct([]);
+ids_wall.description_2d{1}.mobile.unit = {};
+ids_wall.description_2d{1}.vessel.unit = {};
 
 %% TODO 
 ids_description = [];