* In ./src/core_atmosphere//physics/mpas_atmphys_driver_microphysics.F, corrected

  the calculation of the height Above Ground Level (AGL) prior to computing radar
  in subroutine compute_radar_reflectivity.

Author: ldfowler58

src/core_atmosphere/physics/mpas_atmphys_driver_microphysics.F

@@ -705,7 +705,7 @@ subroutine compute_radar_reflectivity(configs,diag_physics,its,ite)
              qs1d(k)  = qs_p(i,k,j)
              qg1d(k)  = qg_p(i,k,j)
              dBZ1d(k) = -35._RKIND
-             zp(k) = z_p(i,k,j) - z_p(i,1,j)+0.5*dz_p(i,1,j) ! height AGL
+             zp(k) = z_p(i,k,j) - z_p(i,1,j) + 0.5*dz_p(i,k,j) ! height AGL
           enddo
 
           call refl10cm_wsm6(qv1d,qr1d,qs1d,qg1d,t1d,p1d,dBZ1d,kts,kte)
@@ -756,7 +756,7 @@ subroutine compute_radar_reflectivity(configs,diag_physics,its,ite)
              qg1d(k)  = qg_p(i,k,j)
              nr1d(k)  = nr_p(i,k,j)
              dBZ1d(k) = -35._RKIND
-             zp(k) = z_p(i,k,j) - z_p(i,1,j)+0.5*dz_p(i,1,j) ! height AGL
+             zp(k) = z_p(i,k,j) - z_p(i,1,j) + 0.5*dz_p(i,k,j) ! height AGL
           enddo
 
           call calc_refl10cm(qv1d,qc1d,qr1d,nr1d,qs1d,qg1d,t1d,p1d,dBZ1d,kts,kte,i,j)