Add a missing term for (-qv*d(rho_zz)/dt) in updating the tendency for theta_m. (MPAS-Dev#32)

Author: syha

src/core_atmosphere/dynamics/mpas_atm_iau.F

@@ -220,7 +220,8 @@ subroutine atm_add_tend_anal_incr (configs, structs, itimestep, dt, tend_ru, ten
       do i = 1, nCellsSolve
          do k = 1, nVertLevels
             tend_th(k,i) = (1. + rvord * scalars(index_qv,k,i))   * tend_th(k,i) &
-                               + rvord * theta(k,i) * tend_scalars(index_qv,k,i)
+                               + rvord * theta(k,i) * tend_scalars(index_qv,k,i) &
+                               - scalars(index_qv,k,i) * tend_rho(k,i)
             tend_rtheta(k,i) = tend_rtheta(k,i) + tend_th(k,i)
          enddo
       enddo