!--------------------------------------------------------------------------
   !---------------------------- compute_pression ----------------------------------
   IF(ap_bp_present) THEN
      IF(offset>0) THEN
         DO l = ll_begin, ll_endp1
            !DIR$ SIMD
            DO ij=ij_begin_ext, ij_end_ext
               p(ij,l) = AP(ij,l) + BP(ij,l) * ps(ij)
            END DO
         END DO
      ELSE
         DO l = ll_begin, ll_endp1
            !DIR$ SIMD
            DO ij=ij_begin, ij_end
               p(ij,l) = AP(ij,l) + BP(ij,l) * ps(ij)
            END DO
         END DO
      END IF
   END IF
   !---------------------------- compute_pression ----------------------------------
   !--------------------------------------------------------------------------