splotA3_IC Subroutine

subroutine splotA3_IC(pname, X, Y1, Y2, append)

Arguments

Type IntentOptional Attributes Name
character(len=*) :: pname
integer, dimension(size(Y1,3)) :: X
complex(kind=8), dimension(:,:,:) :: Y1
complex(kind=8), optional, dimension(size(Y1,1),size(Y1,2),size(Y1,3)) :: Y2
logical, optional :: append

Calls

proc~~splota3_ic~~CallsGraph proc~splota3_ic splotA3_IC dimag dimag proc~splota3_ic->dimag

Source Code

subroutine splotA3_IC(pname,X,Y1,Y2,append)
  integer                                                         :: i,j,k,Ny1,Ny2,Ny3
  character(len=*)                                                :: pname
  complex(8),dimension(:,:,:)                                     :: Y1
  complex(8),dimension(size(Y1,1),size(Y1,2),size(Y1,3)),optional :: Y2
  integer,dimension(size(Y1,3))                                   :: X
  logical,optional                    :: append
  logical                             :: check,rw
  rw=.false.;if(present(append))rw=append
  Ny1=size(Y1,1) ; Ny2=size(Y1,2); Ny3=size(Y1,3)
  open(719,file=adjustl(trim(pname)))
  if(present(Y2))then
     do i=1,Ny1
        do j=1,Ny2
           do k=1,Ny3
              write(719,"(I15,4(F21.12))")X(k),dimag(Y1(i,j,k)),dreal(Y1(i,j,k)),dimag(Y2(i,j,k)),dreal(Y2(i,j,k))
           enddo
           write(719,*)""
        enddo
     enddo
  else
     do i=1,Ny1
        do j=1,Ny2
           do k=1,Ny3
              write(719,*)X(k),dimag(Y1(i,j,k)),dreal(Y1(i,j,k))
           enddo
           write(719,*)""
        enddo
     enddo
  endif
  close(719)
end subroutine splotA3_IC