-
Notifications
You must be signed in to change notification settings - Fork 0
/
sliceplot_params_vlos.m
109 lines (91 loc) · 2.21 KB
/
sliceplot_params_vlos.m
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
69
70
71
72
73
74
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
function sliceplot_params_vlos(Rx,Ry,Rz,data1,data2,data3,data4,title1)
lp=50;
xmin=min(Rx(:)); xmax=max(Rx(:));
ymin=min(Ry(:)); ymax=max(Ry(:));
zmin=min(Rz(:)); zmax=max(Rz(:));
% x=linspace(-50,250,lp);
% y=linspace(-50,400,lp);
% z=linspace(100,400,lp);
x=linspace(-50,150,lp);
y=linspace(-50,250,lp);
z=linspace(100,400,lp);
[X,Y,Z]=meshgrid(x,y,z);
datFN=TriScatteredInterp(Rx(:),Ry(:),Rz(:),data1(:));
datacube1=datFN(X(:),Y(:),Z(:));
datacube1=reshape(datacube1,[lp,lp,lp]);
datFN=TriScatteredInterp(Rx(:),Ry(:),Rz(:),data2(:));
datacube2=datFN(X(:),Y(:),Z(:));
datacube2=reshape(datacube2,[lp,lp,lp]);
datFN=TriScatteredInterp(Rx(:),Ry(:),Rz(:),data3(:));
datacube3=datFN(X(:),Y(:),Z(:));
datacube3=reshape(datacube3,[lp,lp,lp]);
datFN=TriScatteredInterp(Rx(:),Ry(:),Rz(:),data4(:));
datacube4=datFN(X(:),Y(:),Z(:));
datacube4=reshape(datacube4,[lp,lp,lp]);
sx=[];
sy=[];
dz=zmax-zmin;
dz=dz/10;
sz=linspace(zmin,zmax,4);
%sz=[120,150,200,250]; %Selected altitudes for plotting, updated 2017/03/31
sz=[120,150,250,325];
set(gcf,'PaperPosition',[0 0 12 4.5]);
az=-62.5;
%% el=24;
%el=18;
el=15;
nlims=[10.6 11.6];
Tlims=[0 3e3];
subplot(6,3,[1:3]);
title(title1,'FontSize',18);
axis off;
here=find(datacube1<=0);
datacube1(here)=1;
% subplot(221);
subplot(6,4,[5:4:17]);
slice(X,Y,Z,log10(datacube1),sx,sy,sz);
view([az,el]);
shading flat;
axis tight;
caxis(nlims);
xlabel('E [km]');
ylabel('N [km]');
zlabel('alt. [km]');
h=colorbar;
ylabel(h,sprintf('log_{10} n_e [m^{-3}]'))
% subplot(222);
subplot(6,4,[6:4:18]);
slice(X,Y,Z,datacube2,sx,sy,sz);
view([az,el]);
shading flat;
axis tight;
caxis(Tlims);
xlabel('E [km]');
ylabel('N [km]');
%zlabel('alt. [km]');
h=colorbar;
ylabel(h,'T_i [K]')
% subplot(223);
subplot(6,4,[7:4:19]);
slice(X,Y,Z,datacube3,sx,sy,sz);
view([az,el]);
shading flat;
axis tight;
caxis(Tlims);
xlabel('E [km]');
ylabel('N [km]');
%zlabel('alt. [km]');
h=colorbar;
ylabel(h,'T_e [K]')
subplot(6,4,[8:4:20]);
slice(X,Y,Z,datacube4,sx,sy,sz);
view([az,el]);
shading flat;
axis tight;
caxis([-500 500]);
xlabel('E [km]');
ylabel('N [km]');
%zlabel('alt. [km]');
h=colorbar;
ylabel(h,'v_{los} [m/s]')
end