sreadA4_RC Subroutine

subroutine sreadA4_RC(pname, X, Y1)

Arguments

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

Calls

proc~~sreada4_rc~2~~CallsGraph proc~sreada4_rc~2 sreadA4_RC free_unit free_unit proc~sreada4_rc~2->free_unit ioread_control ioread_control proc~sreada4_rc~2->ioread_control reg reg proc~sreada4_rc~2->reg

Source Code

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