11 integer, parameter :: fp = selected_real_kind(13)
12 complex(kind=fp) :: c(:,:)
13 where( abs( aimag( c ) ) < 1.e-10_fp ) &
14 & c = cmplx( real( c , fp ) , 0._fp , fp )
15 where( abs( real( c , fp ) ) < 1.e-10_fp ) &
16 & c = cmplx( 0._fp , aimag( c ) , fp )