function [Ref,Code,Fac,H,D10,D50,D90,Ws,Rhos,Grad,Tflow,Tosc,Tasym,Tspec,hw,t,Urms,Uon,Uoff,r,b,Unet,Ang,Zref,...
    Bedstate,Bedform,Rh,Rl,Mut,Muts,Mc,Mct,Mcts,Qs,N,D1,P1,D2,P2,D3,P3]=readinput()
%--------------------------------------------------------------------------
%%% THE SANTOSS PRACTICAL SAND TRANSPORT MODEL        
%--------------------------------------------------------------------------
%%% will be written down in version 2 paper
%%% 2009, January 7
%%% René Buijsrogge - Add 'calculation only' option. Read data
%%%                                         from file. Version santoss1c01
%--------------------------------------------------------------------------
%%% 
% %Getting the data to use in 'calc'  mode

[num,txt,raw]=xlsread('Transport_input.xls');
no = num(1,1);
for i=1:no
     Ref(i,1)    = txt(10,i)   ; % Not  used here;  7-digit Code 
     Code(i,1)  = txt(11,i)    ;% Not  used here;  7-digit Code  
     Fac(i,1)    = txt(12,i)    ;% Not  used here;  facility
     H(i,1)       = num(5,i)        ; % water depth [m]
     D10(i,1)    = num(6,i)     ; % D10  [m]
     D50(i,1)    = num(7,i)    ; % D50 [m]
     D90(i,1)    = num(8,i)    ; % D90  [m]
     Ws(i,1)     =  num(9,i)    ; % mean settling velocity of the bed sand [m/s]
     Rhos(i,1)   = num(10,i)    ; % sand density [kg/m3]
     Grad(i,1)   = txt(19,i)     ; %  graded sand (g)  or uniform (u) sand 
     Tflow(i,1)  = txt(20,i)    ; % type of flow 
     Tosc(i,1)   = txt(21,i)       ; % type of oscillatoy flow 
     Tasym(i,1) = txt(22,i)        ; % type of oscillatoy flow asymetry
     Tspec(i,1)  = num(15,i)              ; % type of oscillatoy flow/wave spectrum
     hw(i,1)      = num(16,i)            ; % type of oscillatoy flow/wave spectrum
     t(i,1)         = num(17,i)              ; % oscillatoy flow period [s]
     Urms(i,1)    = num(18,i)  ; % root-mean-square orbital velocity [m/s]
     Uon(i,1)     = num(19,i)          ; % positive (onshire) orbital velocity peak [m/s]
     Uoff(i,1)     = num(20,i)   ; % negative (offshore) orbital velocity peak [m/s]
     r(i,1)         = num(21,i)       ; % degree of orbital flow aymetry [-]
     b(i,1)         = num(22,i)         ; % degree of flow acceleration asymetry [-]
     Unet(i,1)    = num(23,i)        ; % magnitude net current velocity at reference level zref [m/s]
     Ang(i,1)     =  num(24,i)       ; % angle between net cur. dir. and pos. (onshore) flow dir. [deg]
     Zref(i,1)     = num(25,i)        ; % reference level net current velocity [m]
     Bedstate(i,1)=txt(34,i)      ; % fixed bed or mobile bed
     Bedform(i,1)= txt(35,i)   ; % type of bed form
     Rh(i,1)       = num(28,i)    ; % ripple height [m]
     Rl(i,1)        = num(29,i)    ; % ripple length [m]
     Mut(i,1)     = txt(38,i)      ;  % intra-wave velocity measurements in the suspension layer
     Muts(i,1)    = txt(39,i)      ;  % intra-wave velocity measurements in the sheet-flow layer
     Mc(i,1)       = txt(40,i)       ;  % wave averaged concentration measurements in the suspension layer
     Mct(i,1)      = txt(41,i)        ;  % intra-wave concentration  measurements in the suspension layer
     Mcts(i,1)    = txt(42,i)      ;  % intra-wave concentration measurements in the sheet-flow layer
     Qs(i,1)       = num(35,i)      ; % MEASURED net transport rate [m2/s]
     D1(i,1)       = num(36,i)  ; % median grain size of fraction 1 [mm]; NaN for uniform sands
     P1(i,1)       = num(37,i)  ; % percentage of fraction 1 [%]; NaN for uniform sands. 	
     D2(i,1)       = num(38,i)  ; % median grain size of fraction 2 [mm]; NaN for uniform sands
     P2(i,1)       = num(39,i)  ; % percentage of fraction 2 [%]; NaN for uniform sands. 	
     D3(i,1)       = num(40,i)  ; % median grain size of fraction 3 [mm]; NaN for uniform sands
     P3(i,1)       = num(41,i)  ; % percentage of fraction 3 [%]; NaN for uniform sands. 	
end   
N=no;

% Ref,Code,Fac,H,D10,D50,D90,Ws,Rhos,Grad,Tflow,Tosc,Tasym,Tspec,hw,t,Urms,Uon,Uoff,r,b,Unet,Ang,Zref,...
%     Bedstate,Bedform,Rh,Rl,Mut,Muts,Mc,Mct,Mcts,Qs,N
                
                
                