Program Listing for File NextsimPhysics.hpp

↰ Return to documentation for file (src/include/NextsimPhysics.hpp)

#ifndef SRC_INCLUDE_NEXTSIMPHYSICS_HPP
#define SRC_INCLUDE_NEXTSIMPHYSICS_HPP
#include <memory>

#include "BaseElementData.hpp"
#include "Configured.hpp"
#include "IPhysics1d.hpp"
#include "constants.hpp"

namespace Nextsim {

class PrognosticData;
class PhysicsData;
class ExternalData;
class UnusedData;

class IIceAlbedo;
class IThermodynamics;
class IIceOceanHeatFlux;
class IConcentrationModel;

template <class Phys> class ElementData;

class NextsimPhysics;

class NextsimPhysics : public BaseElementData,
                       public Configured<NextsimPhysics>,
                       public IPhysics1d {
public:
    NextsimPhysics();

    void configure() override;
    enum {
        DRAGOCEANQ_KEY,
        DRAGOCEANT_KEY,
        DRAGICET_KEY,
        OCEANALBEDO_KEY,
        I0_KEY,
        MINC_KEY,
        MINH_KEY,
    };

    void calculate(const PrognosticData&, const ExternalData&, PhysicsData&) override;

    void newIceFormation(const PrognosticData& prog, const ExternalData& exter, PhysicsData& phys);
    inline double newIce() const { return m_newice; };

    inline double QIceAtmosphere() const { return m_Qia; };
    inline double& QIceOceanHeat() { return m_Qio; };

    inline double sublimationRate() const { return m_subl; };

    inline double QDerivativeWRTTemperature() const { return m_dQ_dT; };

    inline double& totalIceFromSnow() { return m_hifroms; }

    static double minimumIceConcentration() { return minc; };
    static double minimumIceThickness() { return minh; };
    static double i0() { return m_I0; };

    class SpecificHumidity {
    public:
        SpecificHumidity();
        double operator()(const double temperature, const double pressure) const;
        double operator()(
            const double temperature, const double pressure, const double salinity) const;

    protected:
        SpecificHumidity(double, double, double, double, double, double, double);
        double f(const double temperature, const double pressurePa) const;
        double est(const double temperature, const double salinity) const;
        const double m_a;
        const double m_b;
        const double m_c;
        const double m_d;
        const double m_bigA;
        const double m_bigB;
        const double m_bigC;
        const double m_alpha;
        const double m_beta;
    };
    class SpecificHumidityIce : public SpecificHumidity {
    public:
        SpecificHumidityIce();
        double operator()(const double temperature, const double pressure) const;
        double dq_dT(const double temperature, const double pressure) const;
    };

protected:
    void updateSpecificHumidityAir(const ExternalData& exter, PhysicsData& phys) override;
    void updateSpecificHumidityWater(
        const PrognosticData& prog, const ExternalData& exter, PhysicsData& phys) override;
    void updateSpecificHumidityIce(
        const PrognosticData& prog, const ExternalData& exter, PhysicsData& phys) override;
    void updateAirDensity(const ExternalData& exter, PhysicsData& phys) override;
    void updateHeatCapacityWetAir(const ExternalData& exter, PhysicsData& phys) override;

private:
    void massFluxOpenWater(PhysicsData& phys);
    void momentumFluxOpenWater(PhysicsData& phys);
    void heatFluxOpenWater(
        const PrognosticData& prog, const ExternalData& exter, PhysicsData& phys);

    void massFluxIceAtmosphere(const PrognosticData& prog, PhysicsData& phys);
    void heatFluxIceAtmosphere(
        const PrognosticData& prog, const ExternalData& exter, PhysicsData& phys);
    void massFluxIceOcean(const PrognosticData& prog, const ExternalData& exter, PhysicsData& phys);
    void heatFluxIceOcean(const PrognosticData& prog, const ExternalData& exter, PhysicsData& phys);
    void lateralGrowth(const PrognosticData& prog, const ExternalData& exter, PhysicsData& phys);

    // Phase change rates
    double m_evap;
    double m_subl;

    // Open water heat fluxes
    double m_Qow;
    double m_Qlwow;
    double m_Qswow;
    double m_Qlhow;
    double m_Qshow;

    // Ice heat fluxes
    double m_Qi;
    double m_Qlwi;
    double m_Qswi;
    double m_Qlhi;
    double m_Qshi;
    double m_dQ_dT;

    // ice-ocean fluxes
    double m_Qio;

    double m_Qia;

    // Thickness of ice generated from flooding of snow [m]
    double m_hifroms;

    static double dragOcean_q;
    static double dragOcean_m(double windSpeed);
    static double dragOcean_t;
    static double dragIce_t;

    static double m_oceanAlbedo;

    // Ice-ocean heat flux
    static IIceOceanHeatFlux* iceOceanHeatFluxImpl;
    // New ice created by cooling below freezing
    double m_newice;

    static IConcentrationModel* iConcentrationModelImpl;

    static double m_I0;

    static double latentHeatWater(double temperature);
    static double latentHeatIce(double temperature);

    static double minc; // minimum ice concentration
    static double minh; // minimum ice true thickness [m]

    static SpecificHumidity specHumWater;
    static SpecificHumidityIce specHumIce;

    static IIceAlbedo* iIceAlbedoImpl;
    static IThermodynamics* iThermo;
};

} /* namespace Nextsim */

#endif /* SRC_INCLUDE_NEXTSIMPHYSICS_HPP */