main.py: now we first take sparse observations because calculating variance of the noise, as it should be.

Author: ray-chew

RKLM_Python/__main__.py

@@ -120,10 +120,11 @@
 # where are my observations?
 if N > 1:
     obs = dap.load_obs(dap.obs_path)
-    obs_noisy, obs_covar = obs_noiser(obs,dap,rloc)
-    # obs_noisy_interp, obs_mask = sparse_obs_selector(obs_noisy, elem, node, ud, dap)
     # obs_mask, no calculations where entries are True
-    obs_mask = sparse_obs_selector(obs_noisy, elem, node, ud, dap)
+    obs_mask = sparse_obs_selector(obs, elem, node, ud, dap)
+    obs_noisy, obs_covar = obs_noiser(obs,obs_mask,dap,rloc,elem)
+    # obs_noisy_interp, obs_mask = sparse_obs_selector(obs_noisy, elem, node, ud, dap)
+
 
 # add ensemble info to filename
 ud.output_suffix = fn_gen(ud, dap, N)