# === Hydro model # precipitation from precipitation import precipitation_coefficients # distribution for the initial model values from initial import initial # construct Hydro model from hydro import Hydro model = Hydro(precipitation_coefficients, initial, dt=10, integrator='RKSSP', numbify=0) # === SPUX # LIKELIHOOD # construct Particle Filtering likelihood from spux.likelihoods.pf import PF likelihood = PF(particles=[4, 16]) #likelihood = PF (particles=16) # REPLICATES # use Replicates aggregator to combine likelihood with multiple data sets from spux.aggregators.replicates import Replicates replicates = Replicates() # SAMPLER # construct EMCEE sampler