<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 => mesh % zgrid % array
zz => mesh % zz % array
@@ -366,8 +370,16 @@
theta => state % theta % array
rh => diag % rh % array
+ rt_diabatic_tend => tend % rt_diabatic_tend % array
+
qv => 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>