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(is:ie,js:je,ks:ke)
10 
11    type (maxmin_field_type), intent(out) :: max(xp%kts:xp%kte)
12    type (maxmin_field_type), intent(out) :: min(xp%kts:xp%kte)
13 
14    do k = ks, ke
15       max(k)%value = -1.E20
16       min(k)%value =  1.E20
17       do j = js, je
18          do i = is, ie
19             if (field(i,j,k) > max(k)%value) then
20                max(k)%value = field(i,j,k)
21                max(k)%i     = i
22                max(k)%j     = j
23             end if
24             if (field(i,j,k) < min(k)%value) then
25                min(k)%value = field(i,j,k)
26                min(k)%i     = i
27                min(k)%j     = j
28             end if
29          end do
30       end do
31    end do
32 
33 end subroutine da_maxmin_in_field
34 
35