In the IDL version of the code, there's an additional term in the cost function, which constrains the boxes to be within some fraction of their original distribution. I.e. in the numerator, there's another term like:
D_lat @ S_lat^{-1} @ x
and in the denominator/uncertainty:
D_lat.T @ S_lat^{-1} @ D_lat
Typically, S_lat is just 100% of the prior emissions.
In the IDL version of the code, there's an additional term in the cost function, which constrains the boxes to be within some fraction of their original distribution. I.e. in the numerator, there's another term like:
D_lat @ S_lat^{-1} @ x
and in the denominator/uncertainty:
D_lat.T @ S_lat^{-1} @ D_lat
Typically, S_lat is just 100% of the prior emissions.