gravfftgmt - Dalam Talian di Awan

JADUAL:

NAMA

gravfft - Kira daya tarikan graviti permukaan 3-D dalam nombor gelombang (atau
kekerapan) domain

SINOPSIS

gravfft masuk [ ingrid2 ] fail luar [ n/panjang gelombang/min_depth/tbw ] [ ketumpatan|rhogrid ] [
n_syarat ] [ [f[+]|g|v|n|e] ] [ w|b|c|t |k ] [ [f|q|s|nx/ny][+a|d|h
|l][+e|n|m][+tlebar][+w[akhiran]][+z[p]] [ ] [ te/rl/rm/rw[+m] ] [ [tahap] ] [ wd] [
zm[zl] ] [ -fg ]

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

DESCRIPTION

gravfft boleh digunakan dalam tiga mod utama. Mod 1: Hanya mengira geopotensi disebabkan oleh
permukaan yang diberikan dalam fail topo.grd. Memerlukan kontras ketumpatan (-D) dan mungkin a
tahap pemerhatian yang berbeza (-W). Ia akan mengambil FFT hadapan 2-D grid dan digunakan
kaedah penuh Parker sehingga terma yang dipilih. Mod 2: Kirakan geopotensi
tindak balas disebabkan lenturan fail topografi. Ia akan mengambil FFT hadapan 2-D bagi
grid dan gunakan kaedah penuh Parker yang digunakan pada model isostatik yang dipilih. The
model yang tersedia ialah "pemuatan dari atas", atau model plat elastik, dan "pemuatan dari
bawah" yang merangkumi tindak balas plat kepada beban bawah permukaan (sesuai untuk panas
pemodelan spot - jika anda mempercayainya). Dalam kedua-dua kes, parameter model ditetapkan dengan -T
and -Z pilihan. Mod 3: hitung kemasukan atau koheren antara dua grid. Keluaran
ialah purata dalam arah jejari. Secara pilihan, kemasukan model juga mungkin
dikira. Dimensi mendatar grdfiles diandaikan dalam meter.
Grid geografi boleh digunakan dengan menyatakan -fg pilihan yang skala darjah kepada meter.
Jika anda mempunyai grid dengan dimensi dalam km, anda boleh menukar ini kepada meter menggunakan grdedit or
skala output dengan grdmath. Memandangkan bilangan pilihan yang ditawarkan program ini, adalah
sukar untuk menyatakan apakah pilihan dan apakah hujah yang diperlukan. Ia bergantung kepada apa
anda sedang lakukan; lihat contoh untuk panduan selanjutnya.

DIKEHENDAKI HUJAH

masuk Fail grid binari 2-D untuk dikendalikan. (Lihat FORMAT FAIL GRID di bawah). Untuk
operasi spektrum silang, juga memberikan fail grid kedua ingrd2.

-Gfail luar
Nyatakan nama fail grid output atau jadual spektrum 1-D (lihat -E). (Lihat
FORMAT FAIL GRID di bawah).

PILIHAN HUJAH

-Cn/panjang gelombang/min_depth/tbw
Kira hanya keluk kemasukan teori bagi model dan keluar yang dipilih. n
and panjang gelombang digunakan untuk mengira (n * panjang gelombang) jumlah panjang profil dalam
meter. min_depth ialah purata kedalaman air. Tambahkan bendera data (satu atau dua) daripada tbw in
sebarang pesanan. t = gunakan model "dari atas", b = gunakan model "dari bawah". Tentukan secara pilihan
w untuk menulis panjang gelombang dan bukannya frekuensi.

-Dketumpatan|rhogrid
Menetapkan kontras ketumpatan merentasi permukaan. Digunakan, sebagai contoh, untuk mengira graviti
tarikan lapisan air yang kemudiannya boleh digabungkan dengan anomali udara bebas
untuk mendapatkan anomali Bouguer. Dalam kes ini jangan gunakan -T. Ia juga secara tersirat menetapkan
-N+h. Sebagai alternatif, tentukan grid yang didaftarkan bersama dengan kontras ketumpatan jika a
kontras ketumpatan berubah-ubah diperlukan.

-En_syarat
Bilangan istilah yang digunakan dalam pengembangan Parker (had ialah 10, sebaliknya istilah bergantung pada
n akan meletupkan atur cara) [Lalai = 3]

-F[f[+]|g|v|n|e]
Nyatakan medan geopotensi yang diingini: hitung geoid dan bukannya graviti
f = Anomali udara bebas (mGal) [Lalai]. Lampirkan + untuk menambah papak yang tersirat
apabila mengeluarkan nilai min daripada topografi. Ini memerlukan topografi sifar
bermaksud tiada anomali jisim.

g = Anomali geoid (m).

v = Gradien Graviti Menegak (VGG; 1 Eotvos = 0.1 mGal/km).

e = Pesongan timur menegak (mikro-radian).

n = Pesongan utara menegak (mikro-radian).

-Iw|b|c|t |k
Penggunaan ingrd2 and ingrd1 (grid dengan topografi/batimetri) untuk dianggarkan
kemasukan|kepaduan dan tuliskannya ke stdout (-G diabaikan jika ditetapkan). Grid ini sepatutnya
mengandungi graviti atau geoid untuk kawasan yang sama ingrd1. Pengiraan lalai
kemasukan. Output mengandungi 3 atau 4 lajur. Kekerapan (panjang gelombang), kemasukan
(koheren) satu bar ralat sigma dan, secara pilihan, kemasukan teori. Lampirkan
dataflags (satu hingga tiga) daripada w|b|c|t. w menulis panjang gelombang bukannya nombor gelombang,
k memilih km untuk unit panjang gelombang [m], c mengira koheren bukannya kemasukan, b
menulis lajur keempat dengan kemasukan teori "memuatkan dari bawah", dan t
menulis lajur keempat dengan kemasukan teori "plat elastik".

-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.

-Q Menulis grid dengan topografi lentur (dengan z positif ke atas) yang puratanya
telah ditetapkan oleh -Zzm dan parameter model oleh -T (dan output oleh -G). Itulah
"gravimetrik Moho". -Q secara tersirat menetapkan -N+h

-S Kira graviti atau grid geoid yang diramalkan disebabkan oleh beban subplat yang dihasilkan oleh
batimetri semasa dan model teori. Parameter yang diperlukan ditetapkan
dalam -T and -Z pilihan. Bilangan kuasa dalam pengembangan Parker adalah terhad kepada
1. Lihat contoh lebih jauh ke bawah.

-Tte/rl/rm/rw[+m]
Kira pampasan isostatik daripada beban topografi (fail grid input) pada
plat elastik ketebalan te. Tambah juga ketumpatan untuk beban, mantel dan air masuk
unit SI. Berikan purata kedalaman mantel melalui -Z. Jika ketebalan anjal ialah > 1e10 ia
akan ditafsirkan sebagai ketegaran lentur (secara lalai ia dikira daripada te and
Modulus muda). Secara pilihan, tambahkan +m untuk menulis grid dengan geopotensi Moho
kesan (lihat -F) daripada model yang dipilih oleh -T. Jika te = 0 maka respons Airy ialah
kembali. -T+m secara tersirat menetapkan -N+h

-Wwd Tetapkan kedalaman air (atau ketinggian cerapan) berbanding dengan topografi [0]. Lampirkan k kepada
menunjukkan km.

-Zzm[zl]
Moho [dan membengkak] purata kedalaman pampasan. Untuk model "beban dari atas" anda
hanya perlu menyediakan zm, tetapi untuk "pemuatan dari bawah" jangan lupa zl.

-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

Penyelesaian FFT kepada lentur plat elastik memerlukan ketumpatan isian untuk menyamai beban
ketumpatan. Ini biasanya hanya benar di bawah beban; di luar beban yang
infill cenderung menjadi sedimen berketumpatan rendah atau bahkan air (atau udara). Wessel [2001] mencadangkan
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 kesan lapisan air di atas kelawar.grd batimetri menggunakan 2700 dan 1035
untuk ketumpatan kerak dan air dan menulis hasilnya pada water_g.grd (mengira
kepada kuasa keempat batimetri dalam pengembangan Parker):

gmt gravfft bat.grd -D1665 -Gwater_g.grd -E4

Sekarang tolakkannya kepada faa.grd anomali udara bebas anda dan anda akan mendapat anomali Bouguer. awak
mungkin tertanya-tanya mengapa kita menolak dan tidak menambah. Lagipun anomali Bouguer berpura-pura
untuk membetulkan kekurangan jisim yang dibentangkan oleh lapisan air, jadi kita harus menambah kerana
air kurang tumpat daripada batuan di bawah. Jawapannya bergantung pada cara kesan graviti
dikira dengan kaedah Parker dan aspek praktikal menggunakan FFT.

gmt grdmath faa.grd water_g.grd SUB = bouguer.grd

Mahukan anomali MBA? Kira sumbangan mantel kerak dan tambahkannya pada
anomali dasar laut. Dengan mengandaikan kerak tebal 6 km berketumpatan 2700 dan mantel dengan 3300
ketumpatan kita boleh mengulangi arahan yang digunakan untuk mengira anomali lapisan air, menggunakan 600
(3300 - 2700) sebagai kontras ketumpatan. Tetapi kita sekarang mempunyai masalah kerana kita perlu tahu
min kedalaman Moho. Iaitu apabila skala/offset yang boleh dilampirkan pada nama grid
datang dalam tangan. Perhatikan bahawa kami tidak perlu berbuat demikian sebelum ini kerana purata kedalaman air adalah
dikira terus daripada data (perhatikan juga tanda negatif offset disebabkan oleh fakta
Bahawa z adalah positif):

gmt gravfft bat.grd=nf/1/-6000 -D600 -Gmoho_g.grd

Sekarang, tolakkannya kepada anomali dasar laut untuk mendapatkan anomali MBA. Itu dia:

gmt grdmath water_g.grd moho_g.grd SUB = mba.grd

Untuk mengira kesan graviti Moho bagi kelawar plat anjal.grd dengan Te = 7 km, ketumpatan
2700, di atas mantel berketumpatan 3300, pada kedalaman purata 9 km

gmt gravfft bat.grd -Gelastic.grd -T7000/2700/3300/1035+m -Z9000

Jika anda menambah sekarang dasar laut dan kesan Moho, anda akan mendapat tindak balas graviti penuh
model isostatik anda. Kami akan menggunakan di sini hanya penggal pertama dalam pengembangan Parker.

gmt gravfft bat.grd -D1665 -Gwater_g.grd -E1
gmt gravfft bat.grd -Gelastic.grd -T7000/2700/3300/1035+m -Z9000 -E1
gmt grdmath water_g.grd elastik.grd ADD = model.grd

Hasil yang sama boleh diperoleh secara langsung dengan arahan seterusnya. Walau bagaimanapun, PERHATIAN kepada
yang berikut. Saya belum tahu sama ada ia disebabkan oleh pepijat atau disebabkan oleh beberapa had, tetapi
hakikatnya arahan berikut dan sebelumnya hanya memberikan hasil yang sama jika -E1
digunakan. Untuk kuasa batimetri yang lebih tinggi dalam pengembangan Parker, hanya contoh di atas
jahitan untuk memberikan hasil yang betul.

gmt gravfft bat.grd -Gmodel.grd -T7000/2700/3300/1035 -Z9000 -E1

Dan apakah anomali geoid yang dihasilkan oleh beban pada kedalaman 50 km, di bawah kawasan a
yang batimetri diberikan oleh bat.grd, Moho pada kedalaman 9 km dan ketumpatan yang sama seperti
sebelum?

gmt gravfft topo.grd -Gswell_geoid.grd -T7000/2700/3300/1035 -Fg -Z9000/50000 -S -E1

Untuk mengira kemasukan antara batimetri topo.grd dan anomali udara bebas faa.grd
grid menggunakan model plat anjal kerak 6 km ketebalan min dengan 10 km berkesan
ketebalan elastik di kawasan 3 km purata kedalaman air:

gmt gravfft topo.grd faa.grd -Ia -T10000/2700/3300/1035 -Z9000

Untuk mengira kemasukan antara batimetri topo.grd dan grid geoid.grd dengan
model "pemuatan dari bawah" (LFB) dengan beban yang sama seperti di atas dan di bawah permukaan pada 40 km,
tetapi dengan mengandaikan kini grid berada dalam geografi dan kami mahukan panjang gelombang dan bukannya kekerapan:

gmt gravfft topo.grd geoid.grd -Ibw -T10000/2700/3300/1035 -Z9000/40000 -fg

Untuk mengira kemasukan teori graviti LFB sepanjang profil 2000 km menggunakan
parameter yang sama seperti di atas

gmt gravfft -C400/5000/3000/b -T10000/2700/3300/1035 -Z9000/40000

RUJUKAN

Luis, JF dan MC Neves. 2006, Pampasan isostatik Dataran Tinggi Azores: 3D
analisis kemasukan dan koheren. J. Geoterma Volc. Res. Jilid 156, Isu 1-2, Halaman
10-22, http://dx.doi.org/10.1016/j.jvolgeores.2006.03.010 Parker, RL, 1972, The rapid
pengiraan potensi anomali, Geophys. J., 31, 447-455. Wessel. P., 2001, Global
taburan gunung laut yang disimpulkan daripada altimetri Geosat/ERS-1 grid, J. Geophys. Res.,
106(B9), 19,431-19,441, http://dx.doi.org/10.1029/2000JB000083

Gunakan gravfftgmt dalam talian menggunakan perkhidmatan onworks.net