src/base/energy/bedrockThermalUnit.hh

Go to the documentation of this file.

// Copyright (C) 2011 Ed Bueler and Constantine Khroulev
//
// This file is part of PISM.
//
// PISM is free software; you can redistribute it and/or modify it under the
// terms of the GNU General Public License as published by the Free Software
// Foundation; either version 2 of the License, or (at your option) any later
// version.
//
// PISM is distributed in the hope that it will be useful, but WITHOUT ANY
// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
// FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more
// details.
//
// You should have received a copy of the GNU General Public License
// along with PISM; if not, write to the Free Software
// Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301  USA

#ifndef _PISMBEDTHERMALUNIT_H_
#define _PISMBEDTHERMALUNIT_H_

#include "PISMComponent.hh"
#include "iceModelVec.hh"
#include "PISMVars.hh"
#include "basal_resistance.hh"

class IceModelVec3BTU : public IceModelVec3D {
public:
  IceModelVec3BTU() { Lbz = -1.0; }
  virtual ~IceModelVec3BTU() {}

  virtual PetscErrorCode create(IceGrid &mygrid, const char my_short_name[], bool local,
                                int myMbz, PetscReal myLbz, int stencil_width = 1);
                                
  virtual PetscErrorCode get_layer_depth(PetscReal &depth); 
  virtual PetscErrorCode get_spacing(PetscReal &dzb);
private:
  PetscReal Lbz;
  bool good_init();
};



class PISMBedThermalUnit : public PISMComponent_TS {

public:
  PISMBedThermalUnit(IceGrid &g, const NCConfigVariable &conf);

  virtual ~PISMBedThermalUnit() { }

  virtual PetscErrorCode init(PISMVars &vars);

  virtual void add_vars_to_output(string keyword, set<string> &result);
  virtual PetscErrorCode define_variables(set<string> vars, const NCTool &nc, nc_type nctype);  
  virtual PetscErrorCode write_variables(set<string> vars, string filename);

  virtual PetscErrorCode max_timestep(PetscReal /*t_years*/, PetscReal &dt_years);

  virtual PetscErrorCode update(PetscReal t_years, PetscReal dt_years);

  virtual PetscErrorCode get_upward_geothermal_flux(IceModelVec2S &result);
protected:
  virtual PetscErrorCode allocate(int Mbz, double Lbz);

  virtual PetscErrorCode bootstrap();
  virtual PetscErrorCode regrid();

  IceModelVec3BTU  temp;     



  // parameters of the heat equation:  T_t = D T_xx  where D = k / (rho c)
  PetscScalar    bed_rho, bed_c, bed_k, bed_D;
  
  PetscInt Mbz;
  PetscReal Lbz;

  IceModelVec2S *bedtoptemp, 
                *ghf; 
};

#endif /* _PISMBEDTHERMALUNIT_H_ */