subroutine da_print_stats_rad( stats_unit, nchan, rad ) !--------------------------------------------------------------------------- ! 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