__author__ = 'mike' from RockPyV3.Structure.sample import Sample from RockPyV3.Plots.general import Henkel_Plot data = '../test_data/FeNi06Eb006-GC-140624.irm' S = Sample(name='irm_test') M = S.add_measurement(mtype='irm', mfile=data, machine='vsm') M = S.add_measurement(mtype='coe', mfile=data, machine='vsm') Henkel_Plot(sample_list=S)
__author__ = 'mike' from RockPyV3.Structure.sample import Sample # test_data = '/Users/mike/Google Drive/__code/RockPyV3/test_data/TEST.hys' # test_data = '../test_data/TEST.forc' test_data = '../test_data/FeNi06-F-d-24-M001-002.forc' S = Sample(name='test') F = S.add_measurement(mtype='forc', mfile=test_data, machine='VSM') # print F.return_fitting_surface[4] # F.calculate_forc(SF=2) # print F.return_fitting_surface.shape # F.plt_forc() # F.plt_residuals() # F.plt_drift_moment_sat() # F.plt_drift_field_sat() # F.plt_forcs() # F.plt_forc_field_spacing() # F.plt_ha_hb_space() F.plt_backfield() # F.plt_hysteresis()
# ] # x_front_files = ['/Users/mike/Google Drive/__code/RockPyV3/test_data/X-front/LF4c-6c-IRM-aniso-X_front.hys.000', # '/Users/mike/Google Drive/__code/RockPyV3/test_data/X-front/LF4c-6c-IRM-aniso-X_front.hys.001', # '/Users/mike/Google Drive/__code/RockPyV3/test_data/X-front/LF4c-6c-IRM-aniso-X_front.hys.002', # '/Users/mike/Google Drive/__code/RockPyV3/test_data/X-front/LF4c-6c-IRM-aniso-X_front.hys.003', # '/Users/mike/Google Drive/__code/RockPyV3/test_data/X-front/LF4c-6c-IRM-aniso-X_front.hys.004', # '/Users/mike/Google Drive/__code/RockPyV3/test_data/X-front/LF4c-6c-IRM-aniso-X_front.hys.005', # '/Users/mike/Google Drive/__code/RockPyV3/test_data/X-front/LF4c-6c-IRM-aniso-X_front.hys.006', # '/Users/mike/Google Drive/__code/RockPyV3/test_data/X-front/LF4c-6c-IRM-aniso-X_front.hys.007', # '/Users/mike/Google Drive/__code/RockPyV3/test_data/X-front/LF4c-6c-IRM-aniso-X_front.hys.008', # '/Users/mike/Google Drive/__code/RockPyV3/test_data/X-front/LF4c-6c-IRM-aniso-X_front.hys.009', # '/Users/mike/Google Drive/__code/RockPyV3/test_data/X-front/LF4c-6c-IRM-aniso-X_front.hys.010', # '/Users/mike/Google Drive/__code/RockPyV3/test_data/X-front/LF4c-6c-IRM-aniso-X_front.hys.011', # '/Users/mike/Google Drive/__code/RockPyV3/test_data/X-front/LF4c-6c-IRM-aniso-X_front.hys.012', # ] sample = Sample(name='test_sample') for i in x_up_files: hys = sample.add_measurement(mtype='hys', mfile=i, machine='vsm') hys.paramag_slope_correction() hys.limit_fields(1.5) # for i in x_left_files: # hys = sample.add_measurement(mtype='hys', mfile=i, machine='vsm') # hys.paramag_slope_correction() # for i in x_front_files: # hys = sample.add_measurement(mtype='hys', mfile=i, machine='vsm') # hys.paramag_slope_correction() # general.Hys_Fabian2003(sample, plot='show', plt_opt={'ls':':'}) general.Hysteresis(sample, plot='show', plt_opt={'ls':':'}, plt_virgin=True, plt_branches=False)