subroutine da_print_stats_gpspw(stats_unit, ntpw, Gpspw) !----------------------------------------------------------------------- ! 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