GridSorting.f90 File Reference

Go to the source code of this file.

Modules
module	gridsorting

Functions/Subroutines
subroutine, public	gridsorting::quicksortgrid (array, arraySize, idxToGlobal)
subroutine	init ()
logical(lgp) function	lessthan (n, m)
subroutine	swap (a, b)
subroutine	rshift (left, right)

Function/Subroutine Documentation

init()

subroutine quicksortgrid::init

private

Definition at line 23 of file GridSorting.f90.

Here is the caller graph for this function:

lessthan()

logical(lgp) function quicksortgrid::lessthan ( integer(i4b), intent(in)  n, integer(i4b), intent(in)  m  )

private

Definition at line 28 of file GridSorting.f90.

      integer(I4B), intent(in) :: n
      integer(I4B), intent(in) :: m
      logical(LGP) :: isLess
      ! local
      type(GlobalCellType), pointer :: gcn, gcm
      real(DP) :: xn, yn, zn, xm, ym, zm
      real(DP), dimension(:), pointer, contiguous :: dis_top_n, dis_bot_n, &
                                                     dis_top_m, dis_bot_m
      real(DP), dimension(:), pointer, contiguous :: dis_xc_n, dis_yc_n, &
                                                     dis_xc_m, dis_yc_m
      real(DP) :: xorigin_n, yorigin_n, angrot_n, &
                  xorigin_m, yorigin_m, angrot_m
 
      ! get coordinates
      gcn => idxtoglobal(array(n))
      gcm => idxtoglobal(array(m))
 
      ! get model data
      ! for n:
      dis_top_n => gcn%v_model%dis_top%get_array()
      dis_bot_n => gcn%v_model%dis_bot%get_array()
      dis_xc_n => gcn%v_model%dis_xc%get_array()
      dis_yc_n => gcn%v_model%dis_yc%get_array()
      xorigin_n = gcn%v_model%dis_xorigin%get()
      yorigin_n = gcn%v_model%dis_yorigin%get()
      angrot_n = gcn%v_model%dis_angrot%get()
      ! for m:
      dis_top_m => gcm%v_model%dis_top%get_array()
      dis_bot_m => gcm%v_model%dis_bot%get_array()
      dis_xc_m => gcm%v_model%dis_xc%get_array()
      dis_yc_m => gcm%v_model%dis_yc%get_array()
      xorigin_m = gcm%v_model%dis_xorigin%get()
      yorigin_m = gcm%v_model%dis_yorigin%get()
      angrot_m = gcm%v_model%dis_angrot%get()
 
      ! convert coordinates
      call dis_transform_xy(dis_xc_n(gcn%index), dis_yc_n(gcn%index), &
                            xorigin_n, yorigin_n, angrot_n, &
                            xn, yn)
      zn = dhalf * (dis_top_n(gcn%index) + &
                    dis_bot_n(gcn%index))
 
      call dis_transform_xy(dis_xc_m(gcm%index), dis_yc_m(gcm%index), &
                            xorigin_m, yorigin_m, angrot_m, &
                            xm, ym)
      zm = dhalf * (dis_top_m(gcm%index) + &
                    dis_bot_m(gcm%index))
 
      ! compare
      if (.not. is_close(zn, zm, 10 * epsilon(zn))) then
        isless = zn > zm
      else if (.not. is_close(yn, ym, 10 * epsilon(yn))) then
        isless = yn > ym
      else if (.not. is_close(xn, xm, 10 * epsilon(xn))) then
        isless = xn < xm
      else
        isless = .false.
      end if
 

Here is the call graph for this function:

rshift()

subroutine quicksortgrid::rshift ( integer, intent(in)  left, integer, intent(in)  right  )

private

Definition at line 102 of file GridSorting.f90.

      integer, intent(in) :: left, right
      integer :: hold
 
      hold = array(right)
      array(left + 1:right) = array(left:right - 1)
      array(left) = hold
 

swap()

subroutine quicksortgrid::swap ( integer, intent(in)  a, integer, intent(in)  b  )

private

Definition at line 91 of file GridSorting.f90.

      integer, intent(in) :: a, b
      integer :: hold
 
      hold = array(a)
      array(a) = array(b)
      array(b) = hold