Skip to content
Snippets Groups Projects
Commit 51de653b authored by Michele Marin's avatar Michele Marin
Browse files

adding work in progress sections

parent 0580b163
Branches
Tags
1 merge request!184Adding errorbars for ion temperature and impurity density
...@@ -167,10 +167,10 @@ end ...@@ -167,10 +167,10 @@ end
% Need to check if the temperature is local or averaged. Best way is to check the data themselves % Need to check if the temperature is local or averaged. Best way is to check the data themselves
te_local_time = gdat(params_cores_profiles.shot, /'proffit:local_time:teft'); te_local_time = gdat(params_cores_profiles.shot, '\proffit:local_time:teft');
te_avg_time = gdat(params_cores_profiles.shot, /'proffit:avg_time:teft'); te_avg_time = gdat(params_cores_profiles.shot, '\proffit:avg_time:teft');
te_conf = gdat(params_cores_profiles.shot, /'results::conf:te'); te_conf = gdat(params_cores_profiles.shot, '\results::conf:te');
% Getting the call done to get the data in the description % Getting the call done to get the data in the description
...@@ -185,7 +185,9 @@ te_profit_standard_call = te_profit_calls{2}; ...@@ -185,7 +185,9 @@ te_profit_standard_call = te_profit_calls{2};
[aa, bb] = hldsi(shot); [aa, bb] = hldsi(shot);
aa
bb
return
...@@ -416,31 +418,31 @@ if ~isempty(temp_1d.cxrs_rho.data) && ~isempty(temp_1d.cxrs_rho.ti.fit.data) ...@@ -416,31 +418,31 @@ if ~isempty(temp_1d.cxrs_rho.data) && ~isempty(temp_1d.cxrs_rho.ti.fit.data)
end end
time_measurement_slice_method %time_measurement_slice_method
time_measurement_slice_method.name='exact' or 'local' %time_measurement_slice_method.name='exact' or 'local'
time_measurement_slice_method.index=0 %time_measurement_slice_method.index=0
time_measurement_slice_method.description='using Thomson scattering measurement at the given time' %time_measurement_slice_method.description='using Thomson scattering measurement at the given time'
time_measurement_width %time_measurement_width
time_measurement_slice_method %time_measurement_slice_method
time_measurement_slice_method.name='exact' or 'local' %time_measurement_slice_method.name='exact' or 'local'
time_measurement_slice_method.index=0 %time_measurement_slice_method.index=0
time_measurement_slice_method.description='using Thomson scattering measurement at the given time' %time_measurement_slice_method.description='using Thomson scattering measurement at the given time'
local = 1 %local = 1
chi_squared = \tcv_shot::top.results.proffit.avg_time:teft_chi2 %chi_squared = \tcv_shot::top.results.proffit.avg_time:teft_chi2
time_measurement_width(:) = gdat(84927,'\tcv_shot::top.results.proffit.avg_time:time_avg','doplot',1); % nargout = 1 %time_measurement_width(:) = gdat(84927,'\tcv_shot::top.results.proffit.avg_time:time_avg','doplot',1); % nargout = 1
% I am guessing that in gdat params_eff.fit_type = 'avg' takes the data from the avg_time and not local_time proffit nodes? % I am guessing that in gdat params_eff.fit_type = 'avg' takes the data from the avg_time and not local_time proffit nodes?
name='exact' or 'local' %name='exact' or 'local'
index=0 %index=0
description='using Thomson scattering measurement at the given time' %description='using Thomson scattering measurement at the given time'
(if only one TS laser time was used ) %(if only one TS laser time was used )
or may be: %or may be:
name='avg' or 'average' %name='avg' or 'average'
index=-1 %index=-1
description='average with TS at this time and other TS within +-delta_t' (delta_t is provided in the proffit:local tree), conf should know if nodes from local or avg %description='average with TS at this time and other TS within +-delta_t' (delta_t is provided in the proffit:local tree), conf should know if nodes from local or avg
if ~isempty(temp_1d.cxrs_rho.error_bar) if ~isempty(temp_1d.cxrs_rho.error_bar)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment