Skip to content
Merged
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion ensemble/ensemble.py
Original file line number Diff line number Diff line change
Expand Up @@ -260,7 +260,7 @@ def calc_prediction(self, enX=None, save_prediction=None):
en_pred = []
pbar = tqdm(enumerate(enX), total=self.ne, **progbar_settings)
for member_index, state in pbar:
en_pred.append(self.sim.run_fwd_sim(state, member_index))
en_pred.append(deepcopy(self.sim.run_fwd_sim(state, member_index)))

# Parallelization on HPC using SLURM
elif self.sim.input_dict.get('hpc', False): # Run prediction in parallel on hpc
Expand Down
4 changes: 2 additions & 2 deletions pipt/loop/assimilation.py
Original file line number Diff line number Diff line change
Expand Up @@ -217,11 +217,11 @@ def run(self):
# always store posterior forcast and state, unless specifically told not to
if 'nosave' not in self.ensemble.keys_da:
try: # first try to save as npz file
np.savez(f'{self.save_folder}/posterior_state_estimate.npz', **self.ensemble.enX)
np.savez(f'{self.save_folder}/posterior_state_estimate.npz', **entools.matrix_to_dict(self.ensemble.enX, self.ensemble.idX))
np.savez(f'{self.save_folder}/posterior_forecast.npz', **{'pred_data': self.ensemble.pred_data})
except: # If this fails, store as pickle
with open(f'{self.save_folder}/posterior_state_estimate.p', 'wb') as file:
pickle.dump(self.ensemble.enX, file)
pickle.dump(entools.matrix_to_dict(self.ensemble.enX, self.ensemble.idX), file)
with open(f'{self.save_folder}/posterior_forecast.p', 'wb') as file:
pickle.dump(self.ensemble.pred_data, file)

Expand Down
Loading