Skip to content
This repository was archived by the owner on Aug 28, 2024. It is now read-only.

Commit

Permalink
Cleaned up formatting of swiftest_parameters definition
Browse files Browse the repository at this point in the history
  • Loading branch information
@daminton
daminton committed Dec 1, 2022
1 parent b7afb0d commit 3e791e9
Showing 1 changed file with 37 additions and 37 deletions.
74 changes: 37 additions & 37 deletions src/modules/swiftest_classes.f90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -32,43 +32,43 @@ module swiftest_classes
!> User defined parameters that are read in from the parameters input file.
!> Each paramter is initialized to a default values.
type :: swiftest_parameters
character(STRMAX) :: integrator = UNKNOWN_INTEGRATOR !! Symbolic name of the nbody integrator used
character(STRMAX) :: param_file_name = "param.in" !! The default name of the parameter input file
integer(I4B) :: maxid = -1 !! The current maximum particle id number
integer(I4B) :: maxid_collision = 0 !! The current maximum collision id number
real(DP) :: t0 = -1.0_DP !! Integration start time
real(DP) :: t = -1.0_DP !! Integration current time
real(DP) :: tstop = -1.0_DP !! Integration stop time
real(DP) :: dt = -1.0_DP !! Time step
integer(I8B) :: iloop = 0_I8B !! Main loop counter
integer(I8B) :: ioutput = 0_I8B !! Output counter
character(STRMAX) :: incbfile = CB_INFILE !! Name of input file for the central body
character(STRMAX) :: inplfile = PL_INFILE !! Name of input file for massive bodies
character(STRMAX) :: intpfile = TP_INFILE !! Name of input file for test particles
character(STRMAX) :: in_netcdf = NC_INFILE !! Name of system input file for NetCDF input
character(STRMAX) :: in_type = "ASCII" !! Data representation type of input data files
character(STRMAX) :: in_form = "XV" !! Format of input data files ("EL" or "XV")
integer(I4B) :: istep_out = -1 !! Number of time steps between saved outputs
character(STRMAX) :: outfile = NETCDF_OUTFILE !! Name of output binary file
character(STRMAX) :: out_type = "NETCDF_DOUBLE" !! Binary format of output file
character(STRMAX) :: out_form = "XVEL" !! Data to write to output file
character(STRMAX) :: out_stat = 'NEW' !! Open status for output binary file
integer(I4B) :: dump_cadence = 1 !! Number of output steps between dumping simulation data to file
real(DP) :: rmin = -1.0_DP !! Minimum heliocentric radius for test particle
real(DP) :: rmax = -1.0_DP !! Maximum heliocentric radius for test particle
real(DP) :: rmaxu = -1.0_DP !! Maximum unbound heliocentric radius for test particle
real(DP) :: qmin = -1.0_DP !! Minimum pericenter distance for test particle
character(STRMAX) :: qmin_coord = 'HELIO' !! Coordinate frame to use for qmin
real(DP) :: qmin_alo = -1.0_DP !! Minimum semimajor axis for qmin
real(DP) :: qmin_ahi = -1.0_DP !! Maximum semimajor axis for qmin
real(QP) :: MU2KG = -1.0_QP !! Converts mass units to grams
real(QP) :: TU2S = -1.0_QP !! Converts time units to seconds
real(QP) :: DU2M = -1.0_QP !! Converts distance unit to centimeters
real(DP) :: GU = -1.0_DP !! Universal gravitational constant in the system units
real(DP) :: inv_c2 = -1.0_DP !! Inverse speed of light squared in the system units
character(NAMELEN) :: interaction_loops = "ADAPTIVE" !! Method used to compute interaction loops. Options are "TRIANGULAR", "FLAT", or "ADAPTIVE"
character(NAMELEN) :: encounter_check_plpl = "ADAPTIVE" !! Method used to compute pl-pl encounter checks. Options are "TRIANGULAR", "SORTSWEEP", or "ADAPTIVE"
character(NAMELEN) :: encounter_check_pltp = "ADAPTIVE" !! Method used to compute pl-tp encounter checks. Options are "TRIANGULAR", "SORTSWEEP", or "ADAPTIVE"
character(STRMAX) :: integrator = UNKNOWN_INTEGRATOR !! Symbolic name of the nbody integrator used
character(STRMAX) :: param_file_name = "param.in" !! The default name of the parameter input file
integer(I4B) :: maxid = -1 !! The current maximum particle id number
integer(I4B) :: maxid_collision = 0 !! The current maximum collision id number
real(DP) :: t0 = -1.0_DP !! Integration start time
real(DP) :: t = -1.0_DP !! Integration current time
real(DP) :: tstop = -1.0_DP !! Integration stop time
real(DP) :: dt = -1.0_DP !! Time step
integer(I8B) :: iloop = 0_I8B !! Main loop counter
integer(I8B) :: ioutput = 0_I8B !! Output counter
character(STRMAX) :: incbfile = CB_INFILE !! Name of input file for the central body
character(STRMAX) :: inplfile = PL_INFILE !! Name of input file for massive bodies
character(STRMAX) :: intpfile = TP_INFILE !! Name of input file for test particles
character(STRMAX) :: in_netcdf = NC_INFILE !! Name of system input file for NetCDF input
character(STRMAX) :: in_type = "ASCII" !! Data representation type of input data files
character(STRMAX) :: in_form = "XV" !! Format of input data files ("EL" or "XV")
integer(I4B) :: istep_out = -1 !! Number of time steps between saved outputs
character(STRMAX) :: outfile = NETCDF_OUTFILE !! Name of output binary file
character(STRMAX) :: out_type = "NETCDF_DOUBLE" !! Binary format of output file
character(STRMAX) :: out_form = "XVEL" !! Data to write to output file
character(STRMAX) :: out_stat = 'NEW' !! Open status for output binary file
integer(I4B) :: dump_cadence = 1 !! Number of output steps between dumping simulation data to file
real(DP) :: rmin = -1.0_DP !! Minimum heliocentric radius for test particle
real(DP) :: rmax = -1.0_DP !! Maximum heliocentric radius for test particle
real(DP) :: rmaxu = -1.0_DP !! Maximum unbound heliocentric radius for test particle
real(DP) :: qmin = -1.0_DP !! Minimum pericenter distance for test particle
character(STRMAX) :: qmin_coord = 'HELIO' !! Coordinate frame to use for qmin
real(DP) :: qmin_alo = -1.0_DP !! Minimum semimajor axis for qmin
real(DP) :: qmin_ahi = -1.0_DP !! Maximum semimajor axis for qmin
real(QP) :: MU2KG = -1.0_QP !! Converts mass units to grams
real(QP) :: TU2S = -1.0_QP !! Converts time units to seconds
real(QP) :: DU2M = -1.0_QP !! Converts distance unit to centimeters
real(DP) :: GU = -1.0_DP !! Universal gravitational constant in the system units
real(DP) :: inv_c2 = -1.0_DP !! Inverse speed of light squared in the system units
character(NAMELEN) :: interaction_loops = "ADAPTIVE" !! Method used to compute interaction loops. Options are "TRIANGULAR", "FLAT", or "ADAPTIVE"
character(NAMELEN) :: encounter_check_plpl = "ADAPTIVE" !! Method used to compute pl-pl encounter checks. Options are "TRIANGULAR", "SORTSWEEP", or "ADAPTIVE"
character(NAMELEN) :: encounter_check_pltp = "ADAPTIVE" !! Method used to compute pl-tp encounter checks. Options are "TRIANGULAR", "SORTSWEEP", or "ADAPTIVE"

! The following are used internally, and are not set by the user, but instead are determined by the input value of INTERACTION_LOOPS
logical :: lflatten_interactions = .false. !! Use the flattened upper triangular matrix for pl-pl interaction loops
Expand Down

0 comments on commit 3e791e9

Please sign in to comment.