lwang-astro / PeTar

PeTar is a high-performance N-body code for modelling the evolution of star clusters and tidal streams, including the effect of galactic potential, dynamics of binary and hierarchical system, single and binary stellar evolution.
MIT License
71 stars 19 forks source link

Running petar.find.dt with bse enabled #21

Closed JoRdAnB421 closed 2 years ago

JoRdAnB421 commented 2 years ago

I am having some problems understanding the error from the check.perf.test.log file when I run the petar.find.dt tool on my cluster (which I created from mcluster). This is the command that I submitted:

petar.find.dt -o 8 -a "--stellar-evolution 1 --detect-interrupt" ../mcluster/test_10000_run.txt.input

and this fails after a single time step with the error message in check.perf.test.log as follows:

image

I am not sure what I am doing wrong here and so any assistance is much appreciated.

lwang-astro commented 2 years ago

The option used in petar.find.dt is incorrect. --detect-interrupt requires one argument. But you don't need to provide '--stellar-evolution' and '--detect-interrupt' if the default values are used.

JoRdAnB421 commented 2 years ago

Okay I see, thank you for the clarification! Can I check what the default values are for these two parameters?

JoRdAnB421 commented 2 years ago

As an extra comment, I have added the parameter required for --detect-interrupt so that the command line reads:

petar.find.dt -m 2 -o 8 -a "--stellar-evolution 1 --detect-interrupt 1" ../mcluster/test_10000_run.txt.input

This resolves the error I was facing before, however presents a new issue in the log files:

image

Thank you very much for your help.

lwang-astro commented 2 years ago

Which fortran compiler did you use? For gfortran, I did not get this error.

JoRdAnB421 commented 2 years ago

As far as I am aware I am also using gfortran.

lwang-astro commented 2 years ago

Are you using Mac OS?

lwang-astro commented 2 years ago

I did not got this error in Linux. One simple solution is comment out the line 23-25

23      do K=0,14
24         write(*,*) ktype(k,:)
25      end do

in bse-interface/bse/print_const.f

JoRdAnB421 commented 2 years ago

I am running this on a computing cluster which uses Linux, I will try commenting out those lines and see if that helps.

Thanks for your help