src/base/stressbalance/PISMBedSmoother.hh

Go to the documentation of this file.

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

#include <petsc.h>
#include "NCVariable.hh"
#include "grid.hh"
#include "iceModelVec.hh"


class PISMBedSmoother {
public:
  PISMBedSmoother(IceGrid &g, const NCConfigVariable &conf, PetscInt MAX_GHOSTS);
  virtual ~PISMBedSmoother();

  virtual PetscErrorCode preprocess_bed(IceModelVec2S topg, PetscReal n, 
                                        PetscReal lambda);
  virtual PetscErrorCode preprocess_bed(IceModelVec2S topg, PetscReal n,
                                        PetscInt Nx_in, PetscInt Ny_in);

  virtual PetscErrorCode get_smoothing_domain(PetscInt &Nx_out, PetscInt &Ny_out);
  virtual PetscInt       get_max_ghosts() { return maxGHOSTS; }

  virtual PetscErrorCode get_smoothed_thk(IceModelVec2S usurf, IceModelVec2S thk, IceModelVec2Int mask,
                                          PetscInt GHOSTS, IceModelVec2S *thksmooth);
  virtual PetscErrorCode get_theta(IceModelVec2S usurf, PetscReal n,
                                   PetscInt GHOSTS, IceModelVec2S *theta);

  IceModelVec2S topgsmooth;  

protected:
  IceGrid &grid;
  const NCConfigVariable &config;
  IceModelVec2S maxtl,C2,C3,C4;

  PetscInt Nx,Ny;  


  PetscInt maxGHOSTS; 



  PetscErrorCode allocate();
  PetscErrorCode deallocate();

  Vec g2, g2natural;  
  VecScatter scatter; 
  Vec topgp0,         
      topgsmoothp0,   
      maxtlp0,        
      C2p0,C3p0,C4p0;

  PetscErrorCode smooth_the_bed_on_proc0();
  PetscErrorCode compute_coefficients_on_proc0(PetscReal n);
};

#endif  // __PISMBedSmoother_hh