PLOT_GEOMETRY_INFO

 plot information of the geometry
 
 Input:
     1. cmpts_bdys_mat
     2. OUT_cmpts_index
     3. IN_cmpts_index
     4. ECS_cmpts_index
     5. VOL_cmpts
     6. SA_cmpts
     7. SAu_cmpts
     
 Output:
     1. 1 figure with title of "Connections boundary-compartment"
     2. 1 figure with title of "VOL_cmpts"
     3. 1 figure with title of "Surface Area"
     4. 1 figure with title of "SAu_cmpts/VOL_cmpts"

0001 function PLOT_GEOMETRY_INFO(cmpts_bdys_mat,OUT_cmpts_index,IN_cmpts_index,ECS_cmpts_index,...
0002     VOL_cmpts,SA_cmpts,SAu_cmpts)
0003 
0004 % plot information of the geometry
0005 %
0006 % Input:
0007 %     1. cmpts_bdys_mat
0008 %     2. OUT_cmpts_index
0009 %     3. IN_cmpts_index
0010 %     4. ECS_cmpts_index
0011 %     5. VOL_cmpts
0012 %     6. SA_cmpts
0013 %     7. SAu_cmpts
0014 %
0015 % Output:
0016 %     1. 1 figure with title of "Connections boundary-compartment"
0017 %     2. 1 figure with title of "VOL_cmpts"
0018 %     3. 1 figure with title of "Surface Area"
0019 %     4. 1 figure with title of "SAu_cmpts/VOL_cmpts"
0020 
0021 Ncmpt = length(VOL_cmpts);
0022 Nboundary = size(cmpts_bdys_mat,2);
0023 
0024 figure; 
0025 hold all;
0026 cmpts_bdys_mat_plot=zeros(size(cmpts_bdys_mat));
0027 cmpts_bdys_mat_plot(OUT_cmpts_index,:)=cmpts_bdys_mat(OUT_cmpts_index,:);
0028 spy(cmpts_bdys_mat_plot,'ko',10);
0029 
0030 
0031 cmpts_bdys_mat_plot(:,:)=0;
0032 cmpts_bdys_mat_plot(ECS_cmpts_index,:)=cmpts_bdys_mat(ECS_cmpts_index,:);
0033 spy(cmpts_bdys_mat_plot,'ks',10);
0034 
0035 cmpts_bdys_mat_plot(:,:)=0;
0036 cmpts_bdys_mat_plot(IN_cmpts_index,:)=cmpts_bdys_mat(IN_cmpts_index,:);
0037 spy(cmpts_bdys_mat_plot,'kd',10); 
0038 
0039 
0040 xlabel('iboundary'); ylabel('icmpt');
0041 title('Connections boundary-compartment');
0042 set(gca,'Ytick',[1:Ncmpt]);
0043 set(gca,'Xtick',[1:Nboundary]);
0044 grid on;
0045 
0046 figure; 
0047 hold on;
0048 bar(1:Ncmpt,VOL_cmpts,'b');
0049 bar(Ncmpt+1,sum(VOL_cmpts),'r');
0050 title('VOL\_cmpts');
0051 set(gca,'Xtick',[1:Ncmpt+1]);
0052 xlabel('icmpt (last: all cmpts) ');
0053 grid on;
0054 
0055 figure; 
0056 hold on;
0057 bar(1:Ncmpt,SA_cmpts,'b');
0058 bar(Ncmpt+1,sum(SA_cmpts),'r');
0059 title('Surface Area');
0060 set(gca,'Xtick',[1:Ncmpt+1]);
0061 xlabel('icmpt (last: all cmpts)');
0062 grid on;
0063 
0064 figure; 
0065 hold on;
0066 bar(1:Ncmpt,SAu_cmpts./VOL_cmpts','b');
0067 bar(Ncmpt+1,sum(SAu_cmpts)./sum(VOL_cmpts),'r');
0068 title('SAu\_cmpts/VOL\_cmpts');
0069 set(gca,'Xtick',[1:Ncmpt+1]);
0070 xlabel('icmpt (last: all cmpts)');
0071 grid on;