1. Matthew Turk
  2. enzo-3.0-event-handling

Commits

Michael Kuhlen  committed 37de62a

The new configuration now works with CosmologySimulation (ProblemType==30). Try it out with run/Cosmology/Hydro/AMRCosmology/AMRCosmology.cfg. I moved all problem types that have not yet been converted to problem_types/not_yet_converted/. Most of these I've already converted in my new-config branch of enzo-dev, and I simply need to copy those over.

  • Participants
  • Parent commits 9433149
  • Branches active_particles

Comments (0)

Files changed (155)

File run/Cosmology/Hydro/AMRCosmology/AMRCosmology.cfg

View file
+#
+# 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];
+	};
+};
+
+

File src/enzo/Makefile

View file
 #-----------------------------------------------------------------------
 
 .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

File 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

File 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

File src/enzo/control/SetDefaultGlobalValues.C

View file
 /
 ************************************************************************/
  
-// 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;
 }

File 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"
-  "	};\n"
-  "\n"
-  "	TracerParticleDump:\n"
-  "	{\n"
-  "		Name     = \"TracerOutput\";\n"
-  "		Dir      = \"TD\";\n"
-  "		dt       = 0.0;\n"
-  "	};\n"
-  "\n"
-  "	OutputTriggers:\n"
-  "	{\n"
-  "		OutputFirstTimeAtLevel = 0;\n"
-  "		StartDensityOutputs    = -99999;\n"
-  "		OutputOnDensity        = False;\n"
-  "		CurrentDensityOutput   = -99999;\n"
-  "		IncrementDensityOutput = -99999.0;\n"
-  "	};\n"
-  "\n"
-  "	SupplementalFields:\n"
-  "	{\n"
-  "		WritePotential                 = False;\n"
-  "		OutputCoolingTime              = False;\n"
-  "		OutputTemperature              = False;\n"
-  "		OutputSmoothedDarkMatter       = False;\n"
-  "		SmoothedDarkMatterNeighbors    = 32;\n"
-  "		OutputGriddedStarParticle      = False;\n"
-  "		BAnyl                          = False;\n"
-  "		VelAnyl                        = False;\n"
-  "	};\n"
-  "\n"
-  "}; # end OutputControlParameters\n"
-  "\n"
-  "\n"
-  "\n"
-  "\n"
-  "\n"
-  "Physics:\n"
-  "{\n"
-  "\n"
-  "	TopGridGravityBoundary     = 0;\n"
-  "	ParticleBoundaryType       = 3;\n"
-  "	GridVelocity               = [0.0 ,0.0, 0.0]; #(remove)\n"
-  "	AccretionKernel            = 0;\n"
-  "	UsePhysicalUnit            = False;\n"
-  "	UseFloor                   = False;\n"
-  "	SmallRho                   = 1.0e-30;\n"
-  "	SmallP                     = 1.0e-35;\n"
-  "	SmallT                     = 1.0e-10;\n"
-  "	Coordinate                 = 0;   \n"
-  "	PoissonApproximationThreshold = 0.001;\n"
-  "	PoissonBoundaryType        = 0;\n"
-  "\n"
-  "	AngularVelocity              = 0.001;\n"
-  "	VelocityGradient             = 1.0;\n"
-  "	ShearingVelocityDirection    = -1;\n"
-  "	ShearingBoundaryDirection    = -1;\n"
-  "	ShearingBoxProblemType       = 0;\n"
-  "	UseMHD                       = 0;\n"
-  "	\n"
-  "	Hydro:\n"
-  "	{\n"
-  "		#split between ppm and zeus?\n"
-  "	\n"
-  "		CourantSafetyNumber    = 0.6;\n"
-  "		HydroMethod            = 0;   # PPM Direct Euler\n"
-  "		Gamma                  = 1.6666667;\n"
-  "\n"
-  "		PPMFlatteningParameter = False;\n"
-  "		PPMDiffusionParameter  = False;\n"
-  "		PPMSteepeningParameter = False;\n"
-  "\n"
-  "		EOSType                    = 0;\n"
-  "		EOSSoundSpeed              = 2.65e4;\n"
-  "		EOSCriticalDensity         = 1e-13;\n"
-  "		EOSGamma                   = 1.6666667;\n"
-  "		Mu                         = 0.6;\n"
-  "\n"
-  "		PressureFree                   = False;\n"
-  "\n"
-  "		UseMinimumPressureSupport        = False;\n"
-  "		MinimumPressureSupportParameter  = 100.0;\n"
-  "\n"
-  "		FluxCorrection                 = True;\n"
-  "		InterpolationMethod            = 1;  # SecondOrderA\n"
-  "		ConservativeInterpolation      = True;\n"
-  "\n"
-  "		DualEnergyFormalism         = False;\n"
-  "		DualEnergyFormalismEta1     = 0.001;\n"
-  "		DualEnergyFormalismEta2     = 0.1;\n"
-  "\n"
-  "		ZEUSLinearArtificialViscosity    = 0.0;\n"
-  "		ZEUSQuadraticArtificialViscosity = 2.0;\n"
-  "\n"
-  "		ParticleCourantSafetyNumber = 0.5;\n"
-  "		RootGridCourantSafetyNumber = 1.0;\n"
-  "\n"
-  "		UseViscosity               = False;\n"
-  "		ViscosityCoefficient       = 0.0;\n"
-  "		UseAmbipolarDiffusion      = False;\n"
-  "		UseResistivity             = False;\n"
-  "\n"
-  "		UseHydro                   = True;\n"
-  "		Theta_Limiter              = 1.5;\n"
-  "		RiemannSolver              = -99999;\n"
-  "		RiemannSolverFallback      = False;\n"
-  "		IsothermalSoundSpeed       = 1.0;\n"
-  "\n"
-  "		ConservativeReconstruction = True;\n"
-  "		PositiveReconstruction     = False;\n"
-  "		ReconstructionMethod       = -99999;\n"
-  "		RKOrder                    = 2;\n"
-  "	}; # end hydro\n"
-  "\n"
-  "	OtherParticles:\n"
-  "	{\n"
-  "		StarParticleCreation                  = 0;\n"
-  "		StarMetalYield                        = 0.02;\n"
-  "		StarMassEjectionFraction              = 0.25;\n"
-  "		StarEnergyToThermalFeedback           = 1.0e-05;\n"
-  "		StarEnergyToStellarUV                 = 3.0e-06;\n"
-  "		StarEnergyToQuasarUV                  = 5.0e-06;\n"
-  "		StellarWindFeedback                   = 0;\n"
-  "		StellarWindTurnOnMass                 = 0.1;  #(please document)\n"
-  "		MSStellarWindTurnOnMass               = 10.0; #(please document)\n"
-  "		StarParticleFeedback                  = 0;\n"
-  "		NumberOfParticleAttributes            = -99999;\n"
-  "		AddParticleAttributes                 = False;\n"
-  "		S2ParticleSize                        = 3.0;\n"
-  "\n"
-  "		StarFeedbackDistRadius                = 0;\n"
-  "		StarFeedbackDistCellStep	      = 0;\n"
-  "\n"
-  "		BigStar:\n"
-  "		{\n"
-  "			Formation                      = False;\n"
-  "			FormationDone		       = False;\n"
-  "			Separation                     = 0.25;\n"
-  "		};\n"
-  "\n"
-  "		StarMaker:\n"
-  "		{\n"
-  "			# Question: do any of the other star methods use these parameters?\n"
-  "			OverDensityThreshold         = 100.0;\n"
-  "			SHDensityThreshold           = 7.0e-26;\n"
-  "			MassEfficiency               = 1.0;\n"
-  "			MinimumMass                  = 1.0e9;\n"
-  "			MinimumDynamicalTime         = 1.0e6;\n"
-  "		};\n"
-  "\n"
-  "		PopIIIStarParticle:\n"
-  "		{\n"
-  "\n"
-  "			StarMass                        = 100.0;\n"
-  "			InitialMassFunction             = False;\n"
-  "			InitialMassFunctionSeed         = -99999;\n"
-  "			InitialMassFunctionCalls        = 0;\n"
-  "			LowerMassCutoff                 = 1.0;\n"
-  "			UpperMassCutoff                 = 300.0;\n"
-  "			MassRange                       = [1.0, 300.0]; \n"
-  "			InitialMassFunctionSlope        = -1.3;\n"
-  "			HeliumIonization                = False;\n"
-  "			BlackHoles                      = False;\n"
-  "			BHLuminosityEfficiency          = 0.1;\n"
-  "			OverDensityThreshold            = 1.0e6;\n"
-  "			H2CriticalFraction              = 5.0e-4;\n"
-  "			MetalCriticalFraction           = 1.0e-4;\n"
-  "			SupernovaRadius                 = 1.0;\n"
-  "			SupernovaUseColour              = False;\n"
-  "			SupernovaMustRefine             = False;\n"
-  "			SupernovaMustRefineResolution   = 32;\n"
-  "			ColorDensityThreshold           = 1.0e6;\n"
-  "			ColorMass                       = 1.0e6;\n"
-  "		};\n"
-  "\n"
-  "		MBHParticle:\n"
-  "		{	\n"
-  "			Accretion                          = False;\n"
-  "			AccretionRadius                    = 50.0;\n"
-  "			AccretingMassRatio                 = 1.0;\n"
-  "			AccretionFixedTemperature          = 3.0e5;\n"
-  "			AccretionFixedRate                 = 1.0e-3;\n"
-  "			TurnOffStarFormation               = False;\n"
-  "			CombineRadius                      = 50.0;\n"
-  "			MinDynamicalTime                   = 1.0e7;\n"
-  "			MinimumMass                        = 1.0e3;\n"
-  "			Feedback                           = False;\n"
-  "			FeedbackRadiativeEfficiency        = 0.1;\n"
-  "			FeedbackEnergyCoupling             = 0.05;\n"
-  "			FeedbackMassEjectionFraction       = 0.1;\n"
-  "			FeedbackMetalYield                 = 0.02;\n"
-  "			FeedbackThermalRadius              = 50.0;\n"
-  "			FeedbackJetsThresholdMass          = 10.0;\n"
-  "		};			\n"
-  "\n"
-  "		StarClusterParticle:\n"
-  "		{\n"
-  "			UseMetalField              = False;\n"
-  "			UnresolvedModel            = False;\n"
-  "			MinDynamicalTime           = 1.0e7;\n"
-  "			IonizingLuminosity         = 1.0e47;\n"
-  "			HeliumIonization           = False;\n"
-  "			SNEnergy                   = 6.8e48;\n"
-  "			SNRadius                   = 10.0;\n"
-  "			FormEfficiency             = 0.1;\n"
-  "			MinimumMass                = 1000.0;\n"
-  "			CombineRadius              = 10.0;\n"
-  "			RegionLeftEdge             = [0.0, 0.0, 0.0];\n"
-  "			RegionRightEdge            = [1.0, 1.0, 1.0];\n"
-  "		};\n"
-  "\n"
-  "		SinkParticle:\n"
-  "		{	\n"
-  "			MergeDistance              = 1.0e16;\n"
-  "			MergeMass                  = 0.1;\n"
-  "			TotalMass                  = 0.0;\n"
-  "		};\n"
-  "\n"
-  "		TracerParticle:\n"
-  "		{\n"
-  "			CreationSpacing            = -99999.0;\n"
-  "			CreationLeftEdge           = [-99999.0, -99999.0, -99999.0];\n"
-  "			CreationRightEdge          = [-99999.0, -99999.0, -99999.0];\n"
-  "		};\n"
-  "\n"
-  "	}; # end otherParticles\n"
-  "\n"
-  "	Gravity:\n"
-  "	{\n"
-  "		UniformGravity               = False;\n"
-  "		UniformGravityDirection      = 0;\n"
-  "		UniformGravityConstant       = 1.0;\n"
-  "		PointSourceGravity           = False;\n"
-  "		PointSourceGravityPosition   = [0.0, 0.0, 0.0]; \n"
-  "		PointSourceGravityConstant   = 1.0;\n"
-  "		PointSourceGravityCoreRadius = 0.0;\n"
-  "		ExternalGravity              = False;\n"
-  "		SelfGravity                  = True;\n"
-  "		SelfGravityGasOff            = False;\n"
-  "		GravitationalConstant        = 12.566370614359172; # 4*pi\n"
-  "		GravityResolution            = 1.0;\n"
-  "		ComputePotential             = False;  #(remove confirmed)\n"
-  "		PotentialIterations          = 4; #(revisit by gravity solver group)\n"
-  "		BaryonSelfGravityApproximation = True;\n"
-  "		GreensFunctionMaxNumber     = 1;\n"
-  "		GreensFunctionMaxSize       = 1;\n"
-  "		UseConstantAcceleration     = False;\n"
-  "		ConstantAcceleration        = [0.0, 0.0, 0.0];\n"
-  "	};\n"
-  "\n"
-  "\n"
-  "	AtomicPhysics:\n"
-  "	{	\n"
-  "		MetalCooling                   = False;\n"
-  "		MetalCoolingTable              = \"metal_cool.dat\";\n"
-  "		CoolingCutOffDensity1          = 0.0; \n"
-  "		CoolingCutOffDensity2          = 1e10; \n"
-  "		CoolingCutOffTemperature       = 0.0;\n"
-  "		CoolingPowerCutOffDensity1     = 0.0; \n"
-  "		CoolingPowerCutOffDensity2     = 1.0e10;\n"
-  "		UseH2OnDust                    = False;\n"
-  "		RadiativeCooling               = True;\n"
-  "		RadiativeCoolingModel          = 1;\n"
-  "		GadgetEquilibriumCooling       = 0; #(is being fixed)\n"
-  "		MultiSpecies                   = 0;\n"
-  "		PrimordialChemistrySolver      = False;\n"
-  "		CIECooling                     = True;\n"
-  "		H2OpticalDepthApproximation    = True;\n"
-  "		ThreeBodyRate                  = 0;\n"
-  "		CoolDataf0to3                  = 0.1;\n"
-  "		NumberOfTemperatureBins        = 600;\n"
-  "		CoolDataIh2co                  = True;\n"
-  "		CoolDataIpiht                  = True;\n"
-  "\n"
-  "		TemperatureStart               = 1.0;\n"
-  "		TemperatureEnd                 = 1.0e8;\n"
-  "\n"
-  "		CoolDataCompXray               = 0.0;\n"
-  "		CoolDataTempXray               = 0.0;\n"
-  "		RateDataCaseBRecombination     = False;\n"
-  "		PhotoelectricHeating           = False;\n"
-  "		HydrogenFractionByMass         = 0.76;\n"
-  "		DeuteriumToHydrogenRatio       = 6.8e-5;\n"
-  "		MultiMetals                    = 0;\n"
-  "\n"
-  "		CloudyCooling:\n"
-  "		{\n"
-  "			CloudyCoolingGridFile          = \"\";\n"
-  "			IncludeCloudyHeating           = False;\n"
-  "			IncludeCloudyMMW               = False;\n"
-  "			CMBTemperatureFloor            = True;\n"
-  "			CloudyMetallicityNormalization = 0.018477;\n"
-  "			CloudyElectronFractionFactor   = 9.153959e-3;\n"
-  "		};\n"
-  "\n"
-  "	}; # end atomicPhysics\n"
-  "\n"
-  "	RadiationField:\n"
-  "	{\n"
-  "		RadiationRedshiftOn            = 7.0;\n"
-  "		RadiationRedshiftOff           = 0.0;\n"
-  "		RadiationRedshiftFullOn        = 6.0;\n"
-  "		RadiationRedshiftDropOff       = 0.0;\n"
-  "\n"
-  "		RadiationFieldType             = 0;\n"
-  "		TabulatedLWBackground          = False;\n"
-  "		AdjustUVBackground             = True;\n"
-  "		SetUVBAmplitude                = True;\n"
-  "		SetHeIIHeatingScale            = 1.8;\n"
-  "		RadiationFieldLevelRecompute   = False;\n"
-  "		RadiationShield                = False;\n"
-  "		RadiationSpectrumNormalization = 1.0e-21;\n"
-  "		RadiationSpectrumSlope         = 1.5;\n"
-  "\n"
-  "		RadHydroParamfile = \"\";\n"
-  "	};\n"
-  "\n"
-  "\n"
-  "	Cosmology:\n"
-  "	{\n"
-  "		HubbleConstantNow   = 0.701;\n"
-  "		OmegaMatterNow      = 0.279;\n"
-  "		OmegaLambdaNow      = 0.721;\n"
-  "		ComovingBoxSize     = 64.0;\n"
-  "		MaxExpansionRate    = 0.01;\n"
-  "		InitialRedshift     = 99.0 ; \n"
-  "		FinalRedshift       = 0.0; \n"
-  "		ComovingCoordinates = True;\n"
-  "	};\n"
-  "\n"
-  "	RadiationTransfer:\n"
-  "	{\n"
-  "		RadiativeTransfer              = False;\n"
-  "		RadiationXRaySecondaryIon      = False;\n"
-  "		RadiationXRayComptonHeating    = False;\n"
-  "		RadiativeTransferRadiationPressure        = False;\n"
-  "		RadiativeTransferRadiationPressureScale   = 0.0;\n"
-  "		RadiativeTransferSourceRadius             = 0.0;\n"
-  "		RadiativeTransferPropagationSpeedFraction = 0.0;\n"
-  "		RadiativeTransferPropagationDistance      = 0.0;\n"
-  "		RadiativeTransferCoupledRateSolver        = False;\n"
-  "		RadiativeTransferOpticallyThinH2          = False;\n"
-  "		RadiativeTransferFLDCallOnLevel           = 0;\n"
-  "		RadiativeTransferPeriodicBoundary         = False;\n"
-  "		RadiativeTransferSplitPhotonRadius        = 0.0;\n"
-  "		RadiativeTransferFluxBackgroundLimit      = 0.0;\n"
-  "		RadiativeTransferRaysPerCell              = 0.0;\n"
-  "		RadiativeTransferTimestepVelocityLimit    = 0.0;\n"
-  "		RadiativeTransferInitialHEALPixLevel      = 0;\n"
-  "		RadiativeTransferPhotonEscapeRadius       = 0.0;\n"
-  "		RadiativeTransferInterpolateField         = False;\n"
-  "		RadiativeTransferSourceClustering         = False;\n"
-  "		RadiativeTransferPhotonMergeRadius        = 0.0;\n"
-  "		RadiativeTransferSourceBeamAngle          = 0.0;\n"
-  "		RadiativeTransferHIIRestrictedTimestep    = False;\n"
-  "		RadiativeTransferAdaptiveTimestep         = False;\n"
-  "		RadiativeTransferHydrogenOnly             = False;\n"
-  "		RadiativeTransferTraceSpectrum            = False;\n"
-  "		RadiativeTransferTraceSpectrumTable       = \"\";\n"
-  "		dtPhoton                                  = -99999.0;\n"
-  "\n"
-  "		ImplicitProblem                           = False;\n"
-  "		RadiativeTransferFLD                      = False;\n"
-  "\n"
-  "		SimpleQ                                   = 1.0e50;\n"
-  "		SimpleRampTime                            = 0.1;\n"
-  "	};\n"
-  "\n"
-  "	Conduction:\n"
-  "	{\n"
-  "		IsotropicConduction             = False;\n"
-  "		IsotropicConductionSpitzerFraction = 1.0;\n"
-  "		AnisotropicConduction           = False;\n"
-  "		AnisotropicConductionSpitzerFraction = 1.0;\n"
-  "		ConductionCourantSafetyNumber   = 0.5;\n"
-  "	};\n"
-  "\n"