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

Add parallel conductivity

parent a2fcdbf9
No related branches found
No related tags found
1 merge request!137Add quantities to ids for MRE
......@@ -359,6 +359,12 @@ for it=1:length(ids_core_profiles.time)
temp_1d_desc.j_ohmic = current_non_inductive.ohm.ohm_data.cd_dens.label;
end
%% parallel conductivity (neoclassical)
signeo = gdat(params_cores_profiles.shot,'\results::ibs:signeo');
for it=1:length(ids_core_profiles.time)
ids_core_profiles.profiles_1d{it}.conductivity_parallel = signeo.data(:,it);
temp_1d_desc.j_bootstrap = signeo.label;
end
%% add descriptions for profiles_1d
ids_core_profiles_description.profiles_1d = temp_1d_desc;
......
function [ids_equilibrium,ids_equilibrium_description,varargout] = tcv_get_ids_equilibrium(shot,ids_equil_empty, gdat_params,varargin)
function [ids_equilibrium,ids_equilibrium_description,varargout] = ...
tcv_get_ids_equilibrium(shot,ids_equil_empty,gdat_params,varargin)
%
% [ids_equilibrium,ids_equilibrium_description,varargout] = tcv_get_ids_equilibrium(shot,ids_equil_empty,varargin);
% [ids_equilibrium,ids_equilibrium_description,varargout] = ...
% tcv_get_ids_equilibrium(shot,ids_equil_empty,gdat_params,varargin);
%
%
% gdat_params: gdat_data.gdat_params to get all params passed from original call, in particular error_bar and cocos_out options
% gdat_params: gdat_data.gdat_params to get all params passed from original call,
% in particular error_bar and cocos_out options
%
if exist('gdat_params','var')
[ids_equilibrium, params_equilibrium] = tcv_ids_headpart(shot,ids_equil_empty,'equilibrium','gdat_params',gdat_params,varargin{:});
[ids_equilibrium, params_equilibrium] = ...
tcv_ids_headpart(shot,ids_equil_empty,'equilibrium','gdat_params',gdat_params,varargin{:});
else
[ids_equilibrium, params_equilibrium] = tcv_ids_headpart(shot,ids_equil_empty,'equilibrium',varargin{:});
[ids_equilibrium, params_equilibrium] = ...
tcv_ids_headpart(shot,ids_equil_empty,'equilibrium',varargin{:});
aa=gdat_tcv;
gdat_params = aa.gdat_params; % to get default params
end
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment