Reallocate code from Numerical Recipes - problem at run time

[Brooks Moses]

Edmund J Dunlop wrote:
[...]
>         INTEGER, PARAMETER :: SP = KIND(1.0)
>         INTEGER, PARAMETER :: DP = KIND(1.0D0)
>         INTEGER, PARAMETER :: SPC = KIND((1.0,1.0))
>         INTEGER, PARAMETER :: DPC = KIND((1.0D0,1.0D0))
[...]

Completely off the topic of your main question (which appears to have 
been answered), but this is redundant.  There is no separate kind for 
complex numbers and real numbers, and SP and SPC are necessarily 
identical to each other, as are DP and DPC.

- Brooks