da_check_xtoy_adjoint_radar.inc

1 subroutine da_check_xtoy_adjoint_radar(iv, y, adjtest_lhs, pertile_lhs)
2 
3    !-----------------------------------------------------------------------
4    ! Purpose: TBD
5    !-----------------------------------------------------------------------
6 
7    implicit none
8 
9    type (iv_type), intent(in)     :: iv            ! obs. inc. vector (o-b).
10    type (y_type) , intent(inout)  :: y             ! y = h (xa)
11    real          , intent(inout)  :: adjtest_lhs, pertile_lhs
12 
13    integer :: n, k          ! Loop counter.
14 
15    if (trace_use) call da_trace_entry("da_jo_and_grady_radar")
16 
17    do n=1, iv%info(radar)%nlocal
18       if (iv%info(radar)%proc_domain(1,n)) then
19          do k=1, iv%info(radar)%levels(n)
20             adjtest_lhs = adjtest_lhs + &
21                (y%radar(n)%rv(k)/typical_rv_rms)**2 + (y%radar(n)%rf(k)/typical_rf_rms)**2
22          end do
23       end if
24 
25       do k=1, iv%info(radar)%levels(n)
26          pertile_lhs = pertile_lhs + &
27             (y%radar(n)%rv(k)/typical_rv_rms)**2 + (y%radar(n)%rf(k)/typical_rf_rms)**2
28 
29          y%radar(n)%rv(k)= y%radar(n)%rv(k) / typical_rv_rms ** 2
30          y%radar(n)%rf(k)= y%radar(n)%rf(k) / typical_rf_rms ** 2
31       end do
32    end do
33 
34 end subroutine da_check_xtoy_adjoint_radar
35 
36