// Scilab code developed by:
// International Research Institute for Climate and Society
// Columbia University in the City of New York, USA
// Last update: June 26, 2013
// --------------------------------------------------------------------------------------------------------
// Dependancies : (see each function for further details)
//      -- function infected_basic.sci for basic model executions
//      -- function Stick_Param that put the moving and fixed parameters in the right order
// --------------------------------------------------------------------------------------------------------
// Functions :
//      -- Put the parameters in the right order
//      -- Check wether the parameters are in the right ranges
//      -- Call infected_basic for the whole period of time (1 -- t)
// --------------------------------------------------------------------------------------------------------


function[I1]=infected_likelihood(param_move,param_fix,t,XXX)
    
    //Reconstruct the array of parameters
    param=Stick_Param(param_move,param_fix);

    qq=size(param)
    for i=1:qq(1)
        if (Pa(i,4)==1) then
            if (Param_min(i)>param(i)) then
                mprintf('parameter %i equal to %f',i,param(i))
                param(i)=Param_min(i); //This is to be sure that the parameters tested are in the defined randes. If not, the min/max of the parameter is taken instead.
                disp(param);
            end
            if (Param_max(i)<param(i)) then
                mprintf('parameter %i equal to %f',i,param(i))
                param(i)=Param_max(i);
                disp(param);
            end
        end
    end
    
    I1=infected_basic(param,1,t,XXX);

endfunction