<p><b>laura@ucar.edu</b> 2010-12-21 15:28:42 -0700 (Tue, 21 Dec 2010)</p><p>added dminfo in calls to subroutine mpas_init and physics_init for use in the initialization of short- and long-wave radiation codes<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_nhyd_atmos/module_core.F
===================================================================
--- branches/atmos_physics/src/core_nhyd_atmos/module_core.F        2010-12-10 21:21:09 UTC (rev 646)
+++ branches/atmos_physics/src/core_nhyd_atmos/module_core.F        2010-12-21 22:28:42 UTC (rev 647)
@@ -31,7 +31,7 @@
dt = config_dt
block => domain % blocklist
do while (associated(block))
- call mpas_init_block(block, block % mesh, dt)
+ call mpas_init_block(domain % dminfo, block, block % mesh, dt)
block => block % next
end do
@@ -40,7 +40,7 @@
end subroutine mpas_init
- subroutine mpas_init_block(block, mesh, dt)
+ subroutine mpas_init_block(dminfo, block, mesh, dt)
use grid_types
! use advection
@@ -49,6 +49,7 @@
use RBF_interpolation
use vector_reconstruction
#ifdef DO_PHYSICS
+ use module_physics_aquaplanet
use module_physics_control
use module_physics_init
use module_physics_manager
@@ -56,6 +57,7 @@
implicit none
+ type (dm_info), intent(in):: dminfo
type (block_type), intent(inout) :: block
type (mesh_type), intent(inout) :: mesh
real (kind=RKIND), intent(in) :: dt
@@ -86,7 +88,10 @@
!initialization of all physics variables in registry:
call physics_registry_init(config_do_restart, mesh, block % diag_physics, block % tend_physics)
call physics_wrf_interface(mesh)
- call physics_init(mesh, block % state % time_levs(1) % state, block % diag_physics)
+ call physics_init(dminfo, mesh, block % state % time_levs(1) % state, block % diag_physics)
+
+!initialization of some surface variables (temporary):
+ call physics_aquaplanet_init(mesh, block % diag_physics)
endif
#endif
</font>
</pre>