<p><b>laura@ucar.edu</b> 2010-07-23 12:08:25 -0600 (Fri, 23 Jul 2010)</p><p>added changes made to mpas_interface.F for non-hydrostatic core for consistency between the two dynamical cores<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_hyd_atmos/mpas_interface.F
===================================================================
--- branches/atmos_physics/src/core_hyd_atmos/mpas_interface.F        2010-07-23 18:03:08 UTC (rev 391)
+++ branches/atmos_physics/src/core_hyd_atmos/mpas_interface.F        2010-07-23 18:08:25 UTC (rev 392)
@@ -20,6 +20,7 @@
    use RBF_interpolation
    use vector_reconstruction
 #ifdef DO_PHYSICS
+   use module_physics_control
    use module_physics_init
    use module_physics_manager
 #endif
@@ -39,8 +40,19 @@
    call reconstruct(block % time_levs(1) % state, mesh)
 
 #ifdef DO_PHYSICS
-   call physics_wrf_allocate(mesh)
-   call physics_init(mesh, block % time_levs(1) % state)
+   !check that all the physics options are correctly defined and that at least one physics
+   !parameterization is called (using the logical moist_physics):
+   call physics_namelist_check
+
+   !proceed with initialization of physics parameterization if moist_physics is set to true:
+   if(moist_physics) then
+      call physics_wrf_interface(mesh)
+      call physics_init(mesh, block % time_levs(1) % state)
+   endif
+   write(6,*)
+   write(6,*) '--- end physics_init:'
+   stop
+
 #endif   
 
 end subroutine mpas_init
@@ -85,10 +97,10 @@
 #endif
 
 #ifdef DO_PHYSICS
-   l_physics = .false.
-   call physics_timetracker(itimestep,l_physics)
-   if(l_physics) call physics_driver(domain,itimestep)
-   call timestep(domain, dt, itimestep)
+!  l_physics = .false.
+!  call physics_timetracker(itimestep,l_physics)
+!  if(l_physics) call physics_driver(domain,itimestep)
+   call timestep(domain, dt, itimestep, config_ntimesteps)
 #else
    call timestep(domain, dt)
 #endif

</font>
</pre>