def make_thing(): """ makes example PintPars object """ cwd = os.getcwd() main_dir = cwd + '/SPD' try: import new_lj_thellier_gui_spd as tgs gui = tgs.Arai_GUI('/magic_measurements.txt', main_dir) specimens = list(gui.Data.keys()) thing = PintPars(gui.Data, '0238x6011044', 473., 623.) thing.calculate_all_statistics() #new_thing = PintPars(gui.Data, '0238x5721062', 100. + 273., 525. + 273.) #new_thing.calculate_all_statistics() gui2 = tgs.Arai_GUI('/magic_measurements.txt', '/Users/nebula/Desktop/MagIC_experiments/ODP-SBG_1') thing2 = PintPars(gui2.Data, '0335x1031411', 273., 743.) return thing, thing2 except Exception as ex: print('could not make standard specimen objects') print(ex)
import spd # K temps: [0.0, 100.0, 150.0, 200.0, 225.0, 250.0, 275.0, 300.0, 325.0, 350.0, 375.0, 400.0, 425.0, 450.0, 475.0, 500.0, 525.0, 550.0] # C temps: [273, 373.0, 423.0, 473.0, 498.0, 523.0, 548.0, 573.0, 598.0, 623.0, 648.0, 673.0, 698.0, 723.0, 748.0, 773.0, 798.0, 823.0] import new_lj_thellier_gui_spd as tgs cwd = os.getcwd() main_dir = cwd + '/SPD' calculate = ['int_n', 'frac', 'fvds', 'b_sigma', 'b_beta', 'scat', 'g', 'k', 'k_sse', 'z', 'int_mad_anc', 'int_dang', 'int_alpha', 'alpha_prime', 'theta', 'gamma', 'int_ptrm_n', 'ptrm', 'drat', 'mdrat', 'maxdev', 'dpal', 'md', 'tail_drat', 'dtr', 'dac', 'DANG'] #calculate = ['int_n', 'frac', 'f', 'fvds', 'b_sigma', 'b_beta', 'scat', 'g', 'k', 'k_sse', 'z', 'z_md', 'q', 'r_sq', 'coeff_det_sq', 'int_mad', 'int_mad_anc', 'int_dang', 'int_alpha', 'alpha_prime', 'theta', 'int_crm', 'gamma', 'int_ptrm_n', 'ptrm', 'drat', 'drats', 'cdrat', 'mdrat', 'dck', 'maxdev', 'mdev', 'dpal', 'int_ptrm_tail_n', 'md', 'tail_drat', 'dtr', 'dt', 'ac_n', 'dac', 'gmax'] #calculate = ['int_n', 'int_alpha', 'f', 'k', 'drats', 'int_ptrm_tail_n'] #calculate = ['drats'] gui = tgs.Arai_GUI('/magic_measurements.txt', main_dir) specimens = gui.Data.keys() example = spd.PintPars(gui.Data, '0238x6011044', 473., 623., 'magic', calculate) example.calculate_all_statistics() PintPars_example = example def make_specimens(calculate=calculate): for stat in calculate: spec = spd.PintPars(gui.Data, '0238x6011044', 473., 623., 'magic', [stat]) spec.reqd_stats() print '---------' print calculate def many_specimens(calculate=calculate): from itertools import combinations