var/da/da_obs/da_fill_obs_structures_radar.inc

   1 subroutine da_fill_obs_structures_radar(iv, ob)
   2
   3    !----------------------------------------------------------------------------
   4    ! Purpose: Allocates observation structure and fills it from iv.
   5    !----------------------------------------------------------------------------
   6
   7    implicit none
   8
   9    type (iv_type), intent(inout) :: iv   ! Obs and header structure.
  10    type (y_type), intent(out)    :: ob   ! (Smaller) observation structure.
  11
  12    integer :: n, k     ! Loop counters.
  13    real    :: rh_error ! RH obs. error.
  14    real    :: q_error  ! q obs. error.
  15    real    :: geometric_h, geopotential_h
  16    integer :: i,j
  17    logical :: outside
  18
  19    if (trace_use) call da_trace_entry("da_fill_obs_structures_radar")
  20
  21    !---------------------------------------------------------------------------
  22    ! Initialise obs error factors (which will be overwritten in use_obs_errfac)
  23    !---------------------------------------------------------------------------
  24
  25    iv % radar_ef_rv = 1.0
  26    iv % radar_ef_rf = 1.0
  27
  28    !----------------------------------------------------------------------
  29    ! [1.0] Allocate innovation vector and observation structures:
  30    !----------------------------------------------------------------------
  31    call da_allocate_y_radar(iv, ob)
  32
  33    !----------------------------------------------------------------------
  34    ! [2.0] Transfer observations:
  35    !----------------------------------------------------------------------
  36
  37    ! [2.18] Transfer radar obs:
  38
  39    if (iv%info(radar)%nlocal > 0) then
  40       do n = 1, iv%info(radar)%nlocal
  41          do k = 1, iv%info(radar)%levels(n)
  42             ! Copy observation variables:
  43             !ob % radar(n) % rv(k) = iv % radar(n) % rv(k) % inv
  44             !ob % radar(n) % rf(k) = iv % radar(n) % rf(k) % inv
  45             if(iv % radar(n) % rv(k) % qc-missing_data>1.0)  then
  46               ob % radar(n) % rv(k) = iv % radar(n) % rv(k) % inv * radar_rv_rscl
  47             endif
  48             if(iv % radar(n) % rf(k) % qc-missing_data>1.0)  then
  49               ob % radar(n) % rf(k) = iv % radar(n) % rf(k) % inv * radar_rf_rscl
  50             endif
  51          end do
  52       end do
  53    end if
  54
  55    if (trace_use) call da_trace_exit("da_fill_obs_structures_radar")
  56
  57 end subroutine da_fill_obs_structures_radar
  58
  59