master: Reversed Heng's initialization of moist_old in...

master: Reversed Heng's initialization of moist_old in module_big_step_utilities_em.F since v3.8.1 does this from a call in solve_em.

dyn_em/module_big_step_utilities_em.F

@@ -6313,19 +6313,14 @@ SUBROUTINE theta_to_thetam ( t_1 , moist_old ,                &
 
    INTEGER :: i , j , k
 
+   !  First RK loop, this info is from the physics packages.  It is modified immediately after the
+   !  call to the physics schemes, and the remains constant for the remainder of the RK loops.
+
    IF ( rk_step .EQ. 1 ) THEN
       DO j = jts , MIN(jte,jde-1)
          DO k = kts , kte-1
             DO i = its , MIN(ite,ide-1)
-               ! We need to initialize moist_old here for the first RK step at every time step
-               ! (not just the first time step). It should be OK. Later in the first RK step rk_update_scalar
-               ! would do the same assignment to moist_old. - Heng Xiao 10-08-2015
-               moist_old(i,k,j) = moist(i,k,j)
-
-   !  First RK loop, this info is from the physics packages.  It is modified immediately after the
-   !  call to the physics schemes, and the remains constant for the remainder of the RK loops.
-               t_tendf(i,k,j) = (1. + (R_v/R_d) * moist_old(i,k,j))*t_tendf(i,k,j) &
-                              + (R_v/R_d)*(t_1(i,k,j)+T0)*moist_tend(i,k,j)
+               t_tendf(i,k,j) = (1. + (R_v/R_d) * moist_old(i,k,j))*t_tendf(i,k,j) + (R_v/R_d)*(t_1(i,k,j)+T0)*moist_tend(i,k,j)
                h_diabatic(i,k,j) =  (1. + (R_v/R_d) * moist_old(i,k,j))*h_diabatic(i,k,j)
             END DO
          END DO