Newer
Older
Olivier Sauter
committed
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
function [torbeam_std_files_struct] = torbeam_plot_output_singlerun;
%
% [torbeam_std_files_struct] = torbeam_plot_output_singlerun;
%
% reads local topfile, Te.dat, ne.dat input files
% as well as outputs files 't1_LIB.dat', 't1tor_LIB.dat', 'volumes.dat', 't2_new_LIB.dat','cntpr.dat'
% and plot results
%
% read topfile
fid=fopen('topfile','r');
% fscanf(fid,'%s %d %.6f\n',' Number of radial and vertical grid points in ',shot,time_for_torbeam(it));
dummy=fgetl(fid);
dummy=fgetl(fid);
dummyaaa = sscanf(dummy,'%d',2);
zzz.nb_rmesh = dummyaaa(1);
zzz.nb_zmesh = dummyaaa(2);
% fscanf(fid,'%s\n','Inside and Outside radius and psi_sep');
dummy=fgetl(fid);
dummy=fgetl(fid);
dummyaaa = sscanf(dummy,'%f',3);
zzz.rmesh_min = dummyaaa(1);
zzz.rmesh_max = dummyaaa(2);
zzz.psi_edge = dummyaaa(3);
% fscanf(fid,'%s\n','Radial grid coordinates');
dummy=fgetl(fid);
zzz.rmesh = fscanf(fid,'%f',zzz.nb_rmesh);
% fscanf(fid,'%s\n','Vertical grid coordinates');
dummy=fgetl(fid);
if isempty(deblank(dummy)); dummy=fgetl(fid); end
zzz.zmesh = fscanf(fid,'%f',zzz.nb_zmesh);
% fscanf(fid,'%s\n','B_r values');
dummy=fgetl(fid);
if isempty(deblank(dummy)); dummy=fgetl(fid); end
zzz.BR = fscanf(fid,'%f',[zzz.nb_rmesh,zzz.nb_zmesh]);
% fscanf(fid,'%s\n','B_t values');
dummy=fgetl(fid);
if isempty(deblank(dummy)); dummy=fgetl(fid); end
zzz.Bphi = fscanf(fid,'%f',[zzz.nb_rmesh,zzz.nb_zmesh]);
% fscanf(fid,'%s\n','B_z values');
dummy=fgetl(fid);
if isempty(deblank(dummy)); dummy=fgetl(fid); end
zzz.BZ= fscanf(fid,'%f',[zzz.nb_rmesh,zzz.nb_zmesh]);
% fscanf(fid,'%s\n','psi values');
dummy=fgetl(fid);
if isempty(deblank(dummy)); dummy=fgetl(fid); end
zzz.psi= fscanf(fid,'%f',[zzz.nb_rmesh,zzz.nb_zmesh]);
fclose(fid);
zzz.prof_to_cp = {'Te.dat','ne.dat'};
zzz.prof_to_cp_var = lower(strrep(zzz.prof_to_cp,'.dat',''));
for i=1:length(zzz.prof_to_cp_var)
zzz.(zzz.prof_to_cp_var{i}) = [];
try
[zzz.(zzz.prof_to_cp_var{i})(:,1),zzz.(zzz.prof_to_cp_var{i})(:,2)] = textread(zzz.prof_to_cp{i},'%f%f','headerlines',1);
catch
end
end
zzz.tb_out_files_to_cp = {'t1_LIB.dat', 't1tor_LIB.dat', 'volumes.dat', 't2_new_LIB.dat'};
zzz.tb_out_files_to_cp_var = strrep(zzz.tb_out_files_to_cp,'.dat','');
for i=1:length(zzz.tb_out_files_to_cp_var)
zzz.(zzz.tb_out_files_to_cp_var{i}) = [];
try
zzz.(zzz.tb_out_files_to_cp_var{i}) = load(zzz.tb_out_files_to_cp{i});
catch
end
end
try
zzz.tbeg_abs=load('cntpr.dat');
zzz.tbeg_abs(2) = zzz.tbeg_abs(2) + 1; % counted from 0?
catch
warning('no cntpr.dat?')
end
Olivier Sauter
committed
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
% plots
figure
contour(zzz.rmesh,zzz.zmesh,zzz.BR',100);
hold on
contour(zzz.rmesh,zzz.zmesh,zzz.psi',[zzz.psi_edge zzz.psi_edge],'linewidth',2,'linecolor','k');
for i=1:2:size(zzz.t1_LIB,2)
hpl=plot(zzz.t1_LIB(:,i)/1e2,zzz.t1_LIB(:,i+1)/1e2);
plotos(zzz.t1_LIB(zzz.tbeg_abs(2):sum(zzz.tbeg_abs),i)/1e2,zzz.t1_LIB(zzz.tbeg_abs(2):sum(zzz.tbeg_abs),i+1)/1e2,'-', ...
[7 0],[],'k');%get(hpl,'color'));
end
axis equal
colorbar
title('BR')
figure
contour(zzz.rmesh,zzz.zmesh,zzz.BZ',100);
hold on
contour(zzz.rmesh,zzz.zmesh,zzz.psi',[zzz.psi_edge zzz.psi_edge],'linewidth',2,'linecolor','k');
for i=1:2:size(zzz.t1_LIB,2)
hpl=plot(zzz.t1_LIB(:,i)/1e2,zzz.t1_LIB(:,i+1)/1e2);
plotos(zzz.t1_LIB(zzz.tbeg_abs(2):sum(zzz.tbeg_abs),i)/1e2,zzz.t1_LIB(zzz.tbeg_abs(2):sum(zzz.tbeg_abs),i+1)/1e2,'-', ...
[7 0],[],'k');%get(hpl,'color'));
end
axis equal
colorbar
title('BZ')
figure
contour(zzz.rmesh,zzz.zmesh,zzz.Bphi',100);
hold on
contour(zzz.rmesh,zzz.zmesh,zzz.psi',[zzz.psi_edge zzz.psi_edge],'linewidth',2,'linecolor','k');
for i=1:2:size(zzz.t1_LIB,2)
hpl=plot(zzz.t1_LIB(:,i)/1e2,zzz.t1_LIB(:,i+1)/1e2);
ixx= plotos(zzz.t1_LIB(zzz.tbeg_abs(2):sum(zzz.tbeg_abs),i)/1e2,zzz.t1_LIB(zzz.tbeg_abs(2):sum(zzz.tbeg_abs),i+1)/1e2,'-', ...
[7 0],[],'k');%get(hpl,'color'));
end
axis equal
colorbar
title('Bphi')
figure;
subplot(2,1,1)
plot(zzz.te(:,1),zzz.te(:,2))
ylabel('Te')
subplot(2,1,2)
plot(zzz.ne(:,1),zzz.ne(:,2))
ylabel('ne')
Olivier Sauter
committed
figure;
subplot(2,1,1)
try
plot(zzz.t2_new_LIB(:,1),zzz.t2_new_LIB(:,2))
Olivier Sauter
committed
ylabel('Pdens')
subplot(2,1,2)
plot(zzz.t2_new_LIB(:,1),zzz.t2_new_LIB(:,3))
ylabel('jcd')
catch
end
Olivier Sauter
committed
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
figure
theta=linspace(0,2.*pi,300);
Rmintheta = (0.3+zzz.rmesh_min).*cos(theta);
Ymintheta = (0.3+zzz.rmesh_min).*sin(theta);
Rmaxtheta = (0.3+zzz.rmesh_max).*cos(theta);
Ymaxtheta = (0.3+zzz.rmesh_max).*sin(theta);
plot(Rmintheta,Ymintheta,'k')
hold on
plot(Rmaxtheta,Ymaxtheta,'k')
for i=1:2:size(zzz.t1_LIB,2)
hpl=plot(zzz.t1tor_LIB(:,i)/1e2,zzz.t1tor_LIB(:,i+1)/1e2);
plotos(zzz.t1tor_LIB(zzz.tbeg_abs(2):sum(zzz.tbeg_abs),i)/1e2,zzz.t1tor_LIB(zzz.tbeg_abs(2):sum(zzz.tbeg_abs),i+1)/1e2,'-', ...
[7 0],[],'k');%get(hpl,'color'));
end
figure
contour(zzz.rmesh,zzz.zmesh,zzz.psi',100);
hold on
contour(zzz.rmesh,zzz.zmesh,zzz.psi',[zzz.psi_edge zzz.psi_edge],'linewidth',2,'linecolor','k');
for i=1:2:size(zzz.t1_LIB,2)
hpl=plot(zzz.t1_LIB(:,i)/1e2,zzz.t1_LIB(:,i+1)/1e2);
plotos(zzz.t1_LIB(zzz.tbeg_abs(2):sum(zzz.tbeg_abs),i)/1e2,zzz.t1_LIB(zzz.tbeg_abs(2):sum(zzz.tbeg_abs),i+1)/1e2,'-', ...
[7 0],[],'k');%get(hpl,'color'));
end
axis equal
colorbar
title('psi')
torbeam_std_files_struct = zzz;