Skip to content
Snippets Groups Projects
Commit 8598a6e5 authored by Antonia Frank's avatar Antonia Frank
Browse files

Correct core_transport setup such that identifier only in description

parent 470a15a7
No related branches found
No related tags found
1 merge request!137Add quantities to ids for MRE
Hide whitespace changes
Inline Side-by-side
matlab/TCV_IMAS/tcv_get_ids_core_transport.m
+ 11
− 14
View file @ 8598a6e5
...@@ -28,17 +28,15 @@ try params_eff_ref=rmfield(params_eff_ref,'source');catch;end % make sure no sou ...@@ -28,17 +28,15 @@ try params_eff_ref=rmfield(params_eff_ref,'source');catch;end % make sure no sou
ids_core_transport_description = []; ids_core_transport_description = [];
%% %%
last_index = 0;
% setup empty model structs (up to 3 for now) % fill model [name: transport solver, desc: output from transport solver]
ids_core_transport.model(1:3) = ids_core_transport.model(1);
% fill model(3) [name: transport solver, desc: output from transport solver]
% setup model and profiles_1d % setup model and profiles_1d
comment = 'Output from a transport solver'; comment = 'Output from a transport solver';
ids_core_transport.model{3}.comment = comment; ids_core_transport.model{last_index+1}.comment = comment;
ids_core_transport.model{3}.name = comment; ids_core_transport.model{last_index+1}.name = comment;
ids_core_transport.model{3}.identifier.index = 2; ids_core_transport.model{last_index+1}.identifier.index = 2;
ids_core_transport.model{3}.name = 'transport_solver'; ids_core_transport.model{last_index+1}.name = 'transport_solver';
% read data and setup time % read data and setup time
params_eff.data_request='\results::conf:chie'; params_eff.data_request='\results::conf:chie';
...@@ -53,14 +51,13 @@ end ...@@ -53,14 +51,13 @@ end
ids_core_transport.time = temp_1d.chie.t; ids_core_transport.time = temp_1d.chie.t;
ids_core_transport_description.time = ['from node' params_eff.data_request]; ids_core_transport_description.time = ['from node' params_eff.data_request];
ids_core_transport.model{3}.profiles_1d(1:length(ids_core_transport.time)) = ... ids_core_transport.model{last_index+1}.profiles_1d(1:length(ids_core_transport.time)) = ...
ids_core_transport.model{3}.profiles_1d(1); ids_core_transport.model{last_index+1}.profiles_1d(1);
%% fill profiles_1d
% fill profiles_1d
for it=1:length(ids_core_transport.time) for it=1:length(ids_core_transport.time)
ids_core_transport.model{3}.profiles_1d{it}.time = ids_core_transport.time(it); ids_core_transport.model{last_index+1}.profiles_1d{it}.time = ids_core_transport.time(it);
ids_core_transport.model{3}.profiles_1d{it}.electrons.energy.d = temp_1d.chie.data(:,it); ids_core_transport.model{last_index+1}.profiles_1d{it}.electrons.energy.d = temp_1d.chie.data(:,it);
temp_1d_desc.electrons.energy.d = temp_1d.chie.label; temp_1d_desc.electrons.energy.d = temp_1d.chie.label;
end end
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment