!--------------------------------------------------------------------------
   !---------------------------- compute_pmid ----------------------------------
   IF(ap_bp_present) THEN
      IF(offset>0) THEN
         DO l = ll_begin, ll_end
            !DIR$ SIMD
            DO ij=ij_begin_ext, ij_end_ext
               pmid(ij,l) = .5*(AP(ij,l)+AP(ij,l+1) + (BP(ij,l)+BP(ij,l+1)) * ps(ij) )
            END DO
         END DO
      ELSE
         DO l = ll_begin, ll_end
            !DIR$ SIMD
            DO ij=ij_begin, ij_end
               pmid(ij,l) = .5*(AP(ij,l)+AP(ij,l+1) + (BP(ij,l)+BP(ij,l+1)) * ps(ij) )
            END DO
         END DO
      END IF
   END IF
   !---------------------------- compute_pmid ----------------------------------
   !--------------------------------------------------------------------------