<p><b>laura@ucar.edu</b> 2012-11-27 10:07:19 -0700 (Tue, 27 Nov 2012)</p><p>corrected to the calculation of precipitable water. commented out some print statements.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_driver_microphysics.F
===================================================================
--- branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_driver_microphysics.F        2012-11-26 22:42:05 UTC (rev 2321)
+++ branches/atmos_physics/src/core_atmos_physics/mpas_atmphys_driver_microphysics.F        2012-11-27 17:07:19 UTC (rev 2322)
@@ -386,6 +386,7 @@
 
 !local variables:
  integer:: i,j,k
+ real(kind=RKIND):: rho_a
 
 !---------------------------------------------------------------------------------------------
 
@@ -397,8 +398,9 @@
     !precipitable water:
     diag_physics % precipw % array(i) = 0._RKIND
     do k = kts,kte
-       diag_physics % precipw % array(i) = diag_physics % precipw % array(i) &amp;
-                                         + qv_p(i,k,j) *  dz_p(i,k,j)
+       rho_a = rho_p(i,k,j) / (1._RKIND + qv_p(i,k,j))
+       diag_physics % precipw % array(i) = &amp;
+            diag_physics % precipw % array(i) + qv_p(i,k,j) * rho_a * dz_p(i,k,j)
     enddo
 
     !time-step precipitation:
@@ -459,8 +461,8 @@
  real(kind=RKIND),dimension(:),allocatable:: qv1d,qr1d,qs1d,qg1d,t1d,p1d,dBZ1d
 
 !---------------------------------------------------------------------------------------------
- write(0,*)
- write(0,*) '--- enter subroutine COMPUTE_RADAR_REFLECTIVITY:'
+!write(0,*)
+!write(0,*) '--- enter subroutine COMPUTE_RADAR_REFLECTIVITY:'
 
  microp_select: select case(microp_scheme)
 
@@ -501,8 +503,8 @@
 !            write(0,201) i,k,dBZ1d(k)
           enddo
           diag_physics % refl10cm_max % array(i) = maxval(dBZ1d(:))
-          if(diag_physics % refl10cm_max % array(i) .gt. 0.) &amp;
-             write(0,201) j,i,diag_physics % refl10cm_max % array(i)
+!         if(diag_physics % refl10cm_max % array(i) .gt. 0.) &amp;
+!            write(0,201) j,i,diag_physics % refl10cm_max % array(i)
        enddo
        enddo
 
@@ -517,7 +519,7 @@
     case default
 
  end select microp_select
- write(0,*) '--- end subroutine COMPUTE_RADAR_REFLECTIVITY'
+!write(0,*) '--- end subroutine COMPUTE_RADAR_REFLECTIVITY'
 
  201 format(2i6,e15.8)
 

</font>
</pre>