# include -
include("Include.jl")

# load my experimental_data file -
path_to_experimental_data_file = joinpath(_PATH_TO_DATA, "Data-3-5-AMP.txt")
experimental_data_array = load_experimental_data_set(path_to_experimental_data_file)

# load my parameters -> configuration_dictionary contains, two sub-dictionaries 
path_to_configuration_file = joinpath(_PATH_TO_CONFIG, "Configuration.toml")
configuration_dictionary = load_configuration_dictionary(path_to_configuration_file)
v_variable_parameter_dictionary = configuration_dictionary["vvariableparameters"]
kinetic_limit_parameter_dictionary = configuration_dictionary["kineticparameters"]

# lets compute the kinetic limit -
kinetic_limit_value = compute_kinetic_limit(kinetic_limit_parameter_dictionary)

# ok, so let's compute a data driven v-variable -
data_driven_v_variable = Array{Float64,1}()
(number_of_effector_values, number_of_cols) = size(experimental_data_array)
for row_index = 1:number_of_effector_values
    value = experimental_data_array[row_index,2]
    scaled_value = value / kinetic_limit_value
    push!(data_driven_v_variable, scaled_value)
end

# make a plot -
plot(experimental_data_array[:,1],data_driven_v_variable[:,1],legend=false)
xlabel!("Effector 3-5-AMP (mM)",fontstze=16)
ylabel!("v-variable (AU)", fontsize=16)