small changes, working get_core_sources now

matlab/TCV_IMAS/tcv_get_ids_core_sources.m

@@ -32,7 +32,7 @@ ids_core_sources_description = [];
 source_template = ids_core_sources.source{1};
 profiles_template = source_template.profiles_1d{1};
 globals_template = source_template.global_quantities{1};
-%% get time grid and data form gdat & initialize source from template
+%% get time grid and data from gdat & initialize source from template
 ohm_gdat = gdat(shot,'ohm_data');
 % fix time grid to ohm time grid, same as bs
 t_grid = ohm_gdat.t; n_t = numel(t_grid); t_spacing = mean(diff(t_grid));
@@ -54,12 +54,13 @@ ids_core_sources.source{3}.identifier = id_ec;
 ids_core_sources.source{3}.profiles_1d(1:n_t) = {profiles_template};
 ids_core_sources.source{3}.global_quantities(1:n_t) = {globals_template};
 % load data
+%nb_launchers = size(aa.ec.ec_data.rho_max_pow_dens.data,1);
 ec_gdat = gdat(shot,'ec_data'); ec_data = ec_gdat.ec.ec_data;
-ec_time = ec_gdat.ec.ec_data.cd_dens.t;
+ec_time = ec_data.cd_dens.t;
 % interpolate totals on t_grid
-ec_total_pow = interp1(ec_data.p_abs_plasma.t,ec_data.p_abs_plasma.data(12,:),t_grid);
+ec_total_pow = interp1(ec_data.p_abs_plasma.t,ec_data.p_abs_plasma.data(10,:),t_grid); %the index 10 works for shot 56044
 ec_total_pow(isnan(ec_total_pow)) = 0;
-ec_total_cur = interp1(ec_data.p_abs_plasma.t,ec_data.cd_tot.data(12,:),t_grid);
+ec_total_cur = interp1(ec_data.cd_tot.t,ec_data.cd_tot.data(10,:),t_grid);
 ec_total_cur(isnan(ec_total_cur)) = 0;
 
 % bs
@@ -72,13 +73,24 @@ ids_core_sources.source{13}.global_quantities(1:n_t) = {globals_template};
 % load data
 bs_gdat = gdat(shot,'bs_data'); bs_data = bs_gdat.bs.bs_data;
 
-% nbi
+% nbi  
+id_nbi.description = 'Source from Neutral Beam Injection';
+id_nbi.index = 2; id_nbi.name = 'nbi';
 ids_core_sources.source{2} = source_template;
-% to be filled
+ids_core_sources.source{2}.identifier = id_nbi;
+ids_core_sources.source{2}.profiles_1d(1:n_t) = {profiles_template};
+ids_core_sources.source{2}.global_quantities(1:n_t) = {globals_template};
 
-% total
+
+% total  
+id_total.description = 'Total source; combines all sources';
+id_total.index = 1; id_total.name = 'total';
 ids_core_sources.source{1} = source_template;
-% to be filled
+ids_core_sources.source{1}.identifier = id_total;
+ids_core_sources.source{1}.profiles_1d(1:n_t) = {profiles_template};
+ids_core_sources.source{1}.global_quantities(1:n_t) = {globals_template};
+
+
 powers_gdat = gdat(shot,'powers');
 
 %% fill profiles for times n t_grid
@@ -93,7 +105,7 @@ for ii = 1:n_t
   ids_core_sources.source{7}.global_quantities{ii}.power = powers_gdat.ohm.data(ii);
     
   % ec, need to find matching time index\
-  ids_core_sources.source{3}.profiles_1d{ii}.time = t_grid(ii); % profiles time
+  ids_core_sources.source{3}.profiles_1d{ii}.time = t_grid(ii); % profiles time (should profiles time be ec_time?)
   ids_core_sources.source{3}.global_quantities{ii}.time = t_grid(ii); % globals time
   [~, ec_idx] = min(abs(ec_time-t_grid(ii)));
   if ec_time(ec_idx) == t_grid(ii) % same grid