subroutine sreadA4_RC(pname,X,Y1)
integer :: Ny1,Ny2,Ny3,Ny4
integer :: i1,i2,i3,i4
character(len=*) :: pname
complex(8),dimension(:,:,:,:) :: Y1
real(8),dimension(size(Y1,4)) :: X
real(8),dimension(&
size(Y1,1),&
size(Y1,2),&
size(Y1,3),&
size(Y1,4)) :: 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)
!
do i1=1,Ny1
do i2=1,Ny2
do i3=1,Ny3
do i4=1,Ny4
read(unit,*)X(i4),imY(i1,i2,i3,i4),reY(i1,i2,i3,i4)
enddo
enddo
enddo
enddo
close(unit)
Y1=dcmplx(reY,imY)
end subroutine sreadA4_RC