<p><b>duda</b> 2010-06-04 14:08:23 -0600 (Fri, 04 Jun 2010)</p><p>Begin using the new RBF interpolation and vector reconstruction<br>
modules in the hydrostatic atmosphere core.<br>
<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: trunk/mpas/src/core_hyd_atmos/module_time_integration.F
===================================================================
--- trunk/mpas/src/core_hyd_atmos/module_time_integration.F        2010-06-04 17:41:39 UTC (rev 335)
+++ trunk/mpas/src/core_hyd_atmos/module_time_integration.F        2010-06-04 20:08:23 UTC (rev 336)
@@ -4,6 +4,7 @@
use configure
use constants
use dmpar
+ use vector_reconstruction
contains
@@ -314,11 +315,13 @@
! END RK loop
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
- ! compute vertical velocity diagnostic
-
+ !
+ ! Compute full velocity vectors at cell centers, and compute vertical velocity diagnostic
+ !
block => domain % blocklist
do while (associated(block))
- call compute_w( block % time_levs(2) % state, block % time_levs(1) % state, block % mesh, dt )
+ call reconstruct(block % time_levs(2) % state, block % mesh)
+ call compute_w(block % time_levs(2) % state, block % time_levs(1) % state, block % mesh, dt)
block => block % next
end do
Modified: trunk/mpas/src/core_hyd_atmos/mpas_interface.F
===================================================================
--- trunk/mpas/src/core_hyd_atmos/mpas_interface.F        2010-06-04 17:41:39 UTC (rev 335)
+++ trunk/mpas/src/core_hyd_atmos/mpas_interface.F        2010-06-04 20:08:23 UTC (rev 336)
@@ -17,6 +17,8 @@
use grid_types
use advection
use time_integration
+ use RBF_interpolation
+ use vector_reconstruction
implicit none
@@ -28,6 +30,9 @@
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 rbfInterp_initialize(mesh)
+ call init_reconstruct(mesh)
+ call reconstruct(block % time_levs(1) % state, mesh)
end subroutine mpas_init
</font>
</pre>