   p_bottom_pipeline_bara_oil = zeros(1,L);% Matrix Init
% Gisle Otto Eikrem
% 21.10.2002

%Define the size of the matrix initially to reduce simulation time
                                                  
   w_pc = zeros(1,L);
   w_oil_pc = zeros(1,L);
   w_oil_in = zeros(1,L);
   
   w_oil_bottom = zeros(1,L);
   w_gas_pc = zeros(1,L);
   
   w_gas_in = zeros(1,L);
   w_gas_bottom = zeros(1,L);
         
   time_log = zeros(1,L);
   rho_m_pipeline = zeros(1,L);
   
   u_PC_simple = zeros(1,L);
    
   rho_m_riser = zeros(1,L);
   rho_gas_top_pipeline = zeros(1,L);
  
   q_pc = zeros(1,L);
   q_bottom_oil = zeros(1,L);
   
   q_bottom_gas = zeros(1,L);
   p_top_riser_bara = zeros(1,L);
   
   p_top_pipeline_bara = zeros(1,L);
   
   p_gas_flow_noflow = zeros(1,L);
   
   p_bottom_riser_onlyoil_bara = zeros(1,L);
   p_bottom_riser_bara = zeros(1,L);
   
   p_bottom_riser_bara = zeros(1,L);
   p_bottom_pipeline_bara_oil = zeros(1,L);
   p_bottom_pipeline_bara_gas = zeros(1,L);   
   m_gas_riser_min = zeros(1,L);
   m = zeros(4,L);
   G_pipeline_bottom_gas = zeros(1,L);
          
   G_pipeline_bottom_oil = zeros(1,L);
   G_riser_top = zeros(1,L);
   delta_p_pipeline = zeros(1,L);
   delta_p_riser = zeros(1,L);

   delta_p_gas = zeros(1,L);
   
   X_gas=zeros(1,L);
   X_oil=zeros(1,L);
   GOR=zeros(1,L);
   
   w_tot_in=zeros(1,L);
   w_gas_in=zeros(1,L);
   w_oil_in=zeros(1,L);
   Q=zeros(1,L);
   P=zeros(11,L);