Shock code current versions – updated in January 2019

Version Access Note
Shock 1.1
rev 89
(January 2019)
Download Associated paper: Godard et al. 2019  

  • computation of C* and CJ stationary shocks
  • change grain treatments – adsorption / erosion widths
  • add thermal desorption mechanisms
  • update chemical network
  • add photoelectric effect for grain charge
  • update photoreaction rates with those given by Heays et al. (2017)
  • update the computation of secondary photon processes
  • add option to set up an external grid of radiation field
  • include radiation field spectrum and radiative transfer
  • integrate photoreaction using cross sections
  • include dust absorption properties and heat capacities
  • compute the grain temperature through thermal balance
  • add radiative pumping of H2 and CO electronic lines
  • compute H2 and CO self-shielding using the FGK approximation
  • change the way of computing H2 heating or cooling
  • modification of dvode for absolute error control on log variables
  • possibility of sorting reactions when computing the derivatives
  • possibility of enforcing elements conserv for unbalanced network
  • add individual species velocities for S and SH
  • introduce an energy criterium to compute the shock size
  • change input / output file – plug ISM services
Shock 1.0
rev 90
(January 2019)
Download Associated paper: Lesaffre et al. 2013  

  • add variables NH2, NH for self-shielding
  • add individual species velocities for CH
  • update treatment of photo-reactions
  • Include analytical self-shielding
  • Allow the interface with dumses
  • Add python fitting routines: in fit/*py
Citation: If you publish results produced by the Paris-Durham Shock code, for statistic reasons, we would appreciate if you could mention it in your publications. A typical sentence could be:
Shock models published in this paper have been produced with the Paris-Durham shock code (Flower & Pineau des Forêts 2003, Lesaffre et al. 2013, Godard et al. 2019,

The output files of the Paris-Durham Shock code are now in HDF5 format and should be read with the dedicated tools: Extractor and Chemistry analyser included in the source code archive.