grdflexuregmt - Online în cloud

PROGRAM:

NUME

grdflexure - Calculează deformarea la încovoiere a suprafețelor 3-D pentru diverse reologii

REZUMAT

grdflexure topogrd rm/rl[/ri]/rw Te[u] outgrid [ Nx/Ny/Nxy ] [ ppește ] [ ytânăr ] [
nu_a[/Ha/nu_m] ] [ listă ] [ [f|q|s|nx/ny][+a|d|h|l][+e|n|m][+tlățime][+w[sufix]][+z[p,
beta ] [ -Tt0[u][/t1[u]/dt[u]|n][+l] ] [ [nivel] ] [ wd] [ zm] [ -fg ]

Notă: Nu este permis niciun spațiu între indicatorul opțiunii și argumentele asociate.

DESCRIERE

grdflexure calculează răspunsul la încovoiere la sarcini folosind o gamă de valori selectabile de utilizator
reologii. Utilizatorul poate alege dintre elastic, vâscoelastic sau firmoviscos (cu unul sau două
straturi vâscoase). Evoluția temporală poate fi modelată și prin furnizarea de încărcare incrementală
grile și specificarea unui interval de timpi de ieșire a modelului.

NECESARE ARGUMENTE

topogrd
Fișier grilă binar 2-D cu topografia încărcăturii (în metri); Vezi FIȘIER GRID
FORMATELE de mai jos. Dacă -T este folosit, topogrd poate fi un șablon de nume de fișier cu un flotant
va fi setat și încărcat pentru un format punct (sintaxă C) și un alt nume de fișier de încărcare
fiecare pas de timp. Timpii de încărcare coincid astfel cu timpii dați prin -T (dar nu
toate timpurile trebuie să aibă un fișier corespunzător). Alternativ, dă topogrd as
=*flist*, unde flist este un tabel ASCII cu unul topogrd numele fișierului și timpul de încărcare per
record. Acești timpi de încărcare pot fi diferiți de timpii de evaluare furnizați prin -T.
Pentru formatul timpului de încărcare, vezi -T.

-Drm/rl[/ri]/rw
Setează densitatea pentru manta, încărcare, umplutură (opțional, altfel se presupune că este egală cu
densitatea sarcinii) și apă sau aer. Dacă ri difera de rl apoi o aproximativă
se va gasi solutie. Dacă ri nu este dat, atunci este implicit rl.

-ETe Setează grosimea plăcii elastice (în metri); adăuga k pentru km. Dacă elasticul
grosimea depășește 1e10 va fi interpretată ca o rigiditate la încovoiere D (în mod implicit
D este calculat din Te, modulul lui Young și raportul lui Poisson; vedea -C pentru a schimba acestea
valori).

-Goutfile
If -T este setat atunci grdfile trebuie să fie un șablon de nume de fișier care conține un flotant
format punct (sintaxa C). Dacă șablonul nume de fișier conține și %s (pentru
numele unității) sau %c (pentru litera unității), apoi folosim timpul corespunzător (în unități
specificat în -T) pentru a genera numele de fișiere individuale, altfel folosim timpul în
ani fără unitate.

OPTIONAL ARGUMENTE

-ANx/Ny/Nxy
Specificați forțele de compresie sau de extensie în plan în direcțiile x și y, ca
precum și orice forță de forfecare [fără forțe în plan]. Compresia este indicată prin negativ
valori, în timp ce forțele de extensie sunt specificate folosind valori pozitive.

-Cppește
Modificați valoarea curentă a raportului lui Poisson [0.25].

-Cytânăr
Modificați valoarea curentă a modulului lui Young [7.0e10 N/m^2].

-Fnu_a[/Ha/nu_m]
Specificați un model firmoviscous împreună cu o grosime a plăcii elastice
specificat prin -E. Dați doar o vâscozitate (nu_a) pentru o placă elastică peste a
semi-spațiu vâscos, sau adăugați, de asemenea, grosimea astenosferei (Ha) si
vâscozitate mai mică a mantalei (nu_m), prima vâscozitate fiind acum cea a lui
astenosferă. Dați vâscozități în Pa*s. Dacă este folosit, dați grosimea
astenosferă în metru; adăuga k pentru km.

-N[f|q|s|nx/ny][+a|[+d|h|l][+e|n|m][+tlățime][+w[sufix]][+z[p]]
Alegeți sau întrebați despre dimensiunile grilei potrivite pentru FFT și setați opțional
parametrii. Controlați dimensiunea FFT:
-Nf va forța FFT să utilizeze dimensiunile reale ale datelor.

-Nq va întreba despre dimensiuni mai potrivite, le va raporta, apoi va continua.

-Ns va prezenta o listă de dimensiuni opționale, apoi va ieși.

-Nnx/ny va face FFT pe dimensiunea matricei nx/ny (trebuie să fie >= dimensiunea fișierului grilă). Mod implicit
alege dimensiuni >= date care optimizează viteza și acuratețea FFT. Dacă FFT
dimensiuni > dimensiuni fișier grid, datele sunt extinse și reduse la zero.

Controlați reducerea tendinței datelor: adăugați modificatori pentru eliminarea unei tendințe liniare:
+d: Date de reducere a tendinței, adică eliminați tendința liniară cea mai potrivită [Implicit].

+a: eliminați numai valoarea medie.

+h: Eliminați doar valoarea medie, adică 0.5 * (max + min).

+l: Lasă datele în pace.

Controlați extinderea și reducerea treptată a datelor: utilizați modificatori pentru a controla modul în care extensia
și înclinarea trebuie efectuată:
+e extinde grila prin impunerea simetriei margine-punct [Implicit],

+m extinde grila prin impunerea simetriei oglinzii marginilor

+n dezactivează extensia de date.

Înclinarea se realizează de la marginea datelor la marginea grilei FFT [100 %]. Schimbare
acest procent prin +tlățime. Când +n este în vigoare, se aplică înclinarea
în schimb la marjele de date, deoarece nu este disponibilă nicio extensie [0%].

Controlați scrierea rezultatelor temporare: pentru investigații detaliate puteți scrie
grila intermediară fiind transmisă către FFT înainte; probabil că asta a fost
descurcat, extins prin simetrie punctuală de-a lungul tuturor marginilor și conic. Adăuga
+w[sufix] din care vor fi create numele fișierelor de ieșire (de exemplu, ingrid_prefix.ext)
[conic], unde ext este extensia ta de fișier. În cele din urmă, puteți salva grila complexă
produs de către FFT forward prin anexare +z. În mod implicit scriem real și
componente imaginare la ingrid_real.ext si ingrid_imag.ext. Adăuga p pentru a salva
în schimb forma polară a mărimii și fazei la fișiere ingrid_mag.ext si
ingrid_fază.ext.

-Llistă Scrieți numele și timpii de evaluare a tuturor grilelor care au fost create în text
fişier listă. Necesită -T.

-Mtm Specificați un model vâscoelastic împreună cu o grosime a plăcii elastice
specificat prin -E. Adăugați ora Maxwell tm pentru modelul vâscoelastic (în ).

-Sbeta Specificați o fracție de șanț înfometată în intervalul 0-1, unde 1 înseamnă că șanțul de șanț este complet
umplut cu material de densitate ri în timp ce 0 înseamnă că este umplut doar cu material de
densitate rw (adică doar apă) [1].

-Tt0[u][/t1[u]/dt[u]|n][+l]
Specifica t0, t1, și increment de timp (dt) pentru succesiunea calculelor [Implicit este
un pas, fără dependență de timp]. Pentru o singură oră specifică, începeți
timp t0. Unitatea este de ani; adăuga k pentru kyr și M pentru Myr. Pentru un timp logaritmic
scară, anexă +l și specificați n pași în loc de dt. Alternativ, dați un fișier
cu timpii doriti în prima coloană (acești timpi pot avea unități individuale
anexat, altfel presupunem anul). Apoi scriem un fișier de grilă model separat pentru
fiecare pas de timp dat.

-Wwd Setați adâncimea de referință la suprafața îndoită neformată în m [0]. Adăuga k la
indica km.

-Zzm Specificați adâncimea de referință la suprafața îndoită (de exemplu, Moho) în m; adăuga k pentru km. Trebuie sa
fii pozitiv. [0].

-V[nivel] (Mai mult ...)
Selectați nivelul de verbozitate [c].

-fg Grilele geografice (dimensiunile de longitudine, latitudine) vor fi convertite în metri
printr-o aproximare „Pământ plat” folosind parametrii elipsoizi actuali.

-^ or doar -
Imprimă un mesaj scurt despre sintaxa comenzii, apoi iese (NOTĂ: pe Windows
foloseste doar -).

-+ or doar +
Imprimați un mesaj extins de utilizare (ajutor), inclusiv explicația oricăruia
opțiunea specifică modulului (dar nu opțiunile comune GMT), apoi iese.

-? or Nu. argumente
Apoi imprimați un mesaj complet de utilizare (ajutor), inclusiv explicația opțiunilor
iesirile.

--versiune
Tipăriți versiunea GMT și ieșiți.

--show-datadir
Imprimați calea completă către directorul de partajare GMT și ieșiți.

GRID FILE FORMATE

În mod implicit, GMT scrie grila ca flotoare de precizie unică într-un netCDF de reclamație COARDS
tipul fisierului. Cu toate acestea, GMT este capabil să producă fișiere grilă în multe alte grile utilizate în mod obișnuit
formate de fișiere și facilitează, de asemenea, așa-numita „împachetare” a grilelor, scrierea în virgulă mobilă
date ca numere întregi de 1 sau 2 octeți. Pentru a specifica precizia, scara și decalajul, utilizatorul ar trebui
adăugați sufixul =id[/scară/compensa[/nan]], Unde id este un identificator de două litere al grilei
tipul și precizia și scară si compensa sunt factor de scară opțional și offset să fie
aplicat tuturor valorilor grilei și nan este valoarea folosită pentru a indica datele lipsă. In caz
cele două personaje id nu este prevăzut, ca în =/scară Decât a id=nf este asumat. Cand
citind grile, formatul este, în general, recunoscut automat. Dacă nu, același sufix
poate fi adăugat la numele fișierelor grilei de intrare. Vedea grdconvert și Secțiune grid-file-format al
Referințe tehnice GMT și carte de bucate pentru mai multe informații.

Când citește un fișier netCDF care conține mai multe grile, GMT va citi, în mod implicit, fișierul
prima grilă bidimensională care poate fi găsită în acel fișier. Pentru a convinge GMT să citească altul
variabilă multidimensională în fișierul grilă, anexați ?varname la numele fișierului, unde
varname este numele variabilei. Rețineți că poate fi necesar să scăpați de sensul special
of ? în programul shell, punând o bară oblică inversă în fața acestuia sau plasând
nume de fișier și sufix între ghilimele sau ghilimele duble. The ?varname se poate folosi și sufixul
pentru grilele de ieșire să specifice un nume de variabilă diferit de cel implicit: „z”. Vedea
grdconvert și Modificatori de secțiuni pentru CF și format de fișier grilă din Tehnica GMT
Referințe și Carte de bucate pentru mai multe informații, în special despre cum să citiți îmbinările de 3-,
Grile 4 sau 5-dimensionale.

GRID DISTANTA UNITĂȚI

Dacă grila nu are un contor ca unitate orizontală, adăugați +uunitate la fișierul de intrare
numele de convertit din unitatea specificată în metru. Dacă grila dvs. este geografică, convertiți
distante la metri prin alimentare -fg in schimb.

CONSIDERAȚII

Grilele netCDF COARDS vor fi recunoscute automat ca fiind geografice. Pentru alte grile
grilele geografice în care doriți să convertiți grade în metri, selectați -fg. Dacă datele
sunt aproape de oricare dintre poli, ar trebui să luați în considerare proiectarea fișierului grilă pe un dreptunghiular
folosind sistemul de coordonate grdproject.

FARFURIE FLEXURA NOTE

Soluția FFT la îndoirea plăcii necesită ca densitatea de umplere să fie egală cu densitatea sarcinii.
De obicei, acest lucru este valabil numai direct sub sarcină; dincolo de sarcina la care tinde umplutura
fie sedimente cu densitate mai mică sau chiar apă (sau aer). Wessel [2001] a propus un
aproximare care permite specificarea unei densități de umplere diferite de
densitatea sarcinii, permițând totuși o soluție FFT. Practic, flexiunea plăcii este
rezolvat pentru utilizarea densității de umplere ca densitate efectivă a sarcinii, dar amplitudinile sunt
ajustat de factor A = sqrt ((rm - ri)/(rm - rl)), care este diferența teoretică
în amplitudine datorită unei sarcini punctuale folosind cele două densități diferite de sarcină. The
aproximarea este foarte bună, dar se defectează pentru încărcături mari pe plăci slabe, o zână
situație neobișnuită.

EXEMPLE

Pentru a calcula îndoirea plăcii elastice din sarcină topo.nc, pentru o placă groasă de 10 km cu
densități tipice, încercați

gmt grdflexure topo.nc -Gflex.nc -E10k -D2700/3300/1035

Pentru a calcula răspunsul firmoviscous la o serie de încărcări incrementale date de numele fișierului
iar timpul de încărcare în tabelul l.lis la un singur timp 1 Ma folosind reologicul specificat
valori, încercați

gmt grdflexure -T1M =l.lis -D3300/2800/2800/1000 -E5k -Gflx/smt_fv_%03.1f_%s.nc -F2e20 -Nf+a

REFERINȚE

Cathles, LM, 1975, viscozitate of il ale pământului manta, Princeton University Press.

Wessel. P., 2001, Distribuția globală a munților submarin deduse din Geosat/ERS-1 grid
altimetrie, J. Geophys. Res., 106(B9), 19,431-19,441,
http://dx.doi.org/10.1029/2000JB000083

Utilizați grdflexuregmt online folosind serviciile onworks.net