Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
character(len=*) | :: | pname | ||||
real(kind=8), | dimension(size(Y1,3)) | :: | X | |||
complex(kind=8), | dimension(:,:,:) | :: | Y1 |
subroutine sreadA3_RC(pname,X,Y1) integer :: i,j,k,Ny1,Ny2,Ny3 character(len=*) :: pname complex(8),dimension(:,:,:) :: Y1 real(8),dimension(size(Y1,3)) :: X real(8),dimension(size(Y1,1),size(Y1,2),size(Y1,3)) :: 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) ! do i=1,Ny1 do j=1,Ny2 do k=1,Ny3 read(unit,*)X(k),imY(i,j,k),reY(i,j,k) enddo enddo enddo close(unit) Y1=dcmplx(reY,imY) end subroutine sreadA3_RC