<p><b>laura@ucar.edu</b> 2010-06-14 10:31:18 -0600 (Mon, 14 Jun 2010)</p><p>Added initialization of physics variables<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_hyd_phys/module_physics_init.F
===================================================================
--- branches/atmos_physics/src/core_hyd_phys/module_physics_init.F        2010-06-14 16:29:54 UTC (rev 350)
+++ branches/atmos_physics/src/core_hyd_phys/module_physics_init.F        2010-06-14 16:31:18 UTC (rev 351)
@@ -23,8 +23,47 @@
  TYPE(grid_meta),INTENT(in):: grid
  TYPE(grid_state),INTENT(inout):: vars
 
+!LOCAL VARIABLES:
+ INTEGER:: iCell,nCellsSolve,nVertLevels
+ INTEGER:: i,k
+
 !==============================================================================
 
+ nCellsSolve = grid%nCellsSolve
+ nVertLevels = grid%nVertLevels
+
+!INITIALIZATION OF ARRAYS ON THE GEODESIC GRID:
+!IF(.not. config_do_restart) THEN
+ IF(.not. restart) THEN
+    DO iCell = 1, nCellsSolve
+
+       !cloud microphysics:
+        vars%rainnc%array(i)    = 0.
+        vars%snownc%array(i)    = 0.
+        vars%graupelnc%array(i) = 0.
+        DO k = 1, nVertLevels
+           vars%h_diabatic%array(k,i) = 0.
+        ENDDO
+
+        !convection:
+        vars%nca%array(i)       = 0.
+        vars%cubot%array(i)     = 0.
+        vars%cutop%array(i)     = 0.
+        vars%pratec%array(i)    = 0.
+        vars%rainc%array(i)     = 0.
+        DO k = 1, nVertLevels
+           vars%w0avg%array(k,i)      = 0.
+           vars%rthcuten%array(k,i)   = 0.
+           vars%rqvcuten%array(k,i)   = 0.
+           vars%rqccuten%array(k,i)   = 0.
+           vars%rqrcuten%array(k,i)   = 0.
+           vars%rqicuten%array(k,i)   = 0.
+           vars%rqscuten%array(k,i)   = 0.
+        ENDDO
+
+    ENDDO
+ ENDIF
+
 !INITIALIZATION OF PARAMETERIZED CONVECTIVE PROCESSES:
  CALL init_convection(grid,vars)
 
@@ -71,6 +110,16 @@
  p_qi = index_qi
  p_qs = index_qs
 
+!INITIALIZATION OF PHYSICS ARRAYS:
+ DO j = jts, jte
+ DO i = its, ite
+    cutop_phy(i,j)  = kts
+    cubot_phy(i,j)  = kte
+    raincv_phy(i,j) = 0.0
+    pratec_phy(i,j) = 0.0
+ ENDDO
+ ENDDO
+
 !INITIALIZATION OF KAIN-FRITSCH-ETA CONVECTION SCHEME:
  write(6,*)
  write(6,*) '--- enter subroutine KF_ETA_INIT:'
@@ -95,10 +144,23 @@
  SUBROUTINE init_microphysics
 !==============================================================================
 
+!LOCAL VARIABLES:
+!----------------
+ INTEGER:: i,j
+
+!==============================================================================
+
+!INITIALIZATION OF MICROPHYSICS ARRAYS:
+ DO j = jts,jte
+ DO i = its,ite
+    rainncv_phy(i,j) = 0.0
+ ENDDO
+ ENDDO
+
  CALL thompson_init
 
  END SUBROUTINE init_microphysics
 
 !==============================================================================
  END MODULE module_physics_init
-!==============================================================================
\ No newline at end of file
+!==============================================================================

</font>
</pre>