mymedian5.f File Reference

Go to the source code of this file.

Functions/Subroutines
program	mymedian5
	Program to filter Map Data for GMT Plotting. More...
subroutine	getmedian (zs, max, inumlo, inumhi, indx, maxppcell, dist, iiival)

Function/Subroutine Documentation

subroutine getmedian	(	real*8, dimension(max)	zs,
			max,
		integer*4	inumlo,
		integer*4	inumhi,
		integer*4, dimension(max)	indx,
			maxppcell,
		real*8, dimension(max)	dist,
			iiival
	)

Definition at line 1334 of file mymedian5.f.

References indexxd().

Referenced by mymedian5().

      real*8 zs(max),dist(max)
      integer*4 indx(max),indx2(maxppcell)
      integer*4 inumlo,inumhi

      real*8 zz(maxppcell)
  

c     write(6,1)inumlo,inumhi
c   1 format('Inside getmedian',/,
c    *'inumlo,inumhi = ',i8,1x,i8)

c - When the whole set of "zs" data is sorted by
c - "whichcell", then the "indx" values are our
c - key to that sort.  The "zs" data remains
c - in "read order" (1-nkt), but the indx values   
c - (also in "read order" point to the sorted
c - order.

c     write(6,*) ' UNSORTED data: '

      nval = inumhi-inumlo+1
      iq = 0
      do 2 i=inumlo,inumhi
        iq = iq + 1
        zz(iq) = zs(indx(i))
c       write(6,3)i,iq,zz(iq)
    2 continue
    3 format(i10,1x,i10,1x,f20.10)

c - Sort this mini vector
c     call hpsort(nval,maxppcell,zz)

      call indexxd(nval,maxppcell,zz,indx2)

c     write(6,*) ' SORTED data: '

c     do 4 i=1,nval
c       write(6,3)i,indx2(i),zz(indx2(i))
c   4 continue
    5 format(i10,1x,f20.10)

c - True median comes from an odd number of values.  If
c - it is even, take the one which sits closest to the
c - center of the cell.
      if(mod(nval,2).eq.1)then
        ival = (nval+1)/2
      else
        ival1 = nval/2
        ival2 = ival1 + 1
c       write(6,*) ' MedianLo = ',zz(indx2(ival1))
c       write(6,*) ' MedianHi = ',zz(indx2(ival2))
c       write(6,*) ' Dist  Lo = ',dist(indx2(ival1))
c       write(6,*) ' Dist  Hi = ',dist(indx2(ival2))
        if(dist(indx2(ival1)).lt.dist(indx2(ival2)))then
c         write(6,*) ' Median Set to LO = ',zz(indx2(ival1))
          ival = ival1 
        else 
c         write(6,*) ' Median Set to HI = ',zz(indx2(ival2))
          ival = ival2
        endif
      endif

c - Note that "iiival" is the read-order number of the
c - point chosen as median for the cell we're analyzing 
c - in this subroutine.
c       write(6,*) ' Median = ',zz(indx2(ival))
        igood1 = indx2(ival)
c       write(6,*) ' ival = ',ival
c       write(6,*) ' indx2(ival) = ',indx2(ival)
        iival = inumlo + indx2(ival) - 1
c       write(6,*) ' iival = ',iival
        iiival = indx(iival)
c       write(6,*) ' indx(iival) = ',indx(iival)


      return

Here is the call graph for this function:

Here is the caller graph for this function: