da_residual_ssmi.inc

1 subroutine da_residual_ssmi(oi, y, re,                           &
2                                        np_missing, np_bad_data,             &
3                                        np_obs_used, np_available           )
4 
5    !-----------------------------------------------------------------------
6    ! Purpose: TBD
7    !-----------------------------------------------------------------------
8 
9    implicit none
10 
11    type (ob_type), intent(in)        :: oi     ! Obs increment vector (O-B).
12    type (y_type) , intent(in)        :: y      ! y = H (xa)
13    type (y_type) , intent(out)       :: re     ! Residual structure.
14 
15    integer       , intent(inout)     :: np_available, np_obs_used,          &
16                                         np_missing  , np_bad_data
17 
18    if (oi%num_ssmi_tb .gt. 0) then
19 
20       call da_residual_ssmi_tb(oi, y, re,  &
21                                          np_missing, np_bad_data, &
22                                          np_obs_used, np_available)
23 
24    else if (oi%num_ssmi_retrieval .gt. 0) then
25 
26       call da_residual_ssmi_rv(oi, y, re,  &
27                                          np_missing, np_bad_data, &
28                                          np_obs_used, np_available)
29    end if
30 
31 end subroutine da_residual_ssmi
32 
33