da_check_xtoy_adjoint_sound.inc

1 subroutine da_check_xtoy_adjoint_sound(oi, y, adjtest_lhs, pertile_lhs)
2 
3    !-----------------------------------------------------------------------
4    ! Purpose: TBD
5    !-----------------------------------------------------------------------
6 
7    implicit none
8 
9    type (ob_type), intent(in)     :: oi            ! 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 (oi%num_Sound < 1) return
16 
17    do n=1, oi%num_Sound
18       if (oi%sound(n)%loc%proc_domain) then
19          do k=1, oi%sound(n)%info%levels
20             adjtest_lhs = adjtest_lhs + &
21                           (y%Sound(n)%u(k)/typical_u_rms)**2 + &
22                           (y%Sound(n)%v(k)/typical_v_rms)**2 + &
23                           (y%Sound(n)%t(k)/typical_t_rms)**2 + &
24                           (y%Sound(n)%q(k)/typical_q_rms)**2
25          end do
26       end if
27 
28       do k=1, oi%sound(n)%info%levels
29          pertile_lhs = pertile_lhs + &
30                        (y%Sound(n)%u(k)/typical_u_rms)**2 + &
31                        (y%Sound(n)%v(k)/typical_v_rms)**2 + &
32                        (y%Sound(n)%t(k)/typical_t_rms)**2 + &
33                        (y%Sound(n)%q(k)/typical_q_rms)**2
34 
35          y%sound(n)%u(k) = y%sound(n)%u(k) / typical_u_rms ** 2
36          y%sound(n)%v(k) = y%sound(n)%v(k) / typical_v_rms ** 2
37          y%sound(n)%t(k) = y%sound(n)%t(k) / typical_t_rms ** 2
38          y%sound(n)%q(k) = y%sound(n)%q(k) / typical_q_rms ** 2
39       end do
40    end do
41 
42 end subroutine da_check_xtoy_adjoint_sound
43 
44