sreadA5_RC Subroutine

subroutine sreadA5_RC(pname, X, Y1)

Arguments

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

Calls

proc~~sreada5_rc~2~~CallsGraph proc~sreada5_rc~2 sreadA5_RC free_unit free_unit proc~sreada5_rc~2->free_unit ioread_control ioread_control proc~sreada5_rc~2->ioread_control reg reg proc~sreada5_rc~2->reg

Source Code

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