subroutine sreadA5_RC(pname,X,Y1)
integer :: Ny1,Ny2,Ny3,Ny4,Ny5
integer :: i1,i2,i3,i4,i5
character(len=*) :: pname
complex(8),dimension(:,:,:,:,:) :: Y1
real(8),dimension(size(Y1,5)) :: X
real(8),dimension(&
size(Y1,1),&
size(Y1,2),&
size(Y1,3),&
size(Y1,4),&
size(Y1,5)) :: reY,imY
!
call ioread_control(pname,control)
open(free_unit(unit),file=reg(pname))
!
!
Ny1=size(Y1,1)
Ny2=size(Y1,2)
Ny3=size(Y1,3)
Ny4=size(Y1,4)
Ny5=size(Y1,5)
!
do i1=1,Ny1
do i2=1,Ny2
do i3=1,Ny3
do i4=1,Ny4
do i5=1,Ny5
read(unit,*)X(i5),imY(i1,i2,i3,i4,i5),reY(i1,i2,i3,i4,i5)
enddo
enddo
enddo
enddo
enddo
close(unit)
Y1=dcmplx(reY,imY)
end subroutine sreadA5_RC