<p><b>laura@ucar.edu</b> 2012-02-29 12:47:46 -0700 (Wed, 29 Feb 2012)</p><p>added the initialization of arrays i_rainc and i_rainnc. corrected the initialization of arrays nsteps_accum,ndays_accum,tday_accum,and tyear_accum<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_init.F
===================================================================
--- branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_init.F        2012-02-29 19:45:09 UTC (rev 1554)
+++ branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_init.F        2012-02-29 19:47:46 UTC (rev 1555)
@@ -58,18 +58,29 @@
 !edges:
  call init_dirs_forphys(mesh)
 
+!initialization of counters i_rainc and i_rainnc. i_rainc and i_rainnc track the number of
+!times the accumulated convective (rainc) and grid-scale (rainnc) rain exceed the prescribed
+!threshold value:
+ if(.not. config_do_restart) then
+    do iCell = 1, mesh % nCellsSolve
+       diag_physics % i_rainc  % array(iCell) = 0._RKIND
+       diag_physics % i_rainnc % array(iCell) = 0._RKIND
+    enddo
+ endif
+
 !initialization of temperatures needed for updating the deep soil temperature:
- do iCell = 1, mesh % nCellsSolve
-    diag_physics % nsteps_accum % array(iCell) = 0
-    diag_physics % ndays_accum  % array(iCell) = 1
-
-    diag_physics % tday_accum  % array(iCell) = sfc_input % tmn % array(iCell)
-    diag_physics % tyear_mean  % array(iCell) = sfc_input % tmn % array(iCell)
-    diag_physics % tyear_accum % array(iCell) = sfc_input % tmn % array(iCell)
-    do iLag = 1, mesh % nLags
-       diag_physics % tlag % array(iLag,iCell) = sfc_input % tmn % array(iCell)
+ if(.not. config_do_restart) then
+    do iCell = 1, mesh % nCellsSolve
+       diag_physics % nsteps_accum % array(iCell) = 0._RKIND
+       diag_physics % ndays_accum  % array(iCell) = 0._RKIND
+       diag_physics % tday_accum   % array(iCell) = 0._RKIND
+       diag_physics % tyear_accum  % array(iCell) = 0._RKIND
+       diag_physics % tyear_mean   % array(iCell) = sfc_input % tmn % array(iCell)
+       do iLag = 1, mesh % nLags
+          diag_physics % tlag % array(iLag,iCell) = sfc_input % tmn % array(iCell)
+       enddo
     enddo
- enddo
+ endif
 
 !initialization of global surface properties. set here for now, but may be moved when time
 !manager is implemented:

</font>
</pre>