grdflexuregmt — online w chmurze

PROGRAM:

IMIĘ

grdflexure — Obliczanie zginania powierzchni trójwymiarowych dla różnych reologii

STRESZCZENIE

gdflexure topogród rm/rl[/ri]/rw Te[u] przesłona [ Nx/Ny/Nxy ] [ pryba ] [ ymłody ] [
nu_a[/h_a/liczba_m] ] [ podstęp ] [ [f|q|s|nx/ny][+a|d|h|l][+e|n|m][+tszerokość][+w[przyrostek]][+z[p]] [
beta ] [ -Tt0[u][/t1[u]/dt[u]|n][+l] ] [ [poziom] ] [ wd] [ zm] [ -fg ]

Uwaga: Między flagą opcji a powiązanymi argumentami nie może być spacji.

OPIS

gdflexure oblicza reakcję zginania na obciążenia przy użyciu zakresu wybieranych przez użytkownika
reologie. Użytkownik może wybrać elastyczną, wiskoelastyczną lub jędrną lepkość (z jednym lub dwoma
lepkie warstwy). Ewolucję czasową można również modelować, dostarczając obciążenie przyrostowe
siatki i określenie zakresu czasów wyjścia modelu.

WYMAGANE ARGUMENTY

topogród
Plik siatki binarnej 2D z topografią obciążenia (w metrach); Zobacz PLIK SIATKI
FORMATY poniżej. Gdyby -T jest używany, topogród może być szablonem nazwy pliku z pływającym
format punktu (składnia C) i inna nazwa pliku wczytywania zostanie ustawiona i wczytana dla
każdy krok czasu. Czasy ładowania pokrywają się zatem z czasami podanymi przez -T (ale nie
zawsze trzeba mieć odpowiedni plik). Ewentualnie daj topogród as
=*lista*, gdzie flet jest tabelą ASCII z jednym topogród nazwa pliku i czas ładowania per
nagrywać. Te czasy ładowania mogą różnić się od czasów oceny podanych przez -T.
Aby zapoznać się z formatem czasu ładowania, zobacz -T.

-Drm/rl[/ri]/rw
Ustawia gęstość dla płaszcza, obciążenia, wypełnienia (opcjonalnie, w przeciwnym razie zakłada się, że jest równe
gęstość ładunku) oraz wodę lub powietrze. Gdyby ri różni się od rl potem przybliżony
rozwiązanie zostanie znalezione. Gdyby ri nie jest podany, to domyślnie rl.

-ETe Ustawia grubość płyty elastycznej (w metrach); dodać k na km. Jeśli elastyczna
grubość przekracza 1e10 zostanie zinterpretowana jako sztywność zginania D (domyślnie
D jest obliczane z Te, moduł Younga i współczynnik Poissona; zobaczyć -C zmienić to
wartości).

-Gplik wyjściowy
If -T jest wtedy ustawiony plik grd musi być szablonem nazwy pliku, który zawiera pływający
format punktowy (składnia C). Jeśli szablon nazwy pliku zawiera także %s (dla
nazwa jednostki) lub %c (dla litery jednostki), wtedy używamy odpowiedniego czasu (w jednostkach
określone w -T) do generowania poszczególnych nazw plików, w przeciwnym razie używamy czasu w
lata bez jednostki.

OPCJA ARGUMENTY

-ANx/Ny/Nxy
Określ siły ściskające lub rozciągające w płaszczyźnie w kierunkach x i y, jako
jak również wszelkie siły ścinające [brak sił w płaszczyźnie]. Kompresja jest oznaczona jako ujemna
wartości, natomiast siły rozciągające są określane za pomocą wartości dodatnich.

-Kpryba
Zmień bieżącą wartość współczynnika Poissona [0.25].

-Cymłody
Zmień aktualną wartość modułu Younga [7.0e10 N/m^2].

-Fnu_a[/h_a/liczba_m]
Określ model jędrny w połączeniu z elastyczną grubością płyty
określone przez -E. Wystarczy podać jedną lepkość (nu_a) dla płyty elastycznej nad a
lepka półprzestrzeń, a także dołączanie grubości astenosfery (h_a) I
niższa lepkość płaszcza (liczba_m), przy czym pierwsza lepkość jest teraz lepkość
astenosfera. Podaj lepkości w Pa*s. Jeśli jest używany, podaj grubość
astenosfera w metrach; dodać k na km.

-N[f|q|s|nx/ny][+a|[+d|h|l][+e|n|m][+tszerokość][+w[przyrostek]][+z[p]]
Wybierz lub zapytaj o odpowiednie wymiary siatki dla FFT i ustaw opcjonalnie
parametry. Kontroluj wymiar FFT:
-Nf zmusi FFT do użycia rzeczywistych wymiarów danych.

-Nq zapyta o bardziej odpowiednie wymiary, zgłosi je, a następnie przejdzie dalej.

-Ns przedstawi listę opcjonalnych wymiarów, a następnie wyjdzie.

-Nnx/ny wykona FFT na rozmiarze tablicy nx/ny (musi być >= rozmiar pliku siatki). Domyślny
wybiera wymiary >= dane, które optymalizują szybkość i dokładność FFT. Jeśli FFT
wymiary > wymiary pliku siatki, dane są rozszerzane i zwężane do zera.

Kontroluj detrending danych: Dołącz modyfikatory do usuwania trendu liniowego:
+d: Detrend data, tj. usuń najlepiej pasujący trend liniowy [Domyślne].

+a: Usuń tylko wartość średnią.

+h: Usuń tylko wartość średnią, tj. 0.5 * (maks. + min.).

+l: Zostaw dane w spokoju.

Kontroluj rozszerzenie i zwężanie danych: użyj modyfikatorów, aby kontrolować sposób rozszerzenia
i zwężanie należy wykonać:
+e wydłuża siatkę poprzez nałożenie symetrii punktów brzegowych [Domyślnie],

+m rozszerza siatkę poprzez nałożenie symetrii lustrzanej krawędzi

+n wyłącza rozszerzenie danych.

Zbieżność jest wykonywana od krawędzi danych do krawędzi siatki FFT [100%]. Reszta
ten procent przez +tszerokość. Kiedy +n obowiązuje, stosuje się zwężanie
zamiast tego do marginesów danych, ponieważ nie jest dostępne rozszerzenie [0%].

Zapis kontrolny wyników tymczasowych: W celu szczegółowego zbadania można napisać
pośrednia siatka przekazywana do przodu FFT; prawdopodobnie tak było
odkształcony, wydłużony o symetrię punktową wzdłuż wszystkich krawędzi i zwężający się. Dodać
+w[przyrostek], z którego zostaną utworzone nazwy plików wyjściowych (tj. ingrid_prefix.ext)
[stożkowe], gdzie ext to rozszerzenie pliku. Na koniec możesz zapisać złożoną siatkę
produkowane przez FFT do przodu przez dołączenie +z. Domyślnie piszemy prawdziwe i
wyimaginowane komponenty do Ingrid_prawdziwy.ext i Ingrid_obraz.ext. Dodać p zapisać
zamiast tego biegunowa forma wielkości i fazy do plików Ingrid_mag.ext i
Ingrid_faza.ext.

-Lpodstęp Napisz nazwy i czasy oceny wszystkich siatek, które zostały utworzone do tekstu
filet podstęp. Wymaga -T.

-Mtm Określ model lepkosprężysty w połączeniu z elastyczną grubością płyty
określone przez -E. Dołącz czas Maxwella tm dla modelu lepkosprężystego (w ).

-Sbeta Określ zagłodzoną frakcję fosy w zakresie 0-1, gdzie 1 oznacza, że ​​fosa jest w pełni
wypełniony materiałem o gęstości ri podczas gdy 0 oznacza, że ​​jest wypełniony tylko materiałem
gęstość rw (tj. tylko woda) [1].

-Tt0[u] [/t1[ty]/dt[u]|n][+l]
Sprecyzować t0, t1i przyrost czasu (dt) dla sekwencji obliczeń [Domyślnie jest
jeden krok, bez zależności czasowej]. Na jeden konkretny czas, po prostu daj start
czas t0. Jednostką są lata; dodać k dla kyr i M dla Myr. Przez czas logarytmiczny
skalować, dołączać +l i określić n kroki zamiast dt. Ewentualnie podaj plik
z żądanymi czasami w pierwszej kolumnie (czasy te mogą mieć indywidualne jednostki)
w załączeniu, w przeciwnym razie zakładamy rok). Następnie piszemy osobny plik siatki modelu dla
każdy podany krok czasowy.

-Wwd Ustaw głębokość odniesienia dla niezdeformowanej wygiętej powierzchni wm [0]. Dodać k do
wskazać km.

-Zzm Określ głębokość odniesienia do wygiętej powierzchni (np. Moho) wm; dodać k na km. Musi
bądź pozytywnie nastawiony. [0].

-V[poziom] (więcej ...)
Wybierz poziom szczegółowości [c].

-fg Siatki geograficzne (wymiary długości i szerokości geograficznej) zostaną przeliczone na metry
poprzez przybliżenie „Płaskiej Ziemi” przy użyciu bieżących parametrów elipsoidy.

-^ or właśnie -
Wydrukuj krótką wiadomość o składni polecenia, a następnie wyjdź (UWAGA: w systemie Windows
użyj tylko -).

-+ or właśnie +
Wydrukuj obszerny komunikat dotyczący użytkowania (pomocy), w tym wyjaśnienie wszelkich
opcja specyficzna dla modułu (ale nie typowe opcje GMT), a następnie kończy działanie.

-? or Nie argumenty
Wydrukuj pełny komunikat (pomoc) dotyczący użytkowania, w tym wyjaśnienie opcji, a następnie
wyjść.

--wersja
Wydrukuj wersję GMT i wyjdź.

--show-katalog danych
Wydrukuj pełną ścieżkę do katalogu współdzielonego GMT i wyjdź.

GRID FILE FORMATY

Domyślnie GMT zapisuje siatkę jako zmiennoprzecinkową o pojedynczej precyzji w netCDF z reklamacją COARDS
format pliku. Jednak GMT jest w stanie tworzyć pliki grid w wielu innych powszechnie używanych siatkach
formatów plików, a także ułatwia tzw. „pakowanie” siatek, wypisywanie zmiennoprzecinkowych
dane jako 1- lub 2-bajtowe liczby całkowite. Aby określić precyzję, skalę i przesunięcie, użytkownik powinien:
dodaj przyrostek =id[/skala/offset[/nan]], gdzie id jest dwuliterowym identyfikatorem siatki
rodzaj i precyzja oraz skala i offset są opcjonalnym współczynnikiem skali i przesunięciem, które mają być
zastosowane do wszystkich wartości siatki, oraz nan to wartość używana do wskazania brakujących danych. W razie
dwie postacie id nie jest dostarczony, jak w =/skala niż a id=nf zakłada się. Kiedy
czytanie siatek, format jest zazwyczaj automatycznie rozpoznawany. Jeśli nie, ten sam przyrostek
można dodać do nazw plików siatki wejściowej. Widzieć grdconvert i format pliku siatki sekcji
Informacje techniczne GMT i książka kucharska, aby uzyskać więcej informacji.

Podczas odczytywania pliku netCDF, który zawiera wiele siatek, GMT domyślnie odczyta
pierwsza dwuwymiarowa siatka, którą można znaleźć w tym pliku. Aby nakłonić GMT do przeczytania innego
wielowymiarowa zmienna w pliku siatki, dołącz ?nazwa_zmiennej do nazwy pliku, gdzie
nazwa_zmiennej to nazwa zmiennej. Pamiętaj, że być może będziesz musiał uciec od specjalnego znaczenia
of ? w swoim programie powłoki, umieszczając przed nim ukośnik odwrotny lub umieszczając znak
nazwę pliku i przyrostek między cudzysłowami lub cudzysłowami. ten ?nazwa_zmiennej można również użyć sufiksu
dla siatek wyjściowych, aby określić nazwę zmiennej inną niż domyślna: „z”. Widzieć
grdconvert i sekcje modyfikatory-for-CF i format pliku siatki GMT Technical
Odnośnik i książka kucharska, aby uzyskać więcej informacji, w szczególności o tym, jak czytać sploty 3-,
Siatki 4- lub 5-wymiarowe.

GRID DYSTANS JEDNOSTKI

Jeśli siatka nie ma licznika jako jednostki poziomej, dołącz +ujednostka do pliku wejściowego
nazwa do konwersji z określonej jednostki na metr. Jeśli twoja siatka jest geograficzna, przekonwertuj
odległości na metry przez dostarczanie -fg zamiast.

UWAGI

Siatki netCDF COARDS zostaną automatycznie rozpoznane jako geograficzne. Dla innych siatek
siatki geograficzne, w których chcesz zamienić stopnie na metry, wybierz -fg. Jeśli dane
są blisko jednego z biegunów, powinieneś rozważyć rzutowanie pliku siatki na prostokąt
układ współrzędnych za pomocą grdprojekt.

PŁYTA ZGIĘCIE UWAGI

Rozwiązanie FFT do zginania płyty wymaga gęstości wypełnienia równej gęstości obciążenia.
Zwykle dzieje się tak tylko bezpośrednio pod obciążeniem; poza obciążeniem wypełnienie ma tendencję do przemieszczania się
być osadami o mniejszej gęstości lub nawet wodą (lub powietrzem). Wessel [2001] zaproponował
przybliżenie, które pozwala na określenie gęstości wypełnienia innej niż
gęstość obciążenia, jednocześnie pozwalając na rozwiązanie FFT. Zasadniczo wygięcie płyty jest
rozwiązano za pomocą gęstości wypełnienia jako efektywnej gęstości obciążenia, ale amplitudy są
skorygowany przez współczynnik A = sqrt ((rm - ri)/(rm - rl)), czyli różnica teoretyczna
w amplitudzie ze względu na obciążenie punktowe przy użyciu dwóch różnych gęstości obciążenia. ten
przybliżenie jest bardzo dobre, ale psuje się przy dużych obciążeniach na słabych talerzach, wróżka
rzadka sytuacja.

PRZYKŁADY

Aby obliczyć ugięcie płyty sprężystej z obciążenia topo.nc, dla płyty o grubości 10 km z
typowe gęstości, spróbuj

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

Aby obliczyć odpowiedź lepkości jędrnej na serię przyrostowych obciążeń podanych przez nazwę pliku
i czas obciążenia w tabeli l.lis w jednym czasie 1 Ma przy użyciu określonego reologicznego
wartości, spróbuj

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

LITERATURA

Bydło, LM, 1975, Kurs lepkość of dotychczasowy ziemski płaszcz, Wydawnictwo Uniwersytetu Princeton.

Wessela. P., 2001, Globalny rozkład gór podwodnych wywnioskowany z siatki Geosat/ERS-1
wysokościomierz, J. Geophys. Res., 106(B9), 19,431-19,441,
http://dx.doi.org/10.1029/2000JB000083

Korzystaj z grdflexuregmt online za pomocą usług onworks.net