Missing observations after pymc is called

openghg / openghg_inversions

University of Bristol Atmospheric Chemistry Research Group RHIME Inversion code (with openghg dependency)

MIT License

5 stars 0 forks source link

Missing observations after pymc is called #127

Closed hdelongueville closed 1 month ago

hdelongueville commented 1 month ago

Some observation points are missing in the hbmcmc raw outputs. The missing points are still present after model scenario is called, suggesting the suppression of points happens in inversion_pymc.py.

hdelongueville commented 1 month ago

I have saved fp_all at each stage of fixedbasisMCMC in hbmcmc.py (see branch Iss127_missing_obs).

The removal of some observations happens at line 490: fp_data[site] = fp_data[site].dropna("time") # pymc doesn't like NaNs

hdelongueville commented 1 month ago

*line 501 with branch Iss127_missing_obs