<p><b>duda</b> 2010-05-03 14:44:36 -0600 (Mon, 03 May 2010)</p><p>BRANCH COMMIT<br>
<br>
Add calls to vector reconstruction module to interpolate<br>
full velocity vector at cell locations for each time step.<br>
<br>
M    src/core_hyd_atmos/mpas_interface.F<br>
M    src/core_hyd_atmos/module_time_integration.F<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_hyd_atmos/module_time_integration.F
===================================================================
--- branches/atmos_physics/src/core_hyd_atmos/module_time_integration.F        2010-05-03 20:35:52 UTC (rev 236)
+++ branches/atmos_physics/src/core_hyd_atmos/module_time_integration.F        2010-05-03 20:44:36 UTC (rev 237)
@@ -4,6 +4,7 @@
    use configure
    use constants
    use dmpar
+   use vector_reconstruction
 
 
    contains
@@ -314,6 +315,17 @@
       ! END RK loop 
       !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 
 
+
+      !
+      ! Compute full velocity vectors at cell centers
+      !
+      block =&gt; domain % blocklist
+      do while (associated(block))
+         call reconstruct(block % time_levs(2) % state, block % mesh)
+         block =&gt; block % next
+      end do
+
+
       if(debug) write(0,*) ' rk step complete - mass diagnostics '
 
       if(debug .or. debug_mass_conservation) then

Modified: branches/atmos_physics/src/core_hyd_atmos/mpas_interface.F
===================================================================
--- branches/atmos_physics/src/core_hyd_atmos/mpas_interface.F        2010-05-03 20:35:52 UTC (rev 236)
+++ branches/atmos_physics/src/core_hyd_atmos/mpas_interface.F        2010-05-03 20:44:36 UTC (rev 237)
@@ -2,6 +2,7 @@
 
    use grid_types
    use test_cases
+   use vector_reconstruction
 
    implicit none
 
@@ -28,6 +29,8 @@
    call compute_state_diagnostics(block % time_levs(1) % state, mesh)
    call compute_solve_diagnostics(dt, block % time_levs(1) % state, mesh)
    call initialize_advection_rk(mesh)
+   call init_reconstruct(mesh)
+   call reconstruct(block % time_levs(1) % state, mesh)
 
 end subroutine mpas_init
 

</font>
</pre>