nsl_drw.m

function nsl_drw()
%NSL_DRW - Draw the slug test (pulse) type curves of Neuzil (1982)
%
% Syntax: nsl_drw()
%
% Description:
%   Draw a series of type curves of Neuzil (1982)
%
% See also: nsl_dim, nsl_dls
%


figure(1)
clf
t=logspace(-3,3);                                       %Create x-axis (time) from 1E-03 to 1E+03
e=[0,log10(5),log10(50),5,12,40];
e=[0,1,2,3,4,5,6];                                      
for i = 1:6
  s=nsl_dls(10^e(i),t);                                 %nsl_dls(cd, t); cd from 1 to 1E+06 i.e alpha 1/20 to 0.5E-06
  %sap1=1-2/sqrt(pi)/sqrt(10^e(i)).*sqrt(t) % Sageev (1986) assymptotes
  %sap2=1/2./t;
  [td,ds]=ldiff(t,s);
  semilogx(t,s,'-')
  hold on
  semilogx(td,-ds,'--')
end   
text(2e-1,0.4,'C_D=1')
text(1.3,0.45,'5')
text(2.1,0.5,'50')
text(4.5,0.55,'10^5')
text(1e1,0.6,'10^{12}')
text(3e1,0.65,'10^{40}')

legend('drawdown','derivative')
xlabel('t_D/C_D = 2Tt/r_C^2','Fontsize',14)
ylabel('s_D','Fontsize',14)
hold off


figure(2)
clf
t=logspace(-3,3);
for i = 1:6
  s=csl_dls([10^e(i)],t);
  [td,ds]=ldiff(t,s);
  loglog(t,s,'-')
  hold on
  loglog(td,-ds,'--')
end   
text(3e-1,0.3,'C_D=1')
text(1.3,0.45,'5')
text(2.1,0.5,'50')
text(4.5,0.55,'10^5')
text(1e1,0.6,'10^{12}')
text(3e1,0.65,'10^{40}')

legend('drawdown','derivative')
axis([1e-2 1e3 1e-3 1])
xlabel('t_D/C_D = 2Tt/r_C^2','Fontsize',14)
ylabel('s_D','Fontsize',14)

hold off