Skip to content
Snippets Groups Projects
Commit 65892284 authored by Olivier Sauter's avatar Olivier Sauter
Browse files

adapt mhd for aug and gdat_plot for labels

git-svn-id: https://spcsvn.epfl.ch/repos/TCV/gdat/trunk@7564 d63d8f72-b253-0410-a779-e742ad2e26cf
parent 47b7b5c0
No related branches found
No related tags found
No related merge requests found
...@@ -171,7 +171,7 @@ switch lower(data_request) ...@@ -171,7 +171,7 @@ switch lower(data_request)
mapping.expression = [{'FPG'},{'li'}]; mapping.expression = [{'FPG'},{'li'}];
case 'mhd' case 'mhd'
mapping.timedim = 1; mapping.timedim = 1;
mapping.label = 'Odd and Even n'; mapping.label = {'n\_odd','n\_even'};% this superseeds the expression
mapping.method = 'expression'; mapping.method = 'expression';
mapping.expression = ['params_eff = gdat_data.gdat_params;params_eff.data_request={''MOD'',''OddN''}; ' ... mapping.expression = ['params_eff = gdat_data.gdat_params;params_eff.data_request={''MOD'',''OddN''}; ' ...
'gdat_tmp=gdat_aug(shot,params_eff);gdat_tmp.data=reshape(gdat_tmp.data,length(gdat_tmp.data),1 );' ... 'gdat_tmp=gdat_aug(shot,params_eff);gdat_tmp.data=reshape(gdat_tmp.data,length(gdat_tmp.data),1 );' ...
......
...@@ -94,7 +94,9 @@ if all(isfield(gdat_data,{'data','t'})) && ~isempty(gdat_data.data) && ~isempty( ...@@ -94,7 +94,9 @@ if all(isfield(gdat_data,{'data','t'})) && ~isempty(gdat_data.data) && ~isempty(
imagesc(T+tmhdm(1),F/1e3,20*log10(abs(B)));axis xy;colormap jet; imagesc(T+tmhdm(1),F/1e3,20*log10(abs(B)));axis xy;colormap jet;
ylabel('freq') ylabel('freq')
xlabel(gdat_data.dimunits{1}) xlabel(gdat_data.dimunits{1})
title([upper(gdat_data.gdat_params.machine) '#' num2str(gdat_data.shot) ' ' gdat_data.label{i}]) ylabel_eff = gdat_data.label;
if iscell(gdat_data.label) && length(gdat_data.label)>=i; ylabel_eff = gdat_data.label{i}; end
title([upper(gdat_data.gdat_params.machine) '#' num2str(gdat_data.shot) ' ' ylabel_eff])
mhd_sum_data = mhd_sum_data + gdat_data.data(:,i); mhd_sum_data = mhd_sum_data + gdat_data.data(:,i);
end end
[B,F,T]=specgram(mhd_sum_data./size(gdat_data.data,2),nfft,1/mean(diff(gdat_data.t)),hanning(nfft),nfft/2); [B,F,T]=specgram(mhd_sum_data./size(gdat_data.data,2),nfft,1/mean(diff(gdat_data.t)),hanning(nfft),nfft/2);
...@@ -102,7 +104,9 @@ if all(isfield(gdat_data,{'data','t'})) && ~isempty(gdat_data.data) && ~isempty( ...@@ -102,7 +104,9 @@ if all(isfield(gdat_data,{'data','t'})) && ~isempty(gdat_data.data) && ~isempty(
imagesc(T+tmhdm(1),F/1e3,20*log10(abs(B)));axis xy;colormap jet; imagesc(T+tmhdm(1),F/1e3,20*log10(abs(B)));axis xy;colormap jet;
ylabel('freq') ylabel('freq')
xlabel(gdat_data.dimunits{1}) xlabel(gdat_data.dimunits{1})
title([upper(gdat_data.gdat_params.machine) '#' num2str(gdat_data.shot) ' sum of ' gdat_data.label{:}]) ylabel_eff = gdat_data.label;
if iscell(gdat_data.label); ylabel_eff = sprintf('%s ',gdat_data.label{:}); end
title([upper(gdat_data.gdat_params.machine) '#' num2str(gdat_data.shot) ' sum of ' ylabel_eff])
end end
else else
disp('cannot plot gdat_data, has empty data or t field') disp('cannot plot gdat_data, has empty data or t field')
......
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