da_calculate_grady_ssmi_rv.inc

1 subroutine da_calculate_grady_ssmi_rv(iv, re, jo_grad_y)
2 
3    !-------------------------------------------------------------------------
4    ! Purpose: Applies obs inverse on re-vector
5    !-------------------------------------------------------------------------
6 
7    implicit none
8 
9    type (iv_type), intent(in)   :: iv          ! Ob Inc. structure.
10    type (y_type), intent(inout) :: re          ! Residual structure.
11    type (y_type), intent(inout) :: jo_grad_y   ! Grad_y(Jo)
12 
13    integer                      :: n
14 
15    if (trace_use_dull) call da_trace_entry("da_calculate_grady_ssmi_rv")
16 
17    do n=1, iv%info(ssmi_rv)%nlocal
18       if (iv%ssmi_rv(n)%speed%qc < obs_qc_pointer) then
19          re%ssmi_rv(n)%speed = 0.0
20       end if
21       jo_grad_y%ssmi_rv(n)%Speed = - re%ssmi_rv(n)%Speed / (iv%ssmi_rv(n)%Speed%error * iv%ssmi_rv(n)%Speed%error)
22 
23       if (iv%ssmi_rv(n)%tpw%qc < obs_qc_pointer) then
24          re%ssmi_rv(n)%tpw = 0.0
25       end if
26       jo_grad_y%ssmi_rv(n)%tpw = -re%ssmi_rv(n)%tpw / (iv%ssmi_rv(n)%tpw%error * iv%ssmi_rv(n)%tpw%error)
27    end do  
28 
29    if (trace_use_dull) call da_trace_exit("da_calculate_grady_ssmi_rv")
30 
31 end subroutine da_calculate_grady_ssmi_rv
32 
33