<p><b>laura@ucar.edu</b> 2010-10-07 12:55:50 -0600 (Thu, 07 Oct 2010)</p><p>added physics-related calls<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_nhyd_atmos/mpas_interface.F
===================================================================
--- branches/atmos_physics/src/core_nhyd_atmos/mpas_interface.F        2010-10-07 18:52:52 UTC (rev 525)
+++ branches/atmos_physics/src/core_nhyd_atmos/mpas_interface.F        2010-10-07 18:55:50 UTC (rev 526)
@@ -17,6 +17,7 @@
use grid_types
use advection
use time_integration
+ use vector_reconstruction
#ifdef DO_PHYSICS
use module_physics_control
use module_physics_init
@@ -34,6 +35,9 @@
call init_coupled_diagnostics( block % time_levs(1) % state, mesh)
call compute_solve_diagnostics(dt, block % time_levs(1) % state, mesh) ! ok for nonhydrostatic model
call initialize_advection_rk(mesh)
+ call rbfInterp_initialize(mesh)
+ call init_reconstruct(mesh)
+ call reconstruct(block % time_levs(1) % state, mesh)
call initialize_deformation_weights(mesh)
#ifdef DO_PHYSICS
@@ -72,6 +76,7 @@
use grid_types
use time_integration
#ifdef DO_PHYSICS
+ use module_physics_control
use module_physics_driver
use module_physics_manager
#endif
@@ -83,13 +88,13 @@
real (kind=RKIND), intent(in) :: dt
#ifdef DO_PHYSICS
- call physics_timetracker(itimestep)
- if(l_physics) call physics_driver(domain,itimestep)
+ !proceed with physics if moist_physics is set to true:
+ if(moist_physics) then
+ call physics_timetracker(itimestep)
+ if(l_physics) call physics_driver(domain,itimestep)
+ endif
+#endif
call timestep(domain, dt, itimestep)
- write(6,*)
-#else
- call timestep(domain, dt)
-#endif
end subroutine mpas_timestep
</font>
</pre>