src/base/energy/tempSystem.hh

Go to the documentation of this file.

// Copyright (C) 2009-2011 Ed Bueler
//
// 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 __tempSystem_hh
#define __tempSystem_hh

#include <petscsys.h>

#include "columnSystem.hh"
#include "pism_const.hh"

class IceModelVec3;


class tempSystemCtx : public columnSystemCtx {

public:
  tempSystemCtx(PetscInt my_Mz, string my_prefix);
  PetscErrorCode initAllColumns();
  PetscErrorCode setSchemeParamsThisColumn(
                     PismMask my_mask, bool my_isMarginal, PetscScalar my_lambda);  
  PetscErrorCode setSurfaceBoundaryValuesThisColumn(PetscScalar my_Ts);
  PetscErrorCode setBasalBoundaryValuesThisColumn(
                     PetscScalar my_G0, PetscScalar my_Tshelfbase, PetscScalar my_Rb);

  PetscErrorCode solveThisColumn(PetscScalar **x, PetscErrorCode &pivoterrorindex);  

public:
  // constants which should be set before calling initForAllColumns()
  PetscScalar  dx,
               dy,
               dtTemp,
               dzEQ,
               ice_rho,
               ice_c_p,
               ice_k;
  // pointers which should be set before calling initForAllColumns()
  PetscScalar  *T,
               *u,
               *v,
               *w,
               *Sigma;
  IceModelVec3 *T3;

protected: // used internally
  PetscInt    Mz;
  PetscScalar lambda, Ts, G0, Tshelfbase, Rb;
  PismMask    mask;
  bool        isMarginal;
  PetscScalar nuEQ,
              rho_c_I,
              iceK,
              iceR;
  bool        initAllDone,
              indicesValid,
              schemeParamsValid,
              surfBCsValid,
              basalBCsValid;
};

#endif  /* __tempSystem_hh */