Comments (7)
Yes, it seems to work perfectly now.
from plonk.
Hi @cpinte. Thanks very much for raising this bug.
I noticed that the snapshot file is missing the internal energy. Try the following:
>>> import h5py
>>> f = h5py.File('hydrostatic_18790.h5', mode='r')
>>> list(f['particles'])
['divv', 'dt', 'h', 'itype', 'vxyz', 'xyz']
Is it converted from a standard Phantom file? Perhaps the converter has a bug.
from plonk.
Yes, the dump was converted from a standard phatom dump via phantom2hdf5
from plonk.
Does the standard dump contain the internal energy?
from plonk.
Yes, the internal is in the standard dump. I have added it to the dropbox directory.
from plonk.
I have added a bug fix in ddaaa1f that sets the pressure correctly from the internal energy if it's available. This should fix your issue as long as the internal energy is in the snapshot file.
The temperature was incorrect because the adiabatic sound speed, which is sqrt(gamma * pressure / density)
, was incorrect due to the problem in calculating the pressure.
The problem with the missing internal energy in the converted HDF5 file is in the phantom2hdf5 converter not writing the internal energy. Hopefully we can fix this in Phantom soon. See the ongoing pull request at danieljprice/phantom#39 for more.
from plonk.
Hi @cpinte, are you happy if I close this issue?
from plonk.
Related Issues (20)
- Automate building the Conda and pip packages with Travis CI
- Add convenience function to generate movies HOT 6
- Add Snap method to set units to physical quantities
- Read more Phantom header quantities
- Add function to modify snaps and write a new one
- Should the 'id' array be 'particle_type'? HOT 1
- Sink particles are not rotated with `snap.rotate`
- plonk.visualize.plot passes colorbar_kwargs to imshow HOT 2
- Plot dust density for one-fluid dump HOT 2
- Animation is broken HOT 1
- Creating arbitrary cross sections through SPH data HOT 4
- plonk.visualize has no attribute 'render' HOT 4
- Only the first profile in list is plotted HOT 4
- Error in Column density plots when using logarithmic colour bar. HOT 2
- Change from splash.py to sph-interp
- Add support for more Phantom equations of state
- TypeError: expected dtype object, got 'numpy.dtype[float64]' when animating HOT 7
- Shifting to the centre of the disc HOT 3
- Tracking individual particles HOT 2
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from plonk.