https://github.com/geodynamics/citcoms
Revision df821c06952128812763ec5b7a944cd178ca5ed3 authored by Thorsten Becker on 08 December 2010, 18:43:58 UTC, committed by Thorsten Becker on 08 December 2010, 18:43:58 UTC
GMT/ggrd compile before).
1 parent cf7495e
Tip revision: df821c06952128812763ec5b7a944cd178ca5ed3 authored by Thorsten Becker on 08 December 2010, 18:43:58 UTC
Added definition of TRUE (1) and FALSE (0) in case undefined (was defined for
Added definition of TRUE (1) and FALSE (0) in case undefined (was defined for
Tip revision: df821c0
viscosity_descriptions.h
/*
*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
*
*<LicenseText>
*
* CitcomS by Louis Moresi, Shijie Zhong, Lijie Han, Eh Tan,
* Clint Conrad, Michael Gurnis, and Eun-seo Choi.
* Copyright (C) 1994-2005, California Institute of Technology.
*
* This program 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.
*
* This program 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 this program; if not, write to the Free Software
* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
*
*</LicenseText>
*
*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
*/
#define CITCOM_MAX_VISC_LAYER 40
struct VISC_OPT {
int update_allowed; /* determines whether visc field can evolve */
int SMOOTH;
int smooth_cycles;
int allow_anisotropic_viscosity,anisotropic_viscosity_init;
#ifdef CITCOM_ALLOW_ANISOTROPIC_VISC
int anivisc_start_from_iso; /* start from isotropic solution? */
int anisotropic_init; /* 0: isotropic, 1: random, 2: from file */
char anisotropic_init_dir[1000];
int anivisc_layer; /* layer to assign anisotropic viscosity to for mode 2 */
#endif
char STRUCTURE[20]; /* which option to determine viscosity field, one of .... */
int FROM_SYSTEM;
int FROM_FILE;
int FROM_SPECS;
/* System ... */
int RHEOL; /* 1,2 */
int num_mat;
float zcmb; /* old layer specs */
float zlm;
float z410;
float zlith;
float zbase_layer[CITCOM_MAX_VISC_LAYER]; /* new */
/* low viscosity channel and wedge stuff */
int channel;
int wedge;
float lv_min_radius;
float lv_max_radius;
float lv_channel_thickness;
float lv_reduction;
/* viscosity cut-off */
int MAX;
float max_value;
int MIN;
float min_value;
/* non-Newtonian stress dependence */
int SDEPV;
float sdepv_misfit;
int sdepv_normalize, sdepv_visited;
float *sdepv_expt;
/* compositional viscosity */
int CDEPV;
float cdepv_ff[10]; /* flavor factors */
/* "plasticity" law parameters */
int PDEPV;
float *pdepv_a, *pdepv_b, *pdepv_y;
float pdepv_offset;
int pdepv_eff, pdepv_visited;
int psrw;
/* temperature dependence */
int TDEPV;
float *N0, *E, *T, *Z;
float ET_red, T_sol0; /* for viscosity law 8 */
int layer_control;
char layer_file[255];
} viscosity;
Computing file changes ...