45 lines
944 B
Matlab
45 lines
944 B
Matlab
function [ ] = plot_reservoir_3D(dx,dy,dz,nx,ny,nz )
|
|
%UNTITLED 此处显示有关此函数的摘要
|
|
% 此处显示详细说明
|
|
%figure
|
|
% 绘制x-y平面
|
|
for i=1:(nz+1)
|
|
% if i ==nz+1
|
|
face=[0,0,(i-1)*dz;nx*dx,0,(i-1)*dz;nx*dx,ny*dy,(i-1)*dz;0,ny*dy,(i-1)*dz];
|
|
% fill3(face(:,1),face(:,2),face(:,3),[0.5,0.5,0.5]);
|
|
fill3(face(:,1),face(:,2),face(:,3),'w');
|
|
% end
|
|
|
|
hold on
|
|
end
|
|
% [0.5,0.5,0.5]
|
|
|
|
% 绘制x-z平面
|
|
for i=1:(ny+1)
|
|
% if i == 1
|
|
face=[0,(i-1)*dy,0;nx*dx,(i-1)*dy,0;nx*dx,(i-1)*dy,nz*dz;0,(i-1)*dy,nz*dz];
|
|
% fill3(face(:,1),face(:,2),face(:,3),[0.5,0.5,0.5]);
|
|
fill3(face(:,1),face(:,2),face(:,3),'w');
|
|
% end
|
|
hold on
|
|
end
|
|
|
|
% 绘制y-z平面
|
|
for i=1:(nx+1)
|
|
% if i == 1
|
|
face=[(i-1)*dx,0,0;(i-1)*dx,ny*dy,0;(i-1)*dx,ny*dy,nz*dz;(i-1)*dx,0,nz*dz;];
|
|
% fill3(face(:,1),face(:,2),face(:,3),[0.5,0.5,0.5]);
|
|
fill3(face(:,1),face(:,2),face(:,3),'w');
|
|
% end
|
|
hold on
|
|
end
|
|
|
|
xlabel('x,m');
|
|
ylabel('y,m');
|
|
zlabel('z,m');
|
|
alpha(1);
|
|
hold on;
|
|
view(3)
|
|
end
|
|
|