Hi Petra,
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.
I added other recommendations in the emails.
Best Regards
Gary
|