radius definition in the yt sunrise_exporter.py

[mohammad safarzadeh]

I see this line

in this def :

def prepare_star_particles(pf,star_type,pos=None,vel=None, age=None,
                          creation_time=None,initial_mass=None,
                          current_mass=None,metallicity=None,
                          radius = None,
                          fle=[0.,0.,0.],fre=[1.,1.,1.]):

    if radius is None:
        radius = initial_mass*0.0+10.0/1000.0 #10pc radius

Why is it a constant ? doesn't that influence the SUNRISE run as the radius is important there? Is it taken as typical radius of stellar cluster?

Or should one change this line for sure for each simulation if we know the radius of particles depending on their mass etc ? 

[Christopher Moody]

Hi Mohammad, 

Since this is part of the yt codebase, this might be more appropriate on the yt-user or yt-dev lists. 

Sunrise spawns photon packets within a stellar 'radius' instead of from a single point, so you need to give Sunrise something. It's set to a small number, 10pc, but I'd suggest experimenting with star_radius_factor in the sfrhist config to boost this number up and down instead of re-running the exporter and creating a new FITS file every time. I'd recommend you vary the radius around whatever length scale your simulation resolves. As you note, you can change that line and experiment with modifying the radius as a function of the star particle's mass. 

--
You received this message because you are subscribed to the Google Groups "Sunrise" group.
To unsubscribe from this group and stop receiving emails from it, send an email to sunrisemcrx...@googlegroups.com.
To post to this group, send email to sunri...@googlegroups.com.
Visit this group at http://groups.google.com/group/sunrisemcrx?hl=en.
For more options, visit https://groups.google.com/groups/opt_out.