I would like to touch base with someone who would be willing to review a very short program need to calculate some astro information about redshift and c/Ho.
Happy to pay for time, really just need to get over a stumbling block I have with the code.
Really need someone who is pretty switched on to this language, primarily I am flying blind on this one.
No guarantees Petra - it's been about 30 years since I last programmed in Fortran, but it is the one and only proper language. Happy to have a look if you're still stuck.
I've have sent you some emails but one salient issue was that you never
explicitly initialise the variable z.
I could get it to run and produce output by initialising z=0.0 otherwise
with GNU f95 compiler the executable never terminates.
The last item in the first line of output produces a floating point NaN.
There is a division the first time around where the denominator d_L
is zero or pretty nearly to it.
