read_params_simul_domain

 read simulation domain parameters
 
 Input:
         fname_domain
         
 Output:
     1. params_domain_geom is a structure with 3 elements:
         Rratio_IN
         include_ECS
         ECS_gap
         
     2. params_domain_pde is a structure with 8 elements:
         dcoeff_IN
         dcoeff_OUT
         dcoeff_ECS
         ic_IN
         ic_OUT
         ic_ECS
         kappa_IN_OUT
         kappa_OUT_ECS
         
     3. params_domain_femesh is a structure with 2 elements:
         Htetgen
         tetgen_cmd

0001 function [params_domain_geom,params_domain_pde,params_domain_femesh] ...
0002     = read_params_simul_domain(fname_domain)
0003 
0004 
0005 % read simulation domain parameters
0006 %
0007 % Input:
0008 %         fname_domain
0009 %
0010 % Output:
0011 %     1. params_domain_geom is a structure with 3 elements:
0012 %         Rratio_IN
0013 %         include_ECS
0014 %         ECS_gap
0015 %
0016 %     2. params_domain_pde is a structure with 8 elements:
0017 %         dcoeff_IN
0018 %         dcoeff_OUT
0019 %         dcoeff_ECS
0020 %         ic_IN
0021 %         ic_OUT
0022 %         ic_ECS
0023 %         kappa_IN_OUT
0024 %         kappa_OUT_ECS
0025 %
0026 %     3. params_domain_femesh is a structure with 2 elements:
0027 %         Htetgen
0028 %         tetgen_cmd
0029 
0030 ndim = 3;
0031 fid=fopen(fname_domain);
0032 
0033 tline = fgetl(fid);
0034 params_domain_geom.Rratio_IN = sscanf(tline,'%f',1);
0035   
0036 if (params_domain_geom.Rratio_IN < 0 | params_domain_geom.Rratio_IN > 0.99)
0037     params_domain_geom.Rratio_IN = 0;
0038 end
0039 
0040 tline = fgetl(fid);
0041 params_domain_geom.include_ECS = sscanf(tline,'%f',1);
0042 
0043 tline = fgetl(fid);
0044 params_domain_geom.ECS_gap = sscanf(tline,'%f',1); 
0045 
0046 tline = fgetl(fid);
0047 params_domain_pde.dcoeff_IN = sscanf(tline,'%f',1);
0048 
0049 tline = fgetl(fid);
0050 params_domain_pde.dcoeff_OUT = sscanf(tline,'%f',1);
0051 
0052 tline = fgetl(fid);
0053 params_domain_pde.dcoeff_ECS = sscanf(tline,'%f',1);
0054 
0055 % initial conditions
0056 tline = fgetl(fid);
0057 params_domain_pde.ic_IN = sscanf(tline,'%f',1);
0058 
0059 tline = fgetl(fid);
0060 params_domain_pde.ic_OUT = sscanf(tline,'%f',1);
0061 
0062 tline = fgetl(fid);
0063 params_domain_pde.ic_ECS = sscanf(tline,'%f',1);
0064 
0065 tline = fgetl(fid);
0066 params_domain_pde.kappa_IN_OUT = sscanf(tline,'%f',1);
0067 tline = fgetl(fid);
0068 params_domain_pde.kappa_OUT_ECS = sscanf(tline,'%f',1);
0069 
0070 tline = fgetl(fid);
0071 params_domain_femesh.Htetgen = sscanf(tline,'%f',1);
0072 
0073 
0074 tline = fgetl(fid);
0075 [strpos] = regexp(tline,"'");
0076 params_domain_femesh.tetgen_cmd = tline(strpos(1)+1:strpos(2)-1);
0077 fclose(fid);