<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>