__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)