<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 =&gt; 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 =&gt; 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>