I found this expression I do not understand:
fac(nf+2) = real ( ntry, kind = 8 )
I understand, there is an array on the left, but why is there a declaration on the right. ntry is declared previously as integer*4 and now it is kind=8, why?
Should it be understood as: assign ntry to fact(nf+2) element but first convert ntry to real*8? Why don't people then write the code simpler:
fact(nf+2) = DBLE(ntry)
Thanks, - A.