var/da/da_statistics/da_maxmin_in_field.inc

   1 subroutine da_maxmin_in_field(field, max, min)
   2
   3    !-----------------------------------------------------------------------
   4    ! Purpose: TBD
   5    !-----------------------------------------------------------------------
   6
   7    implicit none
   8
   9    real,                    intent(in)  :: field(its:ite,jts:jte,kts:kte)
  10    type(maxmin_field_type), intent(out) :: max(kts:kte)
  11    type(maxmin_field_type), intent(out) :: min(kts:kte)
  12
  13    if (trace_use_dull) call da_trace_entry("da_maxmin_in_field")
  14
  15    do k = kts, kte
  16       max(k)%value = -1.0e20
  17       min(k)%value =  1.0e20
  18       do j = jts, jte
  19          do i = its, ite
  20             if (field(i,j,k) > max(k)%value) then
  21                max(k)%value = field(i,j,k)
  22                max(k)%i     = i
  23                max(k)%j     = j
  24             end if
  25             if (field(i,j,k) < min(k)%value) then
  26                min(k)%value = field(i,j,k)
  27                min(k)%i     = i
  28                min(k)%j     = j
  29             end if
  30          end do
  31       end do
  32    end do
  33
  34    if (trace_use_dull) call da_trace_exit("da_maxmin_in_field")
  35
  36 end subroutine da_maxmin_in_field
  37
  38