subroutine sreadA3_RR(pname,X,Y1,Y2)
integer :: i,j,k,Ny1,Ny2,Ny3
character(len=*) :: pname
real(8),dimension(:,:,:) :: Y1
real(8),dimension(size(Y1,1),size(Y1,2),size(Y1,3)),optional :: Y2
real(8),dimension(size(Y1,3)) :: X
Ny1=size(Y1,1) ; Ny2=size(Y1,2); Ny3=size(Y1,3)
open(719,file=adjustl(trim(pname)))
if(present(Y2))then
do i=1,Ny1
do j=1,Ny2
do k=1,Ny3
read(719,*)X(k),Y1(i,j,k),Y2(i,j,k)
enddo
enddo
enddo
else
do i=1,Ny1
do j=1,Ny2
do k=1,Ny3
read(719,*)X(k),Y1(i,j,k)
enddo
enddo
enddo
endif
close(719)
end subroutine sreadA3_RR