sreadA6_RC Subroutine

subroutine sreadA6_RC(pname, X, Y1)

Arguments

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

Calls

proc~~sreada6_rc~2~~CallsGraph proc~sreada6_rc~2 sreadA6_RC free_unit free_unit proc~sreada6_rc~2->free_unit ioread_control ioread_control proc~sreada6_rc~2->ioread_control reg reg proc~sreada6_rc~2->reg

Source Code

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