Added development version of code.
Modified configure file to add LAPACK and BLAS compiltation.
Modified configure file to change bin directory dependent on compiler.
Modified double precision to 15 decimal places for compatibility with LAPACK.
Switched method for timestamp in params.f90 to use the FORTRAN intrinsic
subroutine 'DATE_AND_TIME' for intel/gnu portability.

Added check_file subroutine in params.f90 to replace multiple INQUIRE
statements.

Added checking if input files exist in rho_prop.f90 using check_file
subroutine.
Replaced GNU intrinsics REALPART and IMAGPART with FORTRAN intrinsics REAL and
AIMAG fields.f90 for intel/gnu portability.
	modified:   src/fields.f90
Added new files for post processing. The program progs/post_proc.f90 finds the
coefficients of frequency components for a given sinusoidal function of time.
Too many additions and changes to mention...
Add program to calculate the linear response of a system over a given frequency
interval.
Too many changes/additions to mention - see diffs.
	deleted:    field.dat
	new file:   lr.in
	renamed:    src/check_dependencies -> src/non-fortran/check_dependencies
	new file:   src/non-fortran/lin_resp.sh
	modified:   DATA/data_2level/gma.txt
	modified:   configure
	modified:   pp.in
	modified:   progs/post_proc.f90
	modified:   rho_prop.in
	modified:   src/freq_analysis.f90
	modified:   src/global_params.f90
	modified:   src/params.f90
	modified:   src/post_proc_params.f90
	modified:   src/print_mat_mod.f90
	deleted:    src/print_output_mod.f90
	modified:   src/print_params_mod.f90
	modified:   src/runge_mod.f90
This was meant to be done in last commit...
	modified:   configure
	modified:   pp.in
	modified:   progs/post_proc.f90
	modified:   rho_prop.in
	modified:   src/freq_analysis.f90
	modified:   src/global_params.f90
	modified:   src/params.f90
	modified:   src/post_proc_params.f90
	modified:   src/print_mat_mod.f90
	deleted:    src/print_output_mod.f90
	modified:   src/print_params_mod.f90
	modified:   src/runge_mod.f90
Output stdout of rho_prop.exe when called in lin_resp
	modified:   src/non-fortran/lin_resp.sh
Realised that the openmp loop in runge_mod for commuting matrices was not
properly implemented - removed. Speeds have been improved by a factor of the
number of cores on the machine as the loop was performing the same calculation
on each core! Also removed openmp from default compilation in configure. Will
try to use BLAS mat-mat multiplication next instead.
	modified:   configure
	modified:   src/runge_mod.f90
No time to mention changes - see diffs.
	modified:   lr.in
	modified:   progs/post_proc.f90
	modified:   progs/rho_prop.f90
	modified:   src/fields.f90
	modified:   src/freq_analysis.f90
	modified:   src/gold_dielectric.f90
	modified:   src/non-fortran/lin_resp.sh
	modified:   src/params.f90
	modified:   src/post_proc_params.f90
	modified:   src/print_params_mod.f90
	modified:   src/runge_mod.f90
	deleted:    src/sqd_mnp_params.f90
Modified parameters in the DATA/data_2level/ folder to scale the same as in the
SQD-MNP paper.
	modified:   DATA/data_2level/en.txt
	modified:   DATA/data_2level/gma.txt
	modified:   DATA/data_2level/mu.txt
Modified input file where theta is no longer input as a complex number: the
complex part is dealt with in the main program, so theta is a real number which
is calculated from the inputs m_h and gamma_g.

	modified:   lr.in
Updated the comment in the input file, replacing gold_max_height with m_h
	modified:   lr.in
	modified:   rho_prop.in
Modified src/params.f90 to take in m_h (instead of theta) as a real number.
theta is then calculated.
src/runge_mod.f90 modified to calculate Q_MNP approximately for any waveform.
The a^3 factor is also taken into account in calculation of the P_MNP.
lin_resp.sh: amongst other things, the final chi(omega) is divided by
eps_eff1=2.67777 to correspond with the paper.
Parameters input/output changed.
Modified configure to copy binaries to folder $INPATH.
Modified how lin_resp is compiled.