% Bradway & Shafik 1977 metabolite study
% Male rats dosed for 3 days by gavage


rat_parameters 	       % rat parameters
TSTOP=120;      			  % Simulation duration
BW = 0.425;   	 			 % kg
FA=1; 

MAXT=0.001;       %maximum time step
MINT=1.0e-6;      %minimum time step


%3-day exposure
%to dose with mid or low dose, change name of subsheet accordingly
DT = xlsRead("Bradway_dosing.xls", "high", "A1:A6");  % Dosing Times (hr)
DORAL=xlsRead("Bradway_dosing.xls", "high", "B1:B6"); % Dose in ug/kg

prepare @clear @all

load @format=ascii @separator=comma @ignoreheader @file=Bradway_metabolites_high_dose.csv Bradway_high
tdata1=Bradway_high
load @format=ascii @separator=comma @ignoreheader @file=Bradway_metabolites_mid_dose.csv Bradway_mid
tdata2=Bradway_mid
load @format=ascii @separator=comma @ignoreheader @file=Bradway_metabolites_low_dose.csv Bradway_low
tdata3=Bradway_low

start@nocallback

%high dose
plot(_time, _amcaex, "-", tdata1(:,1), tdata1(:,7), "*", 'MCA_DCA_high_dose.aps')

%mid dose
%plot(_time, _amcaex, "-", tdata2(:,1), tdata2(:,7), "*", 'MCA_DCA_mid_dose.aps')

%low dose
%plot(_time, _amcaex, "-", tdata3(:,1), tdata3(:,7), "*", 'MCA_DCA_low_dose.aps')
