z_eq_y_div_x.F

c**********************************************************************
#include "author.inc"
c*    $Id: z_eq_y_div_x.F,v 1.6 1996/02/15 17:36:53 turner Exp $
c*
c*    Divide each element in y by each element in x and store the
c*    result in z.
c*
c*    <PARAMETER LIST>
c*
c*     Input:
c*      n - number of elements
c*      x - x-vector
c*      y - y-vector
c*
c*     Output:
c*      z - destination vector
c*      status - return status
c*
#include "copyright.inc"
c**********************************************************************
      subroutine JT_z_eq_y_div_x (n, x, y, z, status)
      implicit none
c
c ... Input:
      integer n
      real x(n), y(n)
c
c ... Output:
      integer status
      real z(n)
c
c ... Local:
      integer i
      real zero
      parameter (zero=0.0d0)
c
      do i=1,n
       if (x(i) .ne. zero) then
        z(i) = y(i) / x(i)
       else
        y(i) = zero
       endif
      enddo
c
      status = 0
      return
      end