<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) &
- + 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) = &
+ 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.) &
- write(0,201) j,i,diag_physics % refl10cm_max % array(i)
+! if(diag_physics % refl10cm_max % array(i) .gt. 0.) &
+! 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>