!--------------------------------------------------------------------------
   !---------------------------- compute_pmid ----------------------------------
   IF(ap_bp_present) THEN
      IF(offset>0) THEN
         !$OMP DO SCHEDULE(STATIC)
         DO ij = 1, primal_num
            !DIR$ SIMD
            DO l = 1, llm
               pmid(l,ij) = .5*(AP(l,ij)+AP(l+1,ij) + (BP(l,ij)+BP(l+1,ij)) * ps(ij) )
            END DO
         END DO
         !$OMP END DO
      ELSE
         !$OMP DO SCHEDULE(STATIC)
         DO ij = 1, primal_num
            !DIR$ SIMD
            DO l = 1, llm
               pmid(l,ij) = .5*(AP(l,ij)+AP(l+1,ij) + (BP(l,ij)+BP(l+1,ij)) * ps(ij) )
            END DO
         END DO
         !$OMP END DO
      END IF
   END IF
   !---------------------------- compute_pmid ----------------------------------
   !--------------------------------------------------------------------------