sreadA3_RC Subroutine

subroutine sreadA3_RC(pname, X, Y1)

Arguments

Type IntentOptional Attributes Name
character(len=*) :: pname
real(kind=8), dimension(size(Y1,3)) :: X
complex(kind=8), dimension(:,:,:) :: Y1

Calls

proc~~sreada3_rc~2~~CallsGraph proc~sreada3_rc~2 sreadA3_RC free_unit free_unit proc~sreada3_rc~2->free_unit ioread_control ioread_control proc~sreada3_rc~2->ioread_control reg reg proc~sreada3_rc~2->reg

Source Code

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