Commit 608ed48a authored by tilman.metz's avatar tilman.metz

some changes

parent e4ddb5b1
......@@ -30,7 +30,7 @@ void model(st_model *mod){
float depth=3;
float pos=NXG/2*DX;
// Parameters of trench "Ettlinger Linie" velocities in km/s and density in kg/m^3
float evs=0.0917;
float evs=0.1000;
float evp=0.2569;
float erho=1.5;
......
......@@ -25,9 +25,11 @@ Y=dh:dh*outy:ny*outy*dh;
Z=dh:dh*outz:nz*outz*dh;
% Slizes of the 3D-Model in m
slx=5;
slx=10;
sly=2;
slz=5;
slz=20;
if strcmp(parameter,'vs')
caxis_value_1=0;%vs
......@@ -67,7 +69,11 @@ else
y_line=[FW ny-FW ny-FW FW FW]*dh;
end
% Lines to mark slice position
x_line_slice=[slx slx];
y_line_slice=[sly sly];
z_line_slice=[slz slz];
%m --> model indices
slizex=int64(slx/dh/outx);
......@@ -152,6 +158,8 @@ end
figure(1)
imagesc(Z,Y,model_truex);
line(x_line,y_line,'LineStyle','--','Color','k')
line([0 nz*dh],y_line_slice,'LineStyle','-','Color','red')
%line(z_line_slice,[0 ny*dh],'LineStyle','-','Color','green')
colb=colorbar;
coll=get(colb,'xlabel');
set(coll,'String',[parameter ' in m/s'],'FontSize',fontsize);
......@@ -183,6 +191,8 @@ hold off
figure(3)
imagesc(X,Z,model_truey);
line(x_line,z_line,'LineStyle','--','Color','k')
line([0 nz*dh],z_line_slice,'LineStyle','-','Color','green')
%line(x_line_slice,[0 nx*dh],'LineStyle','-','Color','m')
colb=colorbar;
coll=get(colb,'xlabel');
set(coll,'String',[parameter ' in m/s'],'FontSize',fontsize);
......@@ -216,15 +226,16 @@ hold off
figure(5)
imagesc(X,Y,model_truez);
line(x_line,y_line,'LineStyle','--','Color','k')
line([0 nz*dh],y_line_slice,'LineStyle','-','Color','red')
%line(x_line_slice,[0 ny*dh],'LineStyle','-','Color','magenta')
colb=colorbar;
coll=get(colb,'xlabel');
set(coll,'String',[parameter ' in m/s'],'FontSize',fontsize);
colormap('jet');
caxis([caxis_value_1 caxis_value_2])
set(gca,'ydir','reverse');
%set(gca,'xdir','reverse');
xlabel('x in m','FontSize',fontsize)
ylabel('y in m','FontSize',fontsize)
ylabel('y in m (depth)','FontSize',fontsize)
title(['true ' parameter '-model. Slize at z=' num2str(slz) 'm'],'FontSize',fontsize)
axis tight
hold on
......@@ -343,6 +354,6 @@ end
%exportfig(2, [fname parameter '.eps'],'bounds','tight', 'color','rgb', ...
% 'preview','none', 'resolution',200, 'lockaxes',1);
%
% exportfig(3, [fname parameter '.eps'],'bounds','tight', 'color','rgb', ...
% 'preview','none', 'resolution',200, 'lockaxes',1);
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment