Skip to content

add an extra erosional base level that prevents ghost nodes to lower … #27

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
EstherHeck wants to merge 2 commits into
base: fastscape_update
Choose a base branch
from
Open

add an extra erosional base level that prevents ghost nodes to lower … #27

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
4569a1f
add an extra erosional base level that prevents ghost nodes to lower …
EstherHeck Mar 15, 2022
0920ed8
indent fastscape cc and h
EstherHeck Mar 15, 2022
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
41 changes: 23 additions & 18 deletions include/aspect/mesh_deformation/fastscape.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -112,12 +112,12 @@ namespace aspect
*/
void parse_parameters (ParameterHandler &prm);

/**
* A function that fills the viscosity derivatives in the
* MaterialModelOutputs object that is handed over, if they exist.
* Does nothing otherwise.
*/
void set_ghost_nodes(double *h, double *vx, double *vy, double *vz, int nx, int ny) const;
/**
* A function that fills the viscosity derivatives in the
* MaterialModelOutputs object that is handed over, if they exist.
* Does nothing otherwise.
*/
void set_ghost_nodes(double *h, double *vx, double *vy, double *vz, int nx, int ny) const;

private:
// Number of FastScape steps per ASPECT timestep.
Expand Down Expand Up @@ -213,7 +213,7 @@ namespace aspect
double kdd;
// Sediment transport coefficient.
double kdsed;

// Orographic parameters
int mmax;
int wb;
Expand All @@ -222,6 +222,11 @@ namespace aspect
double reduc_wb;
bool stackoro;

// Parameters to set an extra erosional base level
// on the ghost nodes that differs from sea level.
bool use_extra_base_level;
double h_extra_base_level;

/**
* Marine parameters
*/
Expand Down Expand Up @@ -259,18 +264,18 @@ namespace aspect
double precision;


/**
* Interval between the generation of graphical output. This parameter
* is read from the input file and consequently is not part of the
* state that needs to be saved and restored.
*/
double output_interval;
/**
* Interval between the generation of graphical output. This parameter
* is read from the input file and consequently is not part of the
* state that needs to be saved and restored.
*/
double output_interval;

/**
* A time (in seconds) at which the last graphical output was supposed
* to be produced. Used to check for the next necessary output time.
*/
mutable double last_output_time;
/**
* A time (in seconds) at which the last graphical output was supposed
* to be produced. Used to check for the next necessary output time.
*/
mutable double last_output_time;
};
}
}
Expand Down
Loading