Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Add an option to continue simulation from saved snapshot #89

Open
pdziekan opened this issue Nov 7, 2019 · 4 comments
Open

Add an option to continue simulation from saved snapshot #89

pdziekan opened this issue Nov 7, 2019 · 4 comments

Comments

@pdziekan
Copy link
Contributor

pdziekan commented Nov 7, 2019

Snapshots would be huge, though.

@trontrytel
Copy link
Contributor

It would be awesome to have that. Do you think it would be much different than what you wrote for the piggybacker? We are reading in the velocities and scalars there anyhow?

@trontrytel
Copy link
Contributor

Also, for the bulk schemes we wouldn't need to save that much? If you tell me how to start I could try adding that

@pdziekan
Copy link
Contributor Author

Right, with bulk microphysics snapshots would not bo so large.

With bulk, all variables are stored as blitz arrays, so IO to file can be done the same way as for velocity fields in piggybacking.
You simply open a text file and use << for output and >> for input.

It would be much better to store snapshots in binary files (hdf5), but that would take some more effort.

How this could be implemented (from the top of my head, probably missing a lot of details):
Each snapshot should contain the following arrays:
u, v (only in 3D), w, th, rv, rc, rr
and simulation time.

When running from the snapshot, the arrays should be read after "intcond" method of the case class has been called. Also, timestep should be set to the one from the snapshot.

@trontrytel
Copy link
Contributor

Thanks!. I'll work on it

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants