function sys = v1_gaslift_model(t,x,u,mode,par)

%% Run model
[x_dot,y] =v1_gaslift_model0(x,u,par);

%% Choose output of model
if isequal(mode,'derivatives')    
    sys=x_dot;    
elseif isequal(mode,'measurements')
    sys=[y(1);y(2:5)/1e5;y(6);y(7:8)];
end


end

