var/da/da_ssmi/da_transform_xtoy_ssmi_rv.inc

   1 subroutine da_transform_xtoy_ssmi_rv(grid, iv, y)
   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. increment vector (o-b).
  11    type (y_type),  intent(inout) :: y        ! y = h (grid%xa)
  12
  13    integer :: n        ! loop counter.
  14
  15    real, allocatable :: tpw(:)
  16    real, allocatable :: speed(:)
  17
  18    if (trace_use) call da_trace_entry("da_transform_xtoy_ssmi_rv")
  19
  20    ! SSMI observation operator y = H(x):
  21
  22    allocate (tpw(iv%info(ssmi_rv)%n1:iv%info(ssmi_rv)%n2))
  23    allocate (speed(iv%info(ssmi_rv)%n1:iv%info(ssmi_rv)%n2))
  24
  25    call da_interp_lin_2d (grid%xa%tpw,   iv%info(ssmi_rv), 1, tpw)
  26    call da_interp_lin_2d (grid%xa%speed, iv%info(ssmi_rv), 1, speed)
  27
  28    do n=iv%info(ssmi_rv)%n1,iv%info(ssmi_rv)%n2
  29       y%ssmi_rv(n)%tpw   = tpw(n)
  30       y%ssmi_rv(n)%speed = speed(n)
  31    end do
  32
  33    deallocate (tpw)
  34    deallocate (speed)
  35
  36    if (trace_use) call da_trace_exit("da_transform_xtoy_ssmi_rv")
  37
  38 end subroutine da_transform_xtoy_ssmi_rv
  39
  40