%%
% REFDIST = 0.0017;           % half size of head in sec for dB ref.
% MAXORDER = 16;              % max order of reflections computed     
% N_DEEP = 4;                 % no. of image rooms deep               
% N_ROOM = 2*N_DEEP+1;        % no. of image rooms each dimension     
% MAXR = 500;                 % max no. of reflections                
% MAXDEL = 0.12;              % max reflection delay                  
% DBCRIT = -25;               % criterion dB down for rejection  
% ABSORB = 1;                 % air absorption -1 dB/sec         
% BOUNCE = 3;                 % figure 3 dB loss per bounce      
% DPR = 57.29578;             % degrees per radian conversion    
% TINY = 1.0E-30;             % no zero dividing!                
% SPM = .0034;                % sec per meter conversion         
% MPS = 294.1176;             % meter per sec conversion         
% SR = 44100;                 % default sampling rate            
% CLENGTH = 4096;             % default time window	    

% % 
prompt = 'What is the order value for image model? ';
n_deep = input(prompt);
n_room = n_deep*2 + 1;

prompt = 'What is the MAX reflection value for image model? ';
maxdel = input(prompt);

basic_values = struct('REFDIST' , 0.0017,'MAXORDER',16,'N_DEEP',n_deep,'N_ROOM',n_room,'MAXR',500,'MAXDEL',maxdel,'DBCRIT',-25,'ABSORB',1,'BOUNCE',3,'DPR',57.29578,'TINY',0.1^30,'SPM',0.0034,'MPS',294.1176,'SR',44100,'CLENGTH',4096);

save basic_values.mat basic_values