Commits

Michael Kuhlen committed 532d9ab Merge

Merged with tip, to get the hg log info back for everything in not_yet_converted/

  • Participants
  • Parent commits c48f290, 4258797

Comments (0)

Files changed (99)

 Make.config.override
 *.mod
 
+src/enzo/headers/auto_defaults_string.h
+
 src/mpgrafic/degraf/Makefile
 src/mpgrafic/degraf/config.h
 src/mpgrafic/degraf/config.log

run/Cosmology/Hydro/AMRCosmology/AMRCosmology.cfg

+#
+# AMR PROBLEM DEFINITION FILE: Cosmology Simulation (amr version)
+#
+
+Initialization:
+{
+	ProblemType = 30;
+};
+
+SimulationControl:
+{
+	Domain:
+	{
+		TopGridRank                = 3;
+		TopGridDimensions          = [16, 16, 16];
+		LeftFaceBoundaryCondition  = [3, 3, 3];
+		RightFaceBoundaryCondition = [3, 3, 3];
+	};
+
+	AMR:
+	{
+		StaticHierarchy                 = False;
+		MaximumRefinementLevel          = 2;
+		RefineBy                        = 2;
+		CellFlaggingMethod              = [2, 4];
+		MinimumOverDensityForRefinement = [4.0, 1.5];
+		MinimumMassForRefinementLevelExponent = [-0.3, 0.0];
+	};
+
+	Optimization:
+	{
+		ParallelRootGridIO = True;
+		MinimumEfficiency  = 0.4;
+	};
+
+	StopCycle = 1000;
+};
+
+Physics:
+{
+	Cosmology:
+	{
+		HubbleConstantNow   = 0.5;
+		OmegaMatterNow      = 0.3;
+		OmegaLambdaNow      = 0.7;
+		ComovingBoxSize     = 16.0;
+		MaxExpansionRate    = 0.015;
+		InitialRedshift     = 30.0 ; 
+		FinalRedshift       = 0.0; 
+		ComovingCoordinates = True;
+	};
+
+	Gravity:
+	{
+		GravitationalConstant   = 1.0;
+		GreensFunctionMaxNumber = 30;
+	};
+
+	Hydro:
+	{
+		Gamma                       = 1.6666667;
+		DualEnergyFormalism         = True;
+		InterpolationMethod         = 1;
+		PPMDiffusionParameter       = False;
+		CourantSafetyNumber         = 0.6;
+		ParticleCourantSafetyNumber = 0.5;
+                ConservativeReconstruction  = False;
+	};	
+
+	AtomicPhysics:
+	{	
+		RadiativeCooling = 0;
+		MultiSpecies     = 0;
+	};
+};
+
+Problem:
+{
+	CosmologySimulation:
+	{
+		OmegaBaryonNow       = 0.04;
+		OmegaCDMNow          = 0.26;
+		DensityName          = "GridDensity";
+		VelocitiesName       = "GridVelocities";
+		ParticlePositionName = "ParticlePositions";
+		ParticleVelocityName = "ParticleVelocities";
+	};
+};
+
+
+OutputControl:
+{
+	HierarchyFileInputFormat     = 1;
+	HierarchyFileOutputFormat    = 2;
+
+	DataDump:
+	{
+		Name = "output_";
+		dt   = 10.0;
+	};
+
+	RedshiftDump:
+	{
+		OutputRedshifts = [10.0, 8.0, 6.0, 4.0, 3.0, 2.0, 1.0, 0.0];
+	};
+};
+
+

run/Hydro/Hydro-3D/CollapseTestTwoSphere/CollapseTestTwoSphere.cfg

+#  AMR PROBLEM DEFINITION FILE: Non-cosmological Collapse test
+#  Description: a sphere collapses until becoming pressure supported.
+
+SimulationControl:
+{
+	StopTime            = 1.0;
+#	StopFirstTimeAtLevel    = 10;
+
+	Units:
+	{
+		Density = 3.82e-20;  # 10^4 cm^-3
+		Length = 5.0e17;     # 1 pc in cm
+		Time = 1.0e12;       # 10^4 yrs
+	};
+
+	Domain:
+	{
+		TopGridRank         = 3;
+		TopGridDimensions   = [64, 64, 64];
+		LeftFaceBoundaryCondition  = [3, 3, 3];
+		RightFaceBoundaryCondition = [3, 3, 3];
+	};
+	
+	AMR:
+	{
+		RefineBy                       = 2;
+		MaximumRefinementLevel         = 1;
+		StaticHierarchy                = False;
+		CellFlaggingMethod             = [6];
+		RefineByJeansLengthSafetyFactor  = 8.0;
+	};
+
+	Optimization:
+	{
+		MinimumEfficiency              = 0.3;
+	};
+
+}; # end SimulationControl
+
+
+OutputControl:
+{
+	DataDump:
+	{
+		Name = "DD";
+		Dir  = "DD";
+		dt   = 0.01;
+	};
+
+	OutputTriggers:
+        {
+#                OutputFirstTimeAtLevel = 4;
+        };
+
+	SupplementalFields:
+	{
+		OutputTemperature              = True;
+	};
+
+}; # end OutputControlParameters
+
+
+Physics:
+{
+
+	TopGridGravityBoundary     = 0;
+	SmallRho                   = 1.0e-8;
+	SmallEint                  = 1.0e-8;
+	
+	Hydro:
+	{
+		CourantSafetyNumber    = 0.3;
+		HydroMethod            = 0;   # PPM Direct Euler
+		Gamma                  = 1.0000001;
+
+		PPMDiffusionParameter  = False;
+
+		Mu                         = 3.0;
+
+		FluxCorrection                 = True;
+		InterpolationMethod            = 1;  # SecondOrderA
+		ConservativeInterpolation      = True;
+
+		DualEnergyFormalism         = True;
+
+		RiemannSolver              = 4;
+	};
+
+	ActiveParticles:
+	{
+		ActiveParticlesEnabled        = ["AccretingParticle"];
+	};
+
+	Gravity:
+	{
+		SelfGravity                  = True;
+
+		GravitationalConstant        = 0.0320327594; # 4*pi*G_{cgs}*DensityUnits*TimeUnits^2
+
+		PotentialIterations          = 100;
+		GreensFunctionMaxNumber      = 10;
+	};
+
+
+	AtomicPhysics:
+	{	
+		RadiativeCooling               = False;
+		MultiSpecies                   = 0;
+	};
+
+	Cosmology:
+	{
+		ComovingCoordinates = False;
+	};
+
+}; # end PhysicsParameters
+
+
+Initialization:
+{
+	ProblemType                    = 27;
+}; # end Initialization
+
+
+Problem: {
+	 CollapseTest:
+	 {
+		RefineAtStart   = True;
+		NumberOfSpheres = 2;
+		UseParticles    = 0;
+		InitialTemperature = 500.0; # temperature of the background gas
+		UniformVelocity = [0.0, 0.1, 0.0];
+
+		Spheres = ["Sphere1", "Sphere2"];
+
+		Sphere1: {
+			 Position       = [0.7, 0.2, 0.5];
+			 Velocity       = [0.0, 0.0, 0.0];
+			 Radius         = 0.1;
+			 CoreRadius     = 0.05;  # only used with sphere type 5
+			 Density        = 50.0;  # sphere density, the background density is 1
+			 Temperature    = 10.0;  # put sphere in pressure equilibrium (rho * T is constant)
+			 RotationPeriod = 0.0;
+			 Type           = 5;     # 1: uniform
+				 	         # 2: r^-2 power-law
+					         # 3: NFW
+					         # 4: Gaussian
+					         # 5: r^-2 power-law with a core
+					         # 11: Burkert & Bodenheimer setup
+		};
+
+		Sphere2: {
+			 Position   = [0.3, 0.2, 0.5];
+			 Velocity   = [0.0, 0.0, 0.0];
+			 Radius     = 0.1;
+			 CoreRadius = 0.05;
+			 Density    = 50.0;
+			 Temperature = 10.0;
+			 RotationPeriod = 0.0;
+			 Type       = 5;
+		};
+
+	 };
+}; # end Problem

run/Hydro/Hydro-3D/CollapseTestTwoSphere/CollapseTestTwoSphere.enzo

+#
+#  AMR PROBLEM DEFINITION FILE: Non-cosmological Collapse test
+#  Description: a sphere collapses until becoming pressure supported.
+#
+#  define problem
+#
+ProblemType                = 27         // Collapse test
+TopGridRank                = 3
+TopGridDimensions          = 64 64 64
+SelfGravity                = 1          // gravity on
+TopGridGravityBoundary     = 0          // periodic
+LeftFaceBoundaryCondition  = 3 3 3      // periodic
+RightFaceBoundaryCondition = 3 3 3
+#ExtraOutputs               = 1 2 3 4 5 6 7 8 9
+#
+# problem parameters
+#
+
+CollapseTestRefineAtStart   = 1         // check refinement before running
+CollapseTestNumberOfSpheres = 2
+CollapseTestUseParticles    = 0
+CollapseTestInitialTemperature = 500    // temperature of the background gas
+CollapseTestUniformVelocity = 0.0 0.1 0.0
+CollapseTestSpherePosition[0]   = 0.7 0.2 0.5
+CollapseTestSphereVelocity[0]   = 0.0 0.0 0.0
+CollapseTestSphereRadius[0]     = 0.1
+CollapseTestSphereCoreRadius[0] = 0.05  // only used with sphere type 5
+CollapseTestSphereDensity[0]    = 50   // sphere density, the background density is 1
+CollapseTestSphereTemperature[0] = 10    // put sphere in pressure equilibrium (rho * T is constant)
+CollapseTestSphereRotationPeriod[0] = 0
+CollapseTestSphereType[0]       = 5     // constant density
+                                        // 1: uniform
+					// 2: r^-2 power-law
+					// 3: NFW
+					// 4: Gaussian
+					// 5: r^-2 power-law with a core
+					// 11: Burkert & Bodenheimer setup
+CollapseTestSpherePosition[1]   = 0.3 0.2 0.5
+CollapseTestSphereVelocity[1]   = 0.0 0.0 0.0
+CollapseTestSphereRadius[1]     = 0.1
+CollapseTestSphereCoreRadius[1] = 0.05  // only used with sphere type 5
+CollapseTestSphereDensity[1]    = 50   // sphere density, the background density is 1
+CollapseTestSphereTemperature[1] = 10    // put sphere in pressure equilibrium (rho * T is constant)
+CollapseTestSphereRotationPeriod[1] = 0
+CollapseTestSphereType[1]       = 5     // constant density
+                                        // 1: uniform
+					// 2: r^-2 power-law
+					// 3: NFW
+					// 4: Gaussian
+					// 5: r^-2 power-law with a core
+					// 11: Burkert & Bodenheimer setup
+
+
+
+#  no cosmology for this run
+#
+ComovingCoordinates   = 0              // Expansion OFF
+#
+#  units
+#
+DensityUnits          = 3.82e-20      // 10^4 cm^-3
+LengthUnits           = 5e17          // 1 pc in cm
+TimeUnits             = 1.0e12        // 10^4 yrs
+GravitationalConstant = 0.0320327594   // 4*pi*G_{cgs}*DensityUnits*TimeUnits^2
+#
+#  set I/O and stop/start parameters
+#
+StopTime          = 1.0
+dtDataDump        = 0.01
+#CycleSkipDataDump = 1
+DataDumpDir       = DD
+DataDumpName      = DD
+OutputTemperature = 1                  // Output temperature field.
+#
+#  set hydro parameters
+#
+Gamma                       = 1.0000001
+PPMDiffusionParameter       = 0        // diffusion off
+DualEnergyFormalism         = 1        // use total & internal energy
+InterpolationMethod         = 1        // SecondOrderA
+CourantSafetyNumber         = 0.3
+FluxCorrection              = 1
+ConservativeInterpolation   = 1
+RiemannSolver = 4
+HydroMethod                 = 0        // PPM
+SmallRho                    = 1e-8
+SmallE			    = 1e-8
+#
+#  chemistry/cooling
+#
+MultiSpecies                = 0        // chemistry off
+RadiativeCooling            = 0        // cooling off
+Mu  			    = 3.0
+#
+#  set grid refinement parameters
+#
+StaticHierarchy           = 0          // dynamic hierarchy
+MaximumRefinementLevel    = 1         // use up to 10 levels
+RefineBy                  = 2          // refinement factor
+CellFlaggingMethod        = 6         // use Truelove criterion for refinement 
+MinimumEfficiency         = 0.3
+#OutputFirstTimeAtLevel    = 4         // output when level 4, 5, 6, etc reached (commented out for now)
+#StopFirstTimeAtLevel      = 10         // stop if/when level 10 reached
+RefineByJeansLengthSafetyFactor = 8    // resolve Jeans length by 4 cells (used with CellFlaggingMethod 6)
+
+
+#
+#  set some global parameters
+#
+GreensFunctionMaxNumber   = 10         // # of greens function at any one time
+PotentialIterations	  = 100
+#
+#  active particle parameters
+#
+AppendActiveParticleType = AccretingParticle

src/enzo/Makefile

 #-----------------------------------------------------------------------
 
 .PHONY: autogen
-autogen: io/auto_show_config.C io/auto_show_flags.C io/auto_show_version.C io/auto_show_compile_options.C control/auto_defaults_string.h
+autogen: io/auto_show_config.C io/auto_show_flags.C io/auto_show_version.C io/auto_show_compile_options.C headers/auto_defaults_string.h
 
 # Force update of auto_show_config.C
 
 
 # Force update of auto_defaults_string.h
 
-.PHONY: control/auto_defaults_string.h
-control/auto_defaults_string.h:
-	-@awk 'BEGIN {print "char defaults_string[] = "}; {gsub(/"/,"\\\""); print "  \""$$0"\\n\""}; END {print "  \"\";"}' < control/defaults.cfg > control/auto_defaults_string.h
+.PHONY: headers/auto_defaults_string.h
+headers/auto_defaults_string.h:
+	-@awk 'BEGIN {print "char defaults_string[] = "}; {gsub(/"/,"\\\""); print "  \""$$0"\\n\""}; END {print "  \"\";"}' < control/defaults.cfg > headers/auto_defaults_string.h
 
 #-----------------------------------------------------------------------
 # Generate dependency file

src/enzo/build/machines/Make.mach.trestles-gnu

+#=======================================================================
+#
+# FILE:        Make.mach.trestles
+#
+# DESCRIPTION: Makefile settings for the Trestles Resource at SDSC/UCSD
+#
+# AUTHOR:      John Wise (jwise@astro.princeton.edu)
+#
+# DATE:        07 Dec 2010
+#
+#
+#=======================================================================
+
+MACH_TEXT  = Trestles
+MACH_VALID = 1
+MACH_FILE  = Make.mach.trestles
+
+MACHINE_NOTES = "MACHINE_NOTES for Trestles at SDSC/UCSD: \
+	Load these modules, \
+	'module add gnu/4.6.1 openmpi'"
+
+#-----------------------------------------------------------------------
+# Compiler settings
+#-----------------------------------------------------------------------
+
+LOCAL_MPI_INSTALL = /home/diag/opt/openmpi/1.4.3/gnu/
+LOCAL_PYTHON_INSTALL = /home/mturk/yt-x86_64/
+#LOCAL_COMPILER_DIR = /opt/pgi/linux86-64/10.5
+LOCAL_COMPILER_DIR = /home/diag/opt/gcc/4.6.1/
+LOCAL_HYPRE_INSTALL = 
+
+# With MPI
+
+MACH_CPP       = cpp
+MACH_CC_MPI    = $(LOCAL_MPI_INSTALL)/bin/mpicc # C compiler when using MPI
+MACH_CXX_MPI   = $(LOCAL_MPI_INSTALL)/bin/mpicxx # C++ compiler when using MPI
+MACH_FC_MPI    = $(LOCAL_MPI_INSTALL)/bin/mpif90 # Fortran 77 compiler when using MPI
+MACH_F90_MPI   = $(LOCAL_MPI_INSTALL)/bin/mpif90 # Fortran 90 compiler when using MPI
+MACH_LD_MPI    = $(LOCAL_MPI_INSTALL)/bin/mpicxx # Linker when using MPI
+
+# Without MPI
+
+MACH_CC_NOMPI  = $(LOCAL_COMPILER_DIR)/bin/gcc # C compiler when not using MPI
+MACH_CXX_NOMPI = $(LOCAL_COMPILER_DIR)/bin/g++ # C++ compiler when not using MPI
+MACH_FC_NOMPI  = $(LOCAL_COMPILER_DIR)/bin/gfortran # Fortran 77 compiler when not using MPI
+MACH_F90_NOMPI = $(LOCAL_COMPILER_DIR)/bin/gfortran # Fortran 90 compiler when not using MPI
+MACH_LD_NOMPI  = $(LOCAL_COMPILER_DIR)/bin/g++ # Linker when not using MPI
+
+#-----------------------------------------------------------------------
+# Machine-dependent defines
+#-----------------------------------------------------------------------
+# Defines for the architecture; e.g. -DSUN, -DLINUX, etc.
+MACH_DEFINES   = -DLINUX -DH5_USE_16_API
+
+#-----------------------------------------------------------------------
+# Compiler flag settings
+#-----------------------------------------------------------------------
+
+
+MACH_CPPFLAGS = -P -traditional 
+MACH_CFLAGS   = 
+MACH_CXXFLAGS =
+MACH_FFLAGS   = 
+MACH_F90FLAGS = 
+MACH_LDFLAGS  = 
+
+#-----------------------------------------------------------------------
+# Precision-related flags
+#-----------------------------------------------------------------------
+
+MACH_FFLAGS_INTEGER_32 = 
+MACH_FFLAGS_INTEGER_64 = -fdefault-integer-8
+MACH_FFLAGS_REAL_32    = 
+MACH_FFLAGS_REAL_64    = -fdefault-real-8 -fdefault-double-8
+
+#-----------------------------------------------------------------------
+# Optimization flags
+#-----------------------------------------------------------------------
+
+MACH_OPT_WARN        = -Wall # Flags for verbose compiler warnings
+MACH_OPT_DEBUG       = -O0 -g # Flags for debugging
+# Flags for high conservative optimization
+#MACH_OPT_HIGH        = -O1 -ftz -mieee-fp -fp-speculation=off -prec-sqrt -prec-div
+MACH_OPT_HIGH        = -O2
+# Note that this breaks determinism, which is why it's commented out!
+#
+MACH_OPT_AGGRESSIVE  = -O3 # Flags for aggressive optimization
+# This is the best we can do, from what I can tell.
+#MACH_OPT_AGGRESSIVE  = -O1 -ftz -mieee-fp -fp-speculation=off -prec-sqrt -prec-div
+
+#-----------------------------------------------------------------------
+# Includes
+#-----------------------------------------------------------------------
+
+LOCAL_INCLUDES_MPI    = 
+LOCAL_INCLUDES_HDF5   = -I/home/mturk/yt-x86_64/include # HDF5 includes
+LOCAL_INCLUDES_HYPRE  = 
+LOCAL_INCLUDES_PAPI   = # PAPI includes
+LOCAL_INCLUDES_PYTHON = -I$(LOCAL_PYTHON_INSTALL)/include/python2.7 \
+                        -I$(LOCAL_PYTHON_INSTALL)/lib/python2.7/site-packages/numpy/core/include
+
+MACH_INCLUDES         = $(LOCAL_INCLUDES_HDF5)
+MACH_INCLUDES_PYTHON  = $(LOCAL_INCLUDES_PYTHON)
+MACH_INCLUDES_MPI     = $(LOCAL_INCLUDES_MPI)
+MACH_INCLUDES_HYPRE   = $(LOCAL_INCLUDES_HYPRE)
+MACH_INCLUDES_PAPI    = $(LOCAL_INCLUDES_PAPI)
+
+#-----------------------------------------------------------------------
+# Libraries
+#-----------------------------------------------------------------------
+
+LOCAL_LIBS_MPI    = 
+LOCAL_LIBS_HDF5   = -L/home/mturk/yt-x86_64/lib -lhdf5 # HDF5 libraries
+LOCAL_LIBS_HYPRE  = 
+LOCAL_LIBS_PAPI   = # PAPI libraries
+LOCAL_LIBS_PYTHON  = -L$(LOCAL_PYTHON_INSTALL)/lib/ -lpython2.7 \
+                     -lreadline -ltermcap -lutil
+
+#LOCAL_LIBS_MACH   = -L$(LOCAL_COMPILER_DIR)/lib \
+#			-lpgf90 -lpgf90_rpm1 -lpgf902 -lpgf90rtl -lpgftnrtl -lrt
+LOCAL_LIBS_MACH  = -lgfortran
+
+
+MACH_LIBS         = $(LOCAL_LIBS_HDF5) $(LOCAL_LIBS_MACH)
+MACH_LIBS_MPI     = $(LOCAL_LIBS_MPI)
+MACH_LIBS_HYPRE   = $(LOCAL_LIBS_HYPRE)
+MACH_LIBS_PAPI    = $(LOCAL_LIBS_PAPI)
+MACH_LIBS_PYTHON  = $(LOCAL_LIBS_PYTHON)

src/enzo/control/ParameterBackends.h

-#ifndef PARAMETERBACKENDMGR_H
-#define PARAMETERBACKENDMGR_H
-
-#include <string>
-#include <vector>
-#include <map>
-
-#define MAX_PARAMETERS 1024
-#define MAX_PARAM_LENGTH 1024
-
-// interpreter abstract base class
-class interpreter
-{
-public:
-	explicit interpreter( std::string fname )
-	{ }
-	
-	virtual ~interpreter()
-	{ }
-	
-	virtual int query( std::string key, std::string& ret ) = 0;
-	virtual int query_list( std::string key, std::vector< std::string >& ret ) = 0;
-	
-	virtual size_t size( std::string key ) = 0;
-	
-	virtual bool dump( char fname[], char header_string[]=NULL ) = 0;
-	
-	virtual bool remove( std::string key ) = 0;
-
-	virtual bool update( std::string input_string, bool from_file=false, bool overwrite=false ) = 0;
-
-	virtual bool set( std::string key, int value ) = 0;
-	virtual bool set( std::string key, long long value ) = 0;
-	virtual bool set( std::string key, double value ) = 0;
-	virtual bool set( std::string key, bool value ) = 0;
-	virtual bool set( std::string key, char value[] ) = 0;
-	
-	virtual bool set_list( std::string key, size_t n, const int* value ) = 0;
-	virtual bool set_list( std::string key, size_t n, const long long* value ) = 0;
-	virtual bool set_list( std::string key, size_t n, const double* value ) = 0;
-	virtual bool set_list( std::string key, size_t n, const bool* value ) = 0;
-	virtual bool set_list( std::string key, size_t n, char* value[] ) = 0;
-	
-
-};
-
-// abstract factory pattern
-struct interpreter_creator
-{
-  virtual interpreter* create( std::string fname, std::string defaults ) const = 0;
-	virtual ~interpreter_creator() { }
-};
-
-std::map< std::string, interpreter_creator*>& get_interpreter();
-
-template< class Derived >
-struct interpreter_creator_concrete : public interpreter_creator
-{
-	interpreter_creator_concrete( const std::string& interpreter_name )
-	{
-		get_interpreter()[ interpreter_name ] = this;
-	}
-	
-	interpreter * create( std::string fname, std::string defaults ) const
-	{
-	  return new Derived( fname, defaults );
-	}
-	
-};
-
-#include <stdexcept>
-
-//... helper stuff
-//! runtime error that is thrown if type conversion fails
-class ErrInvalidConversion : public std::runtime_error{
-public:
-	ErrInvalidConversion( std::string errmsg )
-	: std::runtime_error( errmsg )
-	{}
-};
-
-#endif

src/enzo/control/ParameterControl.h

-#ifndef PARAMETER_CONTROL_H
-#define PARAMETER_CONTROL_H
-
-#include <map>
-#include <cstring>
-#include <stdexcept>
-#include <vector>
-#include <sstream>
-#include <typeinfo>
-#include <cstdarg>
-#include <cstdio>
-
-#include <iostream>
-
-#include "ParameterBackends.h"
-
-/*********************************************************************/
-// SOME HELPER FUNCTIONS
-
-//! converts between different variable types
-/*!
- *  The main purpose of this function is to parse and convert
- *  a string argument into numbers, booleans, etc...
- * @param ival the input value (typically a std::string)
- * @param oval the interpreted/converted value
- */
-template <class in_value, class out_value>
-inline void convert( const in_value & ival, out_value & oval) 
-{
-	std::stringstream ss;
-	ss << ival; //.. insert value into stream
-	ss >> oval; //.. retrieve value from stream
-	
-	if (! ss.eof()) {
-		//.. conversion error
-		//std::cerr << "Error: conversion of \'" << ival << "\' failed." << std::endl;
-		throw std::runtime_error(std::string("invalid conversion to ")+typeid(out_value).name()+'.');
-	}
-}
-
-
-template<>
-inline void convert<std::string,std::string>( const std::string & ival, std::string & oval)
-{
-	oval = ival;
-}
-
-
-
-//... actual configuration class
-
-class Configuration
-{
-        interpreter_creator *ic;
-        interpreter* the_interpreter;
-
-	char *argbuf;
-	
-public:
-	
-	Configuration( void )
-	{
-		ic = NULL;
-		the_interpreter = NULL;
-
-		argbuf = new char[1024];
-	}
-	
-	void Initialize( std::string interpreter_name, std::string input_file, const char defaults_string[] )
-	{
-		if( the_interpreter != NULL ) delete the_interpreter;
-		
-		ic = get_interpreter()[interpreter_name];
-		the_interpreter = ic->create( input_file, std::string(defaults_string) );
-	}
-
-	void Update( const char input_string[] )
-	{
-	  // update the current configuration
-	  if( ! the_interpreter->update( std::string(input_string) ) )
-	    {
-	      fprintf(stderr, "Could not update parameter.\n");
-	      throw std::runtime_error("Error updating parameters.");
-	    }
-	  
-	}
-
-	void Dump( char fname[], char header_string[]=NULL )
-	{
-		// dump current configuration to a file
-	  if( ! the_interpreter->dump( fname, header_string ) )
-		{
-			fprintf(stderr, "Could not write parameters to file \'%s\'.\n", fname );
-			throw std::runtime_error("Error writing parameters to file.");
-		}
-	}
-
-	~Configuration()
-	{
-		delete the_interpreter;
-		delete[] argbuf;
-	}
-
-	// eventually...
-	//int PutScalar(...)
-	//int PutArray(...)
-	
-	void Remove( const char* key, ... ) const
-	{
-		va_list argptr;
-		va_start(argptr, key);
-		va_end(argptr);
-		vsprintf(argbuf,key,argptr);
-		
-		int status = 0;
-		status = the_interpreter->remove(std::string(argbuf));
-		
-		if( status != 1 ) {
-			fprintf(stderr, "Could not remove <%s> from parameter file.\n",argbuf);
-			throw std::runtime_error("Could not remove a parameter.");
-		}	
-	}
-	
-
-	template< typename T >
-	void GetScalar( T& val, const char* key, ... ) const
-	{
-		va_list argptr;
-		va_start(argptr, key);
-		va_end(argptr);
-		vsprintf(argbuf,key,argptr);		
-
-		std::string strval;
-		int status = 0;
-		status = the_interpreter->query(std::string(argbuf),strval);
-		if( status != 1 ) {
-		  
-		  fprintf(stderr, "Could not find parameter <%s>.\n",argbuf);
-		  throw std::runtime_error("parameter not found!");
-		}
-
-		convert<std::string,T>(strval,val);
-	}
-	
-	template< typename T >
-	void SetScalar( const T& val, const char* key, ... )
-	{
-		va_list argptr;
-		va_start(argptr, key);
-		va_end(argptr);
-		vsprintf(argbuf,key,argptr);		
-		
-		std::string strval;
-		int status = 0;
-		status = the_interpreter->set(std::string(argbuf),val);
-		if( status != 1 ) {
-			fprintf(stderr, "Could not set <%s> in parameter file.\n",argbuf);
-			throw std::runtime_error("Could not assign parameter!");	
-		}
-	}
-	
-	
-	template< class T >
-	void GetArray( T* val, const char* key, ... ) const
-	{
-		va_list argptr;
-		va_start(argptr, key);
-		va_end(argptr);
-		vsprintf(argbuf,key,argptr);
-		
-		std::vector<std::string> s;
-		int status = 0;
-		status = the_interpreter->query_list(std::string(argbuf),s);
-		if( status != 1 ) {
-		  
-		  fprintf(stderr, "Could not find parameter <%s>.\n",argbuf);
-		}
-		
-		for( size_t i=0; i<s.size(); ++i )
-		  {
-		    convert<std::string,T>(s[i],*val);
-		    ++val;
-		  }
-	}
-	
-	template< typename T >
-	void SetArray( T* val, size_t nelem, const char* key, ... ) const
-	{
-		va_list argptr;
-		va_start(argptr, key);
-		va_end(argptr);
-		vsprintf(argbuf,key,argptr);
-		
-		int status = 0;
-		status = the_interpreter->set_list( std::string(argbuf), nelem, val);
-		if( status != 1 ) {
-			fprintf(stderr, "Could not set <%s> in parameter file.\n",argbuf);
-			throw std::runtime_error("Could not assign parameter!");	
-		}
-	}
-		
-	size_t Size( const char* key, ... ) const
-	{
-		va_list argptr;
-		va_start(argptr, key);
-		va_end(argptr);
-		vsprintf(argbuf,key,argptr);
-		
-		return the_interpreter->size(std::string(argbuf));
-		
-	}
-	
-};
-
-
-
-template<>
-void Configuration::GetScalar<char*>( char*& str, const char* key, ... ) const;
-
-template<>
-void Configuration::GetArray<char*>( char** val, const char* key, ... ) const;
-
-/*********************************************************************/
-/*********************************************************************/
-
-
-/*********************************************************************/
-/*********************************************************************/
-/*********************************************************************/
-
-
-#endif

src/enzo/control/SetDefaultGlobalValues.C

 /
 ************************************************************************/
  
-// This routine intializes a new simulation based on the parameter file.
-//
+// Set default values for global variables that are not set in
+// defaults.cfg.
  
-#include "preincludes.h"
+#include <string.h>
+#include <stdio.h>
 #include "ErrorExceptions.h"
 #include "macros_and_parameters.h"
 #include "typedefs.h"
  
 char DefaultDimUnits[] = "cm";
 char *DefaultDimLabel[] = {"x", "y", "z"};
- 
-char DefaultRestartName[] = "restart";
-char DefaultDataName[] = "data";
-char DefaultHistoryName[] = "history";
-char DefaultRedshiftName[] = "RedshiftOutput";
-char DefaultNewMovieName[] = "MoviePack";
-char DefaultTracerParticleName[] = "TracerOutput";
-char DefaultExtraName[] = "ExtraDumpXX";
-char DefaultExtraDir[]="ED00"; 
-
-char DefaultRestartDir[] = "RS";
-char DefaultDataDir[] = "DD";
-char DefaultHistoryDir[] = "HD";
-char DefaultRedshiftDir[] = "RD";
-char DefaultTracerParticleDir[] = "TD";
- 
- 
- 
- 
+  
 int SetDefaultGlobalValues(TopGridData &MetaData)
 {
  
   const float Pi = 3.14159;
   int dim, i, j;
  
-  huge_number               = 1.0e+20;
-  tiny_number               = 1.0e-20;
-
   /* set the default MetaData values. */
  
-  MetaData.CycleNumber     = 0;
-  MetaData.SubcycleNumber     = 0;
-  MetaData.Time            = 0.0;
-  MetaData.CPUTime         = 0.0;
- 
-  MetaData.StopTime        = FLOAT_UNDEFINED;  // This must be set be the user
-  MetaData.StopCycle       = 100000;            // 10000 timesteps
-  MetaData.StopSteps       = 10000;            // 10000 timesteps
-  MetaData.StopCPUTime     = 720.0*3600.0;     // 30 days
-  MetaData.ResubmitOn      = FALSE;
-  MetaData.ResubmitCommand = NULL;
- 
-  MetaData.MaximumTopGridTimeStep = huge_number;
-
-  MetaData.TimeLastRestartDump = 0.0;
-  MetaData.dtRestartDump       = FLOAT_UNDEFINED;
-  MetaData.TimeLastDataDump    = FLOAT_UNDEFINED;
-  MetaData.dtDataDump          = 0.0;
-  MetaData.TimeLastHistoryDump = FLOAT_UNDEFINED;
-  MetaData.dtHistoryDump       = 0.0;
-  MetaData.TimeLastTracerParticleDump = FLOAT_UNDEFINED;
-  MetaData.dtTracerParticleDump       = 0.0;
-  MetaData.TimeLastInterpolatedDataDump    = FLOAT_UNDEFINED;
-  MetaData.dtInterpolatedDataDump          = 0.0;
   MetaData.WroteData           = FALSE;
- 
-  MetaData.CycleLastRestartDump = 0;
-  MetaData.CycleSkipRestartDump = 0;
-  MetaData.CycleLastDataDump    = INT_UNDEFINED;
-  MetaData.CycleSkipDataDump    = 0;
-  MetaData.SubcycleLastDataDump    = 0;
-  MetaData.SubcycleSkipDataDump    = 0;
-  MetaData.CycleLastHistoryDump = INT_UNDEFINED;
-  MetaData.CycleSkipHistoryDump = 0;
-  MetaData.CycleSkipGlobalDataDump = 0; //AK
- 
-  MetaData.OutputFirstTimeAtLevel = 0; // zero is off
-  MetaData.StopFirstTimeAtLevel   = 0; // zero is off
- 
-  MetaData.NumberOfOutputsBeforeExit = 0;
-  MetaData.OutputsLeftBeforeExit     = 0;
-
-  MetaData.RestartDumpNumber   = 0;            // starting restart id number
-  MetaData.RestartDumpName     = DefaultRestartName;
-  MetaData.RestartDumpDir      = DefaultRestartDir;
-  MetaData.DataDumpNumber      = 0;
-  MetaData.DataDumpName        = DefaultDataName;
-  MetaData.DataDumpDir         = DefaultDataDir;
-  MetaData.HistoryDumpNumber   = 0;
-  MetaData.HistoryDumpName     = DefaultHistoryName;
-  MetaData.HistoryDumpDir      = DefaultHistoryDir;
-  MetaData.TracerParticleDumpNumber = 0;
-  MetaData.TracerParticleDumpName   = DefaultTracerParticleName;
-  MetaData.TracerParticleDumpDir    = DefaultTracerParticleDir;
-  //MetaData.RedshiftDumpNumber  = 0;
-  MetaData.RedshiftDumpName    = DefaultRedshiftName;
-  MetaData.RedshiftDumpDir     = DefaultRedshiftDir;
-  MetaData.ExtraDumpDir        = DefaultExtraDir;
-  MetaData.ExtraDumpName        = DefaultExtraName;
-
-  MetaData.MetaDataIdentifier    = NULL;
-  MetaData.SimulationUUID        = NULL;
-  MetaData.RestartDatasetUUID    = NULL;
-  MetaData.InitialConditionsUUID = NULL;
-
-  MetaData.LocalDir            = NULL;
-  MetaData.GlobalDir           = NULL;
-
-  LoadBalancing = 1;     //On, memory equalization method
-  LoadBalancingCycleSkip = 10;  // Load balance root grids every 10 cycles
-  ResetLoadBalancing = FALSE;
+   
   CoresPerNode = 1;
   PreviousMaxTask = 0;
-  LoadBalancingMinLevel = 0;     //All Levels
-  LoadBalancingMaxLevel = MAX_DEPTH_OF_HIERARCHY;  //All Levels
-
-  FileDirectedOutput = 1;
-
-  // Default Hierarchy File IO settings (1 = ASCII; 2 = HDF5+ASCII)
-  HierarchyFileInputFormat = 1;
-  HierarchyFileOutputFormat = 2;
 
   for (i = 0;i < MAX_DEPTH_OF_HIERARCHY;i++) {
     RebuildHierarchyCycleSkip[i] = 1;
   }
-
-  for (i = 0; i < MAX_TIME_ACTIONS; i++) {
-    TimeActionType[i]      = 0;
-    TimeActionRedshift[i]  = -1;
-    TimeActionTime[i]      = 0;
-    TimeActionParameter[i] = FLOAT_UNDEFINED;
-  }
  
   for (i = 0; i < MAX_CUBE_DUMPS; i++) {
     CubeDumps[i] = NULL;
   }
- 
-  MetaData.StaticHierarchy     = TRUE;
-  FastSiblingLocatorEntireDomain = TRUE;
- 
-  MetaData.TopGridRank = INT_UNDEFINED;
-  for (dim = 0; dim < MAX_DIMENSION; dim++) {
-    MetaData.TopGridDims[dim]                = INT_UNDEFINED;
-    MetaData.LeftFaceBoundaryCondition[dim]  = reflecting;
-    MetaData.RightFaceBoundaryCondition[dim] = reflecting;
-  }
-  MetaData.BoundaryConditionName = NULL;
- 
-  MetaData.GravityBoundary        = TopGridPeriodic;
-
-#ifdef TRANSFER
-  MetaData.RadHydroParameterFname = NULL;
-#endif
- 
-  MetaData.ParticleBoundaryType   = periodic;  // only one implemented!
-  MetaData.NumberOfParticles      = 0;         // no particles
- 
-  MetaData.CourantSafetyNumber    = 0.6;
-  MetaData.PPMFlatteningParameter = 0;    // off
-  MetaData.PPMDiffusionParameter  = 0;    // off
-  MetaData.PPMSteepeningParameter = 0;    // off
-
+  
   MetaData.FirstTimestepAfterRestart = TRUE;
  
-  /* set the default global data. */
-  CheckpointRestart         = 0;
-
-  // Debug flag set in main
-  ProblemType               = 0;                 // None
-  HydroMethod               = PPM_DirectEuler;   //
-  Gamma                     = 5.0/3.0;           // 5/3
-  PressureFree              = FALSE;             // use pressure (duh)
-  RefineBy                  = 4;                 // Refinement factor
-  MaximumRefinementLevel    = 2;                 // three levels (w/ topgrid)
-  MaximumGravityRefinementLevel = INT_UNDEFINED;
-  MaximumParticleRefinementLevel = -1;            // unused if negative
-  MustRefineRegionMinRefinementLevel = -1;        // unused if negative
-  MetallicityRefinementMinLevel = -1;
-  MetallicityRefinementMinMetallicity = 1.0e-5;
-  MetallicityRefinementMinDensity = FLOAT_UNDEFINED;
-  FluxCorrection            = TRUE;
-  InterpolationMethod       = SecondOrderA;      // ?
-  ConservativeInterpolation = TRUE;              // true for ppm
-  MinimumEfficiency         = 0.2;               // between 0-1, usually ~0.1
-  MinimumSubgridEdge        = 6;                 // min for acceptable subgrid
-  MaximumSubgridSize        = 32768;             // max for acceptable subgrid
-  SubgridSizeAutoAdjust     = TRUE; // true for adjusting maxsize and minedge
-  OptimalSubgridsPerProcessor = 16;    // Subgrids per processor
-  NumberOfBufferZones       = 1;
  
-  for (i = 0; i < MAX_FLAGGING_METHODS; i++) {
-    MinimumSlopeForRefinement[i]= 0.3;
-    SlopeFlaggingFields[i] = INT_UNDEFINED;
-    CellFlaggingMethod[i]       = INT_UNDEFINED;
-    MinimumMassForRefinement[i] = FLOAT_UNDEFINED;   // usually set by:
-    MinimumOverDensityForRefinement[i]       = 1.5;
-    MinimumMassForRefinementLevelExponent[i] = 0;
-  }
  
   for (dim = 0; dim < MAX_DIMENSION; dim++) {
-    DomainLeftEdge[dim]             = 0.0;
-    DomainRightEdge[dim]            = 1.0;
-    GridVelocity[dim]               = 0.0;
     DimLabels[dim]                  = DefaultDimLabel[dim];
     DimUnits[dim]                   = DefaultDimUnits;
-    RefineRegionLeftEdge[dim]       = FLOAT_UNDEFINED;
-    RefineRegionRightEdge[dim]      = FLOAT_UNDEFINED;
-    RefineRegionAutoAdjust          = FALSE;
-    MetaData.NewMovieLeftEdge[dim]  = 0.0;
-    MetaData.NewMovieRightEdge[dim] = 1.0;
-    PointSourceGravityPosition[dim] = 0.0;
-    MustRefineRegionLeftEdge[dim] = 0.0;
-    MustRefineRegionRightEdge[dim] = 1.0;
-  }
- 
-  for (i = 0; i < MAX_STATIC_REGIONS; i++) {
-    StaticRefineRegionLevel[i] = INT_UNDEFINED;
-    AvoidRefineRegionLevel[i]  = INT_UNDEFINED;
-    for (dim = 0; dim < MAX_DIMENSION; dim++) {
-      AvoidRefineRegionLeftEdge[i][dim] = FLOAT_UNDEFINED;
-      AvoidRefineRegionRightEdge[i][dim] = FLOAT_UNDEFINED;
-    }
   }
 
   /* For evolving refinement regions. */
-  RefineRegionFile = NULL;
-  RefineRegionTimeType = -1; /* 0=time bins 1=redshift bins*/
   for (i = 0; i < MAX_REFINE_REGIONS; i++) {
     EvolveRefineRegionTime[i] = FLOAT_UNDEFINED;
     for (j = 0; j < MAX_DIMENSION; j++) {
     }
   }
 
-  DatabaseLocation = NULL;
-
- 
-  ParallelRootGridIO          = FALSE;
-  ParallelParticleIO          = FALSE;
-  Unigrid                     = FALSE;
-  UnigridTranspose            = FALSE;
-  NumberOfRootGridTilesPerDimensionPerProcessor = 1;
-  PartitionNestedGrids        = FALSE;
-  ExtractFieldsOnly           = TRUE;
-
   ExternalBoundaryIO          = FALSE;
   ExternalBoundaryTypeIO      = FALSE;
   ExternalBoundaryValueIO     = FALSE;
   debug1                      = 0;
   debug2                      = 0;
 
-  TracerParticleOn            = 0;
-
-  OutputOnDensity                  = 0;
-  StartDensityOutputs              = 999;
-  CurrentDensityOutput             = 999;
-  IncrementDensityOutput           = 999;
   CurrentMaximumDensity            = -999;
  
-  CubeDumpEnabled             = 0;
-
 #ifdef STAGE_INPUT
   StageInput                  = 0;
 #endif
 
   First_Pass                  = 0;
 
-  MemoryLimit                 = 4000000000L;
+  CopyGravPotential           = FALSE;             // off
  
-  ExternalGravity             = FALSE;             // off
-  ExternalGravityDensity      = 0.0;
-  ExternalGravityRadius       = 0.0;
-
-  UniformGravity              = FALSE;             // off
-  UniformGravityDirection     = 0;                 // x-direction
-  UniformGravityConstant      = 1.0;
- 
-  PointSourceGravity           = FALSE;             // off
-  PointSourceGravityConstant   = 1.0;
-  PointSourceGravityCoreRadius = 0.0;
-
-  SelfGravity                 = FALSE;             // off
-  SelfGravityGasOff           = FALSE;             // off
-  AccretionKernal             = FALSE;             // off
-  CopyGravPotential           = FALSE;             // off
-  PotentialIterations         = 4;                 // ~4 is reasonable
-  GravitationalConstant       = 4*Pi;              // G = 1
-  S2ParticleSize              = 3.0;               // ~3 is reasonable
-  GravityResolution           = 1.0;               // equivalent to grid
-  ComputePotential            = FALSE;
-  WritePotential              = FALSE;
-  BaryonSelfGravityApproximation = TRUE;           // less accurate but faster
-
-  GreensFunctionMaxNumber     = 1;                 // only one at a time
-  GreensFunctionMaxSize       = 1;                 // not used yet
- 
-  DualEnergyFormalism         = FALSE;             // off
-  DualEnergyFormalismEta1     = 0.001;             // typical 0.001
-  DualEnergyFormalismEta2     = 0.1;               // 0.08-0.1
-  ParticleCourantSafetyNumber = 0.5;
-  RootGridCourantSafetyNumber = 1.0;
-  RandomForcing               = FALSE;             // off //AK
-  RandomForcingEdot           = -1.0;              //AK
-  RandomForcingMachNumber     = 0.0;               //AK
-  RadiativeCooling            = FALSE;             // off
-  RadiativeCoolingModel       = 1;                 //1=cool_rates.in table lookup
-                                                   //3=Koyama&Inutsuka 2002
-  GadgetEquilibriumCooling    = FALSE;             // off
-  RadiativeTransfer           = 0;                 // off
-  RadiativeTransferFLD        = 0;                 // off
-  ImplicitProblem             = 0;                 // off
-  StarMakerEmissivityField    = 0;                 // off
-  uv_param                    = 1.0e-5;            // mid-range value from Razoumov Norman 2002
-
-  MultiSpecies                = FALSE;             // off
-  NoMultiSpeciesButColors     = FALSE;             // off
-  ThreeBodyRate               = 0;                 // ABN02
-  CIECooling                  = 1;
-  H2OpticalDepthApproximation = 1;
-  H2FormationOnDust           = FALSE;
   GloverChemistryModel        = 0;                 // 0ff
   GloverRadiationBackground   = 0;
   GloverOpticalDepth          = 0;
-  ShockMethod                 = 0;                 // off
-  ShockTemperatureFloor       = 1.0;               // Set to 1K
-  StorePreShockFields         = 0;
-  RadiationFieldType          = 0;
-  RadiationFieldRedshift      = 0.0;
-  TabulatedLWBackground       = 0;
-  RadiationFieldLevelRecompute = 0;
-  RadiationData.RadiationShield = 0;
-  AdjustUVBackground          = 1;
-  AdjustUVBackgroundHighRedshift = 0;
-  SetUVBAmplitude             = 1.0;
-  SetHeIIHeatingScale         = 1.8;
-  PhotoelectricHeating	      = 0;
-  PhotoelectricHeatingRate    = 8.5e-26;           // ergs cm-3 s-1
-  RadiationXRaySecondaryIon   = 0;
-  RadiationXRayComptonHeating = 0;
 
-  CoolData.alpha0             = 1.5;               // radiation spectral slope
-  CoolData.f3                 = 1.0e-21;           // radiation normalization
-  CoolData.f0to3                    = 0.1;
-  CoolData.RadiationRedshiftOn      = 7.0;
-  CoolData.RadiationRedshiftOff     = 0.0;
-  CoolData.RadiationRedshiftFullOn  = 6.0;
-  CoolData.RadiationRedshiftDropOff = 0.0;
-  CoolData.HydrogenFractionByMass   = 0.76;
-  /* The DToHRatio is by mass in the code, so multiply by 2. */
-  CoolData.DeuteriumToHydrogenRatio = 2.0*3.4e-5; // Burles & Tytler 1998
-  CoolData.SolarMetalFractionByMass = 0.02041;
-  CoolData.NumberOfTemperatureBins = 600;
-  CoolData.ih2co                   = 1;
-  CoolData.ipiht                   = 1;
-  CoolData.TemperatureStart        = 1.0;
-  CoolData.TemperatureEnd          = 1.0e8;
-  CoolData.comp_xray               = 0;
-  CoolData.temp_xray               = 0;
-  RateData.CaseBRecombination      = 0;   // default to case A rates
-  RateData.NumberOfDustTemperatureBins = 250;
-  RateData.DustTemperatureStart    = 1.0;
-  RateData.DustTemperatureEnd      = 1500.0;
+ 
+  StarFeedbackDistTotalCells       = 1;
 
-  CloudyCoolingData.CloudyCoolingGridRank          = 0;
-  CloudyCoolingData.CloudyCoolingGridFile          = "";
-  CloudyCoolingData.IncludeCloudyHeating           = 0;
-  CloudyCoolingData.CMBTemperatureFloor            = 1;         // use CMB floor.
-  CloudyCoolingData.CloudyElectronFractionFactor = 9.153959e-3; // calculated using Cloudy 07.02 abundances
-
-  OutputCoolingTime = FALSE;
-  OutputTemperature = FALSE;
-  OutputDustTemperature = FALSE;
-
-  OutputSmoothedDarkMatter = FALSE;
-  SmoothedDarkMatterNeighbors = 32;
-
-  OutputGriddedStarParticle = FALSE;
-
-  ZEUSLinearArtificialViscosity    = 0.0;
-  ZEUSQuadraticArtificialViscosity = 2.0;
-  UseMinimumPressureSupport        = FALSE;
-  MinimumPressureSupportParameter  = 100.0;
- 
-  //MinimumSlopeForRefinement        = 0.3;          // 30% change in value
-  MinimumShearForRefinement        = 1.0;          //AK
-  MinimumPressureJumpForRefinement = 0.33;         // As in PPM method paper
-  MinimumEnergyRatioForRefinement  = 0.1;          // conservative!
-  RefineByJeansLengthSafetyFactor  = 4.0;
-  JeansRefinementColdTemperature  = -1.0;
-  RefineByResistiveLengthSafetyFactor  = 2.0;
-  ShockwaveRefinementMinMach = 1.3; // Only above M=1.3
-  ShockwaveRefinementMinVelocity = 1.0e7; //1000 km/s
-  ShockwaveRefinementMaxLevel = 0; 
-  MustRefineParticlesRefineToLevel = 0;
-  MustRefineParticlesRefineToLevelAutoAdjust = FALSE;
-  MustRefineParticlesMinimumMass   = 0.0;
-  ComovingCoordinates              = FALSE;        // No comoving coordinates
-  StarParticleCreation             = FALSE;
-  StarParticleFeedback             = FALSE;
-  BigStarFormation                 = FALSE;
-  BigStarFormationDone             = FALSE;
-  BigStarSeparation                = 0.25;
-  SimpleQ                          = 1e50;
-  SimpleRampTime                   = 0.1;
-  StarFormationOncePerRootGridTimeStep = FALSE;
-  StarMakerTypeIaSNe               = FALSE;
-  StarMakerPlanetaryNebulae        = FALSE;
-  StarMakerOverDensityThreshold    = 100;          // times mean total density
-  StarMakerSHDensityThreshold      = 7e-26;        // cgs density for rho_crit in Springel & Hernquist star_maker5
-  StarMakerMassEfficiency          = 1;
-  StarMakerMinimumMass             = 1.0e9;        // in solar masses
-  StarMakerMinimumDynamicalTime    = 1.0e6;        // in years
-  StarMassEjectionFraction         = 0.25;
-  StarMetalYield                   = 0.02;
-  StarEnergyToThermalFeedback      = 1.0e-5;
-  StarEnergyToStellarUV            = 3.0e-6;
-  StarEnergyToQuasarUV             = 5.0e-6;
-  StarFeedbackDistRadius           = 0;
-  StarFeedbackDistCellStep         = 0;
-  StarFeedbackDistTotalCells       = 1;
-  MultiMetals                      = FALSE;
-  NumberOfParticleAttributes       = INT_UNDEFINED;
-  ParticleTypeInFile               = TRUE;
-  ReadGhostZones                   = FALSE;
-  WriteGhostZones                  = FALSE;
-  OutputParticleTypeGrouping       = FALSE;
-
-  IsotropicConduction = FALSE;
-  AnisotropicConduction = FALSE;
-  IsotropicConductionSpitzerFraction = 1.0;
-  AnisotropicConductionSpitzerFraction = 1.0;
-  ConductionCourantSafetyNumber = 0.5;
-
-  PythonTopGridSkip                = 0;
-  PythonSubcycleSkip               = 1;
-  PythonReloadScript               = FALSE;
-  
-  // EnzoTiming Dump Frequency
-  TimingCycleSkip                  = 1;
-
-  InlineHaloFinder                 = FALSE;
-  HaloFinderSubfind                = FALSE;
-  HaloFinderOutputParticleList     = FALSE;
-  HaloFinderRunAfterOutput         = TRUE;
-  HaloFinderMinimumSize            = 50;
-  HaloFinderLinkingLength          = 0.1;
-  HaloFinderCycleSkip              = 3;
-  HaloFinderTimestep               = FLOAT_UNDEFINED;
-  HaloFinderLastTime               = 0.0;
-
-  StarClusterUseMetalField         = FALSE;
-  StarClusterUnresolvedModel       = FALSE;
-  StarClusterHeliumIonization      = FALSE;
-  StarClusterMinDynamicalTime      = 10e6;         // in years
-  StarClusterIonizingLuminosity    = 1e47;         // ph/s / Msun
-  StarClusterSNEnergy              = 6.8e48;       // erg / Msun (Woosley&Weaver86)
-  StarClusterSNRadius              = 10;           // pc
-  StarClusterFormEfficiency        = 0.1;
-  StarClusterMinimumMass           = 1000;         // Msun
-  StarClusterCombineRadius         = 10;           // pc
-  for (dim = 0; dim < MAX_DIMENSION; dim++) {
-    StarClusterRegionLeftEdge[dim] = 0.0;
-    StarClusterRegionRightEdge[dim] = 1.0;
-  }
-
-  PopIIIStarMass                   = 100;
-  PopIIIInitialMassFunction        = FALSE;
-  PopIIIInitialMassFunctionSeed    = INT_UNDEFINED;
-  PopIIIInitialMassFunctionCalls   = 0;
-  PopIIIHeliumIonization           = FALSE;
-  PopIIILowerMassCutoff            = 1.0;
-  PopIIIUpperMassCutoff            = 300.0;
-  PopIIIInitialMassFunctionSlope   = -1.3;         // high mass slope
-  PopIIIBlackHoles                 = FALSE;
-  PopIIIBHLuminosityEfficiency     = 0.1;
-  PopIIIOverDensityThreshold       = 1e6;          // times mean total density
-  PopIIIH2CriticalFraction         = 5e-4;
-  PopIIIMetalCriticalFraction      = 1e-4;
-  PopIIISupernovaRadius            = 1;            // pc
-  PopIIISupernovaUseColour         = FALSE;
-  PopIIISupernovaMustRefine        = FALSE;
-  PopIIISupernovaMustRefineResolution = 32;
-  PopIIIColorDensityThreshold      = 1e6;          // times mean total density
-  PopIIIColorMass                  = 1e6;          // total mass to color
   IMFData                          = NULL;
 
-  MBHAccretion                     = FALSE;        // 1: Bondi rate, 2: fix temperature, 3: fix rate, 4: Bondi with v_rel=0, 5: Bondi with v_rel=0 and vorticity
-  MBHAccretionRadius               = 50;           // pc
-  MBHAccretingMassRatio            = 1.0;          // 100%, check Star_CalculateMassAccretion.C
-  MBHAccretionFixedTemperature     = 3e5;          // K,       for MBHAccretion = 2
-  MBHAccretionFixedRate            = 1e-3;         // Msun/yr, for MBHAccretion = 3
-  MBHTurnOffStarFormation          = FALSE;        // check Grid_StarParticleHandler.C
-  MBHCombineRadius                 = 50;           // pc
-  MBHMinDynamicalTime              = 10e6;         // in years
-  MBHMinimumMass                   = 1e3;          // Msun
-
-  MBHFeedback                      = FALSE;        // 1: isotropic thermal, 2: jet along z, 3: jet along L, 4: jet along L with 10deg noise, 5: jet along random direction
-  MBHFeedbackRadiativeEfficiency   = 0.1;          // Shakura & Sunyaev (1973)
-  MBHFeedbackEnergyCoupling        = 0.05;         // Springel (2005), Di Matteo (2005)
-  MBHFeedbackMassEjectionFraction  = 0.1;          // 10%, check Star_CalculateFeedbackParameters.C
-  MBHFeedbackMetalYield            = 0.02;         // 2%, check Star_CalculateFeedbackParameters.C
-  MBHFeedbackThermalRadius         = 50;           // pc
-  MBHFeedbackJetsThresholdMass     = 10;           // Msun
-
-  /* Star Class MBH Paricle IO (PARTICLE_TYPE_MBH) */
-  MBHParticleIO                    = FALSE;
-  MBHParticleIOFilename            = (char*) "mbh_particle_io.dat";
-  MBHInsertLocationFilename        = (char*) "mbh_insert_location.in";
   OutputWhenJetsHaveNotEjected     = FALSE;
 
-  H2StarMakerEfficiency = 0.01;
-  H2StarMakerNumberDensityThreshold = 0.0;
-  H2StarMakerMinimumMass = 0.0;
-  H2StarMakerMinimumH2FractionForStarFormation = 1e-5;
-  H2StarMakerStochastic = 1;
-  H2StarMakerUseSobolevColumn = 0;
-  H2StarMakerSigmaOverR = 1.0/30.0;
-  H2StarMakerAssumeColdWarmPressureBalance = 0;
-  H2StarMakerH2DissociationFlux_MW = 1.0;
-  H2StarMakerH2FloorInColdGas = 0.0;
-  H2StarMakerColdGasTemperature = 1e4;
-
-  GMCParticleRNGSeed               = INT_UNDEFINED;
-  GMCParticleRNGCalls              = 0;
-
-  NumberOfParticleAttributes       = INT_UNDEFINED;
-  AddParticleAttributes            = FALSE;
   LastSupernovaTime                = FLOAT_UNDEFINED;
 
-  for (i = 0; i<MAX_MOVIE_FIELDS; i++)
-    MovieDataField[i] = INT_UNDEFINED;
-  MovieSkipTimestep = INT_UNDEFINED;
-  NewMovieName = DefaultNewMovieName;
-  NewMovieDumpNumber = 0;
-  NewMovieParticleOn = FALSE;
-  Movie3DVolumes  = FALSE;
-  MovieVertexCentered = FALSE;
-  MetaData.MovieTimestepCounter      = 0;
-
   ran1_init = 0;
   rand_init = 0;
 
-  SinkMergeDistance     = 1e16;
-  SinkMergeMass         = 0.1;
-  TotalSinkMass         = 0.0;
-  StellarWindFeedback   = 0;
-  StellarWindTurnOnMass = 0.1;
-  MSStellarWindTurnOnMass = 10.0;
-
-  UseHydro		     = 1;
-  Coordinate		     = Cartesian;
   NSpecies		     = INT_UNDEFINED;
   NColor		     = INT_UNDEFINED;
-  Theta_Limiter		     = 1.5;
-  RKOrder		     = 2;
-  UsePhysicalUnit	     = 0;
   NEQ_HYDRO		     = 5;
   NEQ_MHD		     = 9;
-  SmallRho		     = 1e-30;
-  SmallP		     = 1e-35;
   SmallEint		     = 1e-30;
-  SmallT		     = 1e-10;
-  MaximumAlvenSpeed	     = 1e30;
-  RiemannSolver		     = INT_UNDEFINED;
-  RiemannSolverFallback      = 1;
-  ReconstructionMethod	     = INT_UNDEFINED;
-  PositiveReconstruction     = FALSE;
-  ConservativeReconstruction = 0;
-  EOSType		     = 0;
-  EOSSoundSpeed		     = 2.65e4;
-  EOSCriticalDensity	     = 1e-13;
-  EOSGamma		     = 1.667;
-  Mu			     = 0.6;
-  DivBDampingLength          = 1.;
-  CoolingCutOffDensity1	     = 0;
-  CoolingCutOffDensity2	     = 1e10;
-  CoolingCutOffTemperature   = 0.0;
-  CoolingPowerCutOffDensity1 = 0;
-  CoolingPowerCutOffDensity2 = 1e10;
-  UseCUDA		     = 0;
-  UseFloor		     = 0;
-  UseViscosity		     = 0;
-  ViscosityCoefficient       = 0.;
-  UseAmbipolarDiffusion	     = 0;
-  UseResistivity	     = 0;
-
-  StringKick = 0;
-  StringKickDimension = 0;
 
   iden	= 0;
   ivx	= 1;
   iEtot = 4;
   iEint = 0;
 
-  UseDivergenceCleaning		   = 0;
-  DivergenceCleaningBoundaryBuffer = 0;
-  DivergenceCleaningThreshold	   = 0.001;
-  PoissonApproximationThreshold	   = 0.001;
-  PoissonBoundaryType	   = 0;
-
-  UseDrivingField   = 0;
-  DrivingEfficiency = 1.0;
-
   /* End of Stanford Hydro additions */
 
   /* test problem values */
 
   TestProblemData.UseMetallicityField = 0;
   TestProblemData.MetallicityField_Fraction = tiny_number;
-  TestProblemData.MetallicitySNIaField_Fraction = tiny_number;
 
   TestProblemData.UseMassInjection = 0;
   TestProblemData.InitialHydrogenMass = tiny_number;
   // It could also be used (cautiously) for other purposes.
   LoadGridDataAtStart = TRUE;
 
-  IsothermalSoundSpeed = 1.0;
-  RefineByJeansLengthUnits = 0;
-
-  MetalCooling = FALSE;
-  MetalCoolingTable = (char*) "metal_cool.dat";
 
 #ifdef USE_PYTHON
   NumberOfPythonCalls = 0;
     dtThisLevelSoFar[i] = dtThisLevel[i] = 0.0;
   }
 
-  /* Shearing Boundary Conditions variables */
-
-  AngularVelocity=0.001;
-  VelocityGradient=1.0;
-  ShearingBoundaryDirection=-1;
-  ShearingVelocityDirection=-1;
-  ShearingBoxProblemType = 0; 
-  useMHD=0;
-
-  for(int i=0; i<MAX_EXTRA_OUTPUTS; i++) ExtraOutputs[i]=INT_UNDEFINED;
-  MoveParticlesBetweenSiblings = TRUE;
-
-  /* Particle Splitter */
-
-  ParticleSplitterIterations = FALSE;
-  ParticleSplitterChildrenParticleSeparation = 1.0;
-
-  /* Magnetic Field Resetter */
-
-  ResetMagneticField = FALSE;
-  for (dim = 0; dim < MAX_DIMENSION; dim++) {
-    ResetMagneticFieldAmplitude[dim] = 0.0;   // in Gauss
-  }  
-
-  VelAnyl                     = 0;
-  BAnyl                     = 0;
-
-  for (i = 0; i < MAX_ACTIVE_PARTICLE_TYPES; i++) EnabledActiveParticles[i] = NULL;
-  EnabledActiveParticlesCount = 0;
-  UnfulfilledStarFormationMass = 0;
-
-  /* Gas drag parameters */
-  UseGasDrag = 0;
-  GasDragCoefficient = 0.;
-
   return SUCCESS;
 }

src/enzo/control/auto_defaults_string.h

-char defaults_string[] = 
-  "### ENZO DEFAULT PARAMETER VALUES ###\n"
-  "\n"
-  "Internal:\n"
-  "{\n"
-  "	Fields = [\"Density\", \"X-velocity\", \"Y-velocity\", \"Z-velocity\", \"TotalEnergy\", \"GasEnergy\"];\n"
-  "\n"
-  "	Density:\n"
-  "	{\n"
-  "		Name = \"Density\";\n"
-  "		cgsConversionFactor = \"\";\n"
-  "	};	\n"
-  "	X-velocity:\n"
-  "	{\n"
-  "		Name = \"x-velocity\";\n"
-  "		cgsConversionFactor = \"\";\n"
-  "	};\n"
-  "	Y-velocity:\n"
-  "	{\n"
-  "		Name = \"y-velocity\";\n"
-  "		cgsConversionFactor = \"\";\n"
-  "	};\n"
-  "	Z-velocity:\n"
-  "	{\n"
-  "		Name = \"z-velocity\";\n"
-  "		cgsConversionFactor = \"\";\n"
-  "	};\n"
-  "	TotalEnergy:\n"
-  "	{\n"
-  "		Name = \"TotalEnergy\";\n"
-  "		cgsConversionFactor = \"\";\n"
-  "	};\n"
-  "	GasEnergy:\n"
-  "	{\n"
-  "		Name = \"GasEnergy\";\n"
-  "		cgsConversionFactor = \"\";\n"
-  "	};\n"
-  "\n"
-  "#	Units:\n"
-  "#	{\n"
-  "#		Mass = 1.0;\n"
-  "#		Density = 1.0;\n"
-  "#		Length = 1.0;\n"
-  "#		Time = 1.0;\n"
-  "#		Temperature = 1.0;\n"
-  "#	};\n"
-  "\n"
-  "	Provenance:\n"
-  "	{\n"
-  "		MetaDataIdentifier      = \"\";\n"
-  "		SimulationUUID          = \"\";\n"
-  "#		DatasetUUID             = \"\";\n"
-  "		RestartDatasetUUID      = \"\";\n"
-  "		InitialConditionsUUID   = \"\";\n"
-  "		VersionNumber           = \"2.0\";\n"
-  "	};\n"
-  "\n"
-  "	OutputLabeling:\n"
-  "	{\n"
-  "		TimeLastRestartDump = 0.0;\n"
-  "		TimeLastDataDump    = -99999.0;\n"
-  "		TimeLastHistoryDump = -99999.0;\n"
-  "		TimeLastTracerParticleDump = -99999.0;\n"
-  "		TimeLastInterpolatedDataDump    = -99999.0;\n"
-  "		MovieTimestepCounter = 0;\n"
-  "		CycleLastRestartDump = 0;\n"
-  "		CycleLastDataDump    = -99999;\n"
-  "		CycleLastHistoryDump = -99999;\n"
-  "		SubcycleLastDataDump    = -99999;\n"
-  "		RestartDumpNumber   = 0;\n"
-  "		DataDumpNumber      = 0;\n"
-  "		HistoryDumpNumber   = 0;\n"
-  "		TracerParticleDumpNumber = 0;\n"
-  "		CheckpointRestart   = False;\n"
-  "	};\n"
-  "\n"
-  "	InitialCycleNumber  = 0;\n"
-  "	SubcycleNumber      = 0;\n"
-  "	InitialTime         = 0.0;\n"
-  "	Initialdt           = 0.01;\n"
-  "	InitialCPUTime      = 0.0;\n"
-  "	NumberOfParticles      = 0;\n"
-  "	Debug1                          = False;\n"
-  "	Debug2                          = False;\n"
-  "	CosmologyCurrentRedshift   = -99999.0;\n"
-  "	BoundaryConditionName      = \"\";\n"
-  "};\n"
-  "\n"
-  "\n"
-  "SimulationControl: # (subgroup refinement for later on)\n"
-  "{\n"
-  "	StopTime            = -99999.0;\n"
-  "	StopCycle           = 100000;\n"
-  "	StopSteps           = 10000;\n"
-  "	StopCPUTime         = 2592000.0;  # 30 days\n"
-  "	ResubmitOn          = False;\n"
-  "	ResubmitCommand     = \"\";\n"
-  "	MaximumTopGridTimeStep = 1e20;\n"
-  "	TracerParticleOn           = False;\n"
-  "	StopFirstTimeAtLevel    = 0;\n"
-  "	UseCUDA                        = False;\n"
-  "	huge_number                    = 1.000000e20;\n"
-  "	tiny_number                    = 1.000000e-20;\n"
-  "\n"
-  "	Domain:\n"
-  "	{\n"
-  "		TopGridRank         = -99999;\n"
-  "		TopGridDimensions   = [-99999, -99999, -99999];\n"
-  "		LeftFaceBoundaryCondition  = [0, 0, 0];\n"
-  "		RightFaceBoundaryCondition = [0, 0, 0];\n"
-  "		RefineRegionLeftEdge   = [0.0, 0.0, 0.0];\n"
-  "		RefineRegionRightEdge  = [1.0, 1.0, 1.0];\n"
-  "		DomainLeftEdge         = [0.0, 0.0, 0.0];\n"
-  "		DomainRightEdge        = [1.0, 1.0, 1.0];\n"
-  "		RefineRegionAutoAdjust = False;\n"
-  "		MustRefineRegionLeftEdge = [0.0, 0.0, 0.0]; \n"
-  "		MustRefineRegionRightEdge  = [1.0, 1.0, 1.0]; \n"
-  "	};\n"
-  "	\n"
-  "	AMR:\n"
-  "	{\n"
-  "		RefineBy                       = 2;\n"
-  "		MaximumRefinementLevel         = 2;\n"
-  "		MaximumGravityRefinementLevel  = -99999;\n"
-  "		MaximumParticleRefinementLevel = -1;\n"
-  "		NumberOfBufferZones            = 1;\n"
-  "		StaticHierarchy                = True;\n"
-  "		CellFlaggingMethod             = [-99999, -99999, -99999, -99999, -99999, -99999, -99999, -99999, -99999];\n"
-  "		RefineByJeansLengthSafetyFactor  = 4.0;\n"
-  "		JeansRefinementColdTemperature  = -1.0;\n"
-  "		RefineByResistiveLengthSafetyFactor  = 2.000000;\n"
-  "		MustRefineParticlesRefineToLevel = False;\n"
-  "		MustRefineParticlesRefineToLevelAutoAdjust = False;\n"
-  "		MustRefineParticlesMinimumMass = 0.0;\n"
-  "		MustRefineRegionMinRefinementLevel  = -1;\n"
-  "		MetallicityRefinementMinLevel       = -1;\n"
-  "		MetallicityRefinementMinMetallicity = 1e-05;\n"
-  "		MetallicityRefinementMinDensity     = -99999.0;\n"
-  "		RefineRegionTimeType   = -1;\n"
-  "		SlopeFlaggingFields = [-99999, -99999, -99999, -99999, -99999, -99999, -99999, -99999, -99999];\n"
-  "		MinimumSlopeForRefinement = [0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3, 0.3];\n"
-  "		MinimumOverDensityForRefinement = [1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5, 1.5];\n"
-  "		MinimumMassForRefinement = [-99999.0, -99999.0, -99999.0, -99999.0, -99999.0, -99999.0, -99999.0, -99999.0, -99999.0];\n"
-  "		MinimumMassForRefinementLevelExponent = [0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0];\n"
-  "		MinimumShearForRefinement             = 1.0;\n"
-  "		MinimumPressureJumpForRefinement      = 0.33;\n"
-  "		MinimumEnergyRatioForRefinement       = 0.1;\n"
-  "		ShockwaveRefinementMinMach            = 1.3;\n"
-  "		ShockwaveRefinementMinVelocity        = 1.0e7;\n"
-  "		ShockwaveRefinementMaxLevel           = 0;\n"
-  "		RefineByJeansLengthUnits              = False;\n"
-  "	\n"
-  "		StaticRefineRegion:\n"
-  "		{\n"
-  "			Regions=[];\n"
-  "		};\n"
-  "\n"
-  "		RefineRegionFile = \"\";\n"
-  "	\n"
-  "	\n"
-  "	}; # end amr\n"
-  "\n"
-  "	Optimization:\n"
-  "	{\n"
-  "\n"
-  "		LoadBalancing          = 1;\n"
-  "		ResetLoadBalancing     = False;\n"
-  "		LoadBalancingCycleSkip = 10;\n"
-  "		LoadBalancingMinLevel  = 0;\n"
-  "		LoadBalancingMaxLevel  = 50;\n"
-  "		RadiativeTransferLoadBalance    = False;\n"
-  "		ParallelRootGridIO              = False;\n"
-  "		ParallelParticleIO              = False;\n"
-  "		Unigrid                         = False;\n"
-  "		UnigridTranspose                = False;\n"
-  "		MinimumEfficiency              = 0.2;\n"
-  "		SubgridSizeAutoAdjust          = True;\n"
-  "		OptimalSubgridsPerProcessor    = 16;\n"
-  "		MinimumSubgridEdge             = 6;\n"
-  "		MaximumSubgridSize             = 32768;\n"
-  "		NumberOfRootGridTilesPerDimensionPerProcessor = 1;\n"
-  "		MoveParticlesBetweenSiblings     = True;\n"
-  "		ParticleSplitterIterations       = False;\n"
-  "		ParticleSplitterChildrenParticleSeparation     = 1.0;\n"
-  "		MemoryLimit                     = 4000000000;\n"
-  "		FastSiblingLocatorEntireDomain      = True;\n"
-  "	};\n"
-  "\n"
-  "	Timeaction:\n"
-  "	{\n"
-  "		Actions = [];\n"
-  "	};\n"
-  "\n"
-  "}; # end SimulationControl\n"
-  "\n"
-  "\n"
-  "\n"
-  "OutputControl:\n"
-  "{\n"
-  "\n"
-  "	dtInterpolatedDataDump       = 0.0;\n"
-  "	FileDirectedOutput           = True;\n"
-  "	WriteBinaryHierarchy         = False;\n"
-  "	GlobalDir                    = \"\";\n"
-  "	LocalDir                     = \"\";\n"
-  "	ParticleTypeInFile           = True;\n"
-  "	OutputParticleTypeGrouping   = False;\n"
-  "	ExtractFieldsOnly            = True;\n"
-  "	CubeDumpEnabled              = False;\n"
-  "\n"
-  "\n"
-  "	DataDump:\n"
-  "	{\n"
-  "		Name = \"data\";\n"
-  "		Dir  = \"DD\";\n"
-  "		dt   = 0.0;\n"
-  "	};\n"
-  "\n"
-  "	RestartDump:\n"
-  "	{\n"
-  "		Name = \"restart\";\n"
-  "		Dir  = \"RS\";\n"
-  "		dt   = -99999.0;\n"
-  "	};\n"
-  "\n"
-  "	MovieDump:\n"
-  "	{\n"
-  "		Name                  = \"MoviePack\";\n"
-  "		LeftEdge              = [0.0, 0.0, 0.0];\n"
-  "		RightEdge             = [1.0, 1.0, 1.0];\n"
-  "		SkipTimestep          = -99999;\n"
-  "		Volumes3D             = False;\n"
-  "		VertexCentered        = False;\n"
-  "		ParticleOn            = False;\n"
-  "		DataField             = [-99999, -99999, -99999, -99999, -99999, -99999];\n"
-  "		DumpNumber            = 0;\n"
-  "	};\n"
-  "\n"
-  "	CycleDump:\n"
-  "	{\n"
-  "		SkipRestartDump = 0;\n"
-  "		SkipDataDump    = 0;\n"
-  "		SkipHistoryDump = 0;\n"
-  "	        SkipGlobalDataDump = 0;\n"
-  "		SubcycleSkipDataDump    = 0;\n"
-  "	};\n"
-  "\n"
-  "	RedshiftDump:\n"
-  "	{\n"
-  "		Name        = \"RedshiftOutput\";\n"
-  "		Dir         = \"RD\";\n"
-  "		OutputRedshifts = [];\n"
-  "		OutputRedshiftNames = [];\n"
-  "	};\n"
-  "\n"
-  "	HistoryDump:\n"
-  "	{\n"
-  "		Name     = \"history\";\n"
-  "		Dir      = \"HD\";\n"
-  "		dt       = 0.0;\n"