var/da/da_obs/da_obs_missing.inc

   1 subroutine da_obs_missing (obs, miss)
   2
   3    !-----------------------------------------------------------------------
   4    ! Purpose: TBD
   5    !-----------------------------------------------------------------------
   6
   7    implicit none
   8
   9    type(multi_level_type), intent(inout) :: obs
  10    logical,                intent(out)   :: miss
  11
  12    integer :: i
  13
  14    if (trace_use_frequent) call da_trace_entry("da_obs_missing")
  15
  16    miss = .true.
  17
  18    ! PW:
  19    if (abs(obs % loc % pw  % inv - missing_r) > 1.0) then
  20       miss = .false.
  21    else
  22       if (obs % loc % pw  % qc /= missing_data) &
  23          obs % loc % pw  % qc = missing_data
  24    end if
  25
  26    ! SLP:
  27    if (abs(obs % loc % slp % inv - missing_r) > 1.0) then
  28       miss = .false.
  29    else
  30       if (obs % loc % slp % qc /= missing_data) &
  31         obs % loc % slp % qc = missing_data
  32    end if
  33
  34    do i = 1, obs % info % levels
  35       ! U:
  36       if (abs(obs % each(i) % u    % inv - missing_r) > 1.0) then
  37          miss = .false.
  38       else
  39          if (obs % each(i) % u    % qc /= missing_data) &
  40             obs % each(i) % u    % qc = missing_data
  41       end if
  42
  43       ! V:
  44       if (abs(obs % each(i) % v    % inv - missing_r) > 1.0) then
  45          miss = .false.
  46       else
  47          if (obs % each(i) % v    % qc /= missing_data) &
  48             obs % each(i) % v    % qc = missing_data
  49       end if
  50
  51       ! P:
  52       if (abs(obs % each(i) % p    % inv - missing_r) > 1.0) then
  53          miss = .false.
  54       else
  55          if (obs % each(i) % p    % qc /= missing_data) &
  56             obs % each(i) % p    % qc = missing_data
  57       end if
  58
  59       ! T:
  60       if (abs(obs % each(i) % t    % inv - missing_r) > 1.0) then
  61          miss = .false.
  62       else
  63          if (obs % each(i) % t    % qc /= missing_data) &
  64             obs % each(i) % t    % qc = missing_data
  65       end if
  66
  67       ! Q:
  68       if (abs(obs % each(i) % q    % inv - missing_r) > 1.0) then
  69          miss = .false.
  70       else
  71          if (obs % each(i) % q    % qc /= missing_data) &
  72             obs % each(i) % q    % qc = missing_data
  73       end if
  74
  75       ! RH:
  76       if (abs(obs % each(i) % rh  % inv - missing_r) > 1.0) then
  77          miss = .false.
  78       else
  79          if (obs % each(i) % rh  % qc /= missing_data) &
  80             obs % each(i) % rh  % qc = missing_data
  81       end if
  82
  83       ! TD:
  84       if (abs(obs % each(i) % td   % inv - missing_r) > 1.0) then
  85          miss = .false.
  86       else
  87          if (obs % each(i) % td   % qc /= missing_data) &
  88             obs % each(i) % td   % qc = missing_data
  89       end if
  90
  91       ! SPEED:
  92       if (abs(obs % each(i) % speed % inv - missing_r) > 1.0) then
  93          miss = .false.
  94       else
  95          if (obs % each(i) % speed % qc /= missing_data) &
  96             obs % each(i) % speed % qc = missing_data
  97       end if
  98    end do
  99
 100    if (trace_use_frequent) call da_trace_exit("da_obs_missing")
 101
 102 end subroutine da_obs_missing
 103
 104