Merge branch 'NSF-MPAS/physics_big_fix_zd' into hotfix-v8.2.1 (PR #1213)

This merge corrects the calculation of height above ground level (AGL) that
is used in the computation of the 1-km simulated radar reflectivity field in
mpas_atmphys_driver_microphysics.F. Specifically, the fixes in this merge
affect the refl10cm_1km and refl10cm_1km_max fields.

* NSF-MPAS/physics_big_fix_zd:
  * In ./src/core_atmosphere//physics/mpas_atmphys_driver_microphysics.F, corrected...

Author: mgduda

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)