da_setup_obs_interp_wts.inc

1 subroutine da_setup_obs_interp_wts (iv)
2 
3    !-----------------------------------------------------------------------
4    ! Purpose: TBD
5    !-----------------------------------------------------------------------
6 
7    implicit none
8 
9    type (iv_type), intent(inout)  :: iv         ! Innovation vector (O-B).
10 
11    integer                        :: i         ! Loop counter.
12 
13    if (trace_use) call da_trace_entry("da_setup_obs_interp_wts")
14 
15    do i=1,num_ob_indexes
16       if (i /= radiance .and. iv%info(i)%nlocal > 0) then ! i=22 is radiance , should be excluded here
17          call da_store_obs_grid_info (iv%info(i))
18       end if
19    end do
20 
21    do i = 1, iv % num_inst
22       if (iv % instid(i) % num_rad < 1) cycle
23       call da_store_obs_grid_info (iv%instid(i)%info)
24    end do
25 
26    if (trace_use) call da_trace_exit("da_setup_obs_interp_wts")
27 
28 end subroutine da_setup_obs_interp_wts
29 
30