# This is a code for 8.08 - 8.S308 recitation
# ballistic deposition model

# for installing required packages
# using Pkg
# Pkg.add("Plots")
# Pkg.add("Random")
# Pkg.add("LsqFit")

# for plotting
using Plots
using Random
using LsqFit


cd(@__DIR__)
include("Rec9_ballistic_deposition_module.jl")

#random number generator
rng = MersenneTwister(1234)

# parameters
L = 600     # number of sites
r = 1       # rate of deposition

# set parameters
param = Rec.setParam(L,r)

# set state
state = Rec.setState(param)

# create movie of positions
t_gap = 2     # gap between frames
n_frame = 100 # number of frames
fps = 10      # fps of gif
get_movie!(state, param, rng, t_gap, n_frame, fps)

# create plot of roughness
n_pts = 20    # number of time points in the roughness plot
t_gap = 2     # time gap between time points
n_seeds = 80  # number of independent simulations to average over
plot_roughness!(state,param,rng,t_gap,n_pts,n_seeds)