var/da/da_mtgirs/da_jo_and_grady_mtgirs.inc

   1 subroutine da_jo_and_grady_mtgirs(iv, re, jo, jo_grad_y)
   2
   3    !-----------------------------------------------------------------------
   4    ! Purpose: TBD
   5    !-----------------------------------------------------------------------
   6
   7    implicit none
   8
   9    type (iv_type), intent(in)    :: iv         ! Innovation vector.
  10    type (y_type),  intent(in)    :: re         ! Residual vector.
  11    type (y_type),  intent(inout) :: jo_grad_y  ! Grad_y(Jo)
  12    type (jo_type), intent(inout) :: jo         ! Obs cost function.
  13
  14    integer                      :: n, k
  15    ! the following "global" objects are used only when testing
  16    type (iv_type) :: iv_glob         ! Global Innovation vector (O-B).
  17    type (y_type)  :: re_glob         ! Global Residual vector (O-A).
  18    type (y_type)  :: jo_grad_y_glob  ! Global Grad_y(Jo)
  19
  20    if (trace_use_dull) call da_trace_entry("da_jo_and_grady_mtgirs")
  21
  22    jo % mtgirs_u = 0.0
  23    jo % mtgirs_v = 0.0
  24    jo % mtgirs_t = 0.0
  25    jo % mtgirs_q = 0.0
  26
  27    if (test_dm_exact) then
  28       if (iv%info(mtgirs)%ntotal == 0) then
  29          if (trace_use_dull) call da_trace_exit("da_jo_and_grady_mtgirs")
  30          return
  31       end if
  32    else
  33       if (iv%info(mtgirs)%nlocal < 1) then
  34          if (trace_use_dull) call da_trace_exit("da_jo_and_grady_mtgirs")
  35          return
  36       end if
  37    end if
  38
  39    do n=1, iv%info(mtgirs)%nlocal
  40        do k=1, iv%info(mtgirs)%levels(n)
  41           jo_grad_y%mtgirs(n)%u(k) = -re%mtgirs(n)%u(k) / (iv%mtgirs(n)%u(k)%error * iv%mtgirs(n)%u(k)%error)
  42           jo_grad_y%mtgirs(n)%v(k) = -re%mtgirs(n)%v(k) / (iv%mtgirs(n)%v(k)%error * iv%mtgirs(n)%v(k)%error)
  43           jo_grad_y%mtgirs(n)%t(k) = -re%mtgirs(n)%t(k) / (iv%mtgirs(n)%t(k)%error * iv%mtgirs(n)%t(k)%error)
  44           jo_grad_y%mtgirs(n)%q(k) = -re%mtgirs(n)%q(k) / (iv%mtgirs(n)%q(k)%error * iv%mtgirs(n)%q(k)%error)
  45       end do
  46    end do
  47
  48    ! Bitwise-exact reduction preserves operation order of serial code for
  49    ! testing, at the cost of much-increased run-time.  Turn it off when not
  50    ! testing.  This will always be .false. for a serial or 1-MPI-process run.
  51    if (test_dm_exact) then
  52       ! collect all obs in serial order and allocate global objects
  53 !      call da_to_global_mtgirs(iv, re, jo_grad_y, iv_glob, re_glob, jo_grad_y_glob)
  54       ! perform remaining computations
  55       call da_jo_mtgirs_uvtq(iv_glob, re_glob, jo_grad_y_glob, jo)
  56       ! free global objects
  57 !      call da_deallocate_global_mtgirs(iv_glob, re_glob, jo_grad_y_glob)
  58    else
  59       ! perform remaining computations
  60       call da_jo_mtgirs_uvtq(iv, re, jo_grad_y, jo)
  61    end if
  62
  63    jo % mtgirs_u = 0.5 * jo % mtgirs_u
  64    jo % mtgirs_v = 0.5 * jo % mtgirs_v
  65    jo % mtgirs_t = 0.5 * jo % mtgirs_t
  66    jo % mtgirs_q = 0.5 * jo % mtgirs_q
  67
  68    if (trace_use_dull) call da_trace_exit("da_jo_and_grady_mtgirs")
  69
  70 end subroutine da_jo_and_grady_mtgirs
  71
  72