<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 =&gt; 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 =&gt; 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>