<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 > 0) then
+
+ block_ptr => domain % blocklist
+ do while (associated(block_ptr))
+ call physics_idealized_init(block_ptr % mesh, block_ptr % sfc_input)
+ block_ptr => 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>