function d = disturbance(t,SimTime) 
 d=[];                      %Disturbance vector
 
 %% Disturbance in Reservoir Pressure
 A = 5;                    % Amplitude
 T=10*60;                   % Period [s]
 f=1/T;                     % Frequency
 d.P_res=A*sin(2*pi*f*t);   % Disturbance 
 
  %% Disturbance in Gas Source Pressure 
 A = 10;                  % Amplitude
 T=40*60;                   % Period [s]
 f=1/T;                     % Frequency
 d.P_gs=A*sin(2*pi*f*t);     % Disturbance 

 %% Disturbance in Mass-Gas-Oil-Rato (GOR)
 A = 0.01;                  % Amplitude
 T=40*60;                   % Period [s]
 f=1/T;                     % Frequency
 d.GOR=A*sin(2*pi*f*t);     % Disturbance 
 
 %% Plot disturbances
 %figure()
 %subplot(2,1,1)
 %plot(d.P_res)
 %subplot(2,1,1)
 %plot(d.GOR)
end

