grdflexuregmt - Dalam Talian di Awan

JADUAL:

NAMA

grdflexure - Kira ubah bentuk lentur permukaan 3-D untuk pelbagai reologi

SINOPSIS

grdflexure topogrd rm/rl[/ri]/rw Te[u] outgrid [ Nx/Ny/Nxy ] [ pikan ] [ yYoung ] [
nu_a[/h_a/nu_m] ] [ senarai ] [ [f|q|s|nx/ny][+a|d|h|l][+e|n|m][+tlebar][+w[akhiran]][+z[p]] [
beta ] [ -Tt0[u][/t1[u]/dt[u]|n][+l] ] [ [tahap] ] [ wd] [ zm] [ -fg ]

Catatan: Tiada ruang dibenarkan antara bendera pilihan dan hujah yang berkaitan.

DESCRIPTION

grdflexure mengira tindak balas lentur kepada beban menggunakan julat yang boleh dipilih pengguna
reologi. Pengguna boleh memilih daripada elastik, viskoelastik atau firmoviscous (dengan satu atau dua
lapisan likat). Evolusi temporal juga boleh dimodelkan dengan menyediakan beban tambahan
grid dan menentukan julat masa keluaran model.

DIKEHENDAKI HUJAH

topogrd
Fail grid binari 2-D dengan topografi beban (dalam meter); Lihat FAIL GRID
FORMAT di bawah. Jika -T digunakan, topogrd mungkin templat nama fail dengan terapung
format titik (sintaks C) dan nama fail beban yang berbeza akan ditetapkan dan dimuatkan untuk
setiap kali langkah. Oleh itu, masa memuatkan bertepatan dengan masa yang diberikan melalui -T (tetapi tidak
setiap masa perlu mempunyai fail yang sepadan). Sebagai alternatif, berikan topogrd as
=*flist*, di mana flist ialah jadual ASCII dengan satu topogrd nama fail dan masa muat per
rekod. Masa muat ini boleh berbeza daripada masa penilaian yang diberikan melalui -T.
Untuk format masa muat, lihat -T.

-Drm/rl[/ri]/rw
Menetapkan ketumpatan untuk mantel, beban, isian (pilihan, jika tidak, ia diandaikan sama
ketumpatan beban), dan air atau udara. Jika ri berbeza daripada rl kemudian anggaran
penyelesaian akan ditemui. Jika ri tidak diberikan maka ia lalai kepada rl.

-ETe Menetapkan ketebalan plat elastik (dalam meter); tambahkan k untuk km. Jika anjal
ketebalan melebihi 1e10 ia akan ditafsirkan sebagai ketegaran lentur D (secara lalai
D dikira daripada Te, modulus Young, dan nisbah Poisson; lihat -C untuk menukar ini
nilai).

-Gfail luar
If -T ditetapkan kemudian grdfile mestilah templat nama fail yang mengandungi terapung
format titik (sintaks C). Jika templat nama fail juga mengandungi sama ada %s (untuk
nama unit) atau %c (untuk huruf unit) maka kita menggunakan masa yang sepadan (dalam unit
dinyatakan dalam -T) untuk menjana nama fail individu, jika tidak, kami menggunakan masa masuk
tahun tanpa unit.

PILIHAN HUJAH

-ANx/Ny/Nxy
Nyatakan daya mampatan atau lanjutan dalam satah dalam arah x dan y, sebagai
serta sebarang daya ricih [tiada daya dalam satah]. Mampatan ditunjukkan dengan negatif
nilai, manakala daya lanjutan ditentukan menggunakan nilai positif.

-Cpikan
Tukar nilai semasa nisbah Poisson [0.25].

-CyYoung
Tukar nilai semasa modulus Young [7.0e10 N/m^2].

-Fnu_a[/h_a/nu_m]
Nyatakan model firmoviscous bersama dengan ketebalan plat elastik
dinyatakan melalui -E. Hanya berikan satu kelikatan (nu_a) untuk plat kenyal di atas a
separuh ruang likat, atau juga menambah ketebalan astenosfera (h_a) Dan
kelikatan mantel yang lebih rendah (nu_m), dengan kelikatan pertama sekarang ialah kelikatan
astenosfera. Beri kelikatan dalam Pa*s. Jika digunakan, berikan ketebalan
astenosfera dalam meter; tambahkan k untuk km.

-N[f|q|s|nx/ny][+a|[+d|h|l][+e|n|m][+tlebar][+w[akhiran][+z[p]]
Pilih atau tanya tentang dimensi grid yang sesuai untuk FFT dan tetapkan pilihan
parameter. Kawal dimensi FFT:
-Nf akan memaksa FFT menggunakan dimensi sebenar data.

-Nq akan bertanya tentang dimensi yang lebih sesuai, melaporkannya, kemudian meneruskan.

-NS akan membentangkan senarai dimensi pilihan, kemudian keluar.

-Nnx/ny akan melakukan FFT pada saiz tatasusunan nx/ny (mesti >= saiz fail grid). lalai
memilih dimensi >= data yang mengoptimumkan kelajuan dan ketepatan FFT. Jika FFT
dimensi > dimensi fail grid, data dilanjutkan dan ditiruskan kepada sifar.

Kawal detrending data: Tambahkan pengubah suai untuk mengalih keluar aliran linear:
+d: Detrend data, iaitu alih keluar aliran linear paling sesuai [Default].

+a: Hanya keluarkan nilai min.

+h: Hanya keluarkan nilai pertengahan, iaitu 0.5 * (maks + min).

+l: Biarkan data sahaja.

Kawal sambungan dan pengecilan data: Gunakan pengubah suai untuk mengawal cara sambungan
dan tirus perlu dilakukan:
+e memanjangkan grid dengan mengenakan simetri titik tepi [Lalai],

+m memanjangkan grid dengan mengenakan simetri cermin tepi

+n mematikan sambungan data.

Tirus dilakukan dari tepi data ke tepi grid FFT [100%]. Ubah
peratusan ini melalui +tlebar. Apabila +n sedang berkuat kuasa, tirus digunakan
sebaliknya ke margin data kerana tiada sambungan tersedia [0%).

Kawalan penulisan keputusan sementara: Untuk penyiasatan terperinci anda boleh menulis
grid perantaraan dihantar ke FFT hadapan; ini mungkin telah berlaku
detrend, dilanjutkan dengan simetri titik di sepanjang semua tepi, dan tirus. Lampirkan
+w[akhiran] dari mana nama fail output akan dibuat (iaitu, ingrid_prefix.ext)
[tirus], di mana ext ialah sambungan fail anda. Akhir sekali, anda boleh menyimpan grid kompleks
dihasilkan oleh FFT hadapan dengan melampirkan +z. Secara lalai kita menulis yang sebenar dan
komponen khayalan ke masuk_sebenar.ext and masuk_imag.ext. Lampirkan p untuk menyelamatkan
sebaliknya bentuk polar magnitud dan fasa kepada fail masuk_mag.ext and
masuk_fasa.ext.

-Lsenarai Tulis nama dan masa penilaian semua grid yang dicipta pada teks
fail senarai. Memerlukan -T.

-Mtm Tentukan model viskoelastik bersama dengan ketebalan plat kenyal
dinyatakan melalui -E. Tambah masa Maxwell tm untuk model viskoelastik (dalam ).

-Sbeta Nyatakan pecahan parit kelaparan dalam julat 0-1, dengan 1 bermakna parit telah penuh
diisi dengan bahan ketumpatan ri manakala 0 bermakna ia hanya diisi dengan bahan daripada
ketumpatan rw (iaitu, hanya air) [1].

-Tt0[u][/t1[u]/dt[u]|n][+l]
Tentukan t0, t1, dan kenaikan masa (dt) untuk urutan pengiraan [Default is
satu langkah, tanpa pergantungan masa]. Untuk satu masa tertentu, berikan sahaja permulaan
masa t0. Unitnya ialah tahun; tambahkan k untuk kyr dan M untuk Myr. Untuk masa logaritma
skala, tambah +l dan nyatakan n langkah bukannya dt. Sebagai alternatif, berikan fail
dengan masa yang dikehendaki dalam lajur pertama (masa ini mungkin mempunyai unit individu
dilampirkan, jika tidak kita anggap tahun). Kami kemudian menulis fail grid model berasingan untuk
setiap langkah masa yang diberikan.

-Wwd Tetapkan kedalaman rujukan kepada permukaan lentur yang tidak cacat dalam m [0]. Lampirkan k kepada
menunjukkan km.

-Zzm Tentukan kedalaman rujukan kepada permukaan lentur (cth., Moho) dalam m; tambahkan k untuk km. mesti
bersikap positif. [0].

-V[tahap] (lagi ...)
Pilih tahap verbositi [c].

-fg Grid geografi (dimensi longitud, latitud) akan ditukar kepada meter
melalui anggaran "Bumi Rata" menggunakan parameter ellipsoid semasa.

-^ or hanya -
Cetak mesej ringkas tentang sintaks arahan, kemudian keluar (NOTA: pada Windows
gunakan hanya -).

-+ or hanya +
Cetak mesej penggunaan (bantuan) yang meluas, termasuk penjelasan sebarang
pilihan khusus modul (tetapi bukan pilihan biasa GMT), kemudian keluar.

-? or tidak hujah
Cetak mesej penggunaan (bantuan) lengkap, termasuk penjelasan pilihan, kemudian
jalan keluar.

--versi
Cetak versi GMT dan keluar.

--show-datadir
Cetak laluan penuh ke direktori kongsi GMT dan keluar.

GRID FAIL FORMATS

Secara lalai GMT menulis grid apabila ketepatan tunggal terapung dalam netCDF aduan COARDS
format fail. Walau bagaimanapun, GMT mampu menghasilkan fail grid dalam banyak grid lain yang biasa digunakan
format fail dan juga memudahkan apa yang dipanggil "pembungkusan" grid, menulis titik terapung
data sebagai integer 1- atau 2-bait. Untuk menentukan ketepatan, skala dan ofset, pengguna harus
tambah imbuhan =id[/skala/mengimbangi[/nan]], di mana id ialah pengecam dua huruf bagi grid
jenis dan ketepatan, dan skala and mengimbangi adalah faktor skala pilihan dan diimbangi menjadi
digunakan pada semua nilai grid, dan nan ialah nilai yang digunakan untuk menunjukkan data yang hilang. Dalam kes
dua watak itu id tidak disediakan, seperti dalam =/skala daripada a id=nf diandaikan. Bila
grid membaca, format biasanya diiktiraf secara automatik. Jika tidak, akhiran yang sama
boleh ditambah pada nama fail grid input. Lihat grdconvert dan Bahagian grid-fail-format bagi
Rujukan Teknikal dan Buku Masakan GMT untuk maklumat lanjut.

Apabila membaca fail netCDF yang mengandungi berbilang grid, GMT akan membaca, secara lalai, fail
grid 2 dimensi pertama yang boleh ditemui dalam fail itu. Untuk memujuk GMT supaya membaca yang lain
pembolehubah berbilang dimensi dalam fail grid, tambahkan ?varname kepada nama fail, di mana
varname ialah nama pembolehubah. Ambil perhatian bahawa anda mungkin perlu melarikan diri daripada makna istimewa
of ? dalam program shell anda dengan meletakkan garis miring ke belakang di hadapannya, atau dengan meletakkan
nama fail dan akhiran antara petikan atau petikan berganda. The ?varname akhiran juga boleh digunakan
untuk grid output untuk menentukan nama pembolehubah yang berbeza daripada lalai: "z". Lihat
grdconvert dan Bahagian pengubahsuai-untuk-CF dan format fail grid bagi Teknikal GMT
Rujukan dan Buku Masakan untuk maklumat lanjut, terutamanya tentang cara membaca sambatan 3-,
Grid 4-, atau 5 dimensi.

GRID DISTANCE UNIT

Jika grid tidak mempunyai meter sebagai unit mendatar, tambahkan +uunit kepada fail input
nama untuk ditukar daripada unit yang ditentukan kepada meter. Jika grid anda adalah geografi, tukar
jarak ke meter dengan membekalkan -fg sebaliknya.

KONSIDERASI

grid netCDF COARDS secara automatik akan diiktiraf sebagai geografi. Untuk grid lain
grid geografi jika anda ingin menukar darjah kepada meter, pilih -fg. Jika data
hampir dengan mana-mana tiang, anda harus mempertimbangkan untuk menayangkan fail grid ke segi empat tepat
sistem koordinat menggunakan grdproject.

PINGGANG FLEKTUR NOTA

Penyelesaian FFT kepada lentur plat memerlukan ketumpatan isian untuk menyamai ketumpatan beban.
Ini biasanya hanya benar di bawah beban; melebihi beban yang cenderung kepada isian
menjadi sedimen berketumpatan rendah atau air (atau udara). Wessel [2001] mencadangkan satu
anggaran yang membolehkan spesifikasi ketumpatan isian berbeza daripada
ketumpatan beban sementara masih membenarkan penyelesaian FFT. Pada asasnya, lentur plat adalah
diselesaikan untuk menggunakan ketumpatan isian sebagai ketumpatan beban berkesan tetapi amplitudnya
diselaraskan oleh faktor A = sqrt ((rm - ri)/(rm - rl)), iaitu perbezaan teori
dalam amplitud disebabkan oleh beban titik menggunakan dua ketumpatan beban yang berbeza. The
anggaran adalah sangat baik tetapi rosak untuk beban besar pada plat lemah, dongeng
keadaan yang tidak biasa.

CONTOH

Untuk mengira lentur plat elastik daripada beban topo.nc, untuk plat setebal 10 km dengan
ketumpatan tipikal, cuba

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

Untuk mengira tindak balas firmoviscous kepada siri beban tambahan yang diberikan oleh nama fail
dan masa muat dalam jadual l.lis pada masa tunggal 1 Ma menggunakan reologi yang ditentukan
nilai, cuba

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

RUJUKAN

Cathles, LM, 1975, . kelikatan of yang milik bumi mantel, Akhbar Universiti Princeton.

Wessel. P., 2001, Taburan global bagi gunung laut yang disimpulkan daripada Geosat/ERS-1 bergrid
altimetri, J. Geophys. Res., 106(B9), 19,431-19,441,
http://dx.doi.org/10.1029/2000JB000083

Gunakan grdflexuregmt dalam talian menggunakan perkhidmatan onworks.net