https://github.com/geodynamics/citcoms
Revision 45be8f2e32748cd56c3154bbdcd61da3fddda4bd authored by Eh Tan on 09 March 2007, 00:25:03 UTC, committed by Eh Tan on 09 March 2007, 00:25:03 UTC
* Seperating composition codes from tracer codes, several struct members in E->trace are moved to E->composition
* Add chemical buoyancy term
* Renamed thermal_buoyancy() to get_buoyancy() to reflect the fact that this function computes buoyancy according to temperature/composition/phase
* Fix and optimize the fix_*() functions
* Remove some static variables, or converte them to const
* Change the type of E->SinCos from float to double, and replace E->trace.DSinCos
* Provide default values for "required" input parameters
* Renamed tracer_restart to tracer_ic_method, also remapped the options
* Removed option for Cartesian tracer input and trace.iwrite_tracers_every
* Replaced trace.itracer_type by composition.ichemical_buoyancy
* Optionally output tracers and composition fields

1 parent 068a4d5
Raw File
Tip revision: 45be8f2e32748cd56c3154bbdcd61da3fddda4bd authored by Eh Tan on 09 March 2007, 00:25:03 UTC
A big patch for tracers and composition:
Tip revision: 45be8f2
solver.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>
 * 
 *~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 */


/* 
 * This struct parameterizes those functions which are different
 * between the full and regional solvers.
 */

struct Solver {

    /* Boundary_conditions.c */
    void (*velocity_boundary_conditions)(struct All_variables *);
    void (*temperature_boundary_conditions)(struct All_variables *);

    /* Geometry_cartesian.c */
    void (*set_2dc_defaults)(struct All_variables *);
    void (*set_2pt5dc_defaults)(struct All_variables *);
    void (*set_3dc_defaults)(struct All_variables *);
    void (*set_3dsphere_defaults)(struct All_variables *);

    /* Lith_age.c */
    void (*lith_age_read_files)(struct All_variables *, int);

    /* Parallel_related.c */
    void (*parallel_processor_setup)(struct All_variables *);
    void (*parallel_domain_decomp0)(struct All_variables *);
    void (*parallel_domain_boundary_nodes)(struct All_variables *);
    void (*parallel_communication_routs_v)(struct All_variables *);
    void (*parallel_communication_routs_s)(struct All_variables *);
    void (*exchange_id_d)(struct All_variables *, double **, int);

    /* Read_input_from_files.c */
    void (*read_input_files_for_timesteps)(struct All_variables *, int, int);

    /* Version_dependent.c */
    void (*global_derived_values)(struct All_variables *);
    void (*node_locations)(struct All_variables *);
    void (*construct_tic_from_input)(struct All_variables *);
    void (*construct_boundary)(struct All_variables *);

} solver;
back to top