!--------------------------------------------------------------------------
   !---------------------------- compute_pression ----------------------------------
   IF(ap_bp_present) THEN
      IF(offset>0) THEN
         !$OMP DO SCHEDULE(STATIC)
         DO ij = 1, primal_num
            !DIR$ SIMD
            DO l = 1, llm+1
               p(l,ij) = AP(l,ij) + BP(l,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+1
               p(l,ij) = AP(l,ij) + BP(l,ij) * ps(ij)
            END DO
         END DO
         !$OMP END DO
      END IF
   END IF
   !---------------------------- compute_pression ----------------------------------
   !--------------------------------------------------------------------------