splotM_RI Subroutine

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