<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 => 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 => block % next
end do
</font>
</pre>