#Calculate DCAPE so long as the parcel is cooler than the environment
    for i in range(len(p_unsat)-1,-1, -1):
        if p_unsat[i] >= startp and  tpath_corr[i] < t_unsat_corr[i]:
            DCAPE += (-g*((tpath_corr[i]-t_unsat_corr[i])/t_unsat_corr[i])*(h_unsat[i]-h_unsat[i-1]))
        if p_unsat[i] >= startp and  tpath_corr[i] > t_unsat_corr[i]:
            DCAPE += (-g*((tpath_corr[i]-t_unsat_corr[i])/t_unsat_corr[i])*(h_unsat[i]-h_unsat[i-1]))