da_check_xtoy_adjoint_airep.inc

1 subroutine da_check_xtoy_adjoint_airep(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_dull) call da_trace_entry("da_check_xtoy_adjoint_airep")
16 
17    do n=1, iv%info(airep)%nlocal
18       do k=1, iv%info(airep)%levels(n)
19          if (iv%info(airep)%proc_domain(k,n)) then
20             adjtest_lhs = adjtest_lhs + &
21                (y%airep(n)%u(k) / typical_u_rms)**2 + &
22                (y%airep(n)%v(k) / typical_v_rms)**2 + &
23                (y%airep(n)%t(k) / typical_t_rms)**2 
24          end if
25       end do
26 
27       do k=1, iv%info(airep)%levels(n)
28          pertile_lhs = pertile_lhs + &
29             (y%airep(n)%u(k) / typical_u_rms)**2 + &
30             (y%airep(n)%v(k) / typical_v_rms)**2 + &
31             (y%airep(n)%t(k) / typical_t_rms)**2
32 
33         y%airep(n)%u(k) = y%airep(n)%u(k) / typical_u_rms ** 2
34         y%airep(n)%v(k) = y%airep(n)%v(k) / typical_v_rms ** 2
35         y%airep(n)%t(k) = y%airep(n)%t(k) / typical_t_rms ** 2
36       end do
37    end do
38    
39    if (trace_use_dull) call da_trace_exit("da_check_xtoy_adjoint_airep")
40 
41 end subroutine da_check_xtoy_adjoint_airep
42 
43