initial_lithostatic_pressure.h

Go to the documentation of this file.

// ------------------------------------------------------
//
// SPDX-License-Identifier: GPL-2.0-or-later
//
// SPDX-FileCopyrightText: Copyright (C) 2011-2026 by the ASPECT authors.
//
// This file is part of ASPECT.
//
// Detailed license information governing the source code and
// contributions can be found in the folder LICENSES and in
// LICENSE and CONTRIBUTING.md at the top level directory.
//
// ------------------------------------------------------


#ifndef _aspect_boundary_traction_initial_lithospheric_pressure_h
#define _aspect_boundary_traction_initial_lithospheric_pressure_h

#include <aspect/boundary_traction/interface.h>
#include <aspect/simulator_access.h>


namespace aspect
{
  namespace BoundaryTraction
  {
    template <int dim>
    class InitialLithostaticPressure : public Interface<dim>, public SimulatorAccess<dim>
    {
      public:

        void initialize () override;


        Tensor<1,dim>
        boundary_traction (const types::boundary_id boundary_indicator,
                           const Point<dim> &position,
                           const Tensor<1,dim> &normal_vector) const override;


        static
        void
        declare_parameters (ParameterHandler &prm);


        void
        parse_parameters (ParameterHandler &prm) override;

      private:

        unsigned int n_points;

        double delta_z;

        /*
         * The user-specified point where to calculate the pressure profile.
         * The vertical coordinate/radius is ignored.
         */
        Point<dim> representative_point;

        std::vector<double> pressure;

        double interpolate_pressure (const Point<dim> &p) const;

        types::boundary_id bottom_boundary_id;

        bool prescribe_constant_pressure_at_bottom_boundary;
    };
  }
}


#endif