prachikhatri95
unread,Jan 18, 2022, 3:42:14 AM1/18/22Sign in to reply to author
Sign in to forward
You do not have permission to delete messages in this group
Either email addresses are anonymous for this group or you need the view member email addresses permission to view the original message
to ramses_users
Hello everyone,
I have a question about the following geometrical criterion in the star_formation routine
-------------------------------------------------------
! Geometrical criterion
if(ivar_refine>0)then
do i=1,ngrid
d=uold(ind_cell(i),ivar_refine)
if(d<=var_cut_refine)ok(i)=.false.
end do
endif
-------------------------------------------------------
As the default value for ivar_refine is -1, do I understand correctly that this criterion is not always needed? In which simulations is it required and if I start a cosmological hydro simulation with unigrid ICs (from MUSIC), should I use the default value?
My issue is that while generating unigrid ICs with MUSIC, the namelist by MUSIC suggests a positive value of ivar_refine but, it does not generate any initial conditions file with the name 'pvar_0000x' (for refinement). So, when I use a positive value of ivar_refine in my run, I do not form any stars because by default, the array corresponding to the refinement variable takes value 0.0 (in my simulation) and the geometrical criterion for star formation is never fulfilled.
I have attached my namelist below in case there is some other mistake that is causing this problem. I compile ramses with 9 variables as I use metals, delayed cooling, and another hydro variable that I have hardcoded to take the 9th uold array.
Thank you in advance
Best regards,
Prachi