src/earth/PISMBedDef.hh

Go to the documentation of this file.

// Copyright (C) 2010, 2011 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 __PISMBedDef_hh
#define __PISMBedDef_hh

#include "PISMComponent.hh"
#include "iceModelVec.hh"
#include "deformation.hh"


class PISMBedDef : public PISMComponent_TS {
public:
  PISMBedDef(IceGrid &g, const NCConfigVariable &conf);
  virtual ~PISMBedDef() {}
  virtual PetscErrorCode init(PISMVars &vars);
  virtual PetscErrorCode update(PetscReal t_years, PetscReal dt_years) = 0;
  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*/);
protected:
  PetscErrorCode pismbeddef_allocate(); // packaged to simplify error checking
  PetscErrorCode compute_uplift(PetscScalar dt_beddef);
  PetscReal t_beddef_last;              

  IceModelVec2S topg_initial;
  IceModelVec2S topg_last;      
  IceModelVec2S *thk,           
    *topg,                      
    *uplift;                    
};

class PBPointwiseIsostasy : public PISMBedDef {
public:
  PBPointwiseIsostasy(IceGrid &g, const NCConfigVariable &conf); 
  virtual ~PBPointwiseIsostasy() {}
  virtual PetscErrorCode init(PISMVars &vars);
  virtual PetscErrorCode update(PetscReal t_years, PetscReal dt_years);
protected:
  PetscErrorCode allocate();
  IceModelVec2S thk_last;       
};

#if (PISM_HAVE_FFTW==1)
#include <fftw3.h>

class PBLingleClark : public PISMBedDef {
public:
  PBLingleClark(IceGrid &g, const NCConfigVariable &conf);
  virtual ~PBLingleClark();

  PetscErrorCode init(PISMVars &vars);
  PetscErrorCode update(PetscReal t_years, PetscReal dt_years);
protected:
  PetscErrorCode correct_topg();
  PetscErrorCode allocate();
  PetscErrorCode deallocate();
  PetscErrorCode transfer_to_proc0(IceModelVec2S *source, Vec result);
  PetscErrorCode transfer_from_proc0(Vec source, IceModelVec2S *result);
  Vec g2, g2natural;  
  VecScatter scatter; 
  // Vecs on processor 0:
  Vec Hp0,                      
    bedp0,                      
    Hstartp0,                   
    bedstartp0,                 
    upliftp0;                   
  BedDeformLC bdLC;
};
#endif  // PISM_HAVE_FFTW

#endif  // __PISMBedDef_hh