<p><b>laura@ucar.edu</b> 2011-06-09 10:39:34 -0600 (Thu, 09 Jun 2011)</p><p>added subroutine to initialized input physics variables to default values when running idealized cases<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_physics/module_physics_control.F
===================================================================
--- branches/atmos_physics/src/core_physics/module_physics_control.F        2011-06-09 16:32:33 UTC (rev 887)
+++ branches/atmos_physics/src/core_physics/module_physics_control.F        2011-06-09 16:39:34 UTC (rev 888)
@@ -22,6 +22,7 @@
  implicit none
  private
  public:: physics_namelist_check, &amp;
+          physics_idealized_init, &amp;
           physics_registry_init
 
  logical,public:: moist_physics
@@ -202,6 +203,62 @@
  end subroutine physics_registry_init
 
 !=============================================================================================
+ subroutine physics_idealized_init(mesh,sfc_input)
+!=============================================================================================
+
+!input and output arguments:
+ type(mesh_type),intent(in):: mesh
+ type(sfc_input_type),intent(inout):: sfc_input
+
+!local variables:
+ integer:: iCell,iMonth,iSoil
+ 
+!---------------------------------------------------------------------------------------------
+
+!initialization of surface input variables that are not needed if we run the current set of
+!idealized test cases:
+
+ do iCell = 1, mesh % nCells
+    !terrain,soil type, and vegetation:
+    sfc_input % ter      % array(iCell) = 0.
+    sfc_input % xice     % array(iCell) = 0.
+    sfc_input % landmask % array(iCell) = 0
+    sfc_input % ivgtyp   % array(iCell) = 0
+    sfc_input % isltyp   % array(iCell) = 0
+    sfc_input % shdmin   % array(iCell) = 0.
+    sfc_input % shdmax   % array(iCell) = 0.
+    sfc_input % vegfra   % array(iCell) = 0.
+
+    !snow coverage:    
+    sfc_input % snow     % array(iCell) = 0.
+    sfc_input % snowc    % array(iCell) = 0.
+    sfc_input % snoalb   % array(iCell) = 0.08 
+    
+    !surface and sea-surface temperatures:
+    sfc_input % skintemp % array(iCell) = 288.0
+    sfc_input % sst      % array(iCell) = 288.0
+    
+    !soil layers:
+    sfc_input % tmn % array(iCell) = 288.0
+    do iSoil = 1, mesh % nSoilLevels
+       sfc_input % tslb % array(iSoil,iCell)   = 288.0
+       sfc_input % smcrel % array(iSoil,iCell) =   0.0
+       sfc_input % sh2o   % array(iSoil,iCell) =   0.0
+       sfc_input % smois  % array(iSoil,iCell) =   0.0
+       sfc_input % dzs    % array(iSoil,iCell) =   0.0
+    enddo
+
+    !monthly climatological surface albedo and greeness fraction:
+    do iMonth = 1, mesh % nMonths
+       sfc_input % albedo12m % array(iMonth,iCell) = 0.08
+       sfc_input % greenfrac % array(iMonth,iCell) = 0.
+    enddo
+
+ enddo   
+
+ end subroutine physics_idealized_init
+
+!=============================================================================================
  end module module_physics_control
 !=============================================================================================
 

</font>
</pre>