myrms.f File Reference

Go to the source code of this file.

Functions/Subroutines
program	myrms
	Part of the NADCON5 build process, generates gmtbat05 More...
subroutine	getrms (zs, max, inumlo, inumhi, indx, maxppcell, xlat, xlon, alat, alon, rval)

Function/Subroutine Documentation

subroutine getrms	(	real*8, dimension(max)	zs,
			max,
		integer*4	inumlo,
		integer*4	inumhi,
		integer*4, dimension(max)	indx,
			maxppcell,
		real*8, dimension(max)	xlat,
		real*8, dimension(max)	xlon,
		real*8	alat,
		real*8	alon,
		real*8	rval
	)

Definition at line 1252 of file myrms.f.

Referenced by myrms().

c - As of 10/28/2015, no longer using average
c - lat/lon, but rather center of cell, to
c - register the RMS vector

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

      real*8 alat,alon,rval
  
c     write(6,1)inumlo,inumhi
    1 format('Inside getrms',/,
     *'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 - Because we are NOT attempting any kind of
c - median filter, there is no need to sort
c - the "mini vector" of data in this cell.
c - Just compute the average lat, average lon
c - and RMS value and return.
c - Change as of 10/28/2015:  Do not
c - use the avelat/avelon.  Instead, 
c - register the RMS vector to the center of 
c - the cell.  (See DRU-11, p. 145)

      rval = 0.d0
      alat = 0.d0
      alon = 0.d0

      
      nval = inumhi-inumlo+1
      iq = 0
      do 4 i=inumlo,inumhi
        iq = iq + 1
        rval = rval + zs(indx(i))**2
        alat = alat + xlat(indx(i))
        alon = alon + xlon(indx(i))
    4 continue
      rval = dsqrt(rval/nval) 
      alat = alat/nval
      alon = alon/nval

      return

Here is the caller graph for this function: