<HTML> <BODY BGCOLOR=#ccccdd LINK=#0000aa VLINK=#0000ff ALINK=#ff0000 ><BASE TARGET="bottom_target"><PRE>
<A NAME='DA_PRINT_STATS_GPSPW'><A href='../../html_code/gpspw/da_print_stats_gpspw.inc.html#DA_PRINT_STATS_GPSPW' TARGET='top_target'><IMG SRC="../../gif/bar_red.gif" border=0></A>
subroutine da_print_stats_gpspw(stats_unit, ntpw, Gpspw) 2,2
!-----------------------------------------------------------------------
! Purpose: TBD
!-----------------------------------------------------------------------
implicit none
integer, intent(in) :: stats_unit
integer, intent(inout) :: ntpw
type (stats_gpspw_type), intent(in) :: gpspw
if (trace_use_dull) call da_trace_entry
("da_residual_gpspw")
if (ntpw > 0) then
if (use_gpspwObs) then
write(unit=stats_unit, fmt='(a/)') ' var tpw(cm) n'
else if (use_gpsztdObs) then
write(unit=stats_unit, fmt='(a/)') ' var ztd(cm) n'
endif
write(unit=stats_unit, fmt='(a,i14)') ' Number: ', ntpw
write(unit=stats_unit, fmt='(a, f12.4,i5)') &
' Minimum(n): ', gpspw%minimum%tpw%value, &
gpspw%minimum%tpw%n , &
' Maximum(n): ', gpspw%maximum%tpw%value, &
gpspw%maximum%tpw%n
write(unit=stats_unit, fmt='(a, f12.4,5x)') &
' Average : ', gpspw%average%tpw/real(ntpw), &
' RMSE : ', sqrt(gpspw%rms_err%tpw/real(ntpw))
end if
if (trace_use_dull) call da_trace_exit("da_residual_gpspw")
end subroutine da_print_stats_gpspw