splotA3_RR Subroutine

subroutine splotA3_RR(pname,X,Y1,Y2,append)
  integer                                                      :: i,j,k,Ny1,Ny2,Ny3
  character(len=*)                                             :: pname
  real(8),dimension(:,:,:)                                     :: Y1
  real(8),dimension(size(Y1,1),size(Y1,2),size(Y1,3)),optional :: Y2
  real(8),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,*)X(k),Y1(i,j,k),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),Y1(i,j,k)
           enddo
           write(719,*)""
        enddo
     enddo
  endif
  close(719)
end subroutine splotA3_RR