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