da_transform_xtoy_ssmi_rv_adj.inc

1 subroutine da_transform_xtoy_ssmi_rv_adj(grid, iv, jo_grad_y, jo_grad_x)
2 
3    !-----------------------------------------------------------------------
4    ! Purpose: TBD
5    !-----------------------------------------------------------------------
6 
7    implicit none
8 
9    type (domain),  intent(in)       :: grid
10    type (iv_type), intent(in)       :: iv          ! obs. inc. vector (o-b).
11    type (y_type) , intent(in)       :: jo_grad_y   ! grad_y(jo)
12    type (x_type) , intent(inout)    :: jo_grad_x   ! grad_x(jo)
13 
14    integer :: n
15 
16    real, allocatable :: tpw(:)
17    real, allocatable :: speed(:)
18 
19    if (trace_use) call da_trace_entry("da_transform_xtoy_ssmi_rv_adj")
20 
21    allocate (tpw(iv%info(ssmi_rv)%n1:iv%info(ssmi_rv)%n2))
22    allocate (speed(iv%info(ssmi_rv)%n1:iv%info(ssmi_rv)%n2))
23 
24    do n=iv%info(ssmi_rv)%n1,iv%info(ssmi_rv)%n2
25       tpw(n)   = jo_grad_y%ssmi_rv(n)%tpw
26       speed(n) = jo_grad_y%ssmi_rv(n)%speed
27    end do
28 
29    call da_interp_lin_2d_adj(jo_grad_x%tpw,   iv%info(ssmi_rv), 1, tpw)
30    call da_interp_lin_2d_adj(jo_grad_x%speed, iv%info(ssmi_rv), 1, speed)
31    
32    deallocate (tpw)
33    deallocate (speed)
34 
35    if (trace_use) call da_trace_exit("da_transform_xtoy_ssmi_rv_adj")
36 
37 end subroutine da_transform_xtoy_ssmi_rv_adj
38 
39