<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>