Back to site

Пакет праграмнага забеспячэння PmagPy

F.3 праграмнага пакета PmagPy

Праграмны пакет PmagPy ўяўляе сабой поўны набор праграм для paleomagnetists і рок magnetists. Яна напісана на Python і чытачу рэкамендуецца прачытаць Дадатак F.1 так што гэта не "чорны скрыню". Праблемы ў канцы кожнай кіраўніка закліканы дапамагчы чытачам даведацца, як выкарыстоўваць пакет PmagPy а таксама напісаць свой ??уласны код. Тут варта кароткае ўвядзенне.

F.3.1 Загрузка і ўстаноўка PmagPy

Вы можаце знайсці самыя апошнія размеркавання па наступнай спасылцы, каб паштэт дома PmagPy праз:

http://magician.ucsd.edu/Software/PmagPy/.

Усталюйце яго наступным напрамках, асабліва пераканацца, што каталог, у якім сцэнарыі знаходзіцца ў вашым шляху.

F.3.2 Агульная характарыстыка праграм PmagPy

PmagPy сцэнарыі працы, называючы іх у камандным радку (гл. раздзел F.1.3 ). Python сцэнарыі павінны быць змешчаныя ў каталог, які знаходзіцца ў Вашым "шлях". Каб пераканацца, што гэта было правільна зроблена, тыпу реж _cart.py-H ў камандным радку, і вы павінны атрымаць паведамленне. Калі вы атрымліваеце "каманда не знойдзена" паведамленне, неабходна выправіць свой шлях, праверыць "Усталяванне Python" старонкі на сайце праграмы. Іншы магчымай прычынай адмовы зьяўляецца тое, што так ці інакш, Python сцэнары больш не выкананы файл. Каб гэта выправіць, змяніць каталогі ў каталог са скрыптамі і ўвядзіце каманду: CHMOD + х *. ру

Для людзей, якія ненавідзяць лініі праграм каманды і аддаюць перавагу графічныя інтэрфейсы карыстальніка з меню і г.д., некаторыя з ключавых праграм для інтэрпрэтацыі палеомагнитных і рок-магнітных дадзеных, упакаваных разам у праграму пад назвай MagIC.py. Гэта можа быць выкліканы, набраўшы MagIC.py у камандным радку. MagIC.py праграма генеруе жаданага каманды для вас, таму вы не павінны навучыцца UNIX або як выкарыстоўваць камандную радок (за выключэннем прыняцця MagIC.py сама праграма). Тым не менш, некаторы ўяўленне пра тое, што адбываецца на самай справе карысна, таму што праграма MagIC.py больш абмежаваны, чым поўны спектр праграм PmagPy. Такім чынам, вось кароткае ўвядзенне, як праграмы работы PmagPy.

Усе PmagPy праграмы друку даведачнае паведамленне, калі вы тыпу: праграмы _name.py-H ў камандным радку. Многія з іх "intereactive" варыянт спрацоўвае, увёўшы праграму _name.py-я. Многія таксама дазваляюць чытаць з стандартнага ўводу і вываду. Даведка паведамленне будзе растлумачыць, як кожны канкрэтны функцый праграмы. Ёсць некаторыя агульныя рысы для опцыі каманднага радка:

  1. Камутатары ад аднаго да трох знакаў, папярэднічае '-'.
  2. '-H' перамыкач заўсёды друкуе паведамленне і '-я "дазваляе весці інтэрактыўны ўваход варыянтаў.
  3. Функцыі для каманднага радка непасрэдна пасля пераключальніка. Напрыклад:-F-F INPUT OUTPUT устанаўлівае ўваходны файл для ўводу і вываду OUTPUT.
  4. Перамыкач для ўваходных файлаў пачынаюцца з-F і-F для выходных файлаў.
  5. -SPC-Сэм-сядзець-сін-Loc з'яўляюцца перамыкачы, звязаныя з узораў, узоры сайтаў, сінтэтыкі і месцы адпаведна.
  6. Капіталізаваць перамыкачы здушыць варыянт (напрыклад,-сродкі не сярэдні, а-ў ??якім выпадку не ў сярэднім).
  7. CRD-[S, G, T] задае сістэму каардынат і FMT-[SVG, PNG, JPG] фармат малюнка па змаўчанні.

PmagPy сцэнары тэлефанаваць на два спецыяльных модуляў, pmag і pmagplotlib модуляў. Яны ўтрымоўваюць большасць разлікаў і пабудовы функцый.

F.3.3 Прыклады таго, як выкарыстоўваць PmagPy праграм

Ва ўсіх прыкладах, падкажуць варта '%' па нейкай каманднага радку ў вас ёсць. Спампаваць пакет, які змяшчае напрыклад файлы дадзеных для гэтай кнігі ад:

magician.ucsd.edu/Праграмнае забеспячэнне/PmagPy/Datafiles.zip

і распакаваць файл. Дадзеныя файлы для наступных прыкладаў можна знайсці ў каталогу прыкладаў.

aarm _magic.py

[Кіраўнік 13 & Magic: гл. Дадатак E ]

Анізатрапіі anhysteretic ці іншых астаткавага могуць быць ператвораныя ў тэнзарнае і выкарыстоўваецца для карэкціроўкі прыродных дадзеных астаткавага для эфекту анізатрапіі набыцця астаткавага. Напрыклад, кірункі можа быць адхіленая ад геамагнітныя поля напрамкаў і інтэнсіўнасці можа быць прадузятым моцнай анізатрапіі ў магнітных тканіны ўзору. Па надання anhysteretic або цеплавой астаткавага ў многіх канкрэтных напрамкаў, анізатрапіі астаткавага намагнічанасці набыцця могуць быць ахарактарызаваны і выкарыстоўваецца для карэкцыі. Мы робім гэта для анізатрапіі anhysteretic астаткавага (aarm) шляхам надання ARM ў 9, 12 ці 15 пазіцый. Кожны ARM павінен папярэднічаць этап размагничивания AF. 15 пасад прыведзены ў дадатку D.1 . За 9 пазіцый схемы, AARMs прышчапляюцца ў пазіцыі 1,2,3, 6,7,8, 11,12,13, напрыклад. Кто-то ласкава вымярэнняў, ператварыла іх у фармат magic_measurements (гл. Дадатак E ) І пасяліў іх у файл: aarm _measurements.txt. [Адзін са спосабаў зрабіць гэта для сябе, для выкарыстання праграмы магнітнага _magic.py.]

Першы распакаваць thellier _example.zip файл _redo і змены ў каталог thellier_redo_example. Затым з дапамогай праграмы магнітнага _magic.py імпартаваць AARM дадзеных у фармаце Magic. Поле пастаяннага току складала 50 ? T, AC поле пік 180 МТЛ, месца было "Bushveld і лабараторыі пратакол быў А. Ф. і анізатрапіі. Наймення выкарыстоўваецца Варыянт № 3 (гл. меню "Даведка").

Затым з дапамогай праграмы aarm _magic.py для вылічэння найлепшага тэнзар і выпісаць Magic табліцы: rmag _anisotropy і rmag _results. Гэтыя файлы могуць выкарыстоўвацца для карэкцыі астаткавага інфармацыі ў фармаце pmag _specimens табліцы (напрыклад інтэнсіўнасць дадзеных) для эфекту астаткавага анізатрапіі (прыкладам, з выкарыстаннем праграмы thellier _magic _redo.py.)

Вось стэнаграма сесіі, якая працуе. Адзначым, што магнітнае _magic.py каманды ўсё ў адной радку. Акрамя таго, не турбавацца аб Увага, гэта ўсяго толькі "асаблівасць" некаторых магнітных файлаў у фармаце:

%mag_magic.py -f bg.arm -loc Bushveld -LP AF:ANI -F aarm_measurements.txt  
    -ncn 3 -ac 180 -dc 50 0 90  
Warning - inconsistency in mag file with lab field d  
           - overriding file with 0  
.....  
% aarm_magic.py -f aarm_measurements.txt  
Processing:  bg2.01  Number of positions:  9  
Processing:  bg2.03  Number of positions:  9  
....  
specimen tensor elements stored in  ./rmag_anisotropy.txt  
specimen statistics and eigenparameters stored in  ./rmag_results.txt

agm_magic.py

[Chapter  5 & MagIC: see Appendix E]

Гэтая праграма імпарту Micromag гістарэзісу файлы ў magic_measurements файлаў у фармаце. Таму што гэтая праграма імпарту дадзеных у базу дадзеных Magic, узоры неабходна таксама мець ўзор/сайтаў/размяшчэнне інфармацыі, якая можа быць прадастаўлена ў камандным радку. Калі гэтая інфармацыя недаступная, напрыклад, калі гэта сінтэтычных ўзору, пакажыце-сін для сінтэтычных у камандным радку.

Кто-то імем Ліма Танга падлічаных сінтэтычных ўзору імя myspec для гістарэзісу і захаваць дадзеныя ў файл з імем AGM _example.dat. Выкарыстоўвайце AGM _magic.py праграма для імпарту дадзеных у magic_measurements фарматаваны выснову файлаў. Гэтыя вымяраліся СГС, так што для ўстаноўкі камутатараў правільна. [Гэта можа быць пабудавана пры выкарыстанні гістарэзісу _magic.py.] Вы таксама можаце імпартаваць набыцця IRM і пастаяннага току размагничивания крывых з дапамогай IRM _magic.py. Гістарэзісу _magic.py таксама вылічае розныя параметры гістарэзісу. ]

%  agm_magic.py -spc myspec -syn -usr "Lima Tango" -f agm_example.dat  
        -F myspec_agm.magic -u cgs  
specimen name put in  er_specimens.txt  
agm_magic_example.dat  appears to be a dos file  
results put in  myspec_agm.magic

AGM _magic.py вызначылі гэта як DOS фарматаванага і разгледжаны яго належным чынам. Гэта выратавала дадзеных у файл з імем myspec _agm.magic, якія могуць быць аб'яднаны з іншымі magic_measurements адфарматаваных файлы, выкарыстоўваючы аб'яднаць _magic.py. Глядзіце таксама, як імпартаваць IRM і заднюю дадзеных у паведамленне.

Адзначым, што новая праграма Micromag мае трохі іншы фармат са старых праграм. Каб імпартаваць новы фармат файлаў дадзеных, выкарыстайце наступны сінтаксіс:

% agm_magic.py -f new_agm_example.dat -spc myspec -u cgs -new

angle.py

[Appendix A.3.2]

Выкарыстоўвайце праграму для разліку кута кут (?) паміж двума кірункамі D = 350 8. 2, I = 26;. 98. 6 D, I = 67. = 3.

% angle.py -i  
Declination 1: [cntrl-D  to quit] 350.2  
Inclination 1: 26.8  
Declination 2: 98.6  
Inclination 2: 67.3  
   72.1  
Declination 1: [cntrl-D  to quit] ^D  
Good bye  [NB: PC users will get a more angry sounding exit message]

анізатрапіі _magic.py

[Кіраўнік 13 & Magic: гл. Дадатак E ]

Узоры былі сабраныя ад усходняй ускраіны дамбы арыентаваны ложках полюса скланенне 110 ? і падзенні 2 ?. Дадзеныя былі імпартаваныя ў rmag_anisotropy адфарматаваную файл з імем дамбы _anisotropy.txt.

Зрабіць ўчастак даных з выкарыстаннем анізатропнага _magic.py. Выкарыстаньне сайта параметрычнага варыянт загрузкі і зямельны ўчастак з пачатковай загрузкі уласных вектараў. Намалюйце на след дамбы.

Гэтыя рэчы робяцца ў гэтай сесіі:

% aniso_magic.py -f dike_anisotropy.txt -gtc 110 2 -par -v  
plotting  6  points  
Doing bootstrap - be patient  
plotting bootstrapped eigenvectors  
Bootstrap Statistics:  
tau_i, V_D, V_I, V_zeta, V_zeta_D, V_zeta_I, V_eta, V_eta_D, V_eta_I  
0.34040    29.6    14.5    29.0    165.5    66.4     7.2   295.4   15.7  
0.33536   166.3    70.5    25.2     20.4    16.7    11.5   287.4    9.9  
0.32423   296.2    12.8    12.4    131.6    76.5     6.6    27.1    3.4  
compare with [d]irection  
 plot [g]reat circle,  change [c]oord. system, change [e]llipse  
    calculation,   s[a]ve plots, [q]uit  

якая вырабіла гэтыя ўчасткі:

EPSfiles/dike.eps

Накапіцельныя размеркавання пачатковай загрузкі уласных Паказана, што засталося 95% даверны мяжы нанесеныя ў выглядзе вертыкальных ліній. Ўзору ўласных пабудаваны ў сярэдзіне дыяграмы з звычайным канвенцыі, што квадраты V 1 напрамках, трыкутнікі V 2-х напрамках і колы V 3-х напрамках. Усе напрамкі будуць выводзіцца на ніжняй паўсферы. Пачатковую загрузку ўласных паказаны з правага боку.

Падобна на тое, што магма рухаецца ў паўночным і трохі ўверх, уздоўж дамбы.

apwp.py

[Кіраўнік 16 ]

Праграма вылічае apwp.py paleolatitude, скланенне, лад ад полюса шыраты і даўгаты на аснове паперы і Бэс Courtillot (2002, гл 16 для вядзення і поўнага абмеркавання). З яго дапамогай можна вылічыць чаканае кірунак на 100 мільёнаў гадоў парод на мясцовасці ў Ла Джолі ков (шырата: 33N, даўгата 117W). Выкажам здагадку, што мы знаходзімся на паўночнаамерыканскім пліты! (Майце на ўвазе, што няма магчымасці для Ціхаакіянскай пліты ў праграме apwp.py, і што Ла Хоя быў ??на паўночнаамерыканскім пласціну да некалькіх мільёнаў гадоў назад (6?).

% apwp.py -i  
 
 Welcome to paleolatitude calculator - CTRL-D to quit  
 
pick a plate: NA, SA, AF, IN, EU, AU, ANT, GL  
 
Plate  
NA  
Site latitude  
33  
 Site longitude  
-117  
 Age  
100  
Age  Paleolat.  Dec.  Inc.  Pole_lat.  Pole_Long.  
   100.0    38.8   352.4    58.1    81.5   198.3  

Адзначым, што як і ў многіх праграмах PmagPy, увод інфармацыі можа быць прачытана з файла і выхад можа быць змешчаны ў файл (гл. опцыю-H для даведкі).

azdip_magic.py

[Кіраўнік 9 і Magic: гл. Дадатак E ]

Многія paleomagnetists захаваць арыентацыю інфармацыі ў файлы ў гэтым фармаце:

Ўзору Азімут Падзенне Strike Падзенне
дзе азімуце і акунь з'яўляюцца скланенне і лад кірунак свідравання і простиранию і падзення з'яўляюцца стасункі пробы адзінку (з апусканнем на права ўдару). Базы дадзеных канвенцыі Magic з'яўляецца выкарыстанне кірунку X каардынаты ўзору сістэмы вымярэння. Каб пераўтварыць файл у фармаце AzDip (example.az) для узораў, узятых з месца, назва "Паўночная Ісландыя" у фармаце Magic і захоўваць інфармацыю ў фармаце ER Magic файл _samples.txt, выкарыстанне azdip _magic.py праграмы:
% azdip_magic.py -f example.az -loc "Northern Iceland"  
Data saved in  er_samples.txt

Звярніце ўвагу, што Ёсць шмат варыянтаў, якія тычацца ўзору імёны імёны сайтаў. Глядзіце меню "Даведка" для падрабязнасьцяў.

У _vdm.py

[Кіраўнік 2 ]

Выкарыстоўвайце праграму б _vdm пераўтварыць ацэначны кошт paleofield з 33 ? T атрымліваецца з лававых патоку на 22 ? з.ш. да эквівалента VDM ў Am 2. Пакладзеце ўводзіць інфармацыю ў файл пад назвай _input.dat VDM і чытаюцца з яго з выкарыстаннем стандартнага ўводу:

% cat > vdm_input.dat  
33 22  
^D  
% b_vdm.py < vdm_input.dat  
 7.145e+22

базавая карта _magic.py

[Magic: гл. Дадатак E ]

Выкарыстаньне праграмы базавай картай _magic.py зрабіць просты сюжэт базавая карта з сайта месцаў у ER _sites.txt Magic адфарматаваную файл з імем базавая _example.txt.

%basemap_magic.py -f basemap_example.txt

што робіць гэты ўчастак:

EPSfiles/basemap.eps

Выкарыстоўвайце кнопкі ў ніжняй частцы ўчастка, каб змяніць памер або захаваць сюжэт у жаданым фармаце.

biplot _magic.py

[Кіраўнік 8 і Magic: гл. Дадатак E ]

Дадзеныя былі атрыманы з асноўных палеогена з 28 ? S для даследавання адносна палеонапряженности. магнітнай успрымальнасці IRM, ARM, і астаткавага інфармацыі былі загружаны ў базу дадзеных Magic. Magic_measurements адфарматаваную файл для гэтага даследавання захоўваецца ў асноўных _measurements.txt.

Выкарыстаньне праграмы biplot _magic.py зрабіць biplot магнітнай успрымальнасці ад ARM. Заўважым, што праграма выкарыстоўвае метад коды Magic (гл. Дадатак E ), Якія LT-IRM для IRM, LT-AF-I для ARM (AF размагничивания, у полі), і LP-X для магнітнай успрымальнасці.

Па-першае, высветліць, якія дадзеныя маюцца, запусціце праграму, як гэта:

% biplot_magic.py -f core\_measurements.txt  
   which responds with this:  
   [’LT-AF-Z’, ’LT-AF-I’, ’LT-IRM’, ’LP-X’]  
   

Да іх адносяцца: А. Ф. размагничивания NRM, ARM, IRM і успрымальнасці вымярэнняў адпаведна. Таму, каб зрабіць сюжэт пра успрымальнасці супраць ARM, мы б зноў запусціце праграму:

% biplot_magic.py -f core_measurements.txt -x LP-X -y LT-AF-I  
LP-X  selected for X axis  
LT-AF-I  selected for Y axis

што робіць сюжэт:

EPSfiles/ARM-x.eps

bootams.py

[Кіраўнік 13 ]

Праграма вылічае bootams.py загрузкі Статыстыка анізатрапіі тэнзар дадзеных у выглядзе:

x11 x22 x33 x12 x23 x13

Ён робіць гэта, абраўшы пункт-набораў дадзеных і разліку eigenparameters Hext сярэднім. Гэта дадатковы параметрычны загрузкі якога ? для набору дадзеных у цэлым выкарыстоўваецца для малявання новы пункт набораў дадзеных. Бутстрепа eigenparameters лічацца Кент размеркаваных і праграма вылічае эліпсы Кент памылкі для кожнага набору вектараў. Ён таксама ацэнкі стандартных адхіленняў пачатковай загрузкі уласных.

Выкарыстоўвайце гэта для разліку пачатковай загрузкі статыстыку памылак для дадзеных у файл _examples.data загрузкі:

% bootams.py -par -f bootams_example.dat  
Doing bootstrap - be patient  
 
tau sigma V_dec V_inc V_eta V_eta_dec V_eta_inc V_zeta V_zeta_dec V_zeta_inc  
 
0.33505 0.00029     5.3   14.7    15.6  267.0   27.3   26.7   123.2   57.3  
0.33334 0.00030   124.5   61.7    19.5  230.1    8.4   25.5   324.3   26.1  
0.33160 0.00025   268.8   23.6    14.7    5.9   17.4   21.9   128.7   60.0  

Кошык _dir.py

[Кіраўнік 2 ]

Выкарыстоўвайце кошык _dir праграма для пераўтварэння гэтых дэкартавай каардынаты геамагнітныя элементаў:




х 1 х 2 х 3



0.3971 -0.1445 0.9063
-0.5722 0.0400 -0.8192



Каб выкарыстоўваць інтэрактыўныя опцыі:

% cart_dir.py -i  
X: [cntrl-D  to quit] 0.3971  
Y: -0.1445  
Z: 0.9063  
  340.0    65.0  1.000e+00  
X: [cntrl-D  to quit] -.5722  
Y: 0.0400  
Z: -0.8192  
  176.0   -55.0  1.000e+00

To read from a file:

% cart_dir.py -f cart_dir_example.dat  
  340.0    65.0  1.000e+00  
  176.0   -55.0  1.000e+00

chi_magic.py

[Кіраўнік 8 & Magic: гл. Дадатак E ]

Часам бывае карысна для вымярэння успрымальнасці у залежнасці ад тэмпературы, прыкладзенай поля і частаты. Вось набор дадзеных, які прыйшоў з паслядоўнасці Тыва Каньён туф (гл. Картэр-Стиглица, 2006). Выкарыстаньне праграмы хі _magic.py пабудаваць дадзеныя ў файл фармату magic_measurements: core _measurements.txt.

% chi_magic.py -f chi_measurements.txt  
Input magic_measurements file name? [magic_measurements.txt]  chi_measurements.txt  
IRM-Kappa-2352 1 out of  2  
 
  **** not enough data to plot ***  
Skipping susceptibitily - AC field plot asa function of temperature  
Return to save files, [q] to quit  
 

вытворчасці гэтага ўчастка:

EPSfiles/хі-magic.eps

Вы можаце бачыць, залежнасць ад тэмпературы, частаты і прыкладзенай поля. Гэтыя дадзеныя пацвярджаюць меркаванне, што існуе моцная суперпарамагнитных кампанента ў гэтых узорах.

аб'яднаць _magic.py

[Magic: гл. Дадатак E ]

Magic табліцах шмат Стоўбцаў толькі некаторыя з іх выкарыстоўваюцца ў канкрэтным выпадку. Так аб'яднання файлаў таго ж тыпу, павінна быць зроблена дбайна, каб гарантаваць, што патрэбныя дадзеныя падпадаюць пад права загалоўкамі. Праграма аб'яднаць _magic.py могуць быць выкарыстаны для аб'яднання любую колькасць файлаў Magic ад дадзенага тыпу. Выкарыстоўвайце гэтую праграму, каб аб'яднаць magic_measurements адфарматаваных файлы, створаныя ў прыклады для class="cmbx-10x-x-109">agm_magic.py a file agm_measurements.txt:

 combine_magic.py -F agm_measurements.txt -f myspec_agm.magic myspec_irm.magic  
 

This file can be plotted and interpreted with the hyseresis_magic.py program.

common_mean.py

[Chapter  12]

Выкарыстаньне праграмы агульнай _mean.py выкарыстоўваць загрузачны тэст для агульнага сярэдняга і праверыць ці два скланенне, лад набораў дадзеных агульнай сярэдняй на 95% узроўні даверу. Наборы дадзеных: common_mean_ex_file1.dat and common_mean_ex_file2.dat і, у роўнай праекцыі вобласці з'яўляюцца:

EPSfiles/агульнай сярэдняй eq.eps

Вось стэнаграма сесіі, што будзе рабіць гэта:

% common_mean.py -f common_mean_ex_file1.dat -f2 common_mean_ex_file2.dat  
Doing first set of directions, please be patient..  
Doing second set of directions, please be patient..  
Return to save plots - <cntl-D> to quit  
2  saved in  CD_Y.svg  
1  saved in  CD_X.svg  
3  saved in  CD_Z.svg

Тры ўчасткі з'яўляюцца:

EPSfiles/агульнага mean.eps

Гэтыя выказаць здагадку, што два набору дадзеных долю агульнага сярэдняга.

Цяпер параўнайце дадзеныя ў агульнай _file1.dat _ex _mean з чаканым напрамку на 5 паўночнай шыраты ?, што гэтыя дадзеныя былі сабраны (снежань = 0, Inc = 9,9).

Для гэтага, па гэтай стэнаграмы:

% common_mean.py -f common_mean_ex_file1.dat -dir 0 9.9  
Doing first set of directions, please be patient..  
Return to save plots - <cntl-D> to quit  
2  saved in  CD_Y.svg  
1  saved in  CD_X.svg  
3  saved in  CD_Z.svg

EPSfiles/агульнай сярэдняй EX.eps

Мабыць дадзеных (функцыя размеркавання) цалкам адпавядаюць чаканым кірунку (пункцірныя лініі дэкартавай каардынаты).

працяг _rot.py

[Кіраўнік 16 і Дадатку A.3.5 .]

Выкарыстоўвайце праграму працяг _rot.py зрабіць артаганальнай праекцыі з шырынёй = -20 ? і даўгата = 0 ? ў цэнтры Афрыкі і Паўднёвай Амерыкі кантынентах рэканструяваны да 180000000 выкарыстаннем Torsvik і інш. (2008) палюсоў канчатковага павароту. У гэтага першага правядзення Афрыцы фіксаванай. Перамяшчэнне выхадзе ўчастка з фіксаванай _africa.svg. Затым зрабіць сюжэт для Афрыкі даводзяць да палеомагнитных сістэме адліку. Зрабіць кантынентальнага выкладаюцца ў чорных ліній і ўсталяваць дазвол на "нізкі".

% cont_rot.py -con af:sam -prj ortho -eye -20 0 -sym ’k-’ 1 -age 180 -res l  
 S[a]ve to save plot, Return to quit:  a  
1  saved in  Cont_rot.svg  
% mv Cont_rot.svg fixed_africa.svg  
% cont_rot.py -con af:sam -prj ortho -eye -20 0 -sym ’k-’ 1 -age 180 \  
     -res l -sac  
 S[a]ve to save plot, Return to quit:  a  
1  saved in  Cont_rot.svg

Гэтыя каманды спароджаных наступныя ўчасткі (першы злева, другі справа):

EPSfiles/controt.eps

curie.py

[Кіраўнік 6 ]

Выкарыстаньне праграмы curie.py інтэрпрэтаваць Кюры дадзеных тэмпературы ў _example.dat Кюры прыклад файла. Выкарыстаньне згладжвання вокны ад 10 ?.

%curie.py -f curie_example.dat -w 10  
second deriative maximum is at T=552

which generates these plots:

EPSfiles/curie-ex.eps

customize_criteria.py

[Magic: гл. Дадатак E ]

Базы дадзеных Magic дазваляе дакументацыі, якія крытэры выкарыстоўваліся пры выбары дадзеных аб ўзору, узор або мясцовым узроўні. Гэтыя выбары захоўваюцца ў табліцы pmag _criteria. Некаторыя праграмы PmagPy выкарыстоўваць файл дадзеных _criteria.txt pmag для выбару дадзеных (напрыклад _magic.py thellier і ўзоры _results _magic.py). Каб наладзіць гэтыя крытэры для яго сваімі ўласнымі наборамі даных, вы можаце выкарыстоўваць праграму наладкі _criteria.py. Гэтая праграма таксама называюць MagIC.py GUI ў меню Utilities. Паспрабуйце яго на pmag _criteria.txt. Гэта "поўны вектар" набор крытэраў, - гэта азначае, што яна мае абодвух напрамках і інтэнсіўнасць сцягі мноства. Змяненне ўзору _alpha95 абразаньня да 180. ад таго, што ў цяперашні час усталяваны. Захаваць вынік у новы файл з імем новага _criteria.txt.

% customize_criteria.py -f pmag_criteria.txt -F new_criteria.txt  
Acceptance criteria read in from  pmag_criteria.txt  
 [0] Use no acceptance criteria?  
 [1] full vector  
 [2] direction only  
 [3] intensity only  
 1  
specimen_mad 5.49  
Enter new criterion (return to keep default)  
specimen_alpha95 5.49  
Enter new criterion (return to keep default) 180.  
.  
.  
site_int_sigma_perc 15  
Enter new criterion (return to keep default)  
Customize criteria again ? 1/[0]  
Criteria saved in pmag_criteria.txt  
 
 Pmag Criteria stored in  new_criteria.txt  

dayplot _magic.py

[Кіраўнік 5 & Magic: гл. Дадатак E ]

Выкарыстоўвайце dayplot _magic.py праграму, каб зрабіць дзень, прастакутнай-коэрцитивности і прастакутнай-Коэрцитивность участкаў астаткавага намагнічанасці ад rmag_hyseresis адфарматаваных дадзеных у dayplot_example.dat.

The session:

  % dayplot_magic.py -f dayplot_example.txt  
 
 S[a]ve to save plots, return to quit:  a  
4  saved in  ./dayplot_example.txt_bcr1-bcr2.svg  
3  saved in  ./dayplot_example.txt_S-Bcr.svg  
1  saved in  ./dayplot_example.txt_day.svg  
2  saved in  ./dayplot_example.txt_S-Bc.svg

дае участкаў:

EPSfiles/dayplot.eps

di_geo.py

[Кіраўнік 9 ]

Выкарыстаньне di_geo.py праграмы для пераўтварэння D = 8. 1, I = 45. 2 у геаграфічным і нахілу рэгулюецца каардынат. Арыентацыі лабараторыі стрэлка на ўзор: азімут = 347; акунуцца = 27. Дзі _geo.py працуе ў звычайным трыма спосабамі (дадзеныя інтэрактыўнага ўваходу, каманднага радка спецыфікацыі файла ці са стандартнага ўваходнага файла. Такім чынам, захоўваць дадзеныя ў называецца di_geo _example.dat, напрыклад, з дапамогай кот функцыі Unix (гл. Дадатак F.2 :

cat > di_geo_example.dat  
8.1 45.2 347 27  
<cntr-D>

Затым увядзіце:

% di_geo.py -f di_geo_example.dat  
    5.3    71.6

які выдае наш адказ Скланенне = 5,3 і нахілу = 71,6.

di_tilt.py

[Кіраўнік 9 ]

Выкарыстоўвайце праграму di_tilt.py да ROATE кірунку Скланенне = 5,3 і схільнасць = 71,6 на "стратиграфических" каардынат. Страйк была 135 і падзенне было 21. Канвенцыі ў гэтай праграме з'яўляецца выкарыстанне кірунак падзення, якое з'яўляецца "права" гэтага ўдару.

Вось сесіі з di_tilt.py выкарыстаннем інтэрактыўных варыянт:

% di_tilt.py -i  
Declination: <cntl-D> to quit 5.3  
Inclination: 71.6  
Dip direction: 225  
Dip: 21  
  285.7    76.6  
Declination: <cntl-D> to quit ^D  
 Good-bye

di_rot.py

[Кіраўнік 12 ]

Стварыць Фішэр размеркаваных набор дадзеных з насельніцтва з сярэднім кірунку D = ??0, I = 42 з выкарыстаннем праграмы fishrot.py. Разлічыць сярэдняе кірунак набору дадзеных, выкарыстоўваючы gofish.py. Зараз з дапамогай праграмы di_rot.py круціцца мноства напрамкаў для сярэдняга напрамкі. Паглядзіце на дадзеныя да і пасля кручэньня з дапамогай eqarea.py.

% fishrot.py -I 42 >fishrot.out  
% gofish.py <fishrot.out  
    1.7    42.4    100    95.3720     21.4     3.1  
% di_rot.py -f fishrot.out -F dirot.out -D 1.7 -I 42.4  
% eqarea.py -f fishrot.out  
 % eqarea.py -f dirot.out

, Які генеруе участкаў, як гэтыя:

EPSfiles/dirot.eps

di_vgp.py

[Кіраўнік 2 ]

Выкарыстоўвайце праграму di_vgp канвертаваць наступныя:





D Я ? S (N) ? S (E)




11 63 55 13
154 -58 45.5 -73




Вось стэнаграма тыповы сеанс з дапамогай опцыі каманднага радка для запісу файла імя:

% di_vgp.py -f di_vgp_example.dat  
  154.7    77.3  
    6.6   -69.6

dipole_pinc.py

[Кіраўнік 2 ]

Разлічыце чаканы нахіл на paleolatitude ад 24 ? С.

% dipole_pinc.py -i  
Paleolat for converting to inclination: <cntl-D> to quit -24  
  -41.7  
Paleolat for converting to inclination: <cntl-D> to quit ^D  
 Good-bye  

dipole_plat.py

Разлічыць paleolatitude для сярэдняга нахілу ад 23 ?.

% dipole_plat.py -i  
Inclination for converting to paleolatitude: <cntl-D> to quit 23  
   12.0  
Inclination for converting to paleolatitude: <cntl-D> to quit ^D
Да спаткання

реж _cart.py

[Кіраўнік 2 ]

Выкарыстоўвайце реж _cart.py праграма для пераўтварэння наступныя дадзеныя з скланенне D, лад I і інтэнсіўнасць М X 1, X 2, X 3.




D Я M (? Am 2)



20 46 1.3
175 -24 4.2



Вы можаце ўвесці D, I, M дадзеныя ў файл дадзеных, затым запусціць праграму, набраўшы, што пасля падказкі (%) [іншыя рэчы кампутара адказы]:

% cat > dir_cart_example.dat  
 20 46 1.3  
175 -24 4.2  
^D  
% dir_cart.py <dir_cart_example.dat  
8.4859e-01 3.0886e-01 9.3514e-01  
-3.8223e+00 3.3441e-01 -1.7083e+00

Or you could use dir_cart.py interactively as in:

 % dir_cart.py -i  
 
Declination: [cntrl-D  to quit]  
 Good-bye  
 
% dir_cart.py -i  
Declination: [cntrl-D  to quit] 20  
Inclination: 46  
Intensity [return for unity]: 1.3  
8.4859e-01 3.0886e-01 9.3514e-01  
Declination: [cntrl-D  to quit] 175  
Inclination: -24  
Intensity [return for unity]: 4.2  
-3.8223e+00 3.3441e-01 -1.7083e+00  
Declination: [cntrl-D  to quit] ^D
Да спаткання

Прыклад 1 dmag _magic.py [Кіраўнік 9 & Magic: гл. Дадатак E ]

Выкарыстоўвайце dmag _magic.py да сюжэту згасанне ўсіх размагничивания пераменным полем эксперыментаў у magic_measurements адфарматаваны вымярэнняў _examples.txt файл. Паўтарыце гэтыя дзеянні для ўсіх цеплавых вымярэнняў, але выключыць усе дадзеныя, атрыманыя падчас эксперыментаў палеонапряженности.

Вось стэнаграма сесіі:

%  dmag_magic.py -f measurements_example.txt -LT AF  
969  records read from  measurements_example.txt  
 S[a]ve to save plot, [q]uit,  Return to continue:  a  
1  saved in  McMurdo_LT-AF-Z.svg  
% dmag_magic.py -f measurements_example.txt -LT T -XLP PI  
969  records read from  measurements_example.txt  
 S[a]ve to save plot, [q]uit,  Return to continue:  a  
1  saved in  McMurdo_LT-T-Z.svg

якая вырабіла гэтыя ўчасткі:

EPSfiles/dmag.eps

Прыклад 2 спампаваць _magic.py [Magic: гл. Дадатак E ]

Гэтая праграма распакоўвае. TXT файлаў, загружаных з базы дадзеных Magic на асобныя файлы для кожнай табліцы (напрыклад, э _locations.txt, магія _measurements.txt, pmag _results.txt і гэтак далей.) У якасці прыкладу, перайдзіце да базы дадзеных магіі http://earthref.org/MAGIC/ і націсніце на спасылку PMAG партал. Выконвайце "Пошук-даведачную базу даных" спасылку і ўвядзіце "Sbarbori" у поле Аўтар. Вы трапіце на старонку для публікацыі аль Sbarbori ET. (2008). Загрузіце апошнюю версію "смартбук". Двойчы пстрыкніце на файл sbarbori-э-аль-2008.gz.zip для распакоўкі архіва. Гэта дазволіць зрабіць папку, у якой вы павінны змяніць каталогі ў (sbarbori-э-аль-2008). Зараз з дапамогай праграмы спампаваць _magic.py для распакавання. Тэкставы файл (zmab0099659tmp01.txt).

% download_magic.py -f zmab0099659tmp01.txt  
working on:  er_locations  
er_locations  data put in  ./er_locations.txt  
working on:  er_sites  
er_sites  data put in  ./er_sites.txt  
....

Вы можаце праверыць дадзеныя, выкарыстоўваючы PmagPy праграм (напрыклад, zeq _magic.py).

Прыклад 3 EI.py [Кіраўнік 14 ]

Выкарыстоўвайце праграму для стварэння tk03.py мноства напрамкаў ад мадэлі TK03 поле ПСВ для экватарыяльных шырот. Разлічыць падаўжэнне і нахілу мадэлявання дадзеных і іх параўнанне з прадказаннямі TK03 мадэлі з выкарыстаннем праграмы EI.py.

 % tk03.py >lat0.tk03  
% EI.py -f lat0.tk03  
doing bootstrap - be patient  
100  
200  
...  
 
Io,Eo,Is[lower],Is[upper],Es[lower],Es[upper]  
    0.3 3.13     0.1     4.1 2.04 5.54  
 S[a]ve to save plot, [q]uit without saving:  

, Якая вырабляе ўчастка, як:

EPSfiles/EI-Ex.eps

eigs _s.py

[Кіраўнік 13 ]

Друк eigenparameters ў файл _s _example.dat eigs, а затым пераўтварыць іх у тэнзар дадзеных. S фармаце (x11, x22, x33, x12, x13, x23).

Гэтая сесія выкарыстоўвае UNIX коткі утыліта для друку дадзеных (гл. Дадатак F.2 ). [Вы можаце выкарыстоўваць-Dos формы тыпу Г-жа ў акне Windows лініі каманды.] Затым яна друкуе тэнзар дадзеных на экран.

% cat eigs_s_example.dat  
0.33127  239.53   44.70 0.33351  126.62   21.47 0.33521  19.03   37.54  
0.33177  281.12    6.18 0.33218  169.79   73.43 0.33603   12.82   15.32  
...  
% eigs_s.py -f eigs_s_example.dat  
0.3341630.332802 0.333034 -0.000166 0.001231 0.001355  
0.335557 0.3319740.332468 0.000856 0.000252 0.000981  
...

endnote_magic.py

[MagIC]

Гэта зручная праграма ператварае вашыя EndNote спасылка базе дадзеных BibTex фарматаваны выснова ў выглядзе табліцы Magic er_citations. Паспрабуйце гэта на файл class="cmti-10x-x-109">EndNoteExport.txt:

% endnote_magic.py  
Citations saved in  er_citations.txt

Вам прыйдзецца змяніць er_citation_name на "Работа" на ключавых спасылкай на даследаванне вы спрабуеце загрузіць.

eqarea.py

[Кіраўнік 2 і Дадатак B.1 ]

Выкарыстаньне праграмы eqarea.py пабудаваць роўныя праекцыі вобласці напрамкі ў di_example.txt. Зрабіце гэта, паказаўшы імя файла ў камандным радку:

% eqarea.py -f di_example.txt

, Якая вырабляе ўчастка:

EPSfiles/eqarea.eps

eqarea_ell.py

[Chapters  11 and  12 ]

Выкарыстаньне праграмы eqarea _ell.py пабудаваць роўныя праекцыі вобласці напрамкі ў ды _example.txt і ўпэўненасць эліпсы ўчастка. Вось прыклад для эліпса Bingham.

% eqarea_ell.py -f di_example.txt -B

eqarea_magic.py

[Дадатак E ]

Выконвайце інструкцыі па загрузцы і распакаванні файлаў дадзеных з базы дадзеных Magic. Тады ўчастак накіраванасці дадзеныя для даследавання з pmag файл _results.txt нараўне з пачатковай загрузкі эліпса даверу.

% eqarea_magic.py -obj loc -crd g -f pmag_results.txt -ell Be  
24  records read from  pmag_results.txt  
VGP: Site sc02 DE-FM:DE-DI    23.8    49.0  
VGP: Site sc04 DE-FM:DE-DI     8.8    59.2  
VGP: Site sc05 DE-FM-LP:DE-DI   345.4    28.5  
VGP: Site sc09 DE-FM:DE-DI   358.4    21.5  
VGP: Site sc12 DE-FM:DE-DI    19.3    53.6  
VGP: Site sc14 DE-FM:DE-DI    17.3    53.8  
VGP: Site sc15 DE-FM-LP:DE-DI   318.8    20.4  
VGP: Site sc17 DE-FM:DE-DI   347.3    28.3  
VGP: Site sc20 DE-FM-LP:DE-DI     2.0    32.3  
mode  1  
     Zdec   116.0  
     Edec   230.4  
     Eta     6.7  
     n        1  
     Einc    35.3  
     Zinc    30.3  
     Zeta    16.3  
     dec   356.8  
     inc    39.9  
 S[a]ve to save plot, [q]uit, Return to continue:  a  
1  saved in  Whole file: Geographic Coordinates_eqarea.svg

робіць гэты ўчастак:

EPSfiles/eqarea-magic.eps

Інфармацыю, надрукаваную ў акно pmag_result_name ў табліцы даных, метад коды (тут "скіраванасці адзнака: Фішэр маю на ўвазе, і кожны сайт сярэдняй скланення нахілу інфармацыі наступнае." Рэжыме 1 "пачатковую загрузку параметраў эліпса.

eq_di.py

[Дадатак B.1 ]

Дадзеныя часта публікуюцца як роўныя праекцыі вобласці і не пералічаных у табліцах дадзеных. Гэтыя дадзеныя могуць быць алічбаваныя, як х, у дадзеных (пры ўмове вонкавага вобада роўны адзінцы) і пераўтворыцца ў набліжаных напрамках з _di.py экв праграмы. Каб выкарыстоўваць гэтую праграму, усталяваць графік дигитайзер (Datathief III ад www.datathief.org працуе на большасці платформаў і GraphClick гэта выдатная праграма для Mac OS X 10.3 або больш позняй версіі.)

Аблічбоўка дадзеных, прадстаўленых для eqarea.py прыклад, які захоўваецца ў файле eqarea.png. Вы павінны працаваць толькі на аднаго паўшар'я у той час (верхняя або ніжняя) і захаваць кожнае паўшар'е у сваім уласным файле. Тады вы можаце канвертаваць X, Y дадзеныя набліжаныя снежня і вкл дадзеных - якасць дадзеных залежыць ад вашай дапамогі ў аблічбоўцы і якасць малюнка відаць, што вы лічбавання.

Паспрабуйце экв _di.py на ваш файл дадзеных, або выкарыстоўваць ўраўненні _di _example.dat якія алічбаваных дадзеных з ніжняй паўсферы і праверыць працу з eqarea.py. Вы павінны атрымаць ніжнія пункту паўшар'я з eqarea.py прыклад.

% eq_di.py -f eq_di_example.dat >tmp  
% eqarea.py -f tmp

NB: каб паказаць, што дадзеныя верхняй паўсферы (адмоўныя схільнасці), выкарыстоўвайце-да выключальніка.

di_eq.py

[Дадатак B.1 ]

Гэтая праграма мае скланенне, лад дадзеныя і пераўтворыць іх у X, Y пар выкарыстаннем роўныя праекцыі вобласці. Праверыць гэта з _eq _example.dat дзі. Гэта ніжняя паўсфера дадзеных class="cmti-10x-x-109">di_example.dat.

% di_eq.py -f di\_eq\_example.dat >tmp  
-0.239410246726 -0.893491204636  
0.436413031037 0.712161340983  
0.0638442163337 0.760300494699  
0.321447089986 0.68621605692  
0.322719926466 0.670562477292  
0.407412231411 0.540654291959  
0.580156197966 0.340375619675  
.  
.  
.

find_EI.py

[Кіраўнік 14 ]

Файл дадзеных быў падрыхтаваны з выкарыстаннем tk03.py для імітацыі накіравання на шыраце 42 ?. Гэта былі "плоскія", выкарыстоўваючы формулу загар Io = f tanIf to simulate inclination error and saved in a data file find_EI_ex.dat. Use the program find_EI.py знайсці уплощение фактару F, якія пры выкарыстанні на "unflatten" дае нахілу дадзеных і адноснае падаўжэнне (стаўленне вялікіх і малых уласных значэнняў матрыцы арыентацыі, см. ў Дадатку A.3.5 ) Найбольш паслядоўна з TK03.GAD paleosecular мадэлі змены.

% find_EI.py -f find_EI_ex.dat  
Bootstrapping.... be patient  
25  out of  500  
50  out of  500  
.  
.  
.  
Io Inc  I_lower, I_upper, Elon, E_lower, E_upper  
   35.8  =>     53.5 _    43.5 ^    64.3:  1.5858 _ 1.3407 ^ 1.8639  
Return to save plots - <cntl-D> to quit  
2  saved in  findEI_ei.svg  
3  saved in  findEI_cdf.svg  
1  saved in  findEI_eq.svg  
4  saved in  findEI_v2.svg  
 

, Якая вырабляе гэтыя ўчасткі:

EPSfiles/findEI.eps

У гэтым прыкладзе, арыгінальны чакаецца нахілу на paleolatitude 42 (61 ?) аднаўляецца на працягу 95% мяжы даверу.

fisher.py

[Кіраўнік 11 ]

Намалюйце набор з 10 кірункаў з размеркавання Фішэра з ? ад 30 выкарыстаннем fisher.py:

% fisher.py -k 30 -n 10  
  233.7    81.4  
  357.0    76.4  
  272.5    62.8  
  137.0    70.0  
   83.7    71.2  
...  
  

Заўважым, што кожны асобнік гэтай праграмы прыцягвае розныя размеркавання, так што твая будзе выглядаць па-рознаму ў дэталях.

fishqq.py

[Кіраўнік 11 ]

Правярае, набор з 100 кропак дадзеных створаны з fisher.py на самай справе Фішэр распаўсюджваецца з дапамогай квантилей-квантиль ўчастка:

% fisher.py -k 30 -n 100 >fishqq_example.txt  
% fishqq.py -f fishqq_example.txt

вырабляе гэтыя ўчасткі:

EPSfiles/fishqq-ex.eps

якія падтрымліваюць размеркаванне Фішэра для гэтых даных.

fishrot.py

[Кіраўнік 11 ]

Намалюйце набор з 5 напрамкаў звяртаецца з размеркаваннем Фішэра з сапраўднага сярэдняга скланенне 33, праўдзівага сярэдняга нахілу 41, і ? 50:

% fishrot.py -n 5 -D 33 -I 41 -k 50  
   35.8    32.8  
   36.2    30.2  
   37.5    41.8  
   28.6    23.9  
   26.1    32.8  
   

Прыклад 4 foldtest.py [Кіраўнік 12 ]

Выкарыстоўвайце foldtest.py для выканання foldtest на дадзеных foldtest_example.dat.

  foldtest.py -f foldtest_example.dat  
doing  1000  iterations...please be patient.....  
0  
50  
100  
...  
 

што дае участкаў:

EPSfiles/foldtest-ex.eps

Відавочна, гэтыя напрамкі былі зачыненыя ў да складваюцца таму што 95% ацэнкі ступені untilting, неабходных для максімальнага канцэнтрацыі дадзеных (максімум у прынцыпе ўласнае) ? 1 з матрыцы арыентацыі [гл. Дадатак A.3.5 ] Уключае ў сябе 100%.

gobing.py

[Кіраўнікі 12 і 14 ]

Намалюйце набор з 20 кропак дадзеных ад дыстрыбутыў TK03.GAD прадказаў для шыраты 42 ? N (гл. 14), у тым ліку развароты. Разлік параметраў з дапамогай Бінгем gobing.py

% tk03.py -n 20 -lat 42 -rev >tmp; gobing.py -f  tmp  
    2.9    62.0     8.1    15.0   -27.5     6.6   282.4    -5.0 20

якія ў адпаведнасці з даведку з gobing.py з'яўляецца:

азначае снежня, сярэдні Inc, Eta, дэталі, IETA, Зэта, ZDEC, цынк, N

gofish.py

[Кіраўнік 11 ]

Намалюйце набор з 5 напрамкаў звяртаецца з размеркаваннем Фішэра з сапраўднага сярэдняга скланенне 33, праўдзівага сярэдняга нахілу 41, і ? 50 і захаваць яго ў файл, а затым выкарыстоўваць gofish.py для разліку Фішэр параметры:

% fishrot.py -n 5 -D 33 -I 41 -k 50 > fishrot.out  
% gofish.py -f  fishrot.out  
   33.9    34.6    5     4.8723     31.3    13.9  
   

якія ў адпаведнасці з даведку з gofish.py-Н: сярэдняя снежня, сярэдні Inc, N, R, K, A95.

gokent.py

[Кіраўнік 12 ]

Намалюйце набор з 20 кропак дадзеных ад дыстрыбутыў TK03.GAD прадказаў для шыраты 42 ? N (гл. 14), без разваротаў. Разлік параметраў з дапамогай Кент gokent.py

% tk03.py -n 20 -lat 42 >tmp; gokent.py -f tmp  
  359.3    55.7     5.0   157.2    32.3     3.2   253.8    10.3 20

якія ў адпаведнасці з даведку з gobing.py з'яўляецца: сярэдняя снежня, сярэдні Inc, Eta, дэталі, IETA, Зэта, ZDEC, цынк, N

goprinc.py

[Кіраўнік 12 ]

Намалюйце набор з 20 кропак дадзеных ад дыстрыбутыў TK03.GAD прадказаў для шыраты 42 ? N (гл. 14), у тым ліку развароты. Разлічыць eigenparameters арыентацыі матрыцы (асноўныя кампаненты) з выкарыстаннем class="cmbx-10x-x-109">goprinc.py

% tk03.py -n 20 -lat -rev > tmp;  goprinc.py -f tmp  
0.93404     1.4    62.2 0.05505   185.9    27.7 0.01091    94.9     1.9 20

якія ў адпаведнасці з даведку з gobing.py is: τ1V 1D,V 1I2V 2DV 2Iτ3V 3DV 3I,N.

gaussian.py

[Кіраўнік 11 ]

Выкарыстоўвайце gaussian.py для стварэння набору з 100 нармальна размеркаваных дадзеных кропках звяртаецца ад насельніцтва, у сярэднім 10,0 і стандартным адхіленнем 30. Захавайце яго ў файл з імем gauss.dat.

% gaussian.py -s 3 -n 100 -m 10. -F gauss.out

Вы можаце праверыць выбарачнае сярэдняе і стандартнае адхіленне з stats.py ці зрабіць гістаграму дадзеных з histplot.py

histplot.py

Зрабіць гістаграму дадзеных, атрыманых з gaussian.py праграмы.

% histplot.py -f gauss.dat

што робіць сюжэт падобны на:

EPSfiles/hist.eps

гістарэзісу _magic.py

[Кіраўнікі 5 , 7 І Прыкладанні C.1 , E ]

Участак magic_measurements адфарматаваны гістарэзісу эксперыментальных дадзеных, створаных
AGM _magic.py. Выкарыстоўвайце праграму hysteresis_magic.py.

% hysteresis_magic.py -f agm_measurements.txt  
 
IS06a-1 1 out of  8  
S[a]ve plots, [s]pecimen name, [q]uit, <return> to continue  
 a  
1  saved in  IS06a-1_hyst.svg  
3  saved in  IS06a-1_DdeltaM.svg  
2  saved in  IS06a-1_deltaM.svg  
4  saved in  IS06a-1_irm.svg

which makes the plots:

EPSfiles/hysteresis-magic.eps

igrf.py

[Кіраўнік 2 ]

Выкарыстаньне праграмы igrf.py для ацэнкі поля на 1 чэрвеня 1995 года ў Амстэрдаме, Netherlands (52.5N, 5E).

% igrf.py -i  
Decimal year: <cntrl-D to quit> 1995  
Elevation [0]  
Latitude (positive north) 52.5  
Longitude (positive east) 5  
  357.9    67.4    48543  
Decimal year: <cntrl-D to quit> ^D  
Good-bye  

incfish.py

[ 10]

Выкарыстаньне праграмы incfish.py для разліку сярэдняй схільнасці да схільнасці толькі дадзеныя мадэлююцца fishrot.py для сярэдняга нахілу 60 ?. Калі вы захаваеце скланенне, лад пар, вы можаце параўнаць incfish.py answer with the Fisher mean. The datafile incfish_example_di.dat has the declination, inclination pairs and incfish_example_inc.dat has just the inclinations.

% incfish.py -f incfish_example_inc.dat  
   66.5    69.1  100     96.7    29.6     1.0  
% gofish.py -f incfish\_example_di.dat  
   1.8    71.2    100   95.1793     20.5    3.2

The output for incfish.py is: [gaussian] mean inc, Fisher inc, N,R,k,α95. You can see that the incfish.py result is much closer to the Fisherian result than the gaussian mean is.

k15_magic.py

[Chapter  13, Appendix D.1, & Appendix E]

Хтосьці ўзяў мноства узораў з дамбы маржа ў Троодос офиолитов і вымераць іх анізатрапіі магнітнай успрымальнасці ад Kappabridge KLY 2,0 інструментам у лабараторыі SIO. 15 вымярэнняў для кожнага ўзору, а таксама інфармацыю арыентацыі і ўзору імя было захавана ў файле class="cmti-10x-x-109">k15_example.dat. Convert these to the MagIC format using the program k15_magic.py:

k15_magic.py -spc 0 -f k15_example.dat -ins SIO-Bruno -loc "Troodos Ophiolite"  
 

Вы можаце ўчастка выхаду гэтага прыкладу (па змаўчанні файл rmag_anisotropy.txt) using the program aniso_magic.py.

k15_s.py

[Кіраўнік 13]

Выкарыстоўвайце K15 _s.py для разліку найлепшага тэнзар элементаў і рэшткавую хібнасць для дадзеных у файл K15 _example.dat _s. Да іх адносяцца: назва ўзору, азімут і акунуцца і простиранию і падзення, пасля чаго 15 вымярэнняў з выкарыстаннем схемы Елінэк у дадатку D.1 .. Вылічыць S дадзеныя ў узоры, геаграфічныя і куту нахілу рэгулюецца каардынаты:

% k15_s.py -f k15_example.dat  
0.331469 0.334139 0.33439023 0.000750  -0.000834  -0.000166  0.000086  
0.333359 0.333359  0.33328149 -0.001555 1 -0.001321  0.001166  0.000171  
...  
% k15_s.py -f k15_example.dat -crd g  
0.3341268  0.332827  0.333045  -0.000152  0.001248  0.001357  0.000086  
0.335563  0.331982  0.332454  0.000872  0.000241  0.000961  0.000171  
...  
% k15_s.py -f k15_example.dat -crd t  
0.334557  0.331926  0.333516  -0.000435  0.000927  0.001050  0.000086  
0.335855  0.331915  0.332229  0.000559  -0.000053  0.000647  0.000171  
...

kly-asc_magic.py

[Chapter  13 & Appendix E]

Некаторыя магнітнай успрымальнасці інструментаў захавання дадзеных у тэкставы файл, як ў:

KLY _ascii _example.txt.

Яны могуць быць імпартаваныя ў фармаце Magic з дапамогай праграмы KLY-ўзрастанні _magic.py наступным чынам:

%kly-asc_magic.py -f kly_ascii_example.txt

Гэтая каманда стварае шэраг файлаў, неабходных у базе дадзеных, магіі, гэта можа быць пабудавана пры выкарыстанні анізатропнага _magic.py або імпартаваць у Magic кансоль для далейшага рэдагавання.

kly4s_magic.py

[Chapter  13 & Appendix E]

Праграма AMSSpin даступна для загрузкі з сайта earthref.org генеруе дадзеныя для Kappabridge KLY4S спінінг магнітнай успрымальнасці інструмент, як апісана Джы і інш. (2008). Выхад файлы ў фармат файла kly4s _example.ams. Адным з варыянтаў з'яўляецца для арыентацыі інфармацыі, які будзе выводзіцца як azdip файл у фармаце (гл. azdip _magic.py.) Файлы дадзеных могуць быць імпартаваныя ў фармаце Magic з kly4s _magic.py наступным чынам:

%kly4s_magic.py -f kly4s_example.ams -fad kly4s_example.azdip

Гэтая каманда стварае шэраг файлаў, неабходных у базе дадзеных, магіі і дадзеных можа быць пабудавана пры выкарыстанні aniso_magic.py.

lnp_magic.py

[Chapter  11 & Appendix C.2.2

Гэтая праграма будзе pmag _specimen адфарматаваных файлы (напрыклад, спароджаная
zeq _magic.py) і зямельны ўчастак дадзеныя па сайце, аб'яднаўшы найбольш прыдатную ліній і найлепшае самалётаў з метадам, апісаным у Дадатку C.2.2 . Паспрабуйце гэта на дадзеных, загружаных з базы дадзеных Magic у спампаваць _magic.py. Перайдзіце ў каталог з "распакавалі" дадзеных і выканайце наступныя дзеянні:

lnp_magic.py -f  lnp_specimens.txt -crd g  
ns070  
Site lines planes  kappa   a95   dec   inc  
ns070 2  4     23.2     15.8    310.1     58.9  5.8706  
% tilt correction:  0  
s[a]ve plot, [q]uit, <return> to continue:  
 

якіх створана гэтая лічба:

EPSfiles/ЛНП-ex.eps

магнітнага _magic.py

[Дадатак E ]

Праграма магнітнага _magic.py дазваляе пераўтварыць SIO. Магнітнага магнитометр файлы фармату ў фармат магіі і ёсць варыянт для LDGO адфарматаваных дадзеных. Гэта дазваляе розным тыпам эксперыменту так што чытайце даведку. Напрыклад, даследаванне было зроблена на мноства узораў з месца "Сокорро", у тым ліку і А. Ф., цеплавых, TRM набыццё, AARM, thellier і Шоў эксперыментальных дадзеных. Яны захоўваюцца ў af_mag_example.dat, thermal_mag_example.dat, trm_mag_example.dat, aarm_mag_example.dat, thellier_mag_example.dat адпаведна. [Лабараторыі поле для thellier эксперыменту было 25 ? T і было прыкладзена ўздоўж восі Z ва ўзоры (Phi = 0, тэта = 90).]

Пераўтварэнне напрыклад файлы ў файлы фармату magic_measurement з такімі імёнамі, як af_measurements.txt, etc. Then combine them together with combine_magic.py:

 
Note:  all of these should actually be on ONE line:  
 
%mag_magic.py -f af_mag_example.dat -F af_measurements.txt  
    -LP AF -spc 1 -loc Socorro  
 
%mag_magic.py -f thellier_mag_example.dat -F  thellier_measurements.txt  
      -LP T -spc 1 -loc Socorro -dc 25 0 90  
 
%mag_magic.py -f trm_mag_example.dat  -F trmaq_measurements.txt  
      -spc 1 -loc Socorro -LP TRM  
 
%combine_magic.py -F magic_measurements.txt -f  thellier_measurements.txt  
       trmaq_measurements.txt af_measurements.txt  
 

Дадзеныя ў гэтыя файлы могуць быць нанесеныя і інтэрпрэтаваць з dmag_magic.py, zeq_magic.py or thellier_magic.py depending on the experiment.

MagIC.py

[Дадатак E ]

Гэтая праграма ўяўляе сабой графічны інтэрфейс карыстальніка (GUI) для многіх Magic функцый. Ён можа быць выкліканы з любога каталога (без прабелаў у шляху) і генеруе праграма выклікае апісаны ў гэтым дадатку. Гэта дазваляе імпартаваць многіх лабараторыі і інструмент фарматаў, пабудова розных дадзеных і рабіць асерадненні і бухгалтарскага ўліку, неабходных для стварэння файлаў гатовыя для загрузкі ў праграму Magic кансолі.

mk_redo.py

[Appendix E]

The programs zeq_magic.py and thellier_magic.py make pmag _specimen адфарматаваных файлы, якія могуць быць выкарыстаны для далейшай апрацоўкі дадзеных альбо змове або спрыяюць сайце сродкаў і г.д. Часам бывае карысна паўтарыць разлік з выкарыстаннем анізатрапіі карэкціроўкі або змены ў сістэмах каардынат і г.д. паўторна рабіць з гэтых узораў Узровень разлікаў апрацоўваецца, напрыклад zeq_magic_redo.py or thellier_magic_redo.py. These programs use magic_measurements фармат файлаў і выканання разлікаў, як гэта дыктуецца "паўтор" файл, які мае назву ўзору, ацэнкі для разліку і, у выпадку інтэрпрэтацыі размагничивания дадзеных, тыпу разліку жаданага (найлепшае ліній з магічны метад код DE-ОФЗ, найлепшае плоскасцяў з магічны метад код DE-БПС, і г.д.). Зрабіць "Вярнуць" файлы з існых pmag_specimen адфарматаваную файл дадзеныя файлы, загружаныя з сайта Magic, як у download_magic.py and examine them as follows:

% mk_redo.py  
% cat zeq_redo  
sc01h1 DE-BFP 0 0.18  
sc02a1 DE-BFL 273 813  
sc02d1 DE-BFL 273 873  
.......  
% cat thellier_redo  
sc02a1 623 858  
sc02d1 623 873  
sc02e2 623 858  
sc02h1 623 873  
.....

Адзначым, што тэмпература крокі ў градусах Кэльвіна і крокі Ф. размагничивания ў Тэсла як гэта патрабуецца ў базе дадзеных Magic.

nrm_specimens_magic.py

[Appendix E]

Пасля таго, NRM вымярэнняў, гэта часта карысна паглядзець на напрамках ў роўнай праекцыі вобласці, каб атрымаць "хуткі погляд" на вынікі, перш чым прыступіць да мудрым крокам размагничивания. Дадзеных у magic_measurements файлы, як правіла, у узоры каардынат - не геаграфічнае, так што нам патрэбны спосаб, каб павярнуць дадзеныя ў геаграфічнай і стратиграфической або каардынаты і захаваць іх у файл у фармаце pmag_specimens для пабудовы з eqarea_magic.py. The program nrm_specimens_magic.py will do this for you.

Get into the directory you made for tshe download_magic.py example. Use
nrm_specimens_magic.py to convert the NRM measurements in
magic_measurements.txt to geographic coordinates saved in a file named nrm_specimens.txt. The orientation data are in the file er_samples.txt. Then plot the specimen directions for the entire study using eqarea_magic.py:

% nrm_specimens_magic.py -crd g  
% eqarea_magic.py -f nrm_specimens.txt -crd g  
 

Першая каманда створана NRM _specimens.txt файл, а другі стварыў роўныя праекцыі вобласці напрамкаў NRM ў геаграфічных каардынатах.

Прыклад 5 арыентацыі _magic.py [Кіраўнік 9 І Дадатак E ]

Існуе дзіўнае колькасць розных спосабаў, якімі paleomagnetists дадзеных дакумента ў полі і ў лабараторыі. Базы дадзеных Magic чакае ўзору напрамках, якія будуць азімуце і акунуцца зыходных кропак стрэлка выкарыстаны для вымярэння (гл. раздзел 9 ) І арыентацыі ложка, каб быць кірунак падзення і ўніз правал. Ёсць таксама шэраг метад кодаў, якія апісваюць адбору і арыентацыі працэдуры (гл. Дадатак E для некалькіх). Для пераўтварэння з блакнота інфармацыі ў фармат Magic, мы ствараем файл з падзельнікамі табуляцыяй з загалоўка, які змяшчае "ўкладку" слова ў першым слупку і ў патрэбнае месца імя ў другой калонцы. Наступная радок мае імёны зменных (напрыклад sample_name, mag_azimuth, field_dip, дата, рашоткі, доўгі, sample_lithology, sample_type, sample_class...). Гэта зручна ў гэты момант на пастаўку литологии, тыпу i класіфікацыя інфармацыйных матэрыялаў патрабуецца Magic (гл. Дадатак E на працягу кароткага спісу і кантраляванага слоўніка спісе пад earthref.org/MAGIC вэб-сайце для атрымання поўнага спісу). Акрамя таго, можна паставіць у стратиграфических вышыні, сонца компас і GPS арыентацыі інфармацыі, сьцяг узору арыентацыі, як падазраюць, дакумент лічбавыя імёны фатаграфію поле, і хто быў звязаны з выбаркай. Ёсць варыянты для розных канвенцый арыентацыі (свердзел напрамку з прыладай Памёр арыентацыі [азімуце свідравання і Хад] па змаўчанні), розныя пагаднення аб імёнах і рашэнне аб тым, каб аўтаматычна разлічыць кошт IGRF для магнітнага скланення карэкцыі, пакажыце Ваш уласны або ігнараваць карэкцыя. Праграма генеруе er_samples.txt, er_sites.txt files. Be warned that existing files with these names will be overwritten.

Try to import the file orientation_example.txt. Гэта інфармацыйнае поле для некалькіх сайтаў. Узоры былі арыентаваны прылады Памёр арыентацыі (па змаўчанні) і пажадана, каб вылічыць магнітнае скланенне з IGRF на час адбору проб (таксама па змаўчанні). Прыклады назваў прытрымлівацца правілу, што ўзор пазначаецца літарай у канцы імя сайта (Канвенцыі № 1 - якая па змаўчанні). Так мы робім гэта наступным чынам:

%  orientation_magic.py -f orient_example.txt

pca.py

[Кіраўнік 11 ]

Гэтая праграма вылічае найбольш прыдатную лініі, плоскасці або Фішэр сярэднім праз выбраных крокаў лячэння. Фармат файла прасторы файл з падзельнікамі з узорамі імя, лячэнне крок, інтэнсіўнасць, скланення і ладу. Вылічыць найбольш прыдатную лінію праз першыя крокі лячэння дзесяць файла дадзеных zeq_example.txt:

 % pca.py -dir L 1 10 -f zeq_example.txt  
eba24a DE-BFL  
0 0.00 339.9 57.9 9.2830e-05  
1 2.50 325.7 49.1 7.5820e-05  
...  
eba24a DE-BFL 10    2.50  70.00    8.8   334.9    51.5

Па даведку, гэта: ўзор імя, тып разліку, N, прашу, канец, MAD, скланення і ладу. Разлік тыпу код метаду Магія для найлепшага ліній (гл. Дадатак E .)

plotdi_a.py

[Кіраўнік 11 ]

Месца наступныя скланення, вугал нахілу ? 95 дадзеных у файл з падзельнікамі прасторы называецца
plotdi_a_example.dat.




Dec Incα95



39.137.5 5.0
30.336.2 15
29.945.6 7
34.628.4 3



Зрабіць ўчастка з гэтых дадзеных з дапамогай plotdi_a.py:

% plotdi_a.py -f plotdi_a_example.dat  
 S[a]ve to save plot, [q]uit, Return to continue:  a  
1  saved in  eq.svg

which makes the plot:

EPSfiles/plotdi-a-example.eps

pmag_results_extract.py

[Дадатак E ]

Гэта праграма здабывае табуляцыяй тэкставы файл з pmag _results файл у фармаце. Гэта дазваляе публікаваць дадзеныя табліц, якія маюць аднолькавыя дадзеныя, каб дадзеныя былі загружаны ў базу дадзеных Magic. Паспрабуйце гэта ў каталог, створаны дляdownload_magic.py example:

% pmag_results_extract.py -f pmag_results.txt  
data saved in intensites.txt and/or directions.txt

Гэта стварае табуляцыяй файлаў, якія могуць быць уключаныя ў дакумент, напрыклад.

qqlot.py

[Дадатак В.1.5 ]

Робіць квантиль-квантиль ўчастка (гл. Дадатак В.1.5 ) З файла ўводу даных ад нармальнага размеркавання. Участак мае сярэдняе, стандартнае адхіленне і статыстыкі D, а таксама з Статыстыка D чакаць ад нармальнага размеркавання. Выкарыстоўвайце qqplot.py для праверкі дадзеных, атрыманых з gaussian.py на самай справе нармальна размеркаванай. (Гэта будзе 95% часу!).

% gaussian.py  -F  gauss.out  
% qqplot.py -f gauss.out  
mean,sigma, d, Dc  
0.02069909 0.849042146783 0.0528626896977 0.0886  
 S[a]ve to save plot, [q]uit without saving:  a  
1  saved in  qq.svg

, Які генеруе гэты ўчастак:

EPSfiles/qq.eps

revtest.py

[Кіраўнік 12 ]

Выкарыстоўвайце revtest.py для праверкі двух рэжымах ў наборы дадзеных di_example.txt are antipodal or not:

% revtest.py -f di_example.txt  
doing first mode, be patient  
doing second mode, be patient  
s[a]ve plots, [q]uit: a  
2  saved in  REV_Y.svg  
1  saved in  REV_X.svg  
3  saved in  REV_Z.svg  

, Якая вырабляе гэты ўчастак:

EPSfiles/revtest.eps

Таму што 95% даверныя мяжы для кожнага кампанента накладваюцца сябар на сябра, два напрамкі істотна не адрозніваліся.

revtest _magic.py

[Кіраўнік 12 ]

Тое ж, але для revtest.py pmag _sites Magic файлаў у фармаце. Паспрабуйце яго на файл дадзеных revtest _sites.txt. Тады паспрабуйце выкарыстоўваць налады _criteria.py змяніць або стварыць файл pmag _criteria.txt, што адпавядае вашым патрэбам і паўтор развароты выпрабаванне з выкарыстаннем толькі асобных участках.

% revtest_magic.py -f revtest_sites.txt  
% revtest_magic.py -f revtest_sites.txt -exc

s_eigs.py

[Chapter  13]

Convert the.s format data in s_eigs_example.dat to eigenvalues and eigenvectors:

% s_eigs.py -f s_eigs_example.dat  
0.33127186 239.53 44.70 0.33351338 126.62 21.47 0.33521473 19.03 37.54  
0.33177859 281.12  6.18 0.33218277 169.79 73.43 0.33603862 12.82 15.32  
...

S _geo.py

[Кіраўнік 13 ]

Друк дадзеных, якія захаваліся ў S _example.dat _geo і паварот каардынат. S дадзеныя ў геаграфічнай:

% cat s_geo_example.dat  
0.331469 0.334139 0.334390 0.00075095 -.00083439 -.00016688 80.00 -46.00  
0.333359 0.333359 0.333281 -.00155521 -.00132193 0.00116641 52.00 -23.00  
...  
% s_geo.py -f s_geo_example.dat  
0.33412680 0.33282733 0.33304587 -0.00015289 0.00124843 0.00135721  
0.33556300 0.33198264 0.33245432 0.00087259 0.00024141 0.00096166  
...

s_hext.py

[Chapter  13]

Take the output from the s_geo.py example and calculate Hext statistics:

% s_geo.py -f s_geo_exmple.dat | s_hext.py  
F =  5.79 F12 =  3.55 F23 =  3.66  
N =  8  sigma =  0.000641809950  
0.33505     5.3    14.7    25.5   124.5    61.7    13.3   268.8    23.6  
0.33334   124.5    61.7    25.1   268.8    23.6    25.5     5.3    14.7  
0.33161   268.8    23.6    13.3     5.3    14.7    25.1   124.5    61.7  
 
Note: for PC users,  try this command:  
% s_geo.py -f s_geo_exmple.dat -F tmp; s_hext.py -f tmp  

s_tilt.py

[Кіраўнік 13 ]

Паварот каардынат. S дадзеных, захаваных у _tilt S _example.dat ў стратиграфических:

s_tilt.py -f s_tilt_example.dat  
0.33455709 0.33192658 0.33351630 -0.00043562 0.00092778 0.00105006  
0.33585501 0.33191565 0.33222935 0.00055959 -0.00005316 0.00064731  
0.33586669 0.33084923 0.33328408 0.00142267 0.00013233 0.00009202  
0.33488664 0.33138493 0.33372843 -0.00056597 -0.00039086 0.00004873  
.  
.

s_magic.py

[MagIC]

Import.s format file output from the s_tilt.py example into an rmag_anisotropy formatted file. Files of the rmag_anisotropy format can be plotted with aniso_magic.py. To see how this works, use the program s_magic.py as follows:

% s_tilt.py -f s_tilt_example.dat > example.s  
% s_magic.py -f example.s

Гэта стварае выхадны файл rmag _anisotropy.txt па змаўчанні, якія могуць быць пабудаваныя з праграмай aniso_magic.py.

% cat s_tilt_example.dat  
0.33412680 0.33282742 0.33304584 -.00015292 0.00124846 0.00135721 204 25  
0.33556300 0.33198264 0.33245432 0.00087260 0.00024138 0.00096167 204 25  
0.33584911 0.33140624 0.33274472 0.00131844 0.00118815 0.00002988 204 25  
...  
% s_tilt.py -f s_tilt_example.dat  
0.33455709 0.33192658 0.33351630 -0.00043562 0.00092778 0.00105006  
0.33585501 0.33191565 0.33222935 0.00055959 -0.00005316 0.00064731  
0.33586669 0.33084923 0.33328408 0.00142267 0.00013233 0.00009202  
...

scalc.py

[Кіраўнік 14 ]

Разлічыць роскід статыстыкі S для мноства VGPs захоўваюцца ў scalc _example.txt. Паўтарыце выкарыстання абразаньня Vandamme зменнай. Тады атрымаеце загрузкі ацэнкі разліку.

% scalc.py -f scalc_example.txt  
100    20.2    180.0  
% scalc.py -f scalc_example.txt -v  
92    14.7     31.5  
% scalc.py -f scalc_example.txt -b  
100    20.2    16.3     24.4   180.0

Выкарыстанне не абразання, роскід VGP склаў 20,2 ?. Сумесна шыраты абразаньня Vandamme склала 31,5 ? які выкінуў 6 ачкоў і даў роскід 14,7 ?.

scalc _magic.py

[Кіраўнік 14 ]

Гэта ж, як scalc.py але працуе на pmag _results файлаў у фармаце. Паспрабуйце яго на pmag файл _results.txt ў каталог, створаны для загрузкі _magic.py прыклад. Выкарыстоўвайце сумесна шыраты абразаньня VGP ад 30 ?.

% scalc_magic.py -f pmag_results.txt -c 30  
8    19.0     30.0

specimens_results_magic.py

[Дадатак E ]

Пасля таго як файл pmag_specimens фармат быў створаны з выкарыстаннем, напрыклад,
thellier_magic_redo.py or zeq_magic_redo.py which take the boundary picks from the thellier_redo and zeq_redo files to calculated best-guess interpretations of the
магія дадзеных _measurements.txt, гэтыя дадзеныя павінны быць асерадніць па ўзоры і/або па сайце і пераўтворыцца ў V [] ПМ і/або VGPs і пакласці ў pmag_results фарматаванага нараўне з месца і ўзросту інфармацыі, якія даступныя. Дадзеныя павінны быць выбраны або адхілена па некаторых крытэрыях на кожным узроўні (напрыклад, MAD кірунку ўзору павінна быць менш некаторага значэння або сайта ? павінна быць больш некаторага неадназначнасьць). Гэта цяжкая задача можа быць выканана з выкарыстаннем праграмы узораў _results _magic.py. Гэтая праграма мае шмат дадатковых магчымасцяў, таму чытачу прапануецца проста паглядзець на дакументацыю і паспрабаваць яго. Ці, гэтая праграма можа быць выкліканы з MagIC.py GUI, спасылаючыся на дадзеныя скарачэння каманду / Загрузіць> Збярыце вынікі. Паспрабуйце гэта шляхам ўстаноўкі каталога, створанага ў спампаваць _magic.py напрыклад, як каталог праекта. Паглядзіце на каманды фарміруецца з выкарыстаннем розных варыянтаў.

stats.py

[Кіраўнік 11 ]

Разлічвае Гаўсавай статыстыкі для набораў дадзеных. Падлічыць сярэдняе значэнне дадзеных, атрыманых у gaussian.py example and saved in gauss.out:

% stats.py -f gauss.out  
100 0.02069909 2.069909 0.849042146783 4101.83320515

якія ў адпаведнасці з дапамогу паведамленне:

      N, mean, sum, sigma, (%) , stderr, 95% conf.  
      where sigma is the standard deviation  
      where % is sigma as percentage of the mean  
      stderr is the standard error and  
      95% conf.=  1.96*sigma/sqrt(N)

strip_magic.py

[Chapter  15 & Appendix E]

Выконвайце інструкцыі для загрузкі _magic.py але пошук Tauxe і Хартл, 1997. Спампаваць смартбуков (апошняя версія) і распакаваць яго ў новы каталог з дапамогай загрузкі _magic.py каманды і файл zmab0094214tmp02.txt ў якасці ўваходнага файла (у Tauxe і Hartl каталог). Першы запуск Газа _magic.py, каб убачыць тое, што даступна для пабудовы, то ўчастак нахілу дадзеных ад глыбіні (POS). Тады ўчастак VGP шыротах ад ўзросту:

% strip_magic.py  
available X plots:  [’age’, ’pos’]  
available Y plots:  [’dec’, ’int’, ’lat’, ’inc’, ’lon’, ’lon’]  
available method codes:  [’LP-PI-IRM’, ’LP-PI-REL’, ’LT-AF-Z’]  
% strip_magic.py -x pos -y inc  
 S[a]ve to save plot, [q]uit without saving:  a  
1  saved in  fig.svg  
 %strip_magic.py -x age -y lat  
 S[a]ve to save plot, [q]uit without saving:  a  
1  saved in  fig.svg

Апошняя каманда зрабілі гэты ўчастак:

EPSfiles/latVage.eps

sundec.py

[Кіраўнік 9 ]

Выкарыстаньне праграмы sundec.py для разліку азімуце напрамкі свердзела. Вы знаходзіцеся ў 35 ? N і 33 ? Е. Па мясцовым часе ў тры гадзіны апярэджвае сусветнае час. Цень кут напрамкі свідравання склаў 68 ? вымераная ў 16:09 23 мая 1994 года.

Праграма sundec.py працуе альбо на інтэрактыўнай ўводу даных або шляхам чытання з файла.

Захавайце наступны ў файл з імем sundec _example.dat:

3 35 33 1994 5 23 16 верасні 1968

, Якія:

? GMT Lat Lon года пн дзень гг мм shadow_angle

Мы можам прааналізаваць гэты файл альбо:

% sundec.py -f sundec_example.dat  
  154.2

or

% sundec.py < sundec_example.dat  
  154.2

or by manual input:

% sundec.py -i  
Time difference between Greenwich Mean Time (hrs to ADD to  
      GMT for local time):  
<cntl-D> to quit 3  
Year:  <cntl-D to quit> 1994  
Month:  5  
Day:  23  
hour:  16  
minute:  9  
Latitude of sampling site (negative in southern hemisphere): 35  
Longitude of sampling site (negative for western hemisphere): 33  
Shadow angle: 68  
  154.2  
Time difference between Greenwich Mean Time (hrs to ADD to  
      GMT for local time):  
<cntl-D> to quit ^D  
 Good-bye  

У любым выпадку, скланенне 154,2 ?.

thellier _magic.py

[Кіраўнік 10 І Дадатак E ]

Выкарыстоўвайце thellier _magic.py праграма для пабудовы thellier эксперыментальных дадзеных, загружаных у спампаваць _magic.py. У каталогу з загружаных даных, вы ўбачыце тлумачэння арыгінальных аўтараў, калі абраць pmag _specimens.txt як ваш ўзор ўваходнага файла. Перайсці да sc12b2 ўзору, вылічыць лепшых схілаў падыходзяць і палеонапряженности выкарыстаннем інтэрвале паміж 475 і 590 ? C (крокі 8 да 21), захавання і інтэрпрэтацыі фігуры і выхаду.

% thellier_magic.py -fsp pmag_specimens.txt  
90-13a1 1 of  158  
0  0   335.4   -45.1 2.340e-07  
1  100   330.8   -37.9 1.990e-07  
....  
Looking up saved interpretation....  
Saved interpretation:  
 
killed by:  
specimen_b_beta  
 
 
specimen Tmin  Tmax  N  lab_field  B_anc  b  q  f(coe)  Fvds  beta  MAD  
         Dang  Drats     Nptrm  Grade  R  MD%  sigma  Z Gmax  
 
90-13a1  250  550 7 25.0  4.6 -0.185   4.6 0.764 0.480 0.125      5.1  
     1.5     4.4 8     B  0.961 8 0.023    -1.0     2.5  
 
  s[a]ve plot, set [b]ounds for calculation, [d]elete current  
         interpretation,  [p]revious,      [s]ample, [q]uit:  
 
Return for next specimen  
s  
Enter desired specimen name (or first part there of): sc12b2  
 
 sc12b2 107 of  158  
0  0   185.0    10.3 3.000e-05  
1  100   187.5    10.5 2.990e-05  
2  200   187.4     9.4 2.930e-05  
....  
Step Temperature  Gamma  
0 100    77.5  
1 200    89.1  
....  
18 580     4.2  
19 585     2.8  
20 590     0.9  
Looking up saved interpretation....  
Saved interpretation:  
 
killed by:  
specimen_drats  
 
 
specimen Tmin  Tmax  N  lab_field  B_anc  b  q  f(coe)  Fvds  beta  MAD  
                Dang  Drats           Nptrm  Grade  R  MD%  sigma  Z Gmax  
 
sc12b2  475  590 14 25.0 42.3 -1.693  23.5 0.733 0.646 0.028      3.9  
              1.3    13.5     9 B  0.995 1 0.048     0.0     0.9  
 
Optimization terminated successfully.  
         Current function value: 0.000000  
         Iterations: 85  
         Function evaluations: 161  
Optimization terminated successfully.  
         Current function value: 0.000000  
         Iterations: 47  
         Function evaluations: 93  
3.38338116815e-05  
Banc=  44.8088483904  
Banc=  44.8088483904  
 
   s[a]ve plot, set [b]ounds for calculation, [d]elete current  
    interpretation,  [p]revious,      [s]ample, [q]uit:  
Return for next specimen  
b  
Enter index of first point for calculation:  [ 8 ]  
return to keep default  
Enter index  of last point for calculation:  [ 21 ]  
return to keep default  
 
killed by:  
specimen_drats  
 
 
specimen Tmin  Tmax  N  lab_field  B_anc  b  q  f(coe)  Fvds  beta  MAD  
               Dang  Drats          Nptrm  Grade  R  MD%  sigma  Z Gmax  
 
sc12b2  475  590 14 25.0 42.3 -1.693  23.5 0.733 0.646 0.028      3.9  
           1.3    13.5        9 B  0.995 1 0.048     0.0     0.9  
 
Optimization terminated successfully.  
         Current function value: 0.000000  
         Iterations: 85  
         Function evaluations: 161  
Optimization terminated successfully.  
         Current function value: 0.000000  
         Iterations: 47  
         Function evaluations: 93  
Banc=  44.8088483904  
Save this interpretation? [y]/n  
 
y  
sc12c2 108 of  158  
0  0   215.4     2.4 1.310e-05  
1  100   213.1     2.2 1.280e-05  
2  200   212.8     2.2 1.180e-05  
...  
Saved interpretation:  
specimen Tmin  Tmax  N  lab_field  B_anc  b  q  f(coe)  Fvds  beta  MAD  
               Dang  Drats         Nptrm  Grade  R  MD%  sigma  Z Gmax  
 
sc12c2  100  525 8 25.0 35.3 -1.414   7.9 0.583 0.534 0.062      4.1  
          3.9     3.3     5 A  0.989 -1 0.087     2.0     3.1  
 
 
    s[a]ve plot, set [b]ounds for calculation, [d]elete current  
     interpretation,  [p]revious,      [s]ample, [q]uit:  
Return for next specimen  
a  
4  saved in  sc12c2_TRM.svg  
3  saved in  sc12c2_arai.svg  
1  saved in  sc12c2_deremag.svg  
2  saved in  sc12c2_zijd.svg  
sc12c2 108 of  158  
0  0   215.4     2.4 1.310e-05  
1  100   213.1     2.2 1.280e-05  
....  
w  
    s[a]ve plot, set [b]ounds for calculation, [d]elete current  
     interpretation,  [p]revious,      [s]ample, [q]uit:  
Return for next specimen  
q  
Good bye  
 

Гэты пратакол створаны файлthellier_specimens.txt and the plots shown in Figure F.4.


EPSfiles/thellier-magic-example.eps

Figure F.4: Output for Example 77.


thellier_magic_redo.py

[Кіраўнік 10 , 13 І Дадатак E ]

Гэтая праграма дазваляе пералік thellier даных з выкарыстаннем анізатрапіі папраўкі ад AARM або ATRM эліпсоіду (пачата ў rmag _anisotropy файлаў фармату (гл. раздзел 13 і прыкладам для aarm _magic.py) і нелінейных TRM карэкціроўкі, калі TRM збору дадзеных даступныя. Каб скарыстацца гэтым, па-першае прытрымлівайцеся інструкцыям у aarm _magic.py напрыклад для стварэння aarm _measurements.txt і rmag _anisotropy.txt файлы, якія вы будзеце мець патрэбу.

Пераўтварэнне двух іншых вымярэнняў файлы ў фармаце Magic. Файл bg.mag мае палеонапряженности ў ім дадзеныя. Лабараторыя поле прыкладваліся ўздоўж восі Z і 20 ? T і лабараторыі пратакол цеплавой эксперыменту. Стварыце выхадны файл з імем thellier _measurements.txt. Файл bg.trm мае TRM збору дадзеных у ёй. Лабараторыя пратакол быў таксама цеплавыя эксперымент, але вы таксама павінны ўсталяваць сцяг-TRM для папярэджання праграмы, што гэты файл мае TRM збору дадзеных. Стварыце выхадны файл з імем trmaq _measurements.txt. І нарэшце, аб'яднаць усе гэтыя розныя файлы дадзеных вымярэнняў у адным файле _measurements.txt магіі з дапамогай каманды аб'яднаць _magic.py.

Тым не менш у thellier каталог _example _redo, стварыце файл thellier _redo выкарыстаннем
МК _redo.py ад pmag _specimens адфарматаваны thellier _specimens.txt файл.

Зараз з дапамогай праграмы thellier _magic _redo.py для стварэння анізатрапіі выпраўленая
pmag _specimen адфарматаваную выхадны файл называецца AC _specimens.txt, нелінейных TRM выпраўлены файл называецца
NLT _specimens.txt і файл з абодвух карэкціроўкі называецца NLT _AC _specimens.txt. Змяшайце ўсе pmag _specimen-файлаў у адзін файл _specimen.txt pmag.

Вось стэнаграма сесіі, што робіць усё гэта. Звярніце ўвагу, што водступ лініі належаць на папярэднюю радок і проста аддзяліць для фарматавання мэтаў:

% mag_magic.py -f bg.trm -loc Bushveld -F trmaq_measurements.txt  
       -LP T -trm -ncn 3  
results put in  trmaq_measurements.txt  
% mag_magic.py -f bg.mag -loc Bushveld -dc 20 0 90 -LP T -F  
        thellier_measurements.txt -ncn 3  
results put in  thellier_measurements.txt  
% mag_magic.py -f bg.arm -loc Bushveld -LP AF:ANI -F aarm_measurements.txt  
        -ncn 3 -ac 180 -dc 50 0 90  
Warning - inconsistency in mag file with lab field - overriding file with 0  
Warning - inconsistency in mag file with lab field - overriding file with 0  
Warning - inconsistency in mag file with lab field - overriding file with 0  
.......  
% combine_magic.py -F magic_measurements.txt -f trmaq_measurements.txt  
          thellier_measurements.txt  aarm_measurements.txt  
File  ./trmaq_measurements.txt  read in with  70  records  
File  ./thellier_measurements.txt  read in with  476  records  
File  ./aarm_measurements.txt  read in with  126  records  
All records stored in  ./magic_measurements.txt  
% mk_redo.py -f thellier_specimens.txt  
%  thellier_magic_redo.py -ANI  
Processing  8  specimens - please wait  
9  
%  thellier_magic_redo.py -NLT  
Processing  8  specimens - please wait  
Optimization terminated successfully.  
         Current function value: 0.000001  
         Iterations: 67  
         Function evaluations: 130  
Optimization terminated successfully.  
         Current function value: 0.000001  
         Iterations: 45  
         Function evaluations: 89  
         ........  
 % thellier_magic_redo.py -usr "Lisa Tauxe" -NLT -ANI  
            -Fac NLT_AC_specimens.txt  
Processing  8  specimens - please wait  
Optimization terminated successfully.  
         Current function value: 0.000001  
         Iterations: 67  
         Function evaluations: 130  
Optimization terminated successfully.  
 
%combine_magic.py -F pmag_specimens.txt -f thellier_specimens.txt  
         NLT_specimens.txt AC_specimens.txt NLT_AC_specimens.txt

tk03.py

[Кіраўнік 14 ]

Намалюйце 200 кірункаў з paleosecular мадэлі змены TK03.GAD чакаецца шыраце 42 ?, уключаючы развароты і захаваць іх у файл з імем tk03.out :

% tk03.py -n 200 -lat 42 -rev > tk03.out

Вы можаце зрабіць сюжэт з гэтых дадзеных з выкарыстаннем eqarea.py ці многіх іншых праграм, якія прымаюць вектарныя дадзеныя (goprinc.py, gobing.py, gofish.py).

uniform.py

Draw 10 кірунках з раўнамернага размеркавання і захавання ў файл з імем unf.out:

% uniform.py -i  
Desired number of uniform directions 10  
Output file for saving? unf.out

upload_magic.py

[Дадатак E ]

Гэтая праграма мае ўсе Magic фармат файлаў і змяшчае іх у файл, які можа быць імпартаваны ў Magic кансолі праграмнае забеспячэнне для загрузкі ў базу дадзеных Magic. У якасці прыкладу, мы можам "ўпакаваць" файл, які выкарыстоўваецца файл загружаецца і апрацоўваецца ў thellier _magic.py прыклад. Для гэтага перайдзіце ў каталог, створаны для гэтага праекта і тып загрузкі _magic.py. Гэта стварае файл загрузкі _dos.txt, якія можна чытаць у чароўны кансоль для канчатковай падрыхтоўкі для загрузкі ў базу дадзеных Magic. Вядома, вы не павінны паўторна загрузіць гэты файл, толькі дадзеныя з апублікаваных прац або тэзісаў належаць ў базе дадзеных і гэты дакумент быў ужо апублікаваны і дадзеных, загружаных аўтараў.

% upload_magic.py  
Removing:  [’citation_label’, ’compilation’,..., ’average_n_planes’]  
./er_expeditions.txt is bad or non-existent - skipping  
file  ./er_locations.txt  successfully read in  
er_locations written to  upload.txt  
file  ./er_samples.txt  successfully read in  
only first orientation record from er_samples.txt read in  
 
er_samples written to  upload.txt  
file  ./er_specimens.txt  successfully read in  
only measurements that are used for interpretations  
er_specimens written to  upload.txt  
.....  
now converting to dos file ’upload_dos.txt’  
Finished preparing upload file

vgp_di.py

[Кіраўнік 2 ]

Выкарыстаньне праграмы VGP _di.py канвертаваць наступныя:





λp ϕpλs ϕs




6819133243




Пакладзеце дадзеныя ў файл VGP _di _example.dat напрыклад, з дапамогай кот на аперацыйнай сістэме * Nix. Вось стэнаграма адзін спосаб выкарыстання праграмы, якая выплёўвае скланенне, лад:

% vgp_di.py -f vgp_di_example.dat  
  335.6    62.9

watsonsF.py

[Кіраўнік 11 ]

Першы стварэння двух файлаў дадзеных з fishrot.py з ? = 15, N = 10 і I = 42, з D = 10. для першага і 20 для другога:

% fishrot.py -k 15 -n 10 -I 42 -D 10 > file1  
% fishrot.py -k 15 -n 10 -I 42 -D 20 > file2

Каб параўнаць гэтыя два файла, watsonsF.py:

% WatsonsF.py-F file1 file2-F2
0.293622135833 3.32

Першае лік F Статыстыка Ўотсана для гэтых двух файлаў (гл. раздзел 11 ) І другі нумар біць для двух файлаў, якія будуць вынятыя з той жа размеркаванне Фішэра (доля агульнага сярэдняга). У дадзеным выпадку гэта не дзіўна, што дадзеныя прайсці гэты тэст (F менш неабходнай колькасці).

watsonsV.py

[Кіраўнік 11 ]

Выкарыстанне двух дадзеных файлаў, створаных у прыклад для watsonsF.py і паўтарыць выпрабаванне з выкарыстаннем V W Ўотсана статыстыкі.

% watsonsV.py -f file1 -f2 file2  
Doing  500  simulations  
50  
100  
150  
...  
Watson’s V,  Vcrit:  
       0.6        6.8  
 S[a]ve to save plot, [q]uit without saving:  a  
1  saved in  cdf.svg  
       

, Які генеруе ўчастку:

EPSfiles/watsonsV-example.eps

Два файла не значна адрозніваецца тым, што ў V Уотсан (0,6 ў дадзеным прыкладзе) менш, чым значэнне крытычнага V ацэньваецца з выкарыстаннем метаду Монте-Карла (6,8).

Звярніце ўвагу, што вашы вынікі могуць адрознівацца ў падрабязнасці, таму што кожны асобнік fishrot.py спараджае розныя выпадкова звяртаецца набору дадзеных.

zeq.py

[Кіраўнік 9 ]

а) Выкарыстаньне праграмы zeq.py да 1) ўчастак дыяграмы Zijderveld дадзеных у
zeq _example.txt. б) разлік найлепшага лініі ад 15 да 90 МТЛ. в) Паварот дадзеныя, такія, што найлепшае лініі праецыюецца на гарызантальнай восі (а не па змаўчанні, поўнач). г) Calcualte найлепшага плоскасці ад 5 да 80 МТЛ. Захаваць гэтыя ўчасткі.

% zeq.py -f zeq_example.dat -u mT

Пры выбары "B", вы можаце выбраць мяжы і выберыце 'L' для найлепшага лініі або "Р" для найлепшага плоскасці. Вы можаце круціць XY восяў, выбраўшы 'H' і настройка X восі да 312. Нарэшце, вы можаце захоўваць свае ўчасткі з магчымасцю ". Вы павінны захаваць нешта накшталт гэтых участкаў:

EPSfiles/zeq.eps

zeq_magic.py

[Кіраўнік 9 І Дадатак E ]

Участак размагничивания Ф. дадзеныя ў файл, які вы атрымалі ў download_magic.py example using zeq_magic.py. Выкарыстаньне геаграфічных каардынат, дзе арыентацыі даступныя.

%zeq_magic.py  -fsp pmag_specimens.txt -crd g

zeq_magic_redo.py

[Кіраўнік 9 І Дадатак E ]

У той жа каталог, створаны для download_magic.py and used in the
zeq_magic.py Напрыклад, можна выкарыстоўваць праграмуzeq_magic_redo.py to create a pmag_specimens адфарматаваную файл з дадзенымі ў геаграфічных каардынатах. Мяркуючы, што ўзор арыентацыі ў файл з імем er_samples.txt, выкарыстоўвайце МК _redo спачатку стварыць файл з імем zeq _redo. Затым з дапамогай zeq _magic _redo.py стварыць два _specimen pmag файлаў у фармаце: адзін у узоры каардынаты zeq _s.txt _specimens і адзін у геаграфічных каардынатах zeq _specimens _g.txt. Аб'яднайце гэтыя ў адзін файл з імем pmag _specimens.txt.

Звярніце ўвагу, што водступ лініі належаць з лініі вышэй у выглядзе адной радкі.

%mk_redo.py -f pmag_specimens.txt  
% zeq_magic_redo.py -f magic_measurements.txt  -F zeq_specimens_s.txt  -crd s  
% zeq_magic_redo.py -f magic_measurements.txt  -F zeq_specimens_g.txt  -crd g  
% combine_magic.py -F pmag_specimens.txt  
     -f zeq_specimens_s.txt zeq_specimens_g.txt
Popular Links
Published (Last edited): Apr 5 , source: http://magician.ucsd.edu/Essentials/WebBookse128.html