<p><b>duda</b> 2011-02-08 16:13:39 -0700 (Tue, 08 Feb 2011)</p><p>BRANCH COMMIT<br>
<br>
Add missing iCell index to rho array in hydrostatic balance loop.<br>
<br>
<br>
M    src/core_init_nhyd_atmos/module_test_cases.F<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_init_nhyd_atmos/module_test_cases.F
===================================================================
--- branches/atmos_physics/src/core_init_nhyd_atmos/module_test_cases.F        2011-02-08 20:38:50 UTC (rev 719)
+++ branches/atmos_physics/src/core_init_nhyd_atmos/module_test_cases.F        2011-02-08 23:13:39 UTC (rev 720)
@@ -3239,7 +3239,7 @@
                p_check = pp(k,iCell)
                dz = (zgrid(k,iCell) - zgrid(k-1,iCell))
                pp(k,iCell) = pp(k-1,iCell) - 0.5 * (rr(k,iCell) + rr(k-1,iCell))*gravity*dz &amp;
-                                           - 0.5 * (rho(k,iCell)*scalars(state % index_qv,k,iCell) + rho(k-1)*scalars(state % index_qv,k-1,iCell))*gravity*dz
+                                           - 0.5 * (rho(k,iCell)*scalars(state % index_qv,k,iCell) + rho(k-1,iCell)*scalars(state % index_qv,k-1,iCell))*gravity*dz
                diag % pressure % array(k,iCell) = pp(k,iCell) + ppb(k,iCell)
                p(k,iCell) = (diag % pressure % array(k,iCell) / p0) ** (rgas / cp)
                rho(k,iCell) = diag % pressure % array(k,iCell) / rgas / (p(k,iCell)*t(k,iCell)*(1.0 + 1.61*scalars(state % index_qv,k,iCell)))

</font>
</pre>