<p><b>duda</b> 2010-08-09 16:21:42 -0600 (Mon, 09 Aug 2010)</p><p>BRANCH COMMIT<br>
<br>
1) Remove extra dts factor in flux used for omega update<br>
<br>
2) Initialize coftz at bottom and top levels to 0 in compute_vert_imp_coefs()<br>
<br>
<br>
M    src/core_nhyd_atmos/module_time_integration.F<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_nonhydrostatic/src/core_nhyd_atmos/module_time_integration.F
===================================================================
--- branches/atmos_nonhydrostatic/src/core_nhyd_atmos/module_time_integration.F        2010-08-06 21:51:07 UTC (rev 470)
+++ branches/atmos_nonhydrostatic/src/core_nhyd_atmos/module_time_integration.F        2010-08-09 22:21:42 UTC (rev 471)
@@ -570,11 +570,13 @@
         do k=2,nVertLevels
           cofwr(k,i) =.5*dtseps*gravity*(fzm(k)*zz(k,i)+fzp(k)*zz(k-1,i))
         end do
+        coftz(1,i) = 0.0
         do k=2,nVertLevels
            cofwz(k,i) = dtseps*c2*(fzm(k)*zz(k,i)+fzp(k)*zz(k-1,i))  &amp;
                 *rdzu(k)*cqw(k,i)*(fzm(k)*p (k,i)+fzp(k)*p (k-1,i))
            coftz(k,i) = dtseps*   (fzm(k)*t (k,i)+fzp(k)*t (k-1,i))
         end do
+        coftz(nVertLevels+1,i) = 0.0
         do k=1,nVertLevels
 
           qtot = 0.
@@ -801,7 +803,7 @@
           end do
 
           do k=2,nVertLevels
-            flux =  dts*0.5*dvEdge(iEdge)*((zgrid(k,cell2)-zgrid(k,cell1))*(fzm(k)*du(k)+fzp(k)*du(k-1))  )
+            flux =  0.5*dvEdge(iEdge)*((zgrid(k,cell2)-zgrid(k,cell1))*(fzm(k)*du(k)+fzp(k)*du(k-1))  )
             flux2 =  - (fzm(k)*zz(k  ,cell2) +fzp(k)*zz(k-1,cell2))*flux/AreaCell(cell2)
             flux1 =  - (fzm(k)*zz(k  ,cell1) +fzp(k)*zz(k-1,cell1))*flux/AreaCell(cell1)
             ws(k,cell2) = ws(k,cell2) + flux2

</font>
</pre>