From: jingluan <jin...@ca...> - 2010-06-29 20:55:40
|
Hello :-) I run a one-solar-mass star and print the s% R(k)/(7*10**10), (s% mlt_mixing_length(k))/(7*10**10), (s% mlt_mixing_length(k))/(s% scale_height(k)), here is a sample of the output at the beginning of the run 8.0738873839706695 0.44332401385765779 138958140368.74124 8.0718427816721459 0.44463429001342880 138957723368.84402 8.0697921319668957 0.44594916505956783 138957729700.33823 8.0677354208661853 0.44726719350617850 138957765906.98047 8.0656726249463055 0.44858946730054233 138957436889.61722 7*10**10 cm is the solar radius. 1, Why could R(k) be larger than Rsun please? 2, Why the mixing_length_alpha so large, it should be order of unity. I checked star_defaults.dek, mixing_length_alpha is set to be 2. Thank you very much :-) -- Sincerely Jing Ph.D candidate at physics.caltech email: jin...@ca... address: MC350-17,Caltech,1200E.California Blvd, Pasadena, CA 91125 |