<p><b>laura@ucar.edu</b> 2011-06-09 10:21:12 -0600 (Thu, 09 Jun 2011)</p><p>Added initialization of input physics variables to default values when running MPAS with idealized cases, in particular on the IBM.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_nhyd_atmos/Registry
===================================================================
--- branches/atmos_physics/src/core_nhyd_atmos/Registry        2011-06-08 18:47:07 UTC (rev 884)
+++ branches/atmos_physics/src/core_nhyd_atmos/Registry        2011-06-09 16:21:12 UTC (rev 885)
@@ -56,7 +56,6 @@
 dim R3 3
 dim nVertLevels namelist:config_nvertlevels
 dim nVertLevelsP1 nVertLevels+1
-dim nMonths 12
 
 #
 # var  type  name_in_file  ( dims )  iro-  name_in_code super-array array_class
@@ -273,7 +272,10 @@
 namelist character physics  input_landuse_data          USGS
 namelist character physics  input_soil_data             STAS
 namelist integer   physics  num_soil_layers             4
+namelist integer   physics  num_months                  12
 
+dim nMonths namelist:num_months
+
 #... DIMENSION NEEDED FOR LONGWAVE AND SHORTWAVE RADIATION FLUXES TO INCLUDE AN ADDITIONAL LAYER
 #... BETWEEN THE TOP OF THE MODEL AND THE TOP OF THE ATMOSPHERE
 dim nVertLevelsP2 nVertLevels+2

Modified: branches/atmos_physics/src/core_nhyd_atmos/module_test_cases.F
===================================================================
--- branches/atmos_physics/src/core_nhyd_atmos/module_test_cases.F        2011-06-08 18:47:07 UTC (rev 884)
+++ branches/atmos_physics/src/core_nhyd_atmos/module_test_cases.F        2011-06-09 16:21:12 UTC (rev 885)
@@ -5,6 +5,7 @@
    use constants
    use dmpar
    use advection
+   use module_physics_control
 
 
    contains
@@ -93,6 +94,18 @@
          stop
       end if
 
+      !initialization of surface input variables technically not needed to run our current set of
+      !idealized test cases:
+      if (config_test_case &gt; 0)  then
+
+         block_ptr =&gt; domain % blocklist
+         do while (associated(block_ptr))
+            call physics_idealized_init(block_ptr % mesh, block_ptr % sfc_input)
+            block_ptr =&gt; block_ptr % next
+         end do
+
+      endif
+
    end subroutine setup_nhyd_test_case
 
 !----------------------------------------------------------------------------------------------------------
@@ -1287,6 +1300,7 @@
             rb (k,i) = pb(k,i)**(1./rcv)/((rgas/p0)*tb(k,i)*zz(k,i))
             rtb(k,i) = rb(k,i)*tb(k,i)
             rr (k,i) = p (k,i)**(1./rcv)/((rgas/p0)*t (k,i)*zz(k,i))-rb(k,i)
+            ppb(k,i) = p0*(zz(k,i)*rgas*rtb(k,i)/p0)**(cp/cv)
          end do
       end do
 

</font>
</pre>