function [Ps, Is, Ts] = paramsSIP()
% Parameter file for SIP model
% tbh - May 2015

%% Parameters
Ps={%   Value       Name        Description
		0,          'Rr',       'Restoration of reds';
        2,          'Rg',       'Restoration of greens';
        0.32,       'fr',       'fecundity red';
        0.32,       'fg',       'fecundity greens';
        0.15,       'ur',       'mortality reds';
        0.15,       'ug',       'mortality greens';
        0.005,      'c',        'density dependence';
        0.003,      'beta',     'transmission coefficient';
        0.2,        'alpha_r',  'disease on reds';
        0.1,        'alpha_g',  'disease on greens';
        500.0,      'sigma_r',  'shedding from reds';
        50.0,       'sigma_g',  'shedding from greens';
        50.0,       'delta',    'parasite inactivation';
        };

%% Initial Values
Is={%   Value	Name	Description
		10.0,  'Sr_0', 'Initial susceptible red';
        10.0,  'Sg_0', 'Initial susceptible green';
        1.0,	'Ir_0', 'Initial infected red';
        1.0,    'Ig_0','Initial infected green';
        0.1,    'P_0', 'Initial RLO'; 
};

%% Time
Ts={%   Value	Name	Description
		0.0,	'T0',	'Start Time (years)';
		0.0027,	'Inc',	'Time Increment (day)';
		100.0,   'T_F',	'Finish Time (years)'
};