ivar_refine and the geometrical criterion in star_formation routine

26 views
Skip to first unread message

prachikhatri95

unread,
Jan 18, 2022, 3:42:14 AM1/18/22
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
ramses_lvl8.nml
Reply all
Reply to author
Forward
0 new messages