<p><b>duda</b> 2011-01-05 17:23:40 -0700 (Wed, 05 Jan 2011)</p><p>BRANCH COMMIT<br>
<br>
Compute w from rw and rho in test case setup, since rw is re-computed from w <br>
and rho in the init_coupled_diagnostics() routine at the start of a simulation.<br>
<br>
M    src/core_nhyd_atmos/module_test_cases.F<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_nhyd_atmos/module_test_cases.F
===================================================================
--- branches/atmos_physics/src/core_nhyd_atmos/module_test_cases.F        2011-01-05 02:16:18 UTC (rev 675)
+++ branches/atmos_physics/src/core_nhyd_atmos/module_test_cases.F        2011-01-06 00:23:40 UTC (rev 676)
@@ -803,7 +803,15 @@
 
       end do
 
+      ! Compute w from rho and rw
+      do iCell=1,grid%nCells
+         do k=2,grid%nVertLevels
+            state % w % array(k,iCell) = diag % rw % array(k,iCell)     &amp;
+                                       / (fzp(k) * state % rho % array(k-1,iCell) + fzm(k) * state % rho % array(k,iCell))
+         end do
+      end do
 
+
       !
       ! Compute mass fluxes tangential to each edge (i.e., through the faces of dual grid cells)
       !
@@ -1971,7 +1979,15 @@
 
       end do
 
+      ! Compute w from rho and rw
+      do iCell=1,grid%nCells
+         do k=2,grid%nVertLevels
+            state % w % array(k,iCell) = diag % rw % array(k,iCell)     &amp; 
+                                       / (fzp(k) * state % rho % array(k-1,iCell) + fzm(k) * state % rho % array(k,iCell))
+         end do
+      end do
 
+
       do iEdge=1,grid % nEdges
          grid % fEdge % array(iEdge) = 0.
       end do

</font>
</pre>