Skip to content

Utilities for helioseismology in simulations using the Pencil code

Notifications You must be signed in to change notification settings

Kishore96in/pc_komega_utils

Repository files navigation

Introduction

Utilities to plot dispersion relations and fit modes in simulations that use the Pencil code (https://github.com/pencil-code/pencil-code).

  • power.cached contains a replacement for pc.read.power that caches the results in a HDF5 file (allows faster reads and loading only a part of the array into memory)
  • power.decimate.make_decimated_power allows to take the output of Pencil's power_xy subroutine and save only a subset of it into a HDF5 file. power.decimate.m_pxy_decimated is a mixin that allows the classes in read to use the decimated HDF5 output.
  • power.dat2h5 converts the text files generated by Pencil's power_xy subroutine into HDF5 files.
    • powerxy_to_hdf5: converts power*_xy.dat to HDF5 files.
    • read_power: replacement for pc.read.powers.Power that uses the files generated by powerxy_to_hdf5

Repository: https://bitbucket.org/kishore-phy/pc_komega_utils/

Dependencies

  1. numpy
  2. scipy
  3. matplotlib
  4. pencil
  5. h5py

Known issues

  1. Using dill to save dr_base instances is sometimes prevented by the following bug: uqfoundation/dill#332 (fixed in dill 0.3.8)

About

Utilities for helioseismology in simulations using the Pencil code

Resources

Stars

Watchers

Forks

Packages

No packages published

Languages