Ondrea Fall/17 Questions for Christian: which version of matplotlib was used? Whatever is on helix python 2.7 nugridpy from svn section 2: what is hbb_table.txt,massless_grid.txt,hbb_table.txt, isotope_yield_table_MESAonly_4M.txt and isotope_yield_table_MESAonly_4M.txt used for and are the scripts to make them here? hbb_table.txt -> Combonation of yeild tables. Does not make sense to reproduce. isotope_yield_table_*txt -> Uses NuPyCEE script Nugridset.py function to make yeild tables then reads them in massloss_grid.txt -> mass, metallicity,eta parameter produced by plot_spline_fit. section 3: model_properties_xxx.txt files. What are they and how are they made? *ipython Before producing any plots: %pylab In [1]: cd /apod2/NuGrid/scratch/set1ext/paper_figures/section2/ In [2]: ls *.py calculate_massloss.py* hbb_co_ratio2.py* hbb_surface_abundance.py* hdup_kip.py* massloss_agb_set1_5a.py* calc_yields_no_hbb.py* hbb_co_ratio.py* hbb_timescales.py* HDUP_prof.py* plot_spline_fit.py* In [3]: import sys In [4]: sys.path.append('/apod2/NuGrid/scratch/set1ext/paper_figures/') Fig 1: In [5]: ed plot_spline_fit.py In [6]: run massloss_agb_set1_5a.py -produces plot. user must zoom in to desired location and xaxis is inverted relative to what is in the paper *********************************************************************** Fig 2a:In [7]: ed hbb_timescales.py -works as expected Fig 2b: In [8]: ed hbb_co_ratio.py (as does hbb_co_ratio2.py) leads to: /apod2/NuGrid/scratch/set1ext/paper_figures/nugrid_set.pyc in set_plot_CO(self, fig, xaxis, startfirstTP, t0_model, xtime, label, symbs_1, linestyle, markevery, color, end_model, withoutZlabel, savefig, singlefigures) 11017 #age=case.get("star_age")[t0_model[i]:t1_model] 11018 if singlefigures: > 11019 figure(i+int(fig)+1) 11020 plot(x,CO,label=label[i],marker=symbs[i],color=color[i],linestyle=linestyle[i],markevery=markevery) 11021 legend(loc=2) NameError: global name 'figure' is not defined IndexError: list index out of range Fig 2c:In [9]: ed hbb_surface_abundance.py -works as expected ****************************************************************** Fig 3a: In [10]: ed hdup_kip.py -works as expected Fig 3b: In [10]: ed HDUP_prof.py -works as expected. have to zoom to desired bounds ****************************************************************** Fig 4: In [11]: cd /apod2/NuGrid/scratch/set1ext/paper_figures/section3/ In [12]: ls *.py co_intershell_surf.py* kip_conts_agb.py* plot_T_rho.py* co_ratio.py* kipp_conts_massive.py* pre_exp_density_profiles_alltogether.py* evol_T_rho.py* kipp_conts.py* pre_exp_infall_velocity.py* hrds_agb.py* kippenhahn.py* pre_exp_temperature_profiles.py* hrds_massive.py* kipp_final.py* T_CEB_max_mcore.py* hrds.py* model_prop_tab.py* T_rho_center_agb.py* IFMR.py* plot_core_masses.py* T_rho_center_massive.py* kip_cont_1Msun_ingestion.py* plot_DUPlambda.py* T_rho_center.py* Fig 4a:In [13]: run hrds_agb.py leads to.... TypeError Traceback (most recent call last) /apod2/NuGrid/scratch/set1ext/paper_figures/section3/hrds_agb.py in () 14 linestyle=['--','-','-.',':','-','-.','--',':'] 15 symbs=['o','s','D','+','<','*','p','x'] ---> 16 mesaset.set_plot_hrd(fig=0,symbs_1=symbs,linestyle=linestyle,start_model=[800,650],end_model=20*[-1],markevery=10000,single_plot=True) 17 18 /apod2/NuGrid/scratch/set1ext/paper_figures/nugrid_set.py in set_plot_hrd(self, fig, symbs_1, linestyle, markevery, start_model, end_model, single_plot, labelmassonly) 10477 t0_model=case.get("model_number")[0] 10478 logTeff=case.get('log_Teff')[t00_model:(t1_model-t0_model)] > 10479 logL=case.get('log_L')[t00_model:(t1_model-t0_model)] 10480 10481 label=self.run_label[i] TypeError: slice indices must be integers or None or have an __index__ method Fig 4b:'' Fig 4c: '' *********************************************************** Fig 5a -> core_luminosity_relation.py Fig 5b In [13]: ed T_CEB_max_mcore.py -works *********************************************************** CW from top right Fig 6a: In [11]: ed kipp_final.py -works as expected assuming all four panels are made from zoom-ins of this *********************************************************** Fig 7a: In [13]: run co_ratio.py NameError: global name 'figure' is not defined Fig 7b: In [12]: ed co_intershell_surf.py IndexError: only integers, slices (`:`), ellipsis (`...`), numpy.newaxis (`None`) and integer or boolean arrays are valid indices *********************************************************** Fig 8: In [16]: ed plot_DUPlambda.py -works as expected *********************************************************** Fig 9: In [17]: ed IFMR.py -works as expected *********************************************************** Fig 10:In [19]: ed plot_core_masses.py -plots all four panels in fig 10 as expected *********************************************************** Fig 11 : In [32]: ed kipp_conts_massive.py Select either 1e-3 or 1e-4 in the script. Works as expected *********************************************************** Fig 12: In [9]: ed T_rho_center_massive.py NameError: global name 'figure' is not defined *********************************************************** Fig 13: In [35]: ed pre_exp_density_profiles_alltogether.py -works as expected. Creates the 4 panel plot. *********************************************************** Fig 14: In [37] ed pre_ex_infall_velocity.py -works as expected. Creates the 4 panel plot. *********************************************************** Fig 15: In [20]: ed kip _cont_1Msun_ingestion.py -as expected *********************************************************** cd ../section4 *********************************************************** Fig 16: In [50]: ed prodfac_vsZ_all_AGB.py both 16 and 17 reuturn the error below Fig 17:In [24]: ed prodfac_vsZ_low_agb.py TypeError Traceback (most recent call last) /apod2/NuGrid/scratch/set1ext/paper_figures/section4/prodfac_vsZ_low_agb.py in () 50 figname='all_prodfac_vsZ_AGB_low_set1_2_mesa' 51 ---> 52 setse.get_stellar_ejecta(fig=figname,cycles=cycles,isotopes=isotopes,elements=elements,fixed_ini_abu=fixed_ini_abu,yields_output=False,GCE_tables=False,plots_singleZ=True,exp_dir=expdir,exp_type=exp_type,pre_exp=pre_exp,exp_only=exp_only,wind_only=wind_only,plot_masses_only=plot_masses_only,plot_log=plot_log,color=color,marker_type=marker_type,line_style=line_style,markersize=len(isotopes)*[8],line_width=len(isotopes)*[8],title='',label='exp, delay',withlabel=False) 53 54 fig=plt.gcf() /apod2/NuGrid/scratch/set1ext/paper_figures/nugrid_set.pyc in get_stellar_ejecta(self, cycles, isotopes, elements, exclude_species, yields_output, fixed_ini_abu, file_name, GCE_tables, plots, plots_singleZ, plots_singleA, exp_dir, exp_type, fallback_masses, fallback_coords, pre_exp, exp_only, wind_only, plot_masses_only, exp_cycle_inp, plot_log, color, marker_type, line_style, markersize, line_width, title, withlabel, withlabelelement, label, plot_lines, fig, labelZonly, label_isotopes, plot_o16_lines) 7241 case_label_and_first_column.append(str(star_mass)+" Msun, Z="+str(star_z)) 7242 star_mass_array.append(star_mass) -> 7243 prod_factor,isotopes_prod_fac,yields,iniabu,mass_frac_ini,remn_mass =self.weighted_yields(sefiles[k],sefiles_restart[k],isotopes,elements,cycs) 7244 7245 if len(fixed_ini_abu)>0: /apod2/NuGrid/scratch/set1ext/paper_figures/nugrid_set.pyc in weighted_yields(self, sefiles, sefiles_restart, isotopes, elements, cycs) 8495 E=[] 8496 for k in range(len(cycs)): -> 8497 sefiles.average_iso_abund_marco(mass_range=None,cycle=cycs[k],stable=True,i_decay=2,iolevel=0) 8498 E.append(mp.average_mass_frac_decay) 8499 #isotope name - index name list TypeError: average_iso_abund_marco() got an unexpected keyword argument 'iolevel' Problem with all overabundance plots -> first try removing io level from nugridset.py try with svn version *********************************************************** Fig 25: In [52]: ed npr_profiles.py -makes two plots as expected *********************************************************** Fig 26: In [29]: ed yield_ratio_exp_preexp_Zdependence_eles.py Editing... done. Executing edited code... --------------------------------------------------------------------------- IOError Traceback (most recent call last) /apod2/NuGrid/scratch/set1ext/paper_figures/section4/yield_ratio_exp_preexp_Zdependence_eles.py in () 38 39 M=mass ---> 40 tablepre=ry.read_nugrid_yields(file_pre_exp) 41 table=ry.read_nugrid_yields(file_exp) 42 y_pre=tablepre.get(M=M,Z=Z,quantity='Yields') /apod2/NuGrid/scratch/set1ext/paper_figures/section4/read_yields.py in __init__(self, nugridtable) 57 os.system('sudo python delete.py '+nugridtable) 58 else: ---> 59 file1=open(nugridtable) 60 lines=file1.readlines() 61 file1.close() IOError: [Errno 2] No such file or directory: '../../Decay//element_yield_table_MESAonly_fryer12_delay_pre_exp.txt' Just calculates ratios of yield tables 'read_yeilds.py' from NuPyCEE for updated version *********************************************************** SECTION 5: What Fig 32?: In [41]: ed compare_gridres_c13.py Editing... done. Executing edited code... --------------------------------------------------------------------------- IOError Traceback (most recent call last) /apod2/NuGrid/scratch/set1ext/paper_figures/section5/compare_gridres_c13.py in () 3 path='/nfs/rpod2/home/critter/set1extension/AGB_stars/main_sprocess/GridSize_c13pocket/Comparison_yields' 4 ----> 5 y=ry.read_nugrid_yields(path+'/element_yield_table_set1ext.txt') 6 yhigh=ry.read_nugrid_yields(path+'/element_yield_table_highres.txt') 7 /apod2/NuGrid/scratch/set1ext/paper_figures/section4/read_yields.py in __init__(self, nugridtable) 57 os.system('sudo python delete.py '+nugridtable) 58 else: ---> 59 file1=open(nugridtable) 60 lines=file1.readlines() 61 file1.close() IOError: [Errno 2] No such file or directory: '/nfs/rpod2/home/critter/set1extension/AGB_stars/main_sprocess/GridSize_c13pocket/Comparison_yields/element_yield_table_set1ext.txt' *********************************************************** Fig 33: In [43]: ed compare_yields_COissue.py IOError: [Errno 2] No such file or directory: '/nfs/rpod2/home/critter/set1extension/MassiveStars/CO_core_issue/yield_comparison/exp/element_yield_table_MESAonly_fryer12_delay.txt' *********************************************************** Fig 34: In [44]: ed kip_cont.py works