Programs and Data for "Breakdown of the semiclassical approximation during the early stages of preheating", Physical Review D92, 083522 (2015)
Note that two of the Fortran programs use the Numerical Recipies routine odebs and its subroutines to numerically solve the mode equation. This can easily be replaced by a call to any differential equation solver.
Initial data file read in by the fortran program for Figure 1
Fortran program for the plots in Figure 2
Fortran program for the bottom plots in Figure 2
Initial data file read in by the fortran program for the plots on the left in Figure 2
Initial data file read in by the fortran program for the plot on the bottom left in Figure 2
Initial data file read in by the fortran program for the plots on the right in Figure 2
Initial data file read in by the fortran program for the plot on the bottom right in Figure 2
Data file used for the plots on the left side of Figure 2
Data file used for the plot on the bottom left side of Figure 2
Data file used for the plot on the bottom left side of Figure 2
Data file for the plots on the right side of Figure 2
Data file used for the plot on the bottom right side of Figure 2
Data file used for the plot on the bottom right side of Figure 2
