diff --git a/model/fv_mapz.F90 b/model/fv_mapz.F90
index 8952d23f8..fa8c468d4 100644
--- a/model/fv_mapz.F90
+++ b/model/fv_mapz.F90
@@ -2044,6 +2044,7 @@ subroutine cs_profile(qs, a4, delp, km, i1, i2, iv, kord)
        grat = delp(i,km-1) / delp(i,km)
        q(i,km) = ( 3.*(a4(1,i,k-1)+d4(i)*a4(1,i,k)) - grat*qs(i) - q(i,k-1) )/bet
        q(i,km+1) = qs(i)
+       gam(i,km) = d4(i) / bet
     enddo
 
   else ! all others