User Tools

Site Tools


isis:slxfig:hr_diagram

Hertzsprung-Russell diagram

% -------------
% axes, labels:
variable xfig = xfig_plot_new(12,9);
variable xmin = log10(3), xmax = log10(47), ymin = -1., ymax = 7.; % x-axis is in kilo Kelvin
xfig.world(xmax, xmin, ymin, ymax);
xfig.axis(; ticlabels2=0, ticlabels_confine=0);
variable xmajor = log10([5:45:5]);
variable xminor = log10([1:50:1]);
variable x2major = [ymin:ymax:0.1];
variable x2minor = [ymin:ymax:0.02];
xfig.axis(; depth=0, tic_depth=0, ticlabels_confine=0);
xfig.x1axis(; ticlabels=["O","B","A","F","G","K","M"], major=log10([32,18,10,7,5.5,4.25,3.25]), major_width=0);
xfig.xlabel("Spectral class");
xfig.x2axis(; ticlabels=["5","10","15","20","25","30","35","40","45"R], major=xmajor, minor=xminor, format="%.0f");
xfig.x2label("$T_{\mathrm{eff}}\,(10^3\,\mathrm{K})$"R);
xfig.y2axis(; ticlabels=0);
xfig.ylabel("$\log (L/{L}_{\odot})$"R);
% -------------
% --------------------
% evolutionary tracks:
% from Ekstroem et al. 2012, A&A, 537, A146
variable tracks = struct{ m1, m3, m7, m20 };
tracks.m1 = struct{ log_teff, log_l }; % track for 1 solar mass
tracks.m3 = struct{ log_teff, log_l }; % track for 3 solar mass
tracks.m7 = struct{ log_teff, log_l }; % track for 7 solar mass
tracks.m20 = struct{ log_teff, log_l }; % track for 20 solar mass
tracks.m1.log_teff = [3.754330,3.755300,3.755330,3.755400,3.755480,3.755580,3.755700,3.755810,3.755940,3.756070,3.756190,3.756330,3.756470,3.756610,3.756730,3.756860,3.757000,3.757150,3.757270,3.757420,3.757570,3.757700,3.757850,
		      3.757980,3.758140,3.758270,3.758420,3.758580,3.758710,3.758870,3.759000,3.759140,3.759290,3.759430,3.759610,3.759770,3.759930,3.760120,3.760280,3.760440,3.760600,3.760780,3.760940,3.761100,3.761250,3.761410,
		      3.761560,3.761670,3.761800,3.761950,3.762060,3.762180,3.762300,3.762410,3.762520,3.762620,3.762750,3.762820,3.762880,3.762950,3.763020,3.763080,3.763140,3.763190,3.763240,3.763290,3.763330,3.763370,3.763400,
		      3.763430,3.763450,3.763470,3.763480,3.763490,3.763480,3.763470,3.763450,3.763410,3.763360,3.763290,3.763220,3.763140,3.763040,3.762940,3.762850,3.762800,3.762770,3.762750,3.762730,3.762680,3.762650,3.762630,
		      3.762610,3.762560,3.762530,3.762480,3.762460,3.762410,3.762390,3.762350,3.762300,3.762280,3.762230,3.762180,3.762120,3.762070,3.762010,3.761920,3.761860,3.761740,3.761640,3.761580,3.761520,3.761450,3.761440,
		      3.761400,3.761350,3.761300,3.761230,3.761170,3.761100,3.761030,3.760930,3.760850,3.760760,3.760650,3.760490,3.760360,3.760230,3.760080,3.759870,3.759710,3.759530,3.759340,3.759070,3.758850,3.758620,3.758370,
		      3.758020,3.757730,3.757430,3.757110,3.756620,3.756230,3.755810,3.755370,3.754740,3.754240,3.753690,3.753120,3.752290,3.751620,3.750900,3.750130,3.749010,3.748020,3.746900,3.745720,3.743960,3.742490,3.740860,
		      3.739060,3.737080,3.734070,3.731440,3.728450,3.724750,3.718820,3.713870,3.709260,3.703280,3.697400,3.692710,3.688030,3.684480,3.682560,3.680620,3.678870,3.677430,3.676030,3.674470,3.672500,3.670030,3.667120,
		      3.663040,3.656530,3.647770,3.646700,3.626800,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,
		      3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,
		      3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,
		      3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,
		      3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,
		      3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,
		      3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,
		      3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,3.556410,
		      3.556410,3.556410];
tracks.m1.log_l = [-0.138050,-0.129380,-0.126810,-0.124360,-0.122260,-0.119770,-0.117260,-0.114720,-0.112130,-0.109510,-0.107250,-0.104540,-0.101830,-0.099090,-0.096690,-0.093890,-0.091050,-0.088180,-0.085690,-0.082760,-0.079800,
		   -0.077230,-0.074220,-0.071610,-0.068540,-0.065890,-0.062750,-0.059600,-0.056880,-0.053670,-0.050910,-0.048120,-0.044830,-0.042000,-0.038660,-0.035760,-0.032850,-0.029430,-0.026470,-0.023490,-0.020480,-0.016960,
		   -0.013920,-0.010850,-0.007770,-0.004660,-0.001540,0.001590,0.004760,0.008470,0.011680,0.014910,0.018160,0.020890,0.024180,0.027500,0.030830,0.034170,0.037550,0.040950,0.044380,0.047800,0.051270,0.054190,0.057710,
		   0.061270,0.064840,0.068450,0.072080,0.075120,0.078800,0.082480,0.086200,0.089300,0.093030,0.096140,0.099880,0.102980,0.106070,0.109160,0.112220,0.115260,0.118280,0.121310,0.123730,0.124950,0.125560,0.126180,
		   0.126790,0.128040,0.128660,0.129290,0.129920,0.131190,0.131840,0.133140,0.133800,0.135130,0.135810,0.137170,0.138550,0.139250,0.140630,0.142070,0.144260,0.145740,0.148010,0.151080,0.153420,0.157380,0.160590,
		   0.163010,0.165420,0.167860,0.171150,0.173650,0.176170,0.178710,0.182130,0.184720,0.187330,0.189970,0.193520,0.196210,0.198920,0.201650,0.205330,0.208130,0.210940,0.213780,0.217610,0.220520,0.223440,0.226400,
		   0.230380,0.233390,0.236430,0.239500,0.243630,0.246750,0.249900,0.253070,0.257320,0.260540,0.263770,0.267020,0.271370,0.274650,0.277930,0.281210,0.285580,0.288840,0.292080,0.295300,0.299510,0.302540,0.305550,
		   0.308480,0.312160,0.314740,0.317060,0.319060,0.320830,0.322500,0.322870,0.322620,0.321470,0.318760,0.316620,0.315640,0.317440,0.324700,0.337600,0.362110,0.393770,0.418910,0.461060,0.504440,0.550130,0.604380,
		   0.664270,0.729320,0.807500,0.897740,1.007900,1.152030,1.342320,1.390360,1.719880,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,
		   2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,
		   2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,
		   2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,
		   2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,
		   2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,
		   2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,
		   2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,2.638140,
		   2.638140,2.638140,2.638140,2.638140,2.638140,2.638140];
tracks.m3.log_teff = [4.104477,4.101130,4.098784,4.096932,4.095520,4.094152,4.092906,4.091568,4.090532,4.089462,4.088575,4.087603,4.086631,4.085654,4.084584,4.083618,4.082794,4.081863,4.080880,4.079909,4.078954,4.077811,4.076814,
		      4.075977,4.074953,4.073932,4.072905,4.071846,4.070604,4.069515,4.068425,4.067474,4.066353,4.065180,4.063810,4.062611,4.061387,4.060169,4.058875,4.057809,4.056518,4.054984,4.053684,4.052311,4.050895,4.049511,
		      4.048311,4.046844,4.045140,4.043617,4.042110,4.040580,4.038967,4.037586,4.035944,4.034328,4.032409,4.030733,4.028975,4.027469,4.025674,4.023859,4.022054,4.019996,4.018106,4.016439,4.014562,4.012628,4.010719,
		      4.008746,4.006818,4.004866,4.002865,4.000890,3.998964,3.997042,3.995170,3.993303,3.991514,3.989837,3.988263,3.986853,3.985639,3.984805,3.984479,3.984499,3.984547,3.984648,3.984781,3.984972,3.985252,3.985549,
		      3.985906,3.986334,3.986842,3.987464,3.988192,3.989086,3.990155,3.991399,3.992620,3.994326,3.996373,3.998494,4.001379,4.005055,4.009778,4.015627,4.024013,4.035578,4.027620,4.019198,4.014199,4.010884,4.008090,
		      4.005473,4.003593,4.001870,3.999673,3.998636,3.996650,3.995367,3.992839,3.991586,3.990325,3.988684,3.987023,3.985322,3.983579,3.981786,3.981786,3.979712,3.977541,3.975261,3.972866,3.970342,3.970342,3.967673,
		      3.964838,3.961813,3.958571,3.955076,3.955076,3.952064,3.948853,3.941737,3.937805,3.933549,3.928990,3.924054,3.918705,3.912862,3.907195,3.901103,3.894541,3.887346,3.879446,3.861010,3.850280,3.838353,3.824996,
		      3.807248,3.797427,3.777874,3.762034,3.748710,3.736734,3.724115,3.715732,3.707236,3.696032,3.686913,3.681703,3.677099,3.672587,3.668592,3.665801,3.663002,3.659846,3.657180,3.654128,3.651128,3.648095,3.645231,
		      3.642391,3.639711,3.637484,3.634593,3.632145,3.635666,3.641828,3.643987,3.646684,3.649874,3.652920,3.657066,3.660748,3.664139,3.667093,3.671135,3.674116,3.678197,3.680127,3.682554,3.684970,3.687178,3.688573,
		      3.690152,3.693206,3.695129,3.696329,3.697237,3.697691,3.697918,3.697847,3.697883,3.697758,3.697609,3.697464,3.697333,3.697136,3.697183,3.697142,3.697143,3.697211,3.697227,3.697261,3.697284,3.697322,3.697328,
		      3.697378,3.697393,3.697317,3.697438,3.697463,3.697494,3.697490,3.697448,3.697510,3.697517,3.697539,3.697554,3.697560,3.697585,3.697587,3.697635,3.697641,3.697657,3.697695,3.697711,3.697749,3.697767,3.697795,
		      3.697817,3.697833,3.697856,3.697868,3.697910,3.697926,3.697937,3.697953,3.697917,3.697939,3.697946,3.698006,3.697967,3.697968,3.697975,3.697975,3.697966,3.697966,3.698002,3.698044,3.698057,3.698066,3.698132,
		      3.698023,3.698025,3.697961,3.698303,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,
		      3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,
		      3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,3.698114,
		      3.698114,3.698114,3.698114,3.698114,3.698114,3.697860,3.697296,3.697150,3.696984,3.696565,3.696449,3.696182,3.695561,3.695231,3.694880,3.694011,3.693476,3.692530,3.691571,3.690567,3.688184,3.686590,3.684183,
		      3.680467,3.664628];
tracks.m3.log_l = [1.898976,1.898157,1.900368,1.903228,1.906229,1.909473,1.912669,1.915631,1.919331,1.922666,1.926119,1.929846,1.933160,1.936508,1.940077,1.943397,1.946812,1.950456,1.953744,1.956987,1.960515,1.963574,1.966817,1.970432,
		   1.973557,1.976631,1.980111,1.983128,1.986004,1.989179,1.992273,1.995433,1.998425,2.001378,2.004035,2.006924,2.009910,2.012560,2.015303,2.018372,2.020963,2.023431,2.026095,2.028576,2.031044,2.033705,2.036241,2.038570,
		   2.040540,2.043114,2.045241,2.047596,2.049719,2.051901,2.054142,2.056032,2.057709,2.059767,2.061599,2.063530,2.065377,2.067040,2.068619,2.070105,2.071634,2.073346,2.074828,2.076166,2.077634,2.078898,2.080094,2.081460,
		   2.082636,2.083775,2.085083,2.086236,2.087446,2.088730,2.089970,2.091433,2.092853,2.094588,2.096468,2.098816,2.101593,2.102334,2.102991,2.103922,2.104768,2.105640,2.106693,2.107692,2.108753,2.109920,2.111194,2.112641,
		   2.114297,2.116148,2.118172,2.120395,2.122534,2.125387,2.128790,2.132439,2.137256,2.142846,2.150769,2.160729,2.175909,2.197592,2.175125,2.185301,2.197447,2.202487,2.205639,2.208038,2.209410,2.210467,2.211650,2.212162,
		   2.213094,2.213552,2.214353,2.214709,2.215033,2.215437,2.215775,2.216047,2.216275,2.216510,2.216510,2.216610,2.216638,2.216595,2.216464,2.216231,2.216231,2.215881,2.215414,2.214775,2.213947,2.212979,2.212979,2.211950,
		   2.210762,2.207781,2.205965,2.203840,2.201414,2.198580,2.195333,2.191549,2.187744,2.183283,2.178384,2.172654,2.165949,2.149017,2.138421,2.126120,2.111449,2.090319,2.077763,2.047863,2.010163,1.962986,1.910611,1.839891,
		   1.803651,1.793080,1.821048,1.881206,1.940541,1.987798,2.033413,2.077010,2.113085,2.150113,2.187340,2.223158,2.261067,2.296868,2.333229,2.370631,2.407406,2.441670,2.479107,2.513618,2.542172,2.501533,2.418088,2.390079,
		   2.357468,2.324004,2.289329,2.246575,2.195523,2.155987,2.112991,2.068237,2.034009,1.994281,1.971123,1.949987,1.930177,1.912312,1.904220,1.896871,1.889495,1.888052,1.888889,1.890312,1.891323,1.892204,1.892774,1.893772,
		   1.894377,1.895038,1.895777,1.896499,1.897036,1.898205,1.899111,1.900117,1.901103,1.902077,1.902980,1.904079,1.904962,1.905917,1.906815,1.907734,1.908249,1.909652,1.910666,1.911600,1.912941,1.913630,1.914763,1.915684,
		   1.916661,1.917600,1.918530,1.919604,1.920576,1.921636,1.922639,1.923660,1.925384,1.926319,1.927370,1.928294,1.929379,1.930928,1.932756,1.933701,1.934644,1.938850,1.939751,1.940648,1.941545,1.942451,1.943542,1.944417,
		   1.945651,1.946251,1.947151,1.948034,1.948879,1.949817,1.950968,1.952394,1.953985,1.955612,1.957425,1.961586,1.962245,1.963240,1.964008,1.978211,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,
		   1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,
		   1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,
		   1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.977363,1.984200,1.993446,1.997043,2.000126,2.007648,2.010280,2.012891,2.019092,2.020792,2.024098,
		   2.028754,2.030630,2.034344,2.034319,2.035530,2.044368,2.048004,2.053374,2.067782,2.196314];
tracks.m7.log_teff = [4.332096,4.328320,4.326038,4.324511,4.323259,4.322116,4.321212,4.320357,4.319573,4.318899,4.318188,4.317513,4.316760,4.316124,4.315453,4.314796,4.314137,4.313486,4.312861,4.312146,4.311447,4.310766,4.310062,
		      4.309355,4.308616,4.307887,4.307114,4.306342,4.305527,4.304748,4.303917,4.303072,4.302205,4.301346,4.300464,4.299560,4.298623,4.297675,4.296707,4.295747,4.294726,4.293709,4.292675,4.291631,4.290535,4.289442,
		      4.288291,4.287137,4.285987,4.284830,4.283688,4.282377,4.281137,4.279848,4.278558,4.277258,4.275919,4.274649,4.273265,4.271796,4.270365,4.268891,4.267504,4.266029,4.264506,4.262972,4.261453,4.259881,4.258287,
		      4.256736,4.255117,4.253564,4.251932,4.250330,4.248753,4.247167,4.245650,4.244133,4.242741,4.241340,4.240085,4.239005,4.238083,4.237530,4.237449,4.237502,4.237544,4.237651,4.237784,4.237949,4.238114,4.238347,
		      4.238624,4.238950,4.239345,4.239774,4.240287,4.240864,4.241547,4.242350,4.243391,4.244668,4.245955,4.247738,4.249902,4.252660,4.256188,4.261033,4.268576,4.278758,4.276700,4.271091,4.264163,4.257535,4.251551,
		      4.246792,4.242080,4.238544,4.233549,4.230376,4.227304,4.222256,4.218787,4.215385,4.213696,4.210338,4.206665,4.202654,4.198614,4.196575,4.192425,4.188184,4.186028,4.181635,4.177107,4.174761,4.169952,4.167482,
		      4.162382,4.157073,4.154295,4.148542,4.142472,4.139310,4.132646,4.125541,4.121834,4.114001,4.105536,4.101063,4.091519,4.081102,4.075528,4.063547,4.050320,4.043161,4.030129,4.018452,3.998943,3.984378,3.968375,
		      3.943353,3.921980,3.891749,3.865308,3.826374,3.779667,3.725678,3.681198,3.656651,3.641506,3.627405,3.613375,3.607156,3.604405,3.602538,3.601806,3.599727,3.599603,3.600013,3.599306,3.598872,3.597960,3.597543,
		      3.597011,3.596631,3.596427,3.596191,3.595973,3.595834,3.597160,3.598375,3.599965,3.600833,3.600831,3.601515,3.602871,3.605483,3.606549,3.605555,3.607951,3.610090,3.611943,3.613236,3.614876,3.616421,3.617434,
		      3.618779,3.621073,3.622618,3.622776,3.623235,3.623353,3.623821,3.623997,3.624226,3.625528,3.625735,3.625954,3.626158,3.626497,3.626846,3.627171,3.627456,3.627730,3.627993,3.628300,3.628602,3.628912,3.629263,
		      3.630437,3.630933,3.631360,3.631842,3.632390,3.633723,3.634783,3.665681,3.683321,3.703236,3.723511,3.739929,3.752854,3.762691,3.771399,3.778577,3.785004,3.790939,3.796565,3.802501,3.807311,3.811833,3.815278,
		      3.818532,3.876135,3.877868,3.879691,3.881684,3.883706,3.886060,3.887932,3.890105,3.892588,3.894283,3.895841,3.897358,3.898810,3.900199,3.901722,3.903171,3.904857,3.906437,3.908123,3.909673,3.912379,3.913959,
		      3.915249,3.916732,3.918209,3.966178,3.966115,3.965955,3.965722,3.965449,3.965165,3.964824,3.964451,3.964058,3.963271,3.963075,3.962749,3.962339,3.961805,3.961155,3.960438,3.959722,3.959019,3.958136,3.957191,
		      3.956146,3.955251,3.954283,3.953136,3.951916,3.950582,3.949186,3.947747,3.946258,3.944704,3.943043,3.941355,3.939590,3.937691,3.935705,3.933664,3.931502,3.929228,3.926949,3.924556,3.922072,3.919503,3.916736,
		      3.913864,3.910855,3.907689,3.904399,3.900891,3.897316,3.893520,3.889467,3.885164,3.880593,3.875572,3.870274,3.864655,3.858520,3.851880,3.844598,3.836846,3.827031,3.816640,3.804874,3.791261,3.774941,3.754821,
		      3.729229,3.695042,3.659592,3.643731,3.627150,3.626133,3.625057,3.624490,3.623345,3.622775,3.621631,3.621054,3.619903,3.619290,3.618071,3.617455,3.616156,3.614861,3.614229,3.612866,3.612142,3.610937,3.610142,
		      3.610602,3.602978];
tracks.m7.log_l = [3.220316,3.222236,3.225934,3.229981,3.234307,3.238596,3.242995,3.247321,3.251924,3.256357,3.260792,3.265291,3.269632,3.274029,3.278296,3.282821,3.286990,3.291400,3.295676,3.299967,3.304149,3.308462,3.312527,
		   3.316832,3.320878,3.325039,3.329032,3.333184,3.337169,3.341168,3.345029,3.349119,3.352956,3.356905,3.360594,3.364509,3.368194,3.372097,3.375746,3.379510,3.383115,3.386828,3.390284,3.393923,3.397403,3.400997,
		   3.404437,3.407898,3.411295,3.414547,3.418016,3.421233,3.424598,3.427808,3.431126,3.434182,3.437394,3.440545,3.443572,3.446675,3.449643,3.452622,3.455812,3.458679,3.461540,3.464590,3.467340,3.470130,3.472974,
		   3.475755,3.478522,3.481117,3.484004,3.486664,3.489290,3.491944,3.494787,3.497446,3.500011,3.502895,3.505690,3.508442,3.511563,3.514660,3.517971,3.518808,3.519594,3.520310,3.521011,3.521737,3.522363,3.523141,
		   3.523951,3.524827,3.525737,3.526626,3.527718,3.528864,3.529965,3.531167,3.532609,3.534299,3.535914,3.538223,3.540880,3.544190,3.548222,3.554541,3.565152,3.583014,3.576169,3.563020,3.554800,3.556295,3.564870,
		   3.573734,3.582066,3.587452,3.593816,3.597135,3.599927,3.603747,3.605914,3.607810,3.608631,3.610085,3.611448,3.612663,3.613631,3.614050,3.614736,3.615239,3.615435,3.615717,3.615799,3.615779,3.615642,3.615516,
		   3.615084,3.614467,3.614084,3.613218,3.612078,3.611389,3.609764,3.607780,3.606682,3.604167,3.601164,3.599452,3.595528,3.590906,3.588274,3.582331,3.575291,3.571276,3.563641,3.556470,3.543681,3.533634,3.522214,
		   3.503631,3.486943,3.462564,3.440832,3.408274,3.368361,3.318280,3.239048,3.170775,3.203994,3.295598,3.433337,3.496761,3.530118,3.551234,3.561025,3.575027,3.587653,3.598804,3.611776,3.618321,3.631116,3.636021,
		   3.645658,3.653961,3.659961,3.663939,3.667343,3.669468,3.652996,3.633136,3.615651,3.603901,3.591752,3.576569,3.562458,3.550026,3.537842,3.523101,3.509558,3.495461,3.483123,3.474772,3.467494,3.461328,3.455866,
		   3.452826,3.448147,3.447262,3.447461,3.447566,3.447692,3.448052,3.448105,3.448606,3.450289,3.450684,3.451125,3.451498,3.452573,3.453145,3.453837,3.454481,3.455146,3.455889,3.456769,3.457614,3.458545,3.459693,
		   3.463869,3.465136,3.466698,3.468591,3.470875,3.477056,3.480336,3.567194,3.579202,3.588661,3.595847,3.601208,3.605330,3.608748,3.611829,3.614628,3.617675,3.620322,3.623573,3.626370,3.628592,3.631104,3.633412,
		   3.634879,3.667375,3.668497,3.669818,3.671017,3.672650,3.673945,3.675235,3.676565,3.678166,3.679165,3.680143,3.681254,3.682227,3.683366,3.684341,3.685334,3.686696,3.687644,3.689088,3.689945,3.692384,3.693152,
		   3.693934,3.695269,3.696366,3.726207,3.726218,3.726268,3.726353,3.726449,3.726525,3.726563,3.726617,3.726728,3.727225,3.727205,3.727169,3.727039,3.726874,3.726738,3.726605,3.726719,3.726471,3.726221,3.725983,
		   3.725856,3.725869,3.725482,3.725105,3.724742,3.724375,3.723998,3.723633,3.723204,3.722754,3.722301,3.721821,3.721268,3.720674,3.719999,3.719381,3.718731,3.717999,3.717337,3.716520,3.715740,3.714836,3.713875,
		   3.712960,3.711882,3.710743,3.709566,3.708307,3.707170,3.705678,3.704111,3.702612,3.700703,3.698742,3.696672,3.694619,3.692034,3.689234,3.686435,3.682163,3.677869,3.673341,3.667844,3.661392,3.653430,3.643164,
		   3.629464,3.608824,3.573459,3.518706,3.465240,3.465535,3.466360,3.466996,3.468742,3.469831,3.472490,3.474027,3.477530,3.480029,3.484467,3.486669,3.492091,3.497817,3.500731,3.507212,3.510760,3.518050,3.523789,
		   3.538654,3.590785];
tracks.m20.log_teff = [4.553253,4.549346,4.547174,4.545685,4.544487,4.543546,4.542704,4.541930,4.541199,4.540501,4.539833,4.539180,4.538523,4.537867,4.537208,4.536548,4.535864,4.535171,4.534479,4.533819,4.533099,4.532349,4.531685,
		       4.530858,4.530091,4.529258,4.528381,4.527510,4.526607,4.525624,4.524617,4.523586,4.522501,4.521398,4.520290,4.519143,4.517971,4.516772,4.515485,4.514221,4.512837,4.511507,4.510101,4.508667,4.507198,4.505676,
		       4.504136,4.502523,4.500856,4.499166,4.497491,4.495733,4.493921,4.492090,4.490113,4.488147,4.486238,4.484147,4.482163,4.479980,4.477843,4.475607,4.473299,4.471034,4.468654,4.466329,4.463757,4.461328,4.458638,
		       4.456110,4.453542,4.450639,4.447997,4.445199,4.442406,4.439610,4.436652,4.433856,4.430981,4.428301,4.425631,4.423303,4.421079,4.419665,4.417133,4.417171,4.417249,4.417350,4.417494,4.417629,4.417839,4.418064,
		       4.418213,4.418508,4.418812,4.419100,4.419612,4.420245,4.420835,4.421696,4.422720,4.423804,4.424998,4.426548,4.428610,4.431079,4.434245,4.438862,4.446238,4.466819,4.466818,4.466795,4.466869,4.466593,4.465954,
		       4.464978,4.463753,4.460943,4.459266,4.457600,4.455811,4.454226,4.452252,4.450122,4.447809,4.445428,4.440841,4.438668,4.436202,4.433690,4.429020,4.426507,4.423976,4.421396,4.416259,4.413662,4.411091,4.408434,
		       4.405806,4.400451,4.397710,4.394906,4.392073,4.386325,4.383410,4.380459,4.377478,4.371463,4.368414,4.365331,4.362171,4.359369,4.353332,4.350270,4.347179,4.344057,4.337768,4.334598,4.331373,4.328118,4.321521,
		       4.318163,4.314773,4.311343,4.304390,4.300861,4.297283,4.293637,4.289962,4.282533,4.278783,4.274971,4.271118,4.263665,4.261759,4.255988,4.252112,4.246432,4.243446,4.239031,4.233549,4.230085,4.224989,4.220862,
		       4.216878,4.212944,4.209174,4.205576,4.202151,4.198862,4.110521,4.089470,4.074724,4.062205,4.049360,4.034826,4.016791,3.997607,3.977317,3.957115,3.928374,3.903048,3.876210,3.850041,3.808092,3.734719,3.576925,
		       3.571743,3.571057,3.570602,3.570459,3.570285,3.570244,3.570004,3.570037,3.570028,3.569749,3.569502,3.569621,3.569551,3.569749,3.569678,3.569698,3.568869,3.568953,3.568850,3.568921,3.568871,3.568696,3.568741,
		       3.568656,3.568812,3.568712,3.568734,3.568659,3.568854,3.568885,3.568692,3.568926,3.568882,3.568797,3.568945,3.568887,3.569005,3.568945,3.568872,3.569075,3.569540,3.569856,3.569869,3.569758,3.569846,3.569960,
		       3.569608,3.569670,3.569800,3.569511,3.569921,3.570298,3.570322,3.570332,3.570586,3.570489,3.571143,3.571122,3.571005,3.570983,3.571283,3.571207,3.571646,3.572005,3.572294,3.572640,3.572454,3.572985,3.573242,
		       3.573122,3.573409,3.573730,3.573847,3.574047,3.574057,3.574320,3.574750,3.574671,3.574629,3.574749,3.574825,3.575110,3.575288,3.575142,3.575135,3.575430,3.575209,3.575085,3.574964,3.575442,3.575474,3.575504,
		       3.575606,3.575859,3.575655,3.575990,3.576159,3.575925,3.575717,3.576006,3.575427,3.575602,3.575647,3.575675,3.575689,3.575680,3.575821,3.575970,3.575801,3.575997,3.575863,3.575829,3.575950,3.576194,3.576241,
		       3.576257,3.576189,3.576174,3.576303,3.576327,3.576308,3.576258,3.576341,3.576527,3.576310,3.576297,3.576468,3.576551,3.576468,3.576504,3.576346,3.576391,3.576517,3.576198,3.576264,3.576354,3.576364,3.576172,
		       3.576042,3.576061,3.575789,3.575548,3.574428,3.574141,3.574141,3.574141,3.574141,3.574141,3.574141,3.574141,3.574141,3.574141,3.574141,3.574141,3.574141,3.574141,3.574141,3.574141,3.574141,3.574141,3.574141,
		       3.574141,3.572871];
tracks.m20.log_l = [4.614872,4.619018,4.623463,4.628064,4.632737,4.637466,4.642282,4.646902,4.651751,4.656385,4.661088,4.665791,4.670511,4.675113,4.679815,4.684374,4.688887,4.693485,4.698080,4.702542,4.707066,4.711544,4.715935,
		    4.720303,4.724738,4.729042,4.733292,4.737656,4.741863,4.746114,4.750368,4.754505,4.758829,4.762955,4.767060,4.771277,4.775252,4.779310,4.783414,4.787365,4.791604,4.795450,4.799393,4.803348,4.807345,4.811222,
		    4.815075,4.819143,4.822897,4.826832,4.830508,4.834208,4.838174,4.841852,4.845683,4.849399,4.853005,4.856893,4.860380,4.864053,4.867646,4.871420,4.875058,4.878547,4.882123,4.885772,4.889355,4.892642,4.896518,
		    4.899881,4.903175,4.906892,4.910362,4.913800,4.917309,4.920896,4.924358,4.927805,4.931261,4.934769,4.938316,4.941806,4.945487,4.949152,4.953346,4.953955,4.954392,4.954824,4.955324,4.955727,4.956250,4.956724,
		    4.957139,4.957825,4.958402,4.959010,4.959752,4.960451,4.961028,4.961796,4.962637,4.963474,4.964382,4.965509,4.966926,4.968583,4.970716,4.973939,4.979591,5.001445,5.001124,5.001037,5.001511,5.001620,5.001418,
		    5.000994,5.000533,4.999532,4.999009,4.998868,4.998832,4.999701,5.000112,5.000569,5.000940,5.001460,5.002370,5.003414,5.003940,5.004542,5.007079,5.008093,5.009222,5.010330,5.012720,5.013875,5.015108,5.016161,
		    5.017312,5.019459,5.020446,5.021411,5.022304,5.023958,5.024722,5.025430,5.026101,5.027445,5.028046,5.028629,5.029155,5.030731,5.032783,5.033671,5.034464,5.035161,5.036341,5.036911,5.037367,5.037756,5.038404,
		    5.038645,5.038963,5.039237,5.039694,5.039913,5.040099,5.040396,5.040617,5.041060,5.041325,5.041635,5.041972,5.043553,5.043777,5.044529,5.045130,5.046146,5.046747,5.047740,5.049107,5.050077,5.051664,5.053038,
		    5.054470,5.055889,5.057324,5.058754,5.060176,5.061554,5.079041,5.077968,5.076456,5.075491,5.074597,5.073595,5.072307,5.071644,5.070583,5.069726,5.068636,5.068058,5.068567,5.069077,5.068065,5.064140,5.031244,
		    5.035243,5.035701,5.035921,5.036062,5.036010,5.036046,5.036147,5.036309,5.036548,5.036525,5.036667,5.036822,5.036834,5.037020,5.037013,5.037229,5.036882,5.037137,5.037262,5.037391,5.037315,5.037339,5.037681,
		    5.037828,5.037975,5.038057,5.038071,5.038007,5.038149,5.038147,5.038148,5.038245,5.038156,5.038130,5.038445,5.038370,5.038665,5.038515,5.038405,5.038626,5.038808,5.038602,5.038637,5.038637,5.038579,5.038594,
		    5.038617,5.038437,5.038473,5.037165,5.037742,5.039143,5.038086,5.037764,5.037443,5.037571,5.037662,5.037592,5.037355,5.037118,5.036936,5.036902,5.036804,5.036570,5.036718,5.036260,5.035812,5.035607,5.035755,
		    5.034917,5.034217,5.034203,5.033397,5.032968,5.032967,5.033404,5.032948,5.032644,5.032308,5.032250,5.032060,5.032566,5.031949,5.031713,5.031621,5.031498,5.031133,5.031457,5.030324,5.031045,5.030560,5.030498,
		    5.030130,5.030258,5.029916,5.029750,5.030537,5.029141,5.028812,5.029400,5.029178,5.029342,5.029554,5.028975,5.028874,5.028878,5.028918,5.029124,5.028626,5.028719,5.028838,5.028839,5.028714,5.028546,5.028391,
		    5.028626,5.028693,5.028537,5.028599,5.028749,5.028970,5.028942,5.028997,5.029142,5.029036,5.029335,5.029563,5.029759,5.029464,5.030334,5.030595,5.030709,5.031510,5.031447,5.031968,5.032612,5.033531,5.034031,
		    5.034918,5.036393,5.037789,5.040777,5.056702,5.058582,5.058582,5.058582,5.058582,5.058582,5.058582,5.058582,5.058582,5.058582,5.058582,5.058582,5.058582,5.058582,5.058582,5.058582,5.058582,5.058582,5.058582,
		    5.058582,5.177863];
variable mass;
foreach mass (get_struct_field_names(tracks))
{
  variable track = get_struct_field(tracks, mass);
  xfig_plot_text( xfig, sprintf("$%s\\,M_{\\odot}$", strreplace(mass,"m","")), track.log_teff[0]-3, track.log_l[0], .1, 1.; size="small", depth=2, color="black");
  xfig.plot(track.log_teff[[0:108:1]]-3,track.log_l[[0:108:1]]; depth=2, color="black"); % hydrogen core burning => main-sequence phase
  xfig.plot(track.log_teff[[109:369:1]]-3,track.log_l[[109:369:1]]; depth=2, color="gray"); % post main-sequence phase
}
% --------------------
% -----
variable phi = [0:2*PI:#360];
% Sun:
xfig_plot_text(xfig, "$\\odot$", log10(5785)-3, 0, 0, 0; size="large", depth=1, color="goldenrod");
% Dwarfs:
variable angle = -33; % value depends on aspect ratio of figure
variable zams = struct{ log_teff, log_l }; % zero-age main sequence according to Ekstroem et al. 2012, A&A, 537, A146
zams.log_teff = [3.688210,3.724350,3.754330,3.779290,3.803820,3.824326,3.856430,3.909891,3.972774,4.048683,4.104477,4.186443,4.247069,4.332096,4.391638,4.456053,4.501169,4.553253,4.588704,4.623051,4.649728,4.689020,4.713497,4.729690];
zams.log_l = [-0.590900,-0.356030,-0.138050,0.063710,0.295030,0.459878,0.673904,0.913744,1.208899,1.596045,1.898976,2.365394,2.712744,3.220316,3.584321,3.977491,4.266387,4.614872,4.865990,5.125097,5.343182,5.703098,5.980312,6.231065];
variable tams = struct{ log_teff, log_l }; % terminal-age main sequence according to Ekstroem et al. 2012, A&A, 537, A146
tams.log_teff = [3.713490,3.741400,3.762730,3.781080,3.788900,3.803512,3.813527,3.820671,3.860109,3.926968,3.984781,4.077595,4.145771,4.237784,4.295675,4.352039,4.386232,4.417494,4.427468,4.422336,4.397432,4.152124,4.310353,4.433115];
tams.log_l = [-0.258470,-0.046390,0.126790,0.276170,0.493410,0.646927,0.833820,1.074095,1.372272,1.782232,2.104768,2.608506,2.981932,3.521011,3.906572,4.312416,4.603134,4.955324,5.191676,5.421882,5.607502,5.874329,6.106877,6.310356];
array_reverse(tams.log_teff);
array_reverse(tams.log_l);
xfig.shade_region([zams.log_teff,tams.log_teff]-3,[zams.log_l,tams.log_l]; width=0, color="darkgreen", depth=3);
xfig.add_object( xfig_new_text(sprintf("Dwarfs"); rotate=angle, depth=1, color="whitesmoke"), log10(14500)-3, 2.6);
% Subgiants:
xfig.add_object( xfig_new_text(sprintf("Subgiants"); rotate=angle, depth=1, color="black"), log10(11500)-3, 3.15);
% Supergiants:
xfig.shade_region((0.35*cos(phi)+log10(8)),(0.75*sin(phi)+4.75); width=0, depth=3, color="blue");
xfig_plot_text(xfig, "Supergiants", log10(8), 4.75, 0, 0; depth=1, color="whitesmoke");
% Giants:
xfig.shade_region((0.15*cos(phi)+log10(4.5)),(0.75*sin(phi)+2.75); width=0, depth=3, color="red");
xfig_plot_text(xfig, "Giants", log10(4.5), 2.75, 0, 0; depth=1, color="whitesmoke");
% White dwarfs
xfig.shade_region([0.1,0.2,0.4,0.3,0.2],[0.1,0.0,0.0,0.1,0.2]; width=0, depth=3, color="powderblue", world00);
xfig_plot_text(xfig, "White", log10(26), -0.05, 0, 0; rotate=angle-5, depth=1, color="black");
xfig_plot_text(xfig, "dwarfs", log10(28), -0.45, 0, 0; rotate=angle-5, depth=1, color="black");
% -----
xfig.render("HR_diagram.pdf");
isis/slxfig/hr_diagram.txt · Last modified: 2014-08-08 15:14 by Andreas Irrgang