Пример #1
0
__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)
Пример #2
0
__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()
Пример #3
0
# ]
# 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)