diff --git a/matlab/TCV_IMAS/tcv_get_ids_core_transport.m b/matlab/TCV_IMAS/tcv_get_ids_core_transport.m index 7715f615096ce14aca1d247b31e9fe3047c35a55..ed3abe093c893322fa5fdb0b2f493cae6d582661 100644 --- a/matlab/TCV_IMAS/tcv_get_ids_core_transport.m +++ b/matlab/TCV_IMAS/tcv_get_ids_core_transport.m @@ -28,17 +28,15 @@ try params_eff_ref=rmfield(params_eff_ref,'source');catch;end % make sure no sou ids_core_transport_description = []; %% +last_index = 0; -% setup empty model structs (up to 3 for now) -ids_core_transport.model(1:3) = ids_core_transport.model(1); - -% fill model(3) [name: transport solver, desc: output from transport solver] +% fill model [name: transport solver, desc: output from transport solver] % setup model and profiles_1d comment = 'Output from a transport solver'; -ids_core_transport.model{3}.comment = comment; -ids_core_transport.model{3}.name = comment; -ids_core_transport.model{3}.identifier.index = 2; -ids_core_transport.model{3}.name = 'transport_solver'; +ids_core_transport.model{last_index+1}.comment = comment; +ids_core_transport.model{last_index+1}.name = comment; +ids_core_transport.model{last_index+1}.identifier.index = 2; +ids_core_transport.model{last_index+1}.name = 'transport_solver'; % read data and setup time params_eff.data_request='\results::conf:chie'; @@ -53,14 +51,13 @@ end ids_core_transport.time = temp_1d.chie.t; 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{3}.profiles_1d(1); - -%% fill profiles_1d +ids_core_transport.model{last_index+1}.profiles_1d(1:length(ids_core_transport.time)) = ... + ids_core_transport.model{last_index+1}.profiles_1d(1); +% fill profiles_1d 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{3}.profiles_1d{it}.electrons.energy.d = temp_1d.chie.data(:,it); + ids_core_transport.model{last_index+1}.profiles_1d{it}.time = ids_core_transport.time(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; end