Binary simulation

I don't quite understand the meaning of t1,t2 and radiusw in the binary module. Does the radiusw refer to where the wind is emitted in the unit of base cell? That is, if the sim box is 200 AU and I have 256 base grid, then each cell has 0.8 AU and if radiusw=12 then the wind emit at 10AU?

Attachments (2)

Download all attachments as: .zip

Comments

1. Jonathan -- 11 years ago

So t1 and t2 are times used to calculate 'beta' which is used to modify the value of pointgravity%alpha of the secondary particle. Alpha is a parameter that scales the strength of gravity. Look in source/pointgravity.f90 to see exactly how it is used. And radiusw is multiplied by 'sink_dx' which is the resolution of the finest level. So radiusw is in finest level computational zones - so your statement would be true for MaxLevel=0, but if MaxLevel = 1, then the wind would emit at 5 AU etc…

Also - in looking over the module - it appears that several values are read in from the data file at every time step for every patch… This could hurt performance considerably. All variables should be read in during ProblemModuleInit. If you need to use the variables in BeforeStep, then declare the variables global to the module by placing their declarations before the first 'Contains' statement and remove their declarations within the separate subroutines.

2. Zhuo Chen -- 11 years ago

Thank you very much for your explanation and suggestion. I will modify that part when I start to modify the module.