var/da/da_radiance/da_print_stats_rad.inc

   1 subroutine da_print_stats_rad( stats_unit, nchan, rad )
   2
   3    !---------------------------------------------------------------------------
   4    !  Purpose: print out statistics of omb, oma for radiance data.
   5    !
   6    !  METHOD:  print out average, rms, minimum, maximum of iv, re
   7    !---------------------------------------------------------------------------
   8
   9    implicit none
  10
  11    integer,           intent(in)    :: stats_unit, nchan
  12    type (stats_rad_type), intent(in)    :: rad
  13
  14    integer    :: k, n, nmin, nmax
  15    integer    :: used_nchan
  16
  17    if (trace_use) call da_trace_entry("da_print_stats_rad")
  18
  19    used_nchan = 0
  20    do k=1, nchan                  !! loop for channels
  21       if(rad%ichan(k)%num > 0) used_nchan = used_nchan + 1
  22    end do
  23    write(unit=stats_unit, fmt='((a,i5))')  ' used_nchan: ', used_nchan
  24
  25    write(unit=stats_unit, fmt='(6a)') &
  26         ' Channel ', &
  27         ' num  ', &
  28         ' ave  ', &
  29         ' rms  ', &
  30         ' min  ', &
  31         ' max  '
  32
  33    do k=1, nchan                  !! loop for channels
  34       if (rad%ichan(k)%num > 0) then
  35          n    = rad%ichan(k)%num
  36          nmin = rad%ichan(k)%minimum%n
  37          nmax = rad%ichan(k)%maximum%n
  38
  39          write(unit=stats_unit, fmt='((i3,i7,4f8.2))') &
  40             k, rad%ichan(k)%num, rad%ichan(k)%ave/real(n), &
  41             sqrt(rad%ichan(k)%rms/real(n)), &
  42             rad%ichan(k)%minimum%value, rad%ichan(k)%maximum%value
  43       end if
  44    end do
  45
  46    if (trace_use) call da_trace_exit("da_print_stats_rad")
  47
  48 end subroutine da_print_stats_rad
  49
  50