<p><b>laura@ucar.edu</b> 2011-06-09 10:43:46 -0600 (Thu, 09 Jun 2011)</p><p>corrected z_p needed in kessler parameterization<br>
</p><hr noshade><pre><font color="gray">Modified: branches/atmos_physics/src/core_physics/module_physics_interface_nhyd.F
===================================================================
--- branches/atmos_physics/src/core_physics/module_physics_interface_nhyd.F        2011-06-09 16:39:34 UTC (rev 888)
+++ branches/atmos_physics/src/core_physics/module_physics_interface_nhyd.F        2011-06-09 16:43:46 UTC (rev 889)
@@ -221,6 +221,7 @@
  enddo
  enddo
  enddo
+
  do j = jts,jte
  do i = its,ite
     if(pres_p(i,1,j) .lt. pres_p(i,2,j)) then
@@ -327,7 +328,7 @@
  end subroutine MPAS_to_physics
 
 !=============================================================================================
- subroutine microphysics_from_MPAS(mesh,state,diag)
+ subroutine microphysics_from_MPAS(mesh,state,tend,diag)
 !=============================================================================================
 
 !input variables:
@@ -335,6 +336,8 @@
  type(diag_type) ,intent(in):: diag
  type(mesh_type) ,intent(in):: mesh
 
+ type(tend_type),intent(inout):: tend
+
 !local variables:
  integer:: i,k,j
 
@@ -346,12 +349,13 @@
  real(kind=RKIND),dimension(:,:),pointer:: zgrid
  real(kind=RKIND),dimension(:,:),pointer:: zz,exner,pressure_b,rtheta_p,rtheta_b
  real(kind=RKIND),dimension(:,:),pointer:: rho,theta,qv,rh,pressure_p
+ real(kind=RKIND),dimension(:,:),pointer:: rt_diabatic_tend
 
 !---------------------------------------------------------------------------------------------
 
 !initialization:
  write(0,*)
- write(0,*) '--- enter subroutine microphysics_nhyd_to_phys:'
+ write(0,*) '--- enter subroutine microphysics_from_MPAS:'
 
  zgrid =&gt; mesh % zgrid % array
  zz    =&gt; mesh % zz % array
@@ -366,8 +370,16 @@
  theta      =&gt; state % theta % array
  rh         =&gt; diag % rh % array
 
+ rt_diabatic_tend =&gt; tend % rt_diabatic_tend % array
+
  qv =&gt; state % scalars % array(state%index_qv,:,:)
 
+ do k = kts,kte
+ do i = its,ite
+    rt_diabatic_tend(k,i) = 0.
+ enddo
+ enddo
+
 !copy sounding variables from the geodesic grid to the wrf-physics grid:
  do j = jts, jte
  do k = kts, kte
@@ -378,7 +390,7 @@
     pi_p(i,k,j)   = exner(k,i)
     pres_p(i,k,j) = pressure_b(k,i) + pressure_p(k,i)
 
-    z_p(i,k,j)    = 0.5*(zgrid(k+1,i) + zgrid(k,i))
+    z_p(i,k,j)    = zgrid(k,i)
     dz_p(i,k,j)   = zgrid(k+1,i) - zgrid(k,i)
 
 !   qv_p(i,k,j) = max(0.,state % scalars % array(state%index_qv,k,i))
@@ -390,7 +402,7 @@
  enddo
  enddo
  enddo
- 
+
 !additional initialization as function of cloud microphysics scheme:
  microp_select_init: select case(microp_scheme)
 

</font>
</pre>