read_params_cells

 read cells parameters
 
 Input:
     fname_geometry
             
 Output:
     1. geom is a structure with 8 elements:
         cell_shape
         ncell
         Rmin
         Rmax
         dmin
         dmax
         para_deform
         Hcyl
         
     2. fname_geom

0001 function [geom,fname_geom] = read_params_cells(fname_geometry)
0002 
0003 % read cells parameters
0004 %
0005 % Input:
0006 %     fname_geometry
0007 %
0008 % Output:
0009 %     1. geom is a structure with 8 elements:
0010 %         cell_shape
0011 %         ncell
0012 %         Rmin
0013 %         Rmax
0014 %         dmin
0015 %         dmax
0016 %         para_deform
0017 %         Hcyl
0018 %
0019 %     2. fname_geom
0020 
0021 fid=fopen(fname_geometry);
0022 
0023 tline = fgetl(fid);
0024 geom.cell_shape = sscanf(tline,'%f',1);
0025 
0026 tline = fgetl(fid);
0027 [strpos] = regexp(tline,"'");
0028 fname_geom = tline(strpos(1)+1:strpos(2)-1);
0029 
0030 tline = fgetl(fid);
0031 geom.ncell = sscanf(tline,'%f',1);
0032 
0033 tline = fgetl(fid);
0034 geom.Rmin = sscanf(tline,'%f',1);
0035 tline = fgetl(fid);
0036 geom.Rmax = sscanf(tline,'%f',1);
0037 
0038 tline = fgetl(fid);
0039 geom.dmin = sscanf(tline,'%f',1);
0040 tline = fgetl(fid);
0041 geom.dmax = sscanf(tline,'%f',1);
0042 
0043 tline = fgetl(fid);
0044 geom.para_deform = sscanf(tline,'%f',2);
0045 
0046 if (geom.cell_shape == 2)
0047     tline = fgetl(fid);
0048     geom.Hcyl= sscanf(tline,'%f',1);
0049 end
0050 
0051 fclose(fid);
0052     
0053