function [t] = Cavity_many_cells(fileID,num,para,plottype,i_space_x,i_space_y,i_font,i_contour,i_scale,... i_domain_velocity,i_layer_velocity,i_cilia_velocity,i_cilia_scale,i_plot_X2Y2,... i_particle,i_set_gca_xy,i_set_gca_x_axis,i_set_gca_y_axis); clf file1 = 'moving.'; file2 = 'interface.'; fname1 = sprintf('%s%s%04d',file1,fileID,num,'.dat'); fname2 = sprintf('%s%s%04d',file2,fileID,num,'.dat'); display(fname1) %pause fid = fopen(fname1, 'r'); %pause %display(t) t = fscanf(fid,'%f',1); %pause %display(t) mx = fscanf(fid,'%d',1); my = fscanf(fid,'%d',1); xmin = fscanf(fid,'%f',1); ymin = fscanf(fid,'%f',1); xmax = fscanf(fid,'%f',1); ymax = fscanf(fid,'%f',1); A = fscanf(fid,'%f',[3 mx*my]); fclose(fid); hold on; A = A'; x = linspace(xmin,xmax,mx); y = linspace(ymin,ymax,my); [X,Y] = meshgrid(x,y); X = X'; Y = Y'; h = (xmax-xmin)/(mx-1); if (para == 'p') P = reshape(A(:,3),mx,my); if (plottype==1) mesh(X,Y,P); zlabel('Z') view(145,12) % title(['Biofilm development'],'fontsize',i_font); title(['Nutrient & biofilm at day ',num2str(t)],'fontsize',i_font); else axis equal % surf(X,Y,P); contourf(X,Y,P,i_contour,':k*'); % title(['Biofilm development'],'fontsize',i_font); title(['Nutrient & biofilm at day ',num2str(t)],'fontsize',i_font); colormap(flipud(gray)); colorbar % h.limits=[0,0.0005] % h = colorbar; % caxis([0.0 0.0005]) end end; hold on; if (para == 'p') fid = fopen(fname2, 'r'); N_temp = fscanf(fid,'%d',[6 1]); N = N_temp(1); i_cilia = N_temp(2); species1no=N_temp(3); species2no=N_temp(4); species3no=N_temp(5); species4no=N_temp(6); XY = fscanf(fid,'%f',[2 N*i_cilia])'; fclose(fid); % plot(XY(1:N,1),XY(1:N,2),'g-','LineWidth',2); for i=1:i_cilia N_begin = 1+(i-1)*N; N_end = i*N; if (i<=species1no) plot(XY(N_begin:N_end,1),XY(N_begin:N_end,2),'r-','LineWidth',2); %r elseif ((i>species1no) & (i<=species1no+species2no)) plot(XY(N_begin:N_end,1),XY(N_begin:N_end,2),'b-','LineWidth',2); %b elseif ((i>species1no+species2no) & (i<=species1no+species2no+species3no)) plot(XY(N_begin:N_end,1),XY(N_begin:N_end,2),'y-','LineWidth',2); else plot(XY(N_begin:N_end,1),XY(N_begin:N_end,2),'g-','LineWidth',2); end % plot(XY(N_begin:N_end,1)-(xmax-xmin),XY(N_begin:N_end,2),'r-','LineWidth',2); % plot(XY(N_begin:N_end,1)+(xmax-xmin),XY(N_begin:N_end,2),'r-','LineWidth',2); end; end; %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% axis equal axis([xmin xmax ymin ymax]); if (i_set_gca_xy==1) set(gca,'xtick',[xmin:i_set_gca_x_axis:xmax]); set(gca,'ytick',[ymin:i_set_gca_y_axis:ymax]); set(gca,'fontsize',i_font); end box on hold off