<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, &
+ physics_idealized_init, &
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>