diff --git a/physics/maximum_hourly_diagnostics.F90 b/physics/maximum_hourly_diagnostics.F90
index 615c49bed..1486ac027 100644
--- a/physics/maximum_hourly_diagnostics.F90
+++ b/physics/maximum_hourly_diagnostics.F90
@@ -30,7 +30,8 @@ subroutine maximum_hourly_diagnostics_run(im, levs, reset, lradar, imp_physics,
                                              imp_physics_fer_hires,con_g, phil,                    &
                                              gt0, refl_10cm, refdmax, refdmax263k, u10m, v10m,     &
                                              u10max, v10max, spd10max, pgr, t2m, q2m, t02max,      &
-                                             t02min, rh02max, rh02min, errmsg, errflg)
+                                             t02min, rh02max, rh02min, dtp, rain, pratemax,        &
+                                             errmsg, errflg)
 
        ! Interface variables
        integer, intent(in) :: im, levs
@@ -54,6 +55,9 @@ subroutine maximum_hourly_diagnostics_run(im, levs, reset, lradar, imp_physics,
        real(kind_phys), intent(inout) :: t02min(:)
        real(kind_phys), intent(inout) :: rh02max(:)
        real(kind_phys), intent(inout) :: rh02min(:)
+       real(kind_phys), intent(in   ) :: dtp
+       real(kind_phys), intent(in   ) :: rain(im)
+       real(kind_phys), intent(inout) :: pratemax(im)
        character(len=*), intent(out)  :: errmsg
        integer, intent(out)           :: errflg
 
@@ -96,6 +100,7 @@ subroutine maximum_hourly_diagnostics_run(im, levs, reset, lradar, imp_physics,
              t02min(i)   = 999.
              rh02max(i)  = -999.
              rh02min(i)  = 999.
+             pratemax(i) = 0.
           enddo
        endif
        do i=1,im
@@ -119,6 +124,7 @@ subroutine maximum_hourly_diagnostics_run(im, levs, reset, lradar, imp_physics,
           rh02min(i) = min(rh02min(i),rh02)
           t02max(i)  = max(t02max(i),t2m(i))  !<--- hourly max 2m t
           t02min(i)  = min(t02min(i),t2m(i))  !<--- hourly min 2m t
+          pratemax(i) = max(pratemax(i),(3.6E6/dtp)*rain(i))
        enddo
 
    end subroutine maximum_hourly_diagnostics_run
diff --git a/physics/maximum_hourly_diagnostics.meta b/physics/maximum_hourly_diagnostics.meta
index 67b3df039..0f87c86af 100644
--- a/physics/maximum_hourly_diagnostics.meta
+++ b/physics/maximum_hourly_diagnostics.meta
@@ -233,6 +233,33 @@
   kind = kind_phys
   intent = inout
   optional = F
+[dtp]
+  standard_name = time_step_for_physics
+  long_name = physics timestep
+  units = s
+  dimensions = ()
+  type = real
+  kind = kind_phys
+  intent = in
+  optional = F
+[rain]
+  standard_name = lwe_thickness_of_precipitation_amount_on_dynamics_timestep
+  long_name = total rain at this time step
+  units = m
+  dimensions = (horizontal_loop_extent)
+  type = real
+  kind = kind_phys
+  intent = inout
+  optional = F
+[pratemax]
+  standard_name = maximum_precipitation_rate_over_maximum_hourly_time_interval
+  long_name = maximum precipitation rate over maximum hourly time interval
+  units = mm h-1
+  dimensions = (horizontal_loop_extent)
+  type = real
+  kind = kind_phys
+  intent = inout
+  optional = F
 [errmsg]
   standard_name = ccpp_error_message
   long_name = error message for error handling in CCPP