Scribd
Upload
92 views4 pages

Script:: Adika Bagaskara 03411740000020 Geodinamika B

This document appears to be a technical report or paper describing the computation and visualization of velocity fields, divergence of velocity, and time derivatives of density using MATLAB code. It defines model parameters like grid size and resolution, then uses loops to calculate velocity components, partial derivatives, and divergence. Finally it generates 6 subplots to display colormap and vector plots of the computed velocity fields, derivatives, and divergence.

Uploaded by

Adika Bagaskara
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
92 views4 pages

Script:: Adika Bagaskara 03411740000020 Geodinamika B

This document appears to be a technical report or paper describing the computation and visualization of velocity fields, divergence of velocity, and time derivatives of density using MATLAB code. It defines model parameters like grid size and resolution, then uses loops to calculate velocity components, partial derivatives, and divergence. Finally it generates 6 subplots to display colormap and vector plots of the computed velocity fields, derivatives, and divergence.

Uploaded by

Adika Bagaskara
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 4

ADIKA BAGASKARA

03411740000020
GEODINAMIKA B

Script:
% Computation and visualisation of velocity, divergence of velocity and
% time derivatives of density with ‘pcolor’ and ‘quiver’

% Clearing all variables and arrays


clear;
% Clearing figures
clf;

% Defining new figure


figure(1)

% Defining model size


xsize=1000000; % Horizontal size, m
ADIKA BAGASKARA
03411740000020
GEODINAMIKA B
ysize=1500000; % Vertical size, m
% Defining resolution
xnum=36; % Horizontal resolution (nodal points)
ynum=36; % Vertical resolution (nodal points)
% Step between nodal ponts
xstp=xsize/(xnum-1); % Horizontal grid step
ystp=ysize/(ynum-1); % Vertical grid step
% Defining scales for Vx and Vy
vx0=1e-9*xsize/2/ysize; % Scale for horizontal velocity, m/s
vy0=1e-9; % Scale for vertical velocity, m/s

% 2D colormap
% Creating a vector of arguments for two axes
x=0:xstp:xsize;
y=0:ystp:ysize;

% Creating velocity function


% Smallest index in an array is 1 (not 0)
for i=1:1:ynum
for j=1:1:xnum
% Setup velocity field corresponding
% to circulation wit central upwelling
% Defining Vx = horizontal velocity component distribution
vx(i,j)=-vx0*sin(pi*x(j)/xsize*2)*cos(pi*y(i)/ysize);
% Defining Vy = vertical velocity component distribution
vy(i,j)=vy0*sin(pi*y(i)/ysize)*cos(pi*x(j)/xsize*2);
% Computing partial derivatives
% dVx/dx
dvxdx(i,j)=-vx0*pi/xsize*2*cos(pi*x(j)/xsize*2)*cos(pi*y(i)/ysize);
% dVy/dy
dvydy(i,j)=vy0*pi/ysize*cos(pi*y(i)/ysize)*cos(pi*x(j)/xsize*2);
% Computing div(v)=dVx/dx+dVy/dy
divv(i,j)=dvxdx(i,j)+dvydy(i,j);
end
end

figure(1);
% Plotting Vx velocity
subplot(2,3,1); % defining 1st plotting area in the 3x2 figure
pcolor(x/1000,y/1000,vx); % making a colormap
box on; % making a box around the plot
title('Vx colormap, m/s'); % title for the plot
xlabel('x, km'); % title for the horizontal axis
ylabel('y, km'); % title for the vertical axis
shading interp; % making smooth transitions between colors
colorbar; % showing a colorbar for the map
axis ij image ; % directing vertical axis downward, making proper
dimensions

% Plotting Vy velocity
subplot(2,3,2); % defining 2nd plotting area in the 3x2 figure
pcolor(x/1000,y/1000,vy); % making a colormap
box on; % making a box around the plot
title('Vy colormap, m/s'); % title for the plot
xlabel('x, km'); % title for the horizontal axis
ylabel('y, km'); % title for the vertical axis
ADIKA BAGASKARA
03411740000020
GEODINAMIKA B
shading interp; % making smooth transitions between colors
colorbar; % showing a colorbar for the map
axis ij image ; % directing vertical axis downward, making proper
dimensions

% Plotting velocity vector as arrows


subplot(2,3,3); % defining 3d plotting area in the 3x2 figure
quiver(x/1000,y/1000,vx,vy,'k'); % making a field of arrows
box on; % making a box around the plot
title('velocity field'); % title for the plot
xlabel('x, km'); % title for the horizontal axis
ylabel('y, km'); % title for the vertical axis
axis([0 xsize/1000 0 ysize/1000]); % Defining upper and lower limits for the
axes
axis ij image ; % directing vertical axis downward, making proper
dimensions

% Plotting dVx/dx
subplot(2,3,4); % defining 4th plotting area in the 3x2 figure
pcolor(x/1000,y/1000,dvxdx); % making a colormap
box on; % making a box around the plot
title('dVx/dx colormap, 1/s'); % title for the plot
xlabel('x, km'); % title for the horizontal axis
ylabel('y, km'); % title for the vertical axis
shading interp; % making smooth transitions between colors
colorbar; % showing a colorbar for the map
axis ij image ; % directing vertical axis downward, making proper
dimensions

% Plotting dVy/dy
subplot(2,3,5); % defining 5th plotting area in the 3x2 figure
pcolor(x/1000,y/1000,dvydy); % making a colormap
box on; % making box around the plot
title('dVy/dy colormap, 1/s'); % title for the plot
xlabel('x, km'); % title for the horizontal axis
ylabel('y, km'); % title for the vertical axis
shading interp; % making smooth transitions between colors
colorbar; % showing a colorbar for the map
axis ij image ; % directing vertical axis downward, making proper
dimensions

% Plotting div(v)
subplot(2,3,6); % defining 6th plotting area in the 3x2 figure
pcolor(x/1000,y/1000,divv); % making a colormap
box on; % making a box around the plot
title('div(v) colormap, 1/s'); % title for the plot
xlabel('x, km'); % title for the horizontal axis
ylabel('y, km'); % title for the vertical axis
shading interp; % making smooth transitions between colors
colorbar; % showing acolorbar for the map
axis ij image ; % directing vertical axis downward, making proper
dimensions
ADIKA BAGASKARA
03411740000020
GEODINAMIKA B
Output Gambar:

You might also like