<p><b>duda</b> 2010-11-30 11:36:00 -0700 (Tue, 30 Nov 2010)</p><p>BRANCH COMMIT<br>
<br>
Add calls in non-hydrostatic atmosphere core to reconstruct wind vectors at <br>
cell centers, and to compute zonal and meridional winds at cell centers, too.<br>
<br>
<br>
M    src/core_nhyd_atmos/module_mpas_core.F<br>
M    src/core_nhyd_atmos/module_time_integration.F<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_nonhydrostatic/src/core_nhyd_atmos/module_mpas_core.F
===================================================================
--- branches/atmos_nonhydrostatic/src/core_nhyd_atmos/module_mpas_core.F        2010-11-29 20:28:21 UTC (rev 634)
+++ branches/atmos_nonhydrostatic/src/core_nhyd_atmos/module_mpas_core.F        2010-11-30 18:36:00 UTC (rev 635)
@@ -46,6 +46,7 @@
    !   use advection
       use time_integration
       use configure
+      use vector_reconstruction
    
       implicit none
    
@@ -58,6 +59,9 @@
          call init_coupled_diagnostics( block % state % time_levs(1) % state, block % diag, mesh)
          call compute_solve_diagnostics(dt, block % state % time_levs(1) % state, block % diag, mesh)
       end if
+      call rbfInterp_initialize(mesh)
+      call init_reconstruct(mesh)
+      call reconstruct(block % state % time_levs(1) % state, block % diag, mesh)
    
    !
    ! Note: The following initialization calls have been moved to mpas_setup_test_case()

Modified: branches/atmos_nonhydrostatic/src/core_nhyd_atmos/module_time_integration.F
===================================================================
--- branches/atmos_nonhydrostatic/src/core_nhyd_atmos/module_time_integration.F        2010-11-29 20:28:21 UTC (rev 634)
+++ branches/atmos_nonhydrostatic/src/core_nhyd_atmos/module_time_integration.F        2010-11-30 18:36:00 UTC (rev 635)
@@ -4,6 +4,7 @@
    use configure
    use constants
    use dmpar
+   use vector_reconstruction
 
 
    contains
@@ -380,6 +381,7 @@
       if(do_microphysics) then
       block =&gt; domain % blocklist
         do while (associated(block))
+           call reconstruct(block % state % time_levs(2) % state, block % diag, block % mesh)
            call qd_kessler( block % state % time_levs(1) % state, block % state % time_levs(2) % state, block % diag, block % tend, block % mesh, dt )
            block =&gt; block % next
         end do

</font>
</pre>