parameters
    num_units = {num_units};
    termination_time = {termination_time};
    analysis_time = {analysis_time};
end

species
    U I     [C0]={I_c0}    FW={I_FW}
    M A     [C0]={A_c0}    FW={A_FW}
    M B     [C0]={B_c0}    FW={B_FW}

    P PIA
    P PIB

    P PIAA
    P PIAB
    P PIBA
    P PIBB

    P PAAA
    P PAAB
    P PABA
    P PABB

    P PBAA
    P PBAB
    P PBBA
    P PBBB

    P P_AA PIAA|PAAA|PBAA
    P P_AB PIAB|PAAB|PBAB
    P P_BA PIBA|PABA|PBBA
    P P_BB PIBB|PABB|PBBB

    P D
end

rateconstants
    kpAA = {kpAA}
    kpAB = {kpAB}
    kpBA = {kpBA}
    kpBB = {kpBB}
end

reactions
    
	# ---- Initiation ----
    IN I + A -kpAA-> PIA
    IN I + B -kpAB-> PIB

    PR PIA + A -kpAA-> PIAA
    PR PIA + B -kpAB-> PIAB
    PR PIB + A -kpBA-> PIBA
    PR PIB + B -kpBB-> PIBB

    # ---- Propagation ----
    PR P_AA + A -kpAA-> PAAA
    PR P_AA + B -kpAB-> PAAB
    PR P_AB + A -kpBA-> PABA
    PR P_AB + B -kpBB-> PABB

    PR P_BA + A -kpAA-> PBAA
    PR P_BA + B -kpAB-> PBAB
    PR P_BB + A -kpBA-> PBBA
    PR P_BB + B -kpBB-> PBBB
end