Hi,
very new to Fortran but have written the code attached. However it doesn't seem to run with double precision unless I use the old style method of writing 0d0 and implicit double precision at the top of the subroutines. I have tried to use a module and declare everything in it to be double...