var/da/da_obs/da_store_obs_grid_info.inc

   1 subroutine da_store_obs_grid_info (info)
   2
   3    !-----------------------------------------------------------------------
   4    ! Purpose: this is in parallel of da_store_obs_grid_info_rad but with
   5    !          an extra thinned check to decide proc_domain.
   6    !-----------------------------------------------------------------------
   7
   8    implicit none
   9
  10    type(infa_type), intent(inout) :: info
  11
  12    integer :: n
  13
  14    if (trace_use) call da_trace_entry("da_store_obs_grid_info")
  15
  16    info%proc_domain(:,:) = .false.
  17
  18    do n=1,info%nlocal
  19       if (info%i(1,n) >= its .and. info%i(1,n) <= ite .and. info%j(1,n) >= jts .and. info%j(1,n) <= jte) then
  20          if ( .not. info%thinned(1,n) ) then
  21             info%proc_domain(:,n) = .true.
  22          end if
  23       end if
  24    end do
  25
  26    if (trace_use) call da_trace_exit("da_store_obs_grid_info")
  27
  28 end subroutine da_store_obs_grid_info
  29
  30