There is index error in variable rvcuten whose index is set to (k,k) in loops. But it should be (k,i) because similar variable rucuten has the index (k,i). The bug fixed is at Grell-Freitas parameterization block but since it does not provide momentum tendencies, there is no change to simulation results

Author: lo-y-wni

src/core_atmosphere/physics/mpas_atmphys_driver_convection.F

@@ -744,7 +744,7 @@ subroutine convection_from_MPAS(dt_dyn,configs,mesh,sfc_input,diag_physics,tend_
              rqvblten_p(i,k,j)  = rqvblten(k,i)
              rqvdynten_p(i,k,j) = rqvdynten(k,i)
              rucuten_p(i,k,j)   = rucuten(k,i)
-             rvcuten_p(i,k,j)   = rvcuten(k,k)
+             rvcuten_p(i,k,j)   = rvcuten(k,i)
           enddo
        enddo
        enddo
@@ -941,7 +941,7 @@ subroutine convection_to_MPAS(configs,diag_physics,tend_physics,its,ite)
              qc_cu(k,i)    = qccu_p(i,k,j)
              qi_cu(k,i)    = qicu_p(i,k,j)
              rucuten(k,i)  = rucuten_p(i,k,j)
-             rvcuten(k,k)  = rvcuten_p(i,k,j)
+             rvcuten(k,i)  = rvcuten_p(i,k,j)
           enddo
        enddo
        enddo