var/da/da_obs/da_fill_obs_structures_lightning.inc

   1 subroutine da_fill_obs_structures_lightning(iv, ob)
   2
   3    !----------------------------------------------------------------------------
   4    ! Purpose: Allocates observation structure and fills it from iv.
   5    ! Authors: Z Chen (zchen@fjnu.edu.cn), Jenny Sun (NCAR), X Qie (IAP)
   6    !----------------------------------------------------------------------------
   7
   8    implicit none
   9
  10    type (iv_type), intent(inout) :: iv   ! Obs and header structure.
  11    type (y_type), intent(out)    :: ob   ! (Smaller) observation structure.
  12
  13    integer :: n, k     ! Loop counters.
  14    integer :: i,j
  15
  16    if (trace_use) call da_trace_entry("da_fill_obs_structures_lightning")
  17
  18    !---------------------------------------------------------------------------
  19    ! Initialise obs error factors (which will be overwritten in use_obs_errfac)
  20    !---------------------------------------------------------------------------
  21
  22    iv % lightning_ef_w   = 1.0
  23    iv % lightning_ef_div = 1.0
  24    iv % lightning_ef_qv  = 1.0
  25    !----------------------------------------------------------------------
  26    ! [1.0] Allocate innovation vector and observation structures:
  27    !----------------------------------------------------------------------
  28    call da_allocate_y_lightning(iv, ob)
  29
  30    !----------------------------------------------------------------------
  31    ! [2.0] Transfer observations:
  32    !----------------------------------------------------------------------
  33
  34    ! [2.20] Transfer lightning obs:
  35
  36    if (iv%info(lightning)%nlocal > 0) then
  37       do n = 1, iv%info(lightning)%nlocal
  38          do k = 1, iv%info(lightning)%levels(n)
  39             ob%lightning(n)%w(k)  = iv%lightning(n)%w(k)%inv
  40             ob%lightning(n)%div(k)  = iv%lightning(n)%div(k)%inv
  41             ob%lightning(n)%qv(k) = iv%lightning(n)%qv(k)%inv
  42         end do
  43       end do
  44    end if
  45
  46    if (trace_use) call da_trace_exit("da_fill_obs_structures_lightning")
  47
  48 end subroutine da_fill_obs_structures_lightning