var/da/da_tools/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