--- trunk/libf/dyn3d/exner_hyb.f90	2012/08/29 14:47:17	64
+++ trunk/dyn3d/exner_hyb.f	2014/03/06 15:12:00	83
@@ -30,7 +30,7 @@
 
     use dimens_m, only: iim, jjm, llm
     use comconst, only: kappa, cpp
-    use comvert, only: preff
+    use disvert_m, only: preff
     use comgeom, only: aire_2d, apoln, apols
     use filtreg_m, only: filtreg
 
@@ -59,11 +59,10 @@
             / (p(:, :, l) * unpl2k + p(:, :, l+1) * (beta(:, :, l+1) - unpl2k))
     ENDDO
 
-    pk(:, :, 1) = p(:, :, 1) * pks  &
-         / (p(:, :, 1) * (1. + kappa) &
-         + 0.5 * (beta(:, :, 2) - unpl2k) * p(:, :, 2))
+    pk(:, :, 1) = ps * pks &
+         / (ps * (1. + kappa) + 0.5 * (beta(:, :, 2) - unpl2k) * p(:, :, 2))
     DO l = 2, llm
-       pk(:, :, l) = beta(:, :, l) * pk(:, :, l-1)
+       pk(:, :, l) = beta(:, :, l) * pk(:, :, l - 1)
     ENDDO
 
     if (present(pkf)) then