da_check_xtoy_adjoint_pseudo.inc

1 subroutine da_check_xtoy_adjoint_pseudo(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             ! Loop counter.
14    
15    if (trace_use_dull) call da_trace_entry("da_check_xtoy_adjoint_pseudo")
16 
17    do n=1, iv%info(pseudo)%nlocal
18       if (iv%info(synop)%proc_domain(1,n)) then
19          adjtest_lhs = adjtest_lhs + (y%pseudo(n)%u/typical_u_rms)**2 + &
20             (y%pseudo(n)%v/typical_v_rms)**2 + &
21             (y%pseudo(n)%v/typical_t_rms)**2 + &
22             (y%pseudo(n)%v/typical_p_rms)**2 + &
23             (y%pseudo(n)%v/typical_q_rms)**2
24       end if
25 
26       pertile_lhs = pertile_lhs + (y%pseudo(n)%u/typical_u_rms)**2 + &
27          (y%pseudo(n)%v/typical_v_rms)**2 + &
28          (y%pseudo(n)%v/typical_t_rms)**2 + &
29          (y%pseudo(n)%v/typical_p_rms)**2 + &
30          (y%pseudo(n)%v/typical_q_rms)**2
31 
32       y%pseudo(n)%u = y%pseudo(n)%u/typical_u_rms ** 2
33       y%pseudo(n)%v = y%pseudo(n)%v/typical_v_rms ** 2
34       y%pseudo(n)%t = y%pseudo(n)%p/typical_t_rms ** 2
35       y%pseudo(n)%p = y%pseudo(n)%v/typical_p_rms ** 2
36       y%pseudo(n)%q = y%pseudo(n)%p/typical_q_rms ** 2
37    end do
38    
39    if (trace_use_dull) call da_trace_exit("da_check_xtoy_adjoint_pseudo")
40 
41 end subroutine da_check_xtoy_adjoint_pseudo
42 
43