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