htj_rpt.m

function htj_rpt( p, t, s, d, ttle )
%HTJ_RPT - Produces the final figure and results for the Hantush and Jacob (1955) model
%                
% Syntax: hcj_rpt( p, t, s, d, ttle )
%
%   p(1) = a  = slope of the late time straight line
%   p(2) = t0 = intersept of late time straight line
%   p(3) = r/B 
%   t = measured time
%   s = measured drawdown
%   d(1) = Q = Pumping rate
%   d(2) = r = Distance to the pumping well
%   d(3) = e' = Thickness of the aquitard
%   ttle = Title of the figure 
%
% Description:
%   Produces the final figure and results for Hantush (1955) model.
%
% See also: htj_dmo, htj_dim, htj_gss
%


if(nargin==4) % Default value for d if not given by user
  ttle='Diagnostic plot';
end

% Rename the parameters for a more intuitive check of the formulas
Q=d(1); r=d(2); e=d(3);
a=p(1); t0=p(2); rB=p(3);

% Compute the transmissivity, storativity and radius of influence
T=0.1832339*Q/a;
S=2.2458394*T*t0/r^2;
B=r/p(3);
Ka=T*e/B^2;

% Calls an internal script that computes drawdown, derivative, and residuals
func='htj'; rpt_cmp

% Defines the text of the left side of the legend 
lgdl=char('Test data:');
lgdl=char(lgdl, sprintf(' Discharge rate: %3.2e m^3/s',Q));
lgdl=char(lgdl, sprintf(' Radial distance: %0.2g m',r));
lgdl=char(lgdl, sprintf(' Thickness of aquitard:  %0.2g m',e));
lgdl=char(lgdl, ' ');
lgdl=char(lgdl, 'Hydraulic parameters:');
lgdl=char(lgdl, sprintf(' Transmissivity T: %3.1e m^2/s',T));
lgdl=char(lgdl, sprintf(' Storativity S: %3.1e',S));
lgdl=char(lgdl, sprintf(' Aquitard conductivity k: %2.1e m/s',Ka));

% Defines the text of the right side of the legend 
lgdr=char('Hantush and Jacob (1955) Model');
lgdr=char(lgdr, ' ');
lgdr=char(lgdr, 'Fitting parameters:');
lgdr=char(lgdr, sprintf(' slope a: %0.2g m',a));
lgdr=char(lgdr, sprintf(' intercept t0: %0.2g m',t0));
lgdr=char(lgdr, sprintf(' r/B : %3.1e',rB));
lgdr=char(lgdr, sprintf(' mean residual: %0.2g m',mr));
lgdr=char(lgdr, sprintf(' 2 standard deviation: %0.2g m',sr));


% Calls an internal script that place the legend
rpt_lgd

% Calls an internal script that plot the model and the data
rpt_plt