<HTML> <BODY BGCOLOR=#ccccdd LINK=#0000aa VLINK=#0000ff ALINK=#ff0000 ><BASE TARGET="bottom_target"><PRE>
<A NAME='DA_PRINT_STATS_RAD'><A href='../../html_code/radiance/da_print_stats_rad.inc.html#DA_PRINT_STATS_RAD' TARGET='top_target'><IMG SRC="../../gif/bar_red.gif" border=0></A>
subroutine da_print_stats_rad( stats_unit, nchan, rad ) 2,2
!---------------------------------------------------------------------------
! Purpose: print out statistics of omb, oma for radiance data.
!
! METHOD: print out average, rms, minimum, maximum of iv, re
!---------------------------------------------------------------------------
implicit none
integer, intent(in) :: stats_unit, nchan
type (stats_rad_type), intent(in) :: rad
integer :: k, n, nmin, nmax
integer :: used_nchan
if (trace_use) call da_trace_entry
("da_print_stats_rad")
used_nchan = 0
do k=1, nchan !! loop for channels
if(rad%ichan(k)%num > 0) used_nchan = used_nchan + 1
end do
write(unit=stats_unit, fmt='((a,i5))') ' used_nchan: ', used_nchan
write(unit=stats_unit, fmt='(6a)') &
' Channel ', &
' num ', &
' ave ', &
' rms ', &
' min ', &
' max '
do k=1, nchan !! loop for channels
if (rad%ichan(k)%num > 0) then
n = rad%ichan(k)%num
nmin = rad%ichan(k)%minimum%n
nmax = rad%ichan(k)%maximum%n
write(unit=stats_unit, fmt='((i3,i7,4f8.2))') &
k, rad%ichan(k)%num, rad%ichan(k)%ave/real(n), &
sqrt(rad%ichan(k)%rms/real(n)), &
rad%ichan(k)%minimum%value, rad%ichan(k)%maximum%value
end if
end do
if (trace_use) call da_trace_exit("da_print_stats_rad")
end subroutine da_print_stats_rad