|
|
so now I've built SWIFT and tested it with some of the examples scripts.
|
|
|
#Created cooling.c and cooling.h
|
|
|
|
|
|
Can we new include hydroparticles? yes - Tom's already got an example.
|
|
|
Followed example of external potentials. Struct contains cooling data, have a function which reads
|
|
|
parameters into the struct.
|
|
|
|
|
|
cp -r /gpfs/data/tt/Codes/Swift-git/swiftsim.July21/examples/IsothermalPotential//HydroStatic/ .
|
|
|
Wrote function which takes the 'old particle data', the cooling data, and a timestep as arguments, and returns a new internal energy. This function would have to be different depending on the cooling model used. For now we just have a constant cooling rate for all particles.
|
|
|
|
|
|
and try it :) compare to rgb's attempt
|
|
|
Another functions takes the 'new internal energy' and updates the entropy.
|
|
|
|
|
|
/gpfs/data/rgb/Swift/swiftsim/examples/IsothermalPotential/HydroStatic
|
|
|
#Added the function runner_do_cooling to runner.c
|
|
|
|
|
|
Stefan:
|
|
|
Followed the 'external potentials' example. Calls the 'update entropy' unction for a particle when it is active
|
|
|
|
|
|
GravityOnly works fine
|
|
|
|
|
|
Trying Tom's example:
|
|
|
|
|
|
Need to add timestep_mult to the yml file (set to be 0.03 in the GravityOnly example)
|
|
|
|
|
|
Need to add Statistics:delta_time
|
|
|
|
|
|
Got a segmentation fault when I ran it.
|
|
|
|
|
|
Trying Richard's example:
|
|
|
|
|
|
Seems to run fine
|
|
|
#Added cooling to the engine
|
|
|
|
|
|
Added'engine_policy_cooling' to the list of engine policies.
|
|
|
Added a 'cooling task' to the cell struct
|
|
|
The cooling task is unlocked after the kick is performed. It does not unlock anything
|
|
|
|
|
|
#Added a 'C' flag option. |
|
|
\ No newline at end of file |