In [1]:
import mesa as ms
# quantidade de zonas no modelo 500
folder2 = '/home/glauffer/Dropbox/UFRGS/kepler/mesa/alpha2/'
a = ms.mesa_profile(folder2, 500, num_type='nearest_model')


58 in profiles.index file ...
Found and load nearest profile for cycle 500
reading /home/glauffer/Dropbox/UFRGS/kepler/mesa/alpha2//profile16.data ...
 reading ...100% 


In [2]:
a.cols


Out[2]:
{'abar': 29,
 'am_log_D_DSI': 117,
 'am_log_D_ES': 120,
 'am_log_D_GSF': 121,
 'am_log_D_SH': 118,
 'am_log_D_SSI': 119,
 'am_log_D_ST': 122,
 'am_log_nu_ST': 123,
 'am_log_sig': 116,
 'burn_ar': 86,
 'burn_c': 78,
 'burn_ca': 87,
 'burn_cr': 89,
 'burn_fe': 90,
 'burn_mg': 83,
 'burn_n': 79,
 'burn_na': 82,
 'burn_ne': 81,
 'burn_o': 80,
 'burn_s': 85,
 'burn_si': 84,
 'burn_ti': 88,
 'c12': 62,
 'c12_c12': 91,
 'c12_o16': 92,
 'c13': 63,
 'cno': 76,
 'conv_vel_div_csound': 44,
 'csound': 11,
 'dlnT_dt': 55,
 'dlnd_dt': 54,
 'dq_ratio': 18,
 'dynamo_log_B_phi': 125,
 'dynamo_log_B_r': 124,
 'entropy': 9,
 'eps_grav': 34,
 'eps_nuc': 32,
 'eta': 15,
 'f19': 70,
 'fp_rot': 106,
 'free_e': 28,
 'ft_rot': 107,
 'gradL': 51,
 'gradT': 47,
 'gradT_sub_grada': 37,
 'grad_density': 50,
 'grad_temperature': 49,
 'grada': 27,
 'gradr': 48,
 'gradr_sub_grada': 126,
 'h1': 59,
 'he3': 60,
 'he4': 61,
 'i_rot': 104,
 'j_rot': 105,
 'ledoux_stable': 53,
 'logL': 7,
 'logP': 4,
 'logPgas': 26,
 'logR': 5,
 'logRho': 3,
 'logT': 2,
 'log_D_conv': 38,
 'log_D_mix': 42,
 'log_D_ovr': 40,
 'log_D_semi': 39,
 'log_D_thrm': 41,
 'log_J_div_M53': 100,
 'log_J_inside': 101,
 'log_abs_shear': 102,
 'log_am_nu': 108,
 'log_conv_vel': 43,
 'log_j_rot': 99,
 'log_mlt_D_mix': 45,
 'log_omega': 98,
 'log_opacity': 31,
 'log_r_equatorial': 112,
 'log_r_polar': 110,
 'logdq': 17,
 'logtau': 23,
 'logxq': 58,
 'luminosity': 6,
 'mass': 56,
 'mg24': 73,
 'mixing_type': 10,
 'mlt_mixing_length': 35,
 'mlt_mixing_type': 36,
 'mmid': 57,
 'mu': 16,
 'n13': 64,
 'n14': 65,
 'n15': 66,
 'ne20': 71,
 'ne22': 72,
 'non_nuc_neu': 33,
 'o16': 67,
 'o16_o16': 93,
 'o17': 68,
 'o18': 69,
 'omega': 97,
 'omega_crit': 114,
 'omega_div_omega_crit': 115,
 'other': 96,
 'pgas_div_ptotal': 25,
 'photo': 95,
 'pnhe4': 94,
 'pp': 75,
 'pressure': 24,
 'pressure_scale_height': 46,
 'q': 19,
 'r_e_div_r_p': 113,
 'r_equatorial': 111,
 'r_polar': 109,
 'radius': 20,
 'sch_stable': 52,
 'si28': 74,
 'tau': 22,
 'temperature': 21,
 'total_energy': 13,
 'tri_alfa': 77,
 'v_div_csound': 12,
 'v_div_r': 14,
 'v_rot': 103,
 'velocity': 8,
 'ye': 30,
 'zone': 1}

In [3]:
z = a.get('zone')

In [4]:
z


Out[4]:
array([   1.,    2.,    3.,    4.,    5.,    6.,    7.,    8.,    9.,
         10.,   11.,   12.,   13.,   14.,   15.,   16.,   17.,   18.,
         19.,   20.,   21.,   22.,   23.,   24.,   25.,   26.,   27.,
         28.,   29.,   30.,   31.,   32.,   33.,   34.,   35.,   36.,
         37.,   38.,   39.,   40.,   41.,   42.,   43.,   44.,   45.,
         46.,   47.,   48.,   49.,   50.,   51.,   52.,   53.,   54.,
         55.,   56.,   57.,   58.,   59.,   60.,   61.,   62.,   63.,
         64.,   65.,   66.,   67.,   68.,   69.,   70.,   71.,   72.,
         73.,   74.,   75.,   76.,   77.,   78.,   79.,   80.,   81.,
         82.,   83.,   84.,   85.,   86.,   87.,   88.,   89.,   90.,
         91.,   92.,   93.,   94.,   95.,   96.,   97.,   98.,   99.,
        100.,  101.,  102.,  103.,  104.,  105.,  106.,  107.,  108.,
        109.,  110.,  111.,  112.,  113.,  114.,  115.,  116.,  117.,
        118.,  119.,  120.,  121.,  122.,  123.,  124.,  125.,  126.,
        127.,  128.,  129.,  130.,  131.,  132.,  133.,  134.,  135.,
        136.,  137.,  138.,  139.,  140.,  141.,  142.,  143.,  144.,
        145.,  146.,  147.,  148.,  149.,  150.,  151.,  152.,  153.,
        154.,  155.,  156.,  157.,  158.,  159.,  160.,  161.,  162.,
        163.,  164.,  165.,  166.,  167.,  168.,  169.,  170.,  171.,
        172.,  173.,  174.,  175.,  176.,  177.,  178.,  179.,  180.,
        181.,  182.,  183.,  184.,  185.,  186.,  187.,  188.,  189.,
        190.,  191.,  192.,  193.,  194.,  195.,  196.,  197.,  198.,
        199.,  200.,  201.,  202.,  203.,  204.,  205.,  206.,  207.,
        208.,  209.,  210.,  211.,  212.,  213.,  214.,  215.,  216.,
        217.,  218.,  219.,  220.,  221.,  222.,  223.,  224.,  225.,
        226.,  227.,  228.,  229.,  230.,  231.,  232.,  233.,  234.,
        235.,  236.,  237.,  238.,  239.,  240.,  241.,  242.,  243.,
        244.,  245.,  246.,  247.,  248.,  249.,  250.,  251.,  252.,
        253.,  254.,  255.,  256.,  257.,  258.,  259.,  260.,  261.,
        262.,  263.,  264.,  265.,  266.,  267.,  268.,  269.,  270.,
        271.,  272.,  273.,  274.,  275.,  276.,  277.,  278.,  279.,
        280.,  281.,  282.,  283.,  284.,  285.,  286.,  287.,  288.,
        289.,  290.,  291.,  292.,  293.,  294.,  295.,  296.,  297.,
        298.,  299.,  300.,  301.,  302.,  303.,  304.,  305.,  306.,
        307.,  308.,  309.,  310.,  311.,  312.,  313.,  314.,  315.,
        316.,  317.,  318.,  319.,  320.,  321.,  322.,  323.,  324.,
        325.,  326.,  327.,  328.,  329.,  330.,  331.,  332.,  333.,
        334.,  335.,  336.,  337.,  338.,  339.,  340.,  341.,  342.,
        343.,  344.,  345.,  346.,  347.,  348.,  349.,  350.,  351.,
        352.,  353.,  354.,  355.,  356.,  357.,  358.,  359.,  360.,
        361.,  362.,  363.,  364.,  365.,  366.,  367.,  368.,  369.,
        370.,  371.,  372.,  373.,  374.,  375.,  376.,  377.,  378.,
        379.,  380.,  381.,  382.,  383.,  384.,  385.,  386.,  387.,
        388.,  389.,  390.,  391.,  392.,  393.,  394.,  395.,  396.,
        397.,  398.,  399.,  400.,  401.,  402.,  403.,  404.,  405.,
        406.,  407.,  408.,  409.,  410.,  411.,  412.,  413.,  414.,
        415.,  416.,  417.,  418.,  419.,  420.,  421.,  422.,  423.,
        424.,  425.,  426.,  427.,  428.,  429.,  430.,  431.,  432.,
        433.,  434.,  435.,  436.,  437.,  438.,  439.,  440.,  441.,
        442.,  443.,  444.,  445.,  446.,  447.,  448.,  449.,  450.,
        451.,  452.,  453.,  454.,  455.,  456.,  457.,  458.,  459.,
        460.,  461.,  462.,  463.,  464.,  465.,  466.,  467.,  468.,
        469.,  470.,  471.,  472.,  473.,  474.,  475.,  476.,  477.,
        478.,  479.,  480.,  481.,  482.,  483.,  484.,  485.,  486.,
        487.,  488.,  489.,  490.,  491.,  492.,  493.,  494.,  495.,
        496.,  497.,  498.,  499.,  500.,  501.,  502.,  503.,  504.,
        505.,  506.,  507.,  508.,  509.,  510.,  511.,  512.,  513.,
        514.,  515.,  516.,  517.,  518.,  519.,  520.,  521.,  522.,
        523.,  524.,  525.,  526.,  527.,  528.,  529.,  530.,  531.,
        532.,  533.,  534.,  535.,  536.,  537.,  538.,  539.,  540.,
        541.,  542.,  543.,  544.,  545.,  546.,  547.,  548.,  549.,
        550.,  551.,  552.,  553.,  554.,  555.,  556.,  557.,  558.,
        559.,  560.,  561.,  562.,  563.,  564.,  565.,  566.,  567.,
        568.,  569.,  570.,  571.,  572.,  573.,  574.,  575.,  576.,
        577.,  578.,  579.,  580.,  581.,  582.,  583.,  584.,  585.,
        586.,  587.,  588.,  589.,  590.,  591.,  592.,  593.,  594.,
        595.,  596.,  597.,  598.,  599.,  600.,  601.,  602.,  603.,
        604.,  605.,  606.,  607.,  608.,  609.,  610.,  611.,  612.,
        613.,  614.,  615.,  616.,  617.,  618.,  619.,  620.,  621.,
        622.,  623.,  624.,  625.,  626.,  627.,  628.,  629.,  630.,
        631.,  632.,  633.,  634.,  635.,  636.,  637.,  638.,  639.,
        640.,  641.,  642.,  643.,  644.,  645.,  646.,  647.,  648.,
        649.,  650.,  651.,  652.,  653.,  654.,  655.,  656.,  657.,
        658.,  659.,  660.,  661.,  662.,  663.,  664.,  665.,  666.,
        667.,  668.,  669.,  670.,  671.,  672.,  673.,  674.,  675.,
        676.,  677.,  678.,  679.,  680.,  681.,  682.,  683.,  684.,
        685.,  686.,  687.,  688.,  689.,  690.,  691.,  692.,  693.,
        694.,  695.,  696.,  697.,  698.,  699.,  700.,  701.,  702.,
        703.,  704.,  705.,  706.,  707.,  708.,  709.,  710.,  711.,
        712.,  713.,  714.,  715.,  716.,  717.,  718.,  719.,  720.,
        721.,  722.,  723.,  724.,  725.,  726.,  727.,  728.,  729.,
        730.,  731.,  732.,  733.,  734.,  735.,  736.,  737.,  738.,
        739.,  740.,  741.,  742.,  743.,  744.,  745.,  746.,  747.,
        748.,  749.,  750.,  751.,  752.,  753.,  754.,  755.,  756.,
        757.,  758.,  759.,  760.,  761.,  762.,  763.,  764.,  765.,
        766.,  767.,  768.,  769.,  770.,  771.,  772.,  773.,  774.,
        775.,  776.,  777.,  778.,  779.,  780.,  781.,  782.,  783.,
        784.,  785.,  786.,  787.,  788.,  789.,  790.,  791.,  792.,
        793.,  794.,  795.,  796.,  797.,  798.,  799.,  800.,  801.,
        802.,  803.,  804.,  805.,  806.,  807.,  808.,  809.,  810.,
        811.,  812.,  813.,  814.,  815.,  816.,  817.,  818.,  819.,
        820.,  821.,  822.,  823.,  824.,  825.,  826.,  827.,  828.,
        829.,  830.,  831.,  832.,  833.,  834.,  835.,  836.,  837.,
        838.,  839.,  840.,  841.,  842.,  843.,  844.,  845.,  846.,
        847.,  848.,  849.,  850.,  851.,  852.,  853.,  854.,  855.,
        856.,  857.,  858.,  859.,  860.,  861.,  862.,  863.,  864.,
        865.,  866.,  867.,  868.,  869.,  870.,  871.,  872.,  873.,
        874.,  875.,  876.,  877.,  878.,  879.,  880.,  881.,  882.,
        883.,  884.,  885.,  886.,  887.,  888.,  889.,  890.,  891.,
        892.,  893.,  894.,  895.,  896.,  897.,  898.,  899.,  900.,
        901.,  902.,  903.,  904.,  905.,  906.,  907.,  908.,  909.,
        910.,  911.,  912.,  913.,  914.,  915.,  916.,  917.,  918.,
        919.,  920.,  921.,  922.,  923.,  924.,  925.,  926.,  927.,
        928.,  929.,  930.,  931.,  932.,  933.,  934.,  935.,  936.,
        937.,  938.,  939.,  940.,  941.,  942.,  943.,  944.,  945.,
        946.,  947.,  948.,  949.,  950.,  951.,  952.,  953.,  954.,
        955.,  956.,  957.,  958.,  959.,  960.,  961.,  962.,  963.,
        964.,  965.,  966.,  967.,  968.,  969.,  970.,  971.,  972.,
        973.,  974.,  975.,  976.,  977.,  978.,  979.,  980.,  981.,
        982.,  983.,  984.,  985.,  986.,  987.,  988.,  989.,  990.,
        991.,  992.,  993.,  994.])

In [6]:
# quantidade de zonas no modelo 2090
a2090 = ms.mesa_profile(folder2, 2090, num_type='nearest_model')
z2090 = a2090.get('zone')
z2090


58 in profiles.index file ...
Found and load nearest profile for cycle 2100
reading /home/glauffer/Dropbox/UFRGS/kepler/mesa/alpha2//profile48.data ...
 reading ...100% 

Out[6]:
array([   1.,    2.,    3.,    4.,    5.,    6.,    7.,    8.,    9.,
         10.,   11.,   12.,   13.,   14.,   15.,   16.,   17.,   18.,
         19.,   20.,   21.,   22.,   23.,   24.,   25.,   26.,   27.,
         28.,   29.,   30.,   31.,   32.,   33.,   34.,   35.,   36.,
         37.,   38.,   39.,   40.,   41.,   42.,   43.,   44.,   45.,
         46.,   47.,   48.,   49.,   50.,   51.,   52.,   53.,   54.,
         55.,   56.,   57.,   58.,   59.,   60.,   61.,   62.,   63.,
         64.,   65.,   66.,   67.,   68.,   69.,   70.,   71.,   72.,
         73.,   74.,   75.,   76.,   77.,   78.,   79.,   80.,   81.,
         82.,   83.,   84.,   85.,   86.,   87.,   88.,   89.,   90.,
         91.,   92.,   93.,   94.,   95.,   96.,   97.,   98.,   99.,
        100.,  101.,  102.,  103.,  104.,  105.,  106.,  107.,  108.,
        109.,  110.,  111.,  112.,  113.,  114.,  115.,  116.,  117.,
        118.,  119.,  120.,  121.,  122.,  123.,  124.,  125.,  126.,
        127.,  128.,  129.,  130.,  131.,  132.,  133.,  134.,  135.,
        136.,  137.,  138.,  139.,  140.,  141.,  142.,  143.,  144.,
        145.,  146.,  147.,  148.,  149.,  150.,  151.,  152.,  153.,
        154.,  155.,  156.,  157.,  158.,  159.,  160.,  161.,  162.,
        163.,  164.,  165.,  166.,  167.,  168.,  169.,  170.,  171.,
        172.,  173.,  174.,  175.,  176.,  177.,  178.,  179.,  180.,
        181.,  182.,  183.,  184.,  185.,  186.,  187.,  188.,  189.,
        190.,  191.,  192.,  193.,  194.,  195.,  196.,  197.,  198.,
        199.,  200.,  201.,  202.,  203.,  204.,  205.,  206.,  207.,
        208.,  209.,  210.,  211.,  212.,  213.,  214.,  215.,  216.,
        217.,  218.,  219.,  220.,  221.,  222.,  223.,  224.,  225.,
        226.,  227.,  228.,  229.,  230.,  231.,  232.,  233.,  234.,
        235.,  236.,  237.,  238.,  239.,  240.,  241.,  242.,  243.,
        244.,  245.,  246.,  247.,  248.,  249.,  250.,  251.,  252.,
        253.,  254.,  255.,  256.,  257.,  258.,  259.,  260.,  261.,
        262.,  263.,  264.,  265.,  266.,  267.,  268.,  269.,  270.,
        271.,  272.,  273.,  274.,  275.,  276.,  277.,  278.,  279.,
        280.,  281.,  282.,  283.,  284.,  285.,  286.,  287.,  288.,
        289.,  290.,  291.,  292.,  293.,  294.,  295.,  296.,  297.,
        298.,  299.,  300.,  301.,  302.,  303.,  304.,  305.,  306.,
        307.,  308.,  309.,  310.,  311.,  312.,  313.,  314.,  315.,
        316.,  317.,  318.,  319.,  320.,  321.,  322.,  323.,  324.,
        325.,  326.,  327.,  328.,  329.,  330.,  331.,  332.,  333.,
        334.,  335.,  336.,  337.,  338.,  339.,  340.,  341.,  342.,
        343.,  344.,  345.,  346.,  347.,  348.,  349.,  350.,  351.,
        352.,  353.,  354.,  355.,  356.,  357.,  358.,  359.,  360.,
        361.,  362.,  363.,  364.,  365.,  366.,  367.,  368.,  369.,
        370.,  371.,  372.,  373.,  374.,  375.,  376.,  377.,  378.,
        379.,  380.,  381.,  382.,  383.,  384.,  385.,  386.,  387.,
        388.,  389.,  390.,  391.,  392.,  393.,  394.,  395.,  396.,
        397.,  398.,  399.,  400.,  401.,  402.,  403.,  404.,  405.,
        406.,  407.,  408.,  409.,  410.,  411.,  412.,  413.,  414.,
        415.,  416.,  417.,  418.,  419.,  420.,  421.,  422.,  423.,
        424.,  425.,  426.,  427.,  428.,  429.,  430.,  431.,  432.,
        433.,  434.,  435.,  436.,  437.,  438.,  439.,  440.,  441.,
        442.,  443.,  444.,  445.,  446.,  447.,  448.,  449.,  450.,
        451.,  452.,  453.,  454.,  455.,  456.,  457.,  458.,  459.,
        460.,  461.,  462.,  463.,  464.,  465.,  466.,  467.,  468.,
        469.,  470.,  471.,  472.,  473.,  474.,  475.,  476.,  477.,
        478.,  479.,  480.,  481.,  482.,  483.,  484.,  485.,  486.,
        487.,  488.,  489.,  490.,  491.,  492.,  493.,  494.,  495.,
        496.,  497.,  498.,  499.,  500.,  501.,  502.,  503.,  504.,
        505.,  506.,  507.,  508.,  509.,  510.,  511.,  512.,  513.,
        514.,  515.,  516.,  517.,  518.,  519.,  520.,  521.,  522.,
        523.,  524.,  525.,  526.,  527.,  528.,  529.,  530.,  531.,
        532.,  533.,  534.,  535.,  536.,  537.,  538.,  539.,  540.,
        541.,  542.,  543.,  544.,  545.,  546.,  547.,  548.,  549.,
        550.,  551.,  552.,  553.,  554.,  555.,  556.,  557.,  558.,
        559.,  560.,  561.,  562.,  563.,  564.,  565.,  566.,  567.,
        568.,  569.,  570.,  571.,  572.,  573.,  574.,  575.,  576.,
        577.,  578.,  579.,  580.,  581.,  582.,  583.,  584.,  585.,
        586.,  587.,  588.,  589.,  590.,  591.,  592.,  593.,  594.,
        595.,  596.,  597.,  598.,  599.,  600.,  601.,  602.,  603.,
        604.,  605.,  606.,  607.,  608.,  609.,  610.,  611.,  612.,
        613.,  614.,  615.,  616.,  617.,  618.,  619.,  620.,  621.,
        622.,  623.,  624.,  625.,  626.,  627.,  628.,  629.,  630.,
        631.,  632.,  633.,  634.,  635.,  636.,  637.,  638.,  639.,
        640.,  641.,  642.,  643.,  644.,  645.,  646.,  647.,  648.,
        649.,  650.,  651.,  652.,  653.,  654.,  655.,  656.,  657.,
        658.,  659.,  660.,  661.,  662.,  663.,  664.,  665.,  666.,
        667.,  668.,  669.,  670.,  671.,  672.,  673.,  674.,  675.,
        676.,  677.,  678.,  679.,  680.,  681.,  682.,  683.,  684.,
        685.,  686.,  687.,  688.,  689.,  690.,  691.,  692.,  693.,
        694.,  695.,  696.,  697.,  698.,  699.,  700.,  701.,  702.,
        703.,  704.,  705.,  706.,  707.,  708.,  709.,  710.,  711.,
        712.,  713.,  714.,  715.,  716.,  717.,  718.,  719.,  720.,
        721.,  722.,  723.,  724.,  725.,  726.,  727.,  728.,  729.,
        730.,  731.,  732.,  733.,  734.,  735.,  736.,  737.,  738.,
        739.,  740.,  741.,  742.,  743.,  744.,  745.,  746.,  747.,
        748.,  749.,  750.,  751.,  752.,  753.,  754.,  755.,  756.,
        757.,  758.,  759.,  760.,  761.,  762.,  763.,  764.,  765.,
        766.,  767.,  768.,  769.,  770.,  771.,  772.,  773.,  774.,
        775.,  776.,  777.,  778.,  779.,  780.,  781.,  782.,  783.,
        784.,  785.,  786.,  787.,  788.,  789.,  790.,  791.,  792.,
        793.,  794.,  795.,  796.,  797.,  798.,  799.,  800.,  801.,
        802.,  803.,  804.,  805.,  806.,  807.,  808.,  809.,  810.,
        811.,  812.,  813.,  814.,  815.,  816.,  817.,  818.,  819.,
        820.,  821.,  822.,  823.,  824.,  825.,  826.,  827.,  828.,
        829.,  830.,  831.,  832.,  833.,  834.,  835.,  836.,  837.,
        838.,  839.,  840.,  841.,  842.,  843.,  844.,  845.,  846.,
        847.,  848.,  849.,  850.,  851.,  852.,  853.,  854.,  855.,
        856.,  857.,  858.,  859.,  860.,  861.,  862.,  863.,  864.,
        865.,  866.,  867.,  868.,  869.,  870.,  871.,  872.,  873.,
        874.,  875.,  876.,  877.,  878.,  879.,  880.,  881.,  882.,
        883.,  884.,  885.,  886.,  887.,  888.,  889.,  890.,  891.,
        892.,  893.,  894.,  895.,  896.,  897.,  898.,  899.,  900.,
        901.,  902.,  903.,  904.,  905.,  906.,  907.,  908.,  909.,
        910.,  911.,  912.,  913.,  914.,  915.,  916.,  917.,  918.,
        919.,  920.,  921.,  922.,  923.,  924.,  925.,  926.,  927.,
        928.,  929.,  930.,  931.,  932.,  933.,  934.,  935.,  936.,
        937.,  938.,  939.,  940.,  941.,  942.,  943.,  944.,  945.,
        946.,  947.,  948.,  949.,  950.,  951.,  952.,  953.,  954.,
        955.,  956.,  957.,  958.,  959.,  960.,  961.,  962.,  963.,
        964.,  965.,  966.,  967.,  968.,  969.])

In [8]:
t2090 = a2090.get('logT')

In [9]:
t2090


Out[9]:
array([ 3.78719491,  3.78719492,  3.78719493,  3.78719495,  3.78719496,
        3.78719499,  3.78719503,  3.78719509,  3.78719516,  3.78719525,
        3.78719536,  3.78719554,  3.7871958 ,  3.78719611,  3.78719649,
        3.78719704,  3.78719775,  3.78719894,  3.78720058,  3.78720294,
        3.787206  ,  3.78721017,  3.78721544,  3.78722275,  3.78723212,
        3.78724485,  3.78726095,  3.78728216,  3.78730851,  3.78734232,
        3.78738365,  3.78743523,  3.78749714,  3.78757283,  3.78766243,
        3.78777307,  3.78790501,  3.7880742 ,  3.78828133,  3.78855846,
        3.78890771,  3.78939223,  3.79001915,  3.79089798,  3.79205266,
        3.79363935,  3.79573115,  3.7985719 ,  3.80218789,  3.80680435,
        3.81243658,  3.82021253,  3.83063209,  3.84737008,  3.87571331,
        3.91077412,  3.94239392,  3.9590003 ,  3.97420115,  3.98988484,
        4.00849633,  4.0335859 ,  4.05981097,  4.07882079,  4.09312834,
        4.10457603,  4.11414465,  4.1223866 ,  4.12963993,  4.13613103,
        4.14202205,  4.14743291,  4.15158667,  4.15472422,  4.15773667,
        4.16063751,  4.16343825,  4.16614879,  4.16877767,  4.17133231,
        4.1738192 ,  4.17624405,  4.17979131,  4.18432609,  4.18868756,
        4.19445698,  4.20141835,  4.20805629,  4.21441495,  4.22052493,
        4.22640821,  4.23208141,  4.24026815,  4.25048489,  4.26004018,
        4.26899897,  4.27866584,  4.28996337,  4.29859173,  4.30410354,
        4.31204103,  4.3207097 ,  4.32777506,  4.33340573,  4.33994392,
        4.3482184 ,  4.35562725,  4.36327802,  4.37287358,  4.38064963,
        4.38816457,  4.39651374,  4.40339965,  4.40894661,  4.41543073,
        4.42373638,  4.43120717,  4.43743934,  4.44440392,  4.4537991 ,
        4.46406688,  4.47389712,  4.48332449,  4.49013269,  4.49677337,
        4.50590158,  4.51635712,  4.52796339,  4.53813279,  4.54517266,
        4.55069148,  4.55746565,  4.56675056,  4.5771423 ,  4.58686957,
        4.59448904,  4.60053824,  4.60733503,  4.61484357,  4.62296388,
        4.63307797,  4.64509223,  4.65650487,  4.66496642,  4.67324551,
        4.68453845,  4.69674507,  4.70529312,  4.71070653,  4.71607252,
        4.72138055,  4.72792904,  4.73567148,  4.74203314,  4.74706753,
        4.75453318,  4.76435872,  4.77266723,  4.77947108,  4.78681711,
        4.79660864,  4.80618392,  4.81186067,  4.81561493,  4.82120165,
        4.82858709,  4.83965971,  4.8507882 ,  4.85821651,  4.86466841,
        4.87015833,  4.87803805,  4.88573523,  4.8908242 ,  4.89818774,
        4.90779258,  4.91617667,  4.92335857,  4.93166687,  4.93963371,
        4.94607878,  4.95335697,  4.96146803,  4.96900416,  4.97596983,
        4.98383558,  4.99260894,  5.00063302,  5.00788574,  5.01530633,
        5.02289132,  5.03074098,  5.03684122,  5.04088531,  5.04584338,
        5.05171253,  5.06019856,  5.06927575,  5.07626055,  5.08344753,
        5.09083824,  5.09842874,  5.10607307,  5.11322981,  5.12005093,
        5.12807656,  5.13731549,  5.14499742,  5.15111339,  5.16023824,
        5.17007832,  5.17752113,  5.18495386,  5.19237851,  5.19984273,
        5.20734899,  5.21485255,  5.22151748,  5.22912721,  5.23661172,
        5.24210783,  5.24721696,  5.25193416,  5.2580185 ,  5.26549551,
        5.27230146,  5.27842274,  5.28336687,  5.28711467,  5.2919323 ,
        5.29783809,  5.30290742,  5.30712513,  5.31274788,  5.3198061 ,
        5.32698432,  5.3324727 ,  5.33855058,  5.3473101 ,  5.35386636,
        5.35835414,  5.36325801,  5.3706116 ,  5.37804673,  5.38300289,
        5.38758666,  5.39179018,  5.39739949,  5.40445199,  5.41051907,
        5.41557169,  5.42001007,  5.42382084,  5.42889805,  5.4363717 ,
        5.44477437,  5.4529798 ,  5.46002696,  5.46576742,  5.47230524,
        5.47801775,  5.48207639,  5.4870916 ,  5.49303587,  5.49826301,
        5.50279638,  5.50761447,  5.5127093 ,  5.51754976,  5.52214075,
        5.52690041,  5.53182515,  5.53630306,  5.5403428 ,  5.54475168,
        5.54951986,  5.5542433 ,  5.55832701,  5.56178517,  5.56578608,
        5.56989192,  5.57499026,  5.57966888,  5.58287751,  5.58641233,
        5.59165562,  5.5971668 ,  5.60124834,  5.60528377,  5.60927294,
        5.6131152 ,  5.61696348,  5.62052022,  5.62403856,  5.62790073,
        5.63133477,  5.63510693,  5.63882876,  5.64214062,  5.64541385,
        5.64900996,  5.6543213 ,  5.65952887,  5.6636297 ,  5.6670003 ,
        5.6709966 ,  5.67590562,  5.68017003,  5.68452737,  5.68900695,
        5.69398616,  5.69924234,  5.70384563,  5.7087963 ,  5.71418588,
        5.71940938,  5.72421913,  5.72863652,  5.73267692,  5.736967  ,
        5.74088976,  5.74476315,  5.74863841,  5.75292643,  5.75762839,
        5.76180208,  5.76547308,  5.76910011,  5.77445542,  5.78005269,
        5.78520591,  5.79060194,  5.79622856,  5.80061699,  5.80398923,
        5.80992702,  5.81620727,  5.82192674,  5.8275335 ,  5.83194795,
        5.83629896,  5.84058833,  5.8453434 ,  5.85105499,  5.85564442,
        5.8596711 ,  5.86413391,  5.86853355,  5.8723883 ,  5.87572372,
        5.87973997,  5.88536325,  5.89111804,  5.89452003,  5.89676847,
        5.90011216,  5.9034162 ,  5.90668837,  5.91099519,  5.9162977 ,
        5.92253999,  5.9286606 ,  5.93367086,  5.9376252 ,  5.94153155,
        5.94539114,  5.95125939,  5.95943989,  5.96574559,  5.96988551,
        5.97600224,  5.98215656,  5.98679795,  5.99222476,  5.99865706,
        6.00558776,  6.01237819,  6.01903397,  6.0255606 ,  6.0319635 ,
        6.03824791,  6.04480158,  6.05160681,  6.06157986,  6.07126596,
        6.07758847,  6.08380492,  6.08992077,  6.09792915,  6.1077154 ,
        6.117275  ,  6.1243064 ,  6.12893138,  6.13577724,  6.14474068,
        6.15352528,  6.16213836,  6.16975001,  6.17640877,  6.18297082,
        6.18944028,  6.19897602,  6.21140651,  6.220537  ,  6.22653162,
        6.23538809,  6.24696682,  6.25548682,  6.26108531,  6.26662065,
        6.2720932 ,  6.28018189,  6.29074822,  6.29851063,  6.30360784,
        6.31113607,  6.32095077,  6.3305228 ,  6.33986169,  6.34671878,
        6.35122532,  6.35789011,  6.36660118,  6.37512646,  6.38347515,
        6.39165516,  6.39967329,  6.40753554,  6.41715293,  6.42561481,
        6.43116126,  6.43844044,  6.44735732,  6.4560842 ,  6.46463058,
        6.47134442,  6.47631138,  6.4836552 ,  6.49246762,  6.50033366,
        6.50806986,  6.51568408,  6.5239261 ,  6.5327671 ,  6.5407389 ,
        6.54610435,  6.54964882,  6.55491713,  6.56185692,  6.56937879,
        6.57745599,  6.58474611,  6.59128048,  6.59772933,  6.60409528,
        6.61224859,  6.62210762,  6.63058208,  6.63773409,  6.64711237,
        6.65575939,  6.66144119,  6.66984176,  6.68083222,  6.68945559,
        6.69582248,  6.70210516,  6.70830447,  6.71442134,  6.72045692,
        6.72935919,  6.73808847,  6.74381333,  6.75226196,  6.76055057,
        6.76598973,  6.7740212 ,  6.7843529 ,  6.79429249,  6.80124345,
        6.80535754,  6.81144887,  6.81941722,  6.82721711,  6.83485214,
        6.84232571,  6.84964132,  6.85680197,  6.86438898,  6.87236023,
        6.88067931,  6.88929967,  6.89662896,  6.90276352,  6.9087605 ,
        6.91510741,  6.9217668 ,  6.9310913 ,  6.94097536,  6.94795586,
        6.95391528,  6.96258967,  6.97087961,  6.97891555,  6.98775823,
        6.99465582,  7.00129306,  7.01013979,  7.01865516,  7.02487347,
        7.03183662,  7.03824452,  7.04460233,  7.05232023,  7.05959162,
        7.06964763,  7.07947473,  7.08668327,  7.09543542,  7.10284034,
        7.10825546,  7.11486983,  7.12261384,  7.13035983,  7.13643637,
        7.14095483,  7.14761073,  7.15738663,  7.16801183,  7.17642127,
        7.18378169,  7.19298257,  7.20156806,  7.20864627,  7.21555645,
        7.22443154,  7.23467084,  7.24219615,  7.24784547,  7.255633  ,
        7.26542871,  7.27303778,  7.27836645,  7.28594347,  7.29824791,
        7.30922296,  7.31511397,  7.32250458,  7.33342139,  7.34138407,
        7.34585954,  7.35092084,  7.35654208,  7.36372663,  7.37254393,
        7.38022687,  7.38668759,  7.39316309,  7.40051113,  7.40990219,
        7.41929807,  7.42688938,  7.43508206,  7.44313496,  7.45031635,
        7.4578655 ,  7.46746592,  7.47756634,  7.48370706,  7.4874483 ,
        7.492513  ,  7.49886981,  7.50850097,  7.51776678,  7.52316824,
        7.53172334,  7.54205446,  7.55073389,  7.55851694,  7.56582171,
        7.57566918,  7.58528857,  7.59247373,  7.60003667,  7.60688611,
        7.61259656,  7.6175382 ,  7.6217182 ,  7.62755095,  7.6350236 ,
        7.6420222 ,  7.65049319,  7.65885893,  7.6651919 ,  7.67153369,
        7.68223233,  7.69354102,  7.70096636,  7.71119433,  7.72166419,
        7.72946345,  7.73853651,  7.74588051,  7.75014497,  7.75600104,
        7.76346947,  7.77086631,  7.77819326,  7.78611127,  7.79463718,
        7.80208552,  7.80843173,  7.81600824,  7.8240302 ,  7.831236  ,
        7.84121128,  7.85392944,  7.8634568 ,  7.86981198,  7.88087084,
        7.8924036 ,  7.89932809,  7.9092709 ,  7.92143414,  7.92991895,
        7.93559293,  7.94475106,  7.95626866,  7.96701489,  7.97561174,
        7.9826188 ,  7.99063272,  7.99836523,  8.00456623,  8.0088216 ,
        8.01178358,  8.01376214,  8.0157469 ,  8.01773797,  8.0197354 ,
        8.02173929,  8.0237497 ,  8.02576673,  8.02779046,  8.02982096,
        8.03185833,  8.03492121,  8.03789725,  8.03973018,  8.04156886,
        8.04341335,  8.04618594,  8.04898243,  8.05085083,  8.05272537,
        8.05460614,  8.05649319,  8.0583866 ,  8.06123326,  8.06410584,
        8.06554633,  8.06650781,  8.06800561,  8.0700475 ,  8.07209709,
        8.07415447,  8.07621973,  8.07829298,  8.07985492,  8.08089762,
        8.08246381,  8.08456162,  8.08666781,  8.0887825 ,  8.09090579,
        8.09303779,  8.0951786 ,  8.09732835,  8.09948713,  8.10168329,
        8.10391757,  8.10616173,  8.1084159 ,  8.1106802 ,  8.11295477,
        8.11466961,  8.1158147 ,  8.11753516,  8.11984126,  8.12215811,
        8.12448585,  8.12682458,  8.12917442,  8.13153549,  8.1339079 ,
        8.13629174,  8.1387879 ,  8.14139911,  8.14336936,  8.14468532,
        8.14666295,  8.14931564,  8.15198222,  8.15399415,  8.15533789,
        8.15735719,  8.16026896,  8.16262241,  8.16419466,  8.16655798,
        8.16893982,  8.170531  ,  8.17212701,  8.17372799,  8.17613483,
        8.17856136,  8.18018258,  8.18261968,  8.18507627,  8.18671733,
        8.18836318,  8.19001382,  8.19249455,  8.1949937 ,  8.19666251,
        8.19916954,  8.20169285,  8.20337683,  8.205905  ,  8.20844581,
        8.21013997,  8.21175651,  8.21329443,  8.21559966,  8.21867434,
        8.22097809,  8.2225103 ,  8.22480231,  8.22708702,  8.22860414,
        8.23086992,  8.23312073,  8.23461266,  8.23668532,  8.23932082,
        8.24192418,  8.2444921 ,  8.24702134,  8.24950965,  8.25195603,
        8.25435958,  8.25671431,  8.25901951,  8.26240484,  8.26566247,
        8.26778669,  8.26986682,  8.27190402,  8.27389947,  8.2758541 ,
        8.27729094,  8.27823851,  8.27964456,  8.28148219,  8.28283316,
        8.28372416,  8.28504638,  8.28677468,  8.28804336,  8.28887736,
        8.29011525,  8.29173361,  8.29331959,  8.29487371,  8.29639647,
        8.29788835,  8.29934982,  8.3007813 ,  8.30218324,  8.30355605,
        8.30490014,  8.30621595,  8.30750388,  8.30876436,  8.3101328 ,
        8.31159725,  8.31273124,  8.31355933,  8.31478151,  8.31636141,
        8.31789257,  8.31937688,  8.32081633,  8.32221309,  8.32325391,
        8.32396743,  8.32502141,  8.32638735,  8.32738477,  8.32804098,
        8.32868917,  8.3293298 ,  8.32996331,  8.33059019,  8.33105579,
        8.33136468,  8.33182584,  8.33243557,  8.33304051,  8.33357566,
        8.33404284,  8.33474015,  8.33543218,  8.33589182,  8.33635009,
        8.33680715,  8.33726319,  8.33771839,  8.33817291,  8.33862691,
        8.33908054,  8.33953396,  8.34021379,  8.34089431,  8.34134815,
        8.34202933,  8.34271246,  8.34316842,  8.34362511,  8.34408261,
        8.34476996,  8.34569042,  8.34638401,  8.34684734,  8.34748003,
        8.34828455,  8.34902814,  8.34970973,  8.35056454,  8.35133828,
        8.35185527,  8.3524599 ,  8.35315326,  8.35384888,  8.35437235,
        8.35472179,  8.35524666,  8.35594872,  8.35665281,  8.35735887,
        8.35788996,  8.35824442,  8.35877669,  8.3598268 ,  8.36087417,
        8.36157368,  8.36239141,  8.36332863,  8.36426826,  8.36497488,
        8.36544639,  8.36615425,  8.36725785,  8.36820684,  8.36884006,
        8.36978359,  8.37166697,  8.37473933,  8.37759139,  8.37981087,
        8.38298341,  8.38656088,  8.38933979,  8.39187193,  8.39416327,
        8.39731269,  8.40098673,  8.40365493,  8.4064646 ,  8.40918854,
        8.41171228,  8.41509272,  8.41819752,  8.42075984,  8.42416452,
        8.4269814 ,  8.42977954,  8.43283413,  8.43575787,  8.43965928,
        8.44366426,  8.44692593,  8.44996215,  8.45260319,  8.45552121,
        8.45964085,  8.46329958,  8.4671525 ,  8.47136216,  8.47446316,
        8.47731144,  8.47990979,  8.48266523,  8.48654149,  8.49005795,
        8.49275728,  8.49564274,  8.49933896,  8.50327706,  8.50670804,
        8.5102793 ,  8.51345518,  8.51658413,  8.51951518,  8.52140966,
        8.52415193,  8.5282429 ,  8.53223711,  8.53482789,  8.53739072,
        8.54093134,  8.54404029,  8.5461984 ,  8.54874319,  8.55150881,
        8.55419245,  8.55655847,  8.55917065,  8.56216913,  8.56399057,
        8.56610841,  8.56914105,  8.57044859,  8.57337184,  8.57614926,
        8.57846313,  8.58116662,  8.58365063,  8.58551063,  8.58724599,
        8.5893088 ,  8.59121019,  8.59319254,  8.59497477,  8.59582172,
        8.5963644 ,  8.59715853,  8.59775269,  8.59812504,  8.59870478,
        8.5994718 ,  8.6001113 ,  8.60064926,  8.60100281,  8.60121302,
        8.60133854,  8.60141277,  8.60145565,  8.60148758,  8.60151977,
        8.60154455,  8.6015601 ,  8.60156988,  8.60157603,  8.6015799 ,
        8.60158233,  8.60158387,  8.60158483,  8.60158565])

teste triplicando o numero de camadas

teste alpha2_3x1_camadas


In [2]:
import mesa as ms
folder = '/media/glauffer/Data/mesa/alpha2_3x1_camadas/LOGS'
s = ms.history_data(folder)


No history.datasa file found, create new one from history.data
 reading ...100% 

Closing history.data  tool ...

In [3]:
# extraindo os modelos, temperatura central e quantidade de zonas
model = s.get('model_number')
zones = s.get('num_zones')
ct = s.get('log_center_T')

In [12]:
import matplotlib.pyplot as plt
%matplotlib inline
plt.rcParams["figure.figsize"] = (15., 10.)

# plot do HRD
t = s.get('log_Teff')
L = s.get('log_L')

plt.plot(t, L)
plt.plot(t[154], L[154], 'r.')
plt.plot(t[2149], L[2149], 'go')
plt.gca().invert_xaxis()
plt.title(r'HRD $\alpha_{MLT}=2$', fontsize=18)
plt.xlabel(r'$\log T_{Eff}$', fontsize=16)
plt.ylabel(r'$\log L$', fontsize=16)


Out[12]:
<matplotlib.text.Text at 0x7f0b7815e190>

In [14]:
# grafico de zonas vs model

plt.plot(model, zones, 'b.')
plt.plot(model[154], zones[154], 'ro')
plt.plot(model[2149], zones[2149], 'ro')


Out[14]:
[<matplotlib.lines.Line2D at 0x7f0b788243d0>]

In [17]:
# zoom do HRD
plt.plot(t, L)
plt.plot(t[154], L[154], 'r.')
plt.plot(t[2149], L[2149], 'go')
plt.gca().invert_xaxis()
plt.title(r'HRD $\alpha_{MLT}=2$', fontsize=18)
plt.xlabel(r'$\log T_{Eff}$', fontsize=16)
plt.ylabel(r'$\log L$', fontsize=16)
plt.xlim(4, 3.5)
#plt.ylim(3.5,4.5)
plt.ylim(4,4.4)


Out[17]:
(4, 4.4)

In [ ]: