da_transform_xtoy_gpsref_adj.inc

1 subroutine da_transform_xtoy_gpsref_adj(iv, xp, jo_grad_y, jo_grad_x)
2 
3    !-------------------------------------------------------------------------
4    ! Purpose: TBD
5    !-------------------------------------------------------------------------
6 
7    implicit none
8 
9    type (ob_type), intent(in)    :: iv          ! obs. inc vector (o-b).
10    type (xpose_type), intent(in) :: xp          ! Domain decomposition vars.
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        ! Loop counter.
15    integer                       :: i, j     ! Index dimension.
16    real                          :: dx, dxm  !
17    real                          :: dy, dym  !
18    integer                       :: num_levs ! obs vertical levels
19 
20    if (iv%num_gpsref < 1) return
21 
22    do n=iv%ob_numb(iv%current_ob_time-1)%gpsref + 1, &
23       iv%ob_numb(iv%current_ob_time)%gpsref
24 
25       num_levs = iv % gpsref(n) % info % levels
26 
27       i = iv%gpsref(n)%loc%i
28       j = iv%gpsref(n)%loc%j
29       dx = iv%gpsref(n)%loc%dx
30       dy = iv%gpsref(n)%loc%dy
31       dxm = iv%gpsref(n)%loc%dxm
32       dym = iv%gpsref(n)%loc%dym
33 
34       call da_interp_lin_3d_adj(jo_grad_x %ref, xp, i, j, dx, dy, dxm, dym, &
35                              jo_grad_y%gpsref(n)%ref, num_levs, &
36                              iv%gpsref(n)%zk, num_levs)
37    end do
38 
39 end subroutine da_transform_xtoy_gpsref_adj
40 
41