function data=gen_data
%   Function defining case spesific data for the two_fluid.m
%
%   Written by Espen Storkaas, Jan 20th 2003

%   Fluid properties
data.Mw_G=20.6;
data.Z=1.2;
data.rho_L=850;
data.my_L=0.00014;
data.my_G=1.3e-5;
data.T=273+42;

%   Physical Constants 
data.R=8.314;
data.g=9.81;

%   Choke parameters
data.Ac=pi/4*0.075^2;
data.Cd=0.9;
data.C_vG=1.75e3;
data.C_pG=2.15e3;
data.k=data.C_pG/data.C_vG;
data.C_L=1.75e3;

%   Friction related constants
data.Ann_lim=0.55;
data.Bub_lim=0.65;
data.Db=0.0005;
data.eps=2.8e-5;
data.Gparam=0.01;
data.Lparam=1;
data.iparam=0.5;

%   Discretizaion and geomerty

%25/24
% data.dx_m=[500;500;250;250;250;250;250;250;250;250;250;250;...
%         200;150;100;100;100;50;50;50;50;50;100;100;100];
% data.dx_i=[500;375;250;250;250;250;250;250;250;250;250;225;...
%         175;125;100;100;75;50;50;50;50;75;100;100];
% data.theta=[0;0;0;0;0;0;-0.0137;-0.0274;-0.0274;-0.0274;-0.0274;...
%         -0.0274;-0.0274;-0.0274;-0.0274;-0.0274;-0.0274;-0.0274;-0.0274;...
%         0.7717;1.571;1.571;1.571;0.7854];
%13/12
%Alt. 1
%  data.dx_m=[1000;1000;1000;500;250;250;100;100;50;50;150;150;100];
%  data.dx_i=[1000;1000;750;375;250;175;100;75;50;100;150;125];
%  data.theta=[0;-0.00875;-0.0175;-0.0175;-0.0175;-0.0175;-0.0175;-0.0175;-0.0175;0.7766;pi/2;0.7854];
%Alt. 2

data.dx_m=[1500;1000;1500;100;100;50;40;20;40;100;125;50;75];
data.dx_i=[1250;1250;800;100;75;45;30;30;70;112.5;87.5;62.5];
data.theta=[0;-0.0175;-0.0175;-0.0175;-0.0175;-0.0175;-0.0175;pi/2;pi/2;...
         pi/2;pi/2;0];
data.D=0.12;
data.A=pi/4*data.D^2;

%   Index vectors

data.I_mL=(1:length(data.dx_m))';
data.I_mG=data.I_mL+max(data.I_mL);
data.I_i=(1:length(data.dx_i))';
data.I_iL=data.I_i+max(data.I_mG);
data.I_iG=data.I_i+max(data.I_iL);


