this is my Headline in a IDL-source code:
pro gamow,t_plasma,z1=z1,z2=z2,a1=a1,a2=a2
; displays gamow peak for input value of t (in K)
; default values for protons
if not keyword_set(z1) then z1=1.
if not keyword_set(z2) then z2=1.
if not keyword_set(a1) then a1=1.
if not keyword_set(a2) then a2=1.
I am executing this in the terminal/console, with, for example:
gamow, 1d8
This works, since then z1=z2=a1=a2 = 1.0. And 1d8 means 100 Million. But, this doesn't work:
gamow, 1d8, 2, 2, 4, 4
why?
Best regards