var/da/da_ssmi/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