var/da/da_ssmi/da_jo_and_grady_ssmi_tb.inc

   1 subroutine da_jo_and_grady_ssmi_tb(iv, re, jo, jo_grad_y)
   2
   3    !-----------------------------------------------------------------------
   4    ! Purpose: TBD
   5    !-----------------------------------------------------------------------
   6
   7    implicit none
   8
   9    type (iv_type), intent(in)    :: iv          ! Ob Inc. structure.
  10    type (y_type),  intent(in)    :: re          ! Residual structure.
  11    type (y_type),  intent(inout) :: jo_grad_y   ! Grad_y(Jo)
  12    type (jo_type), intent(out)   :: jo          ! Obs cost function.
  13
  14    integer :: n
  15
  16    if (trace_use) call da_trace_entry("da_jo_and_grady_ssmi_tb")
  17
  18    jo % ssmi_tb19v = 0.0
  19    jo % ssmi_tb19h = 0.0
  20    jo % ssmi_tb22v = 0.0
  21    jo % ssmi_tb37v = 0.0
  22    jo % ssmi_tb37h = 0.0
  23    jo % ssmi_tb85v = 0.0
  24    jo % ssmi_tb85h = 0.0
  25
  26    do n=1, iv%info(ssmi_tb)%nlocal
  27
  28       jo_grad_y%ssmi_tb(n)%tb19v = - re%ssmi_tb(n)%tb19v / &
  29          (iv%ssmi_tb(n)%tb19v%error * iv%ssmi_tb(n)%tb19v%error)
  30
  31       jo_grad_y%ssmi_tb(n)%tb19h = - re%ssmi_tb(n)%tb19h / &
  32           (iv%ssmi_tb(n)%tb19h%error * iv%ssmi_tb(n)%tb19h%error)
  33
  34       jo_grad_y%ssmi_tb(n)%tb22v = - re%ssmi_tb(n)%tb22v / &
  35          (iv%ssmi_tb(n)%tb22v%error * iv%ssmi_tb(n)%tb22v%error)
  36
  37       jo_grad_y%ssmi_tb(n)%tb37v = - re%ssmi_tb(n)%tb37v / &
  38          (iv%ssmi_tb(n)%tb37v%error * iv%ssmi_tb(n)%tb37v%error)
  39
  40       jo_grad_y%ssmi_tb(n)%tb37h = - re%ssmi_tb(n)%tb37h / &
  41          (iv%ssmi_tb(n)%tb37h%error * iv%ssmi_tb(n)%tb37h%error)
  42
  43       jo_grad_y%ssmi_tb(n)%tb85v = - re%ssmi_tb(n)%tb85v / &
  44          (iv%ssmi_tb(n)%tb85v%error * iv%ssmi_tb(n)%tb85v%error)
  45
  46       jo_grad_y%ssmi_tb(n)%tb85h = - re%ssmi_tb(n)%tb85h / &
  47          (iv%ssmi_tb(n)%tb85h%error * iv%ssmi_tb(n)%tb85h%error)
  48
  49       if (iv%info(ssmi_tb)%proc_domain(1,n)) then
  50          jo % ssmi_tb19v = jo % ssmi_tb19v - re%ssmi_tb(n)%tb19v * jo_grad_y%ssmi_tb(n)%tb19v
  51          jo % ssmi_tb19h = jo % ssmi_tb19h - re%ssmi_tb(n)%tb19h * jo_grad_y%ssmi_tb(n)%tb19h
  52          jo % ssmi_tb22v = jo % ssmi_tb22v - re%ssmi_tb(n)%tb22v * jo_grad_y%ssmi_tb(n)%tb22v
  53          jo % ssmi_tb37v = jo % ssmi_tb37v - re%ssmi_tb(n)%tb37v * jo_grad_y%ssmi_tb(n)%tb37v
  54          jo % ssmi_tb37h = jo % ssmi_tb37h - re%ssmi_tb(n)%tb37h * jo_grad_y%ssmi_tb(n)%tb37h
  55          jo % ssmi_tb85v = jo % ssmi_tb85v - re%ssmi_tb(n)%tb85v * jo_grad_y%ssmi_tb(n)%tb85v
  56          jo % ssmi_tb85h = jo % ssmi_tb85h - re%ssmi_tb(n)%tb85h * jo_grad_y%ssmi_tb(n)%tb85h
  57       end if
  58    end do
  59
  60    jo % ssmi_tb19h = 0.5 * jo % ssmi_tb19h
  61    jo % ssmi_tb19v = 0.5 * jo % ssmi_tb19v
  62    jo % ssmi_tb22v = 0.5 * jo % ssmi_tb22v
  63    jo % ssmi_tb37h = 0.5 * jo % ssmi_tb37h
  64    jo % ssmi_tb37v = 0.5 * jo % ssmi_tb37v
  65    jo % ssmi_tb85h = 0.5 * jo % ssmi_tb85h
  66    jo % ssmi_tb85v = 0.5 * jo % ssmi_tb85v
  67
  68    if (trace_use) call da_trace_exit("da_jo_and_grady_ssmi_tb")
  69
  70 end subroutine da_jo_and_grady_ssmi_tb
  71
  72