var/da/da_radiance/da_calculate_grady_rad.inc

   1
   2 subroutine da_calculate_grady_rad( iv, re, jo_grad_y )
   3 !------------------------------------------------------------------------------
   4 !  Purpose: Calculate Gradient_y for radiance data.
   5 !
   6 !  METHOD:  grad_y = -R^-1 (d - H delta_x ).
   7 !
   8 !  HISTORY: 08/2005 - Creation            Zhiquan Liu
   9 !
  10 !------------------------------------------------------------------------------
  11    implicit none
  12
  13    type (iv_type), intent(in)    :: iv          ! Innovation vector.
  14    type (y_type) , intent(inout) :: re          ! Residual vector.
  15    type (y_type) , intent(inout) :: jo_grad_y   ! Grad_y(Jo)
  16
  17    integer                       :: n, k, i
  18
  19    if (trace_use_dull) call da_trace_entry("da_calculate_grady_rad")
  20    do i =1, iv%num_inst
  21       if ( iv%instid(i)%num_rad < 1 .or. iv%instid(i)%rad_monitoring == monitor_on ) cycle
  22       do n=1, iv%instid(i)%num_rad
  23          do k=1, iv%instid(i)%nchan
  24             if ( iv%instid(i)%tb_qc(k,n) < obs_qc_pointer) re%instid(i)%tb(k,n) = 0.0
  25
  26             jo_grad_y%instid(i)%tb(k,n) = -re%instid(i)%tb(k,n) / (iv%instid(i)%tb_error(k,n) * iv%instid(i)%tb_error(k,n))
  27
  28
  29          end do
  30       end do
  31    end do
  32
  33    if (trace_use_dull) call da_trace_exit("da_calculate_grady_rad")
  34
  35 end subroutine da_calculate_grady_rad
  36