COMPASS-DOE / grossMethane

analysis code for paired data sets of laboratory-incubations of gross methane flux and associated field net methane and CO2 flux rates
1 stars 1 forks source link

Allow frac_P and frac_k to vary? #32

Closed bpbond closed 1 year ago

bpbond commented 1 year ago

@kendalynnm

Here are our eight 'problem' samples:

Screenshot 2022-12-30 at 8 51 41 AM Screenshot 2022-12-30 at 8 52 02 AM

I tried allowing frac_P and frac_k to vary, i.e. for the optimizer to fit their values within constraints (0.001-0.1 and 0.9-0.99 respectively). This greatly improves the model fit for five of our eight problem samples.

Screenshot 2022-12-30 at 8 57 23 AM

Dashed line is old model (not allowing the fractionation constants to vary) and solid line is the new one:

Screenshot 2022-12-30 at 8 57 15 AM Screenshot 2022-12-30 at 8 59 49 AM Screenshot 2022-12-30 at 8 57 33 AM
bpbond commented 1 year ago

So, interesting I thought; let's talk about this next week!

kendalynnm commented 1 year ago

Very interesting, tomorrow let's look at this and the dropped-first-timestep fits. I'm curious how divergent the final P and k values are for the three fit types (variable fractionation, dropped-first-timestep, and original).

bpbond commented 1 year ago

Outdated!