hydro_parameters.h

Go to the documentation of this file.

/*******************************************************************************
 * This file is part of SWIFT.
 * Copyright (c) 2019 Josh Borrow (joshua.borrow@durham.ac.uk)
 *
 *
 * This program is free software: you can redistribute it and/or modify
 * it under the terms of the GNU Lesser General Public License as published
 * by the Free Software Foundation, either version 3 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 Lesser General Public License
 * along with this program.  If not, see <http://www.gnu.org/licenses/>.
 *
 ******************************************************************************/
 
#ifndef SWIFT_MINIMAL_HYDRO_PARAMETERS_H
#define SWIFT_MINIMAL_HYDRO_PARAMETERS_H
 
/* Configuration file */
#include <config.h>
 
/* Global headers */
#if defined(HAVE_HDF5)
#include <hdf5.h>
#endif
 
/* Local headers */
//#include "common_io.h"
//#include "error.h"
#include "inline.h"
 
/* Viscosity parameters -- FIXED -- MUST BE DEFINED AT COMPILE-TIME */
 
/* Cosmology default beta=3.0.
 * Alpha can be set in the parameter file.
 * Beta is defined as in e.g. Price (2010) Eqn (103) */
#define const_viscosity_beta 3.0f
 
/* The viscosity that the particles are reset to after being hit by a
 * feedback event. This should be set to the same value as the
 * hydro_props_default_viscosity_alpha in fixed schemes, and likely
 * to hydro_props_default_viscosity_alpha_max in variable schemes. */
#define hydro_props_default_viscosity_alpha_feedback_reset 0.8f
 
/* Viscosity paramaters -- Defaults; can be changed at run-time */
 
/* The "initial" hydro viscosity, or the fixed value for non-variable
 * schemes. This usually takes the value 0.8. */
#define hydro_props_default_viscosity_alpha 0.8f
 
/* Structs that store the relevant variables */
 
struct viscosity_global_data {
  float alpha;
};
 
struct diffusion_global_data {};
 
/* Functions for reading from parameter file */
 
/* Forward declartions */
struct swift_params;
struct phys_const;
struct unit_system;
 
/* Viscosity */
 
static INLINE void viscosity_init(struct swift_params* params,
                                  const struct unit_system* us,
                                  const struct phys_const* phys_const,
                                  struct viscosity_global_data* viscosity) {
 
  /* Read the artificial viscosity parameters from the file, if they exist,
   * otherwise set them to the defaults defined above. */
 
  viscosity->alpha = (hydro_props_default_viscosity_alpha);
}
 
static INLINE void viscosity_init_no_hydro(
    struct viscosity_global_data* viscosity) {
  viscosity->alpha = hydro_props_default_viscosity_alpha;
}
 
static INLINE void viscosity_print(
    const struct viscosity_global_data* viscosity) {
//  printf("Artificial viscosity parameters set to alpha: %.3f\n",
//          viscosity->alpha);
}
 
#if defined(HAVE_HDF5)
static INLINE void viscosity_print_snapshot(
    hid_t h_grpsph, const struct viscosity_global_data* viscosity) {
 
  io_write_attribute_f(h_grpsph, "Alpha viscosity", viscosity->alpha);
  io_write_attribute_f(h_grpsph, "Beta viscosity", const_viscosity_beta);
}
#endif
 
/* Diffusion */
 
static INLINE void diffusion_init(struct swift_params* params,
                                  const struct unit_system* us,
                                  const struct phys_const* phys_const,
                                  struct diffusion_global_data* diffusion) {}
 
static INLINE void diffusion_init_no_hydro(
    struct diffusion_global_data* diffusion) {}
 
static INLINE void diffusion_print(
    const struct diffusion_global_data* diffusion) {}
 
#ifdef HAVE_HDF5
static INLINE void diffusion_print_snapshot(
    hid_t h_grpsph, const struct diffusion_global_data* diffusion) {}
#endif
 
#endif /* SWIFT_MINIMAL_HYDRO_PARAMETERS_H */