tConstantPIK.hh - pism - [fork] customized build of PISM, the parallel ice sheet model (tillflux branch)
  HTML git clone git://src.adamsgaard.dk/pism
   DIR Log
   DIR Files
   DIR Refs
   DIR LICENSE
       ---
       tConstantPIK.hh (2044B)
       ---
            1 // Copyright (C) 2011, 2012, 2013, 2014, 2015, 2016, 2017, 2018, 2019 PISM Authors
            2 //
            3 // This file is part of PISM.
            4 //
            5 // PISM is free software; you can redistribute it and/or modify it under the
            6 // terms of the GNU General Public License as published by the Free Software
            7 // Foundation; either version 3 of the License, or (at your option) any later
            8 // version.
            9 //
           10 // PISM is distributed in the hope that it will be useful, but WITHOUT ANY
           11 // WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
           12 // FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more
           13 // details.
           14 //
           15 // You should have received a copy of the GNU General Public License
           16 // along with PISM; if not, write to the Free Software
           17 // Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301  USA
           18 
           19 #ifndef _POCONSTANTPIK_H_
           20 #define _POCONSTANTPIK_H_
           21 
           22 #include "CompleteOceanModel.hh"
           23 
           24 namespace pism {
           25 namespace ocean {
           26 //! \brief Implements the ocean model used in [@ref Martinetal2011].
           27 //!
           28 //! Uses a parameterization of sub-shelf melting with respect to
           29 //! sub-shelf heat flux like in [@ref BeckmannGoosse2003].
           30 //!
           31 //! Models heat flux into the base of the shelf as
           32 //!
           33 //! @f[ Q_{\text{heat}} = \rho_{o} c_{p_{o}} \gamma_{T} (T_{o} - T_{f}), @f]
           34 //!
           35 //! where @f$\rho_{o}@f$ is the density of ocean water, @f$c_{p_{o}}@f$ and
           36 //! @f$T_{o}@f$ are the heat capacity and temperature of the ocean mixed
           37 //! layer, @f$T_{f}@f$ is the freezing temperature of ocean water at the
           38 //! shelf bottom.
           39 class PIK : public CompleteOceanModel {
           40 public:
           41   PIK(IceGrid::ConstPtr g);
           42   virtual ~PIK();
           43 
           44 private:
           45   MaxTimestep max_timestep_impl(double t) const;
           46   void update_impl(const Geometry &geometry, double my_t, double my_dt);
           47   void init_impl(const Geometry &geometry);
           48 
           49   void melting_point_temperature(const IceModelVec2S &depth, IceModelVec2S &result) const;
           50   void mass_flux(const IceModelVec2S &ice_thickness, IceModelVec2S &result) const;
           51 };
           52 
           53 } // end of namespace ocean
           54 } // end of namespace pism
           55 #endif /* _POCONSTANTPIK_H_ */