search

pypest / pyemu

python modules for model-independent uncertainty analyses, data-worth analyses, and interfacing with PEST(++)
BSD 3-Clause "New" or "Revised" License
174 stars 95 forks source link

JCO2JCO functionality in pyemu #431

Closed Paszka1 closed 1 year ago

Paszka1 commented 1 year ago

Is it possible in pyemu to perform a jco2jco type operation based on a Jacobian matrix in jcb file format? The pyemu.Jco class doesn't seem to provide such a functionality. What I would like to achieve is to modify the Jacobian matrix as a result of changing observation weights in the control file.

mnfienen commented 1 year ago

Hey @Paszka1 - You can definitely change values in jco files within pyemu but the jco or jcb shouldn't have information about weights in it anyway, so sounds like you might not need to ...

Paszka1 commented 1 year ago

Hi @mnfienen - Indeed, that's just what I was thinking as well. What made me to ask this question though, was that in the original PEST manual by John Doherty, there is a list of possible reasons for using JCO2JCO and one of them is altering observation weights in the control file. Is there a particular explanation for that item in the list, or is it a dummy one?

briochh commented 1 year ago

@Paszka1, This might be one for John. I can't think of a reason, off the top of my head.