% Food exposure 
% Once daily oral exposure, dose from cpfs SS F13_49 foodonly exp distribution.csv 90th percentile  
% Female adult, Body weight set to 72.9 kg
% Cecilia Tan
% August 26, 2015

prepare @clear @all
% Setting model parameters
start @NoCallback
use human_parameters_MRP
use preg_female_parameters

BWSW=1;		% Set to 1 if BWT is provided
BWST=72.9;	% Initial BW in kg
PREG=0;		% Not pregnant
TSTOP= 1008; 	% Exposure duration in hrs

%Single oral daily exposure
DT = xlsRead("Food.xls", "90p", "A1:A43"); % Dosing Times
DORAL=xlsRead("Food.xls", "90p", "B1:B43"); % Dose in ug/kg 
prepare time rbcce cv urinetcpy
start@nocallback

mytable = [_time _rbcce _urinetcpy _cv*350.6];	% unit of conc = ug/L
save mytable @file='food_90.txt' @format=ascii

min(_rbcce)		% display peak RBC AChE inhibition  