<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 => domain % blocklist
+ do while (associated(block))
+ call reconstruct(block % time_levs(2) % state, block % mesh)
+ block => 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>