da_residual_new.inc

1 subroutine da_residual_new(yy, qc, inv, rr)
2 
3    !-----------------------------------------------------------------------
4    ! Purpose: TBD
5    !-----------------------------------------------------------------------
6 
7    implicit none
8 
9    integer, intent(in)  :: qc(:,:)
10    real,    intent(in)  :: inv(:,:)
11    real,    intent(in)  :: yy(:,:)
12    real,    intent(out) :: rr(:,:)
13 
14    if (trace_use) call da_trace_entry("da_residual_new")
15   
16    where (qc(:,:) >= obs_qc_pointer)
17       rr(:,:) = inv(:,:) - yy(:,:)
18    elsewhere
19       rr(:,:) = 0.0
20    endwhere
21 
22    if (trace_use) call da_trace_exit("da_residual_new")
23 
24 end subroutine da_residual_new
25 
26