% read_opacity 
% there is script for reading binary file with opacity coefficient.
global T_net Q_net X_net OpaCoef
filename='howdek-opal95-alexander94.ote.bin';
fid = fopen(filename,'r'); % open file, check file name
A = fread(fid, 1, 'int32');
X_Num=10; T_Num=247; Q_Num=207; % dimensions of net
T_net = fread(fid, T_Num, 'float64');
Q_net = fread(fid, Q_Num, 'float64');
X_net = fread(fid, X_Num, 'float64');
Tcol=X_Num*(T_Num-6+2)*(Q_Num-6+2);     % total size of coef array
OpaCoef = fread(fid,Tcol,'float32');
OpaCoef=reshape(OpaCoef,X_Num,T_Num-6+2,Q_Num-6+2);
fclose(fid);
% saving of coefficient is needed for debug purpose only
save opac_coef;