create_geom

 create cells (canonical configuration)
 
 Input:
     fname_params_cells
     
 Output: 
     1. params_cells is a structure with
         a. 7 elements for spheres (cell_shape = 1):
             cell_shape
             ncell
             Rmin
             Rmax
             dmin
             dmax
             para_deform 

         b. 8 elements for cylinders (cell_shape = 2):
             cell_shape
             ncell
             Rmin
             Rmax
             dmin
             dmax
             para_deform
             Hcyl   
     2. fname_cells

0001 function [params_cells,fname_cells] = create_geom(fname_params_cells)
0002 
0003 % create cells (canonical configuration)
0004 %
0005 % Input:
0006 %     fname_params_cells
0007 %
0008 % Output:
0009 %     1. params_cells is a structure with
0010 %         a. 7 elements for spheres (cell_shape = 1):
0011 %             cell_shape
0012 %             ncell
0013 %             Rmin
0014 %             Rmax
0015 %             dmin
0016 %             dmax
0017 %             para_deform
0018 %
0019 %         b. 8 elements for cylinders (cell_shape = 2):
0020 %             cell_shape
0021 %             ncell
0022 %             Rmin
0023 %             Rmax
0024 %             dmin
0025 %             dmax
0026 %             para_deform
0027 %             Hcyl
0028 %     2. fname_cells
0029 
0030 [params_cells,fname_cells] = read_params_cells(fname_params_cells);
0031 
0032 if (params_cells.cell_shape == 1)
0033     create_ellipses_inputfile(params_cells.ncell,params_cells.Rmin,params_cells.Rmax,params_cells.dmin,params_cells.dmax,fname_cells);
0034 elseif (params_cells.cell_shape == 2)
0035     create_cylinders_inputfile(params_cells.ncell,params_cells.Rmin,params_cells.Rmax,params_cells.dmin,params_cells.dmax,fname_cells,params_cells.Hcyl);
0036 end