da_residual.inc

1 function da_residual(n, k, yy, ox, n_bad) result (rr)
2 
3    !-----------------------------------------------------------------------
4    ! Purpose: TBD
5    !-----------------------------------------------------------------------
6 
7    implicit none
8 
9    integer          ,    intent(in) :: n, k
10    type (field_type),    intent(in) :: ox
11    real             ,    intent(in) :: yy
12    type (bad_info_type)  , intent(inout) :: n_bad
13    real                             :: rr
14 
15    ! WHY
16    if (k ==0 .OR. n==0) then
17    end if
18 
19    rr = 0.0
20 
21    if (ox % qc >= obs_qc_pointer) then
22       ! WHY
23       ! n_bad % num % use = n_bad % num % use + 1
24       rr = ox % inv - yy
25       ! else
26       !    if (ox % qc /= -88) then
27       !       n_bad % num % bad = n_bad % num % bad + 1
28       !       n_bad % nn(n_bad % num % bad) = n
29       !       n_bad % kk(n_bad % num % bad) = k
30       !     else
31       !       n_bad % num % miss = n_bad % num % miss + 1
32       !    end if
33    end if
34 
35 end function da_residual
36 
37