gravfftgmt - Bulutta Çevrimiçi

Program:

ADI

gravfft - Dalga sayısında (veya
frekans) etki alanı

SİNOPSİS

yerçekimi ingrid [ ingrid2 ] dış dosya [ n/dalga boyu/ortalama_derinlik/tbw ] [ yoğunluk|rogrid ] [
n_terms ] [ [f[+]|g|v|n|e] ] [ w|b|c|t |k ] [ [f|q|s|nx/ny][+a|d|h
|l][+e|n|m][+tgenişlik][+w[sonek]][+z[p]] [ ] [ te/rl/rm/rw[+m] ] [seviye] ] [ wd] [
zm[zl] ] [ -fg ]

Not: Seçenek bayrağı ve ilişkili bağımsız değişkenler arasında boşluk bırakılamaz.

TANIM

yerçekimi üç ana modda kullanılabilir. Mod 1: Basitçe, aşağıdakilerden kaynaklanan jeopotansiyeli hesaplayın
topo.grd dosyasında verilen yüzey. Yoğunluk kontrastı gerektirir (-D) ve muhtemelen bir
farklı gözlem seviyesi (-W). Izgaranın 2 boyutlu ileri FFT'sini alacak ve
seçilen koşullara kadar tam Parker yöntemi. Mod 2: Jeopotansiyelin hesaplanması
topografya dosyasının bükülmesinden kaynaklanan tepki. 2 boyutlu ileri FFT'yi alacaktır.
ızgara yapın ve seçilen izostatik modele uygulanan tam Parker yöntemini kullanın. NS
mevcut modeller "üstten yükleme" veya elastik plaka modeli ve "yukarıdan yükleme" modelidir.
altında", plakanın bir yüzey altı yüküne tepkisini açıklar (sıcak için uygun
nokta modelleme - onlara inanıyorsanız). Her iki durumda da model parametreleri şu şekilde ayarlanır: -T
ve -Z seçenekler. Mod 3: iki ızgara arasındaki kabulü veya tutarlılığı hesaplayın. Çıktı
radyal yönde ortalamadır. İsteğe bağlı olarak model kabul de olabilir.
hesaplandı. grd dosyalarının yatay boyutlarının metre cinsinden olduğu varsayılır.
Coğrafi ızgaralar, aşağıdakiler belirtilerek kullanılabilir: -fg dereceyi metreye ölçekleyen seçenek.
Boyutları km cinsinden olan ızgaralarınız varsa, bunu kullanarak bunu metre olarak değiştirebilirsiniz. grededit or
çıktıyı ölçeklendir matematik. Bu programın sunduğu seçeneklerin sayısı göz önüne alındığında,
seçeneklerin ve gerekli argümanların ne olduğunu belirtmek zor. Neye bağlı
sen yapıyorsun; daha fazla rehberlik için örneklere bakın.

GEREKLİ ARGÜMANLAR

ingrid Çalıştırılacak 2 boyutlu ikili ızgara dosyası. (Aşağıdaki IZGARA DOSYA FORMATLARI'na bakın). İçin
çapraz spektral işlemler, ayrıca ikinci ızgara dosyasını verir inrd2.

-Gdış dosya
Çıktı ızgara dosyasının veya 1-D spektrum tablosunun adını belirtin (bkz. -E). (Görmek
IZGARA DOSYA FORMATLARI aşağıdadır).

İSTEĞE BAĞLI ARGÜMANLAR

-Cn/dalga boyu/ortalama_derinlik/tbw
Seçilen modelin sadece teorik kabul eğrilerini hesaplayın ve çıkın. n
ve dalga boyu toplam profil uzunluğunu (n * dalga boyu) hesaplamak için kullanılır
metre. ortalama_derinlik ortalama su derinliğidir. Veri bayraklarını (bir veya iki) ekleyin TVS in
herhangi bir sipariş. t = "üstten" modeli kullanın, b = "aşağıdan" modelini kullanın. İsteğe bağlı olarak belirtin
w frekans yerine dalga boyu yazmak için.

-Dyoğunluk|rogrid
Yüzey boyunca yoğunluk kontrastını ayarlar. Örneğin, yerçekimini hesaplamak için kullanılır
daha sonra serbest hava anomalisi ile birleştirilebilen su tabakasının çekiciliği
Bouguer anomalisini elde etmek için. Bu durumda kullanmayın -T. Ayrıca dolaylı olarak ayarlar
-N+sa. Alternatif olarak, eğer bir
değişken yoğunluklu kontrast gereklidir.

-En_terms
Parker genişletmesinde kullanılan terim sayısı (sınır 10'dur, aksi takdirde terimler
n programı patlatacak) [Varsayılan = 3]

-F[f[+]|g|v|n|e]
İstenen jeopotansiyel alanı belirtin: yerçekimi yerine jeoidi hesaplayın
f = Serbest hava anormallikleri (mGal) [Varsayılan]. Ekle + ima edilen levhaya eklemek
ortalama değeri topografyadan çıkarırken. Bu sıfır topografya gerektirir
kütle anomalisi olmaması anlamına gelir.

g = Jeoid anomalileri (m).

v = Dikey Yerçekimi Gradyanı (VGG; 1 Eotvos = 0.1 mGal/km).

e = Dikeyin doğu sapmaları (mikro radyan).

n = Dikeyin kuzey sapmaları (mikro radyan).

-Iw|b|c|t |k
kullanım inrd2 ve inrd1 (topografi/batimetri içeren bir ızgara) tahmin etmek için
kabul|tutarlılık ve stdout'a yaz (-G ayarlanmışsa yoksayılır). Bu ızgara gerekir
aynı bölge için yerçekimi veya jeoid içerir inrd1. Varsayılan hesaplamalar
kabul. Çıktı 3 veya 4 sütun içerir. Frekans (dalga boyu), kabul
(tutarlılık) bir sigma hata çubuğu ve isteğe bağlı olarak teorik bir kabul. Ekle
veri bayrakları (birden üçe kadar) w|b|c|t. w dalga sayısı yerine dalga boyu yazar,
k dalga boyu birimi [m] için km'yi seçer, c kabul yerine tutarlılığı hesaplar, b
"alttan yükleme" teorik kabulü ile dördüncü bir sütun yazar ve t
"elastik plaka" teorik kabulü ile dördüncü bir sütun yazar.

-N[f|q|s|nx/ny][+a|[+d|h|l][+e|n|m][+tgenişlik][+w[sonek]][+z[p]]
FFT için uygun ızgara boyutlarını seçin veya sorgulayın ve isteğe bağlı olarak ayarlayın
parametreler. FFT boyutunu kontrol edin:
-NF FFT'yi verilerin gerçek boyutlarını kullanmaya zorlar.

-Nq daha uygun boyutlar hakkında bilgi alacak, bunları rapor edecek ve devam edecek.

-N'ler isteğe bağlı boyutların bir listesini sunacak ve ardından çıkacaktır.

-Nnx/ny dizi boyutunda FFT yapacak nx/ny (>= ızgara dosya boyutu olmalıdır). Varsayılan
FFT'nin hızını ve doğruluğunu optimize eden boyutları >= verileri seçer. FFT ise
boyutlar > ızgara dosyası boyutları, veriler uzatılır ve sıfıra doğru incelir.

Verilerin trendden arındırılmasını kontrol edin: Doğrusal bir trendi kaldırmak için değiştiriciler ekleyin:
+d: Trend verileri, yani en uygun doğrusal trendi kaldırın [Varsayılan].

+a: Yalnızca ortalama değeri kaldırın.

+h: Yalnızca orta değeri kaldırın, yani 0.5 * (maks + min).

+l: Verileri rahat bırakın.

Verilerin genişletilmesini ve daraltılmasını kontrol edin: Uzantının nasıl kontrol edileceğini kontrol etmek için değiştiricileri kullanın.
ve daraltma yapılacaktır:
+e kenar noktası simetrisi uygulayarak ızgarayı genişletir [Varsayılan],

+m kenar ayna simetrisi uygulayarak ızgarayı genişletir

+n veri uzantısını kapatır.

Daraltma, veri kenarından FFT ızgara kenarına [%100] gerçekleştirilir. Değiştirmek
yoluyla bu yüzde +tgenişlik. Ne zaman +n yürürlüktedir, konik uygulanır
hiçbir uzantı mevcut olmadığından veri marjları yerine [0%].

Geçici sonuçların yazılmasını kontrol edin: Ayrıntılı araştırma için aşağıdakileri yazabilirsiniz.
ileri FFT'ye geçirilen ara ızgara; bu büyük ihtimalle olmuştur
Eğilimi bozulmuş, tüm kenarlar boyunca nokta simetrisi ile uzatılmış ve konik. Ekle
+w[sonek] hangi çıktı dosya adlarının oluşturulacağı (yani, ingrid_prefix.ext)
[konik], nerede ext dosya uzantınızdır. Son olarak, karmaşık ızgarayı kaydedebilirsiniz.
ekleyerek ileri FFT tarafından üretilen +z. Varsayılan olarak gerçeği yazıyoruz ve
hayali bileşenler ingrid_gerçek.ext ve ingrid_resim.ext. Ekle p kaydetmek için
bunun yerine dosyalara büyüklük ve fazın kutupsal biçimi ingrid_mag.ext ve
ingrid_faz.ext.

-Q Ortalaması eğilme topografyasına sahip (z pozitif yukarı ile) bir ızgara yazar.
tarafından kuruldu -Zzm ve model parametreleri tarafından -T (ve çıktı -G). bu
"gravimetrik Moho". -Q örtük olarak ayarlar -N+sa

-S tarafından üretilen bir alt plaka yükü nedeniyle tahmin edilen yerçekimi veya jeoit ızgarasını hesaplar.
mevcut batimetri ve teorik model. Gerekli parametreler ayarlandı
içinde -T ve -Z seçenekler. Parker genişletmesindeki güçlerin sayısı şunlarla sınırlıdır:
1. Aşağıdaki bir örneğe bakın.

-Tte/rl/rm/rw[+m]
Topografi yükünden (giriş ızgara dosyası) izostatik telafiyi hesaplayın.
elastik levha kalınlık te. Ayrıca yük, manto ve su için yoğunlukları ekleyin.
SI birimleri. üzerinden ortalama manto derinliği verin -Z. Elastik kalınlık > 1e10 ise
eğilme rijitliği olarak yorumlanacaktır (varsayılan olarak te ve
Genç modül). İsteğe bağlı olarak, ekle +m Moho'nun jeopotansiyeliyle bir ızgara yazmak için
etkisi (bkz. -F) tarafından seçilen modelden -T. Eğer te = 0 ise, Airy yanıtı
iade. -T+m örtük olarak ayarlar -N+sa

-Wwd Topografyaya göre su derinliğini (veya gözlem yüksekliğini) ayarlayın [0]. Ekle k için
km'yi gösterir.

-Zzm[zl]
Moho [ve kabarma] ortalama telafi derinlikleri. "Üstten yükle" modeli için
sadece sağlamak zorunda zm, ancak "aşağıdan yükleme" için unutmayın zl.

-V[seviye] (Daha ...)
Ayrıntı düzeyini seçin [c].

-fg Coğrafi ızgaralar (boylam, enlem boyutları) metreye dönüştürülecek
mevcut elipsoid parametreleri kullanılarak bir "Düz Dünya" yaklaşımı yoluyla.

-^ or sadece -
Komutun sözdizimi hakkında kısa bir mesaj yazdırın, ardından çıkar (NOT: Windows'ta
sadece kullan -).

-+ or sadece +
Açıklama da dahil olmak üzere kapsamlı bir kullanım (yardım) mesajı yazdırın.
modüle özgü seçenek (ancak GMT ortak seçenekleri değil), ardından çıkar.

-? or yok hayır argümanlar
Seçeneklerin açıklamasını da içeren eksiksiz bir kullanım (yardım) mesajı yazdırın, ardından
çıkışlar.

--versiyon
GMT sürümünü yazdırın ve çıkın.

--gösteri-veri dizini
GMT paylaşım dizininin tam yolunu yazdırın ve çıkın.

GRID DOSYA BİÇİMLER

Varsayılan olarak GMT, ızgarayı COARDS şikayeti netCDF'de tek duyarlıklı yüzer olarak yazar
dosya formatı. Bununla birlikte, GMT, yaygın olarak kullanılan diğer birçok ızgarada ızgara dosyaları üretebilir.
dosya biçimlerini destekler ve ayrıca ızgaraların "paketlenmesini" kolaylaştırır, kayan noktayı yazar
1 veya 2 baytlık tamsayılar olarak veri. Kesinlik, ölçek ve ofseti belirtmek için kullanıcı şunları yapmalıdır:
son eki ekle =id[/ölçek/ofset[/nan]], nerede id ızgaranın iki harfli tanımlayıcısıdır
tip ve hassasiyet ve ölçek ve ofset isteğe bağlı ölçek faktörü ve ofset
tüm ızgara değerlerine uygulanır ve nan eksik verileri belirtmek için kullanılan değerdir. Durumunda
iki karakter id olarak sağlanmamaktadır. =/ölçek a'dan daha id=nf varsayılır. Ne zaman
ızgaraları okurken, format genellikle otomatik olarak tanınır. Değilse, aynı son ek
ızgara dosya adlarını girmek için eklenebilir. Görmek grddönüştürme ve Bölüm ızgara-dosya formatı
Daha fazla bilgi için GMT Teknik Referans ve Yemek Kitabı.

Birden çok ızgara içeren bir netCDF dosyasını okurken, GMT varsayılan olarak şunları okuyacaktır:
o dosyada bulabilen ilk 2 boyutlu ızgara. GMT'yi başka bir okumaya ikna etmek için
ızgara dosyasındaki çok boyutlu değişken, ekleme ?Değişkenadı dosya adına, nerede
Değişkenadı değişkenin adıdır. Özel anlamdan kaçmanız gerekebileceğini unutmayın.
of ? Shell programınızda önüne ters eğik çizgi koyarak veya
tırnaklar veya çift tırnaklar arasında dosya adı ve sonek. NS ?Değişkenadı eki de kullanılabilir
çıkış ızgaralarının varsayılandan farklı bir değişken adı belirtmesi için: "z". Görmek
grddönüştürme ve GMT Teknik bölümünün CF için değiştiricileri ve ızgara dosya formatı
Daha fazla bilgi için Referans ve Yemek Kitabı, özellikle 3-,
4 veya 5 boyutlu ızgaralar.

GRID MESAFE BİRİMLERİ

Izgara yatay birim olarak metreye sahip değilse, +ubirim giriş dosyasına
Belirtilen birimden metreye dönüştürülecek isim. Izgaranız coğrafiyse, dönüştürün
sağlayarak metrelere olan mesafeleri -fg yerine.

İLGİLİ HUSUSLAR

netCDF COARDS ızgaraları otomatik olarak coğrafi olarak tanınacaktır. Diğer ızgaralar için
dereceleri metreye dönüştürmek istediğiniz coğrafi ızgaraları seçin, -fg. eğer veri
her iki kutba da yakınsanız, ızgara dosyasını bir dikdörtgen üzerine yansıtmayı düşünmelisiniz.
kullanarak koordinat sistemi grd projesi.

PLAKA ESNEK

Elastik plaka bükülmesine yönelik FFT çözümü, dolgu yoğunluğunun yüke eşit olmasını gerektirir
yoğunluk. Bu tipik olarak yalnızca doğrudan yükün altında geçerlidir; yükün ötesinde
dolgu, düşük yoğunluklu çökeltiler veya hatta su (veya hava) olma eğilimindedir. Wessel [2001] önerdi
farklı bir dolgu yoğunluğunun belirtilmesine izin veren bir yaklaşımdır.
hala bir FFT çözümüne izin verirken yük yoğunluğu. Temel olarak, plaka bükülmesi
Etkili yük yoğunluğu olarak dolgu yoğunluğunu kullanmak için çözüldü, ancak genlikler
bir faktör tarafından ayarlanmış A = sqrt ((rm - ri)/(rm - rl)) teorik farktır
iki farklı yük yoğunluğunu kullanan bir nokta yükü nedeniyle genlikte. NS
yaklaşıklık çok iyidir ancak zayıf plakalardaki büyük yükler için bozulur, bir peri
alışılmadık durum.

ÖRNEKLER

2700 ve 1035 kullanarak bat.grd batimetrisinin üzerindeki su tabakasının etkisini hesaplamak
kabuk ve su yoğunlukları için ve sonucu water_g.grd'ye yazmak (hesaplama
Parker genişlemesindeki dördüncü batimetri gücüne):

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

Şimdi bunu serbest hava anomaliniz faa.grd'den çıkarın ve Bouguer anomalisini elde edeceksiniz. Sen
neden çıkardığımızı ve eklemediğimizi merak edebilir. Bouguer anomalisi gibi davrandıktan sonra
su tabakası tarafından sunulan kütle eksikliğini düzeltmek için, bu yüzden eklemeliyiz çünkü
su aşağıdaki kayalardan daha az yoğundur. Cevap, yerçekimi etkilerinin nasıl olduğuna bağlıdır.
Parker'ın yöntemi ve FFT'yi kullanmanın pratik yönleriyle hesaplanır.

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

MBA anomalisi ister misiniz? Kabuk manto katkısını iyi hesaplayın ve
deniz dibi anomalisi. Yoğunluğu 6 olan 2700 km kalınlığında bir kabuk ve 3300 olan bir manto varsayıldığında
yoğunluk, 600 kullanarak su tabakası anomalisini hesaplamak için kullanılan komutu tekrarlayabiliriz.
(3300 - 2700) yoğunluk kontrastı olarak. Ama şimdi bir sorunumuz var çünkü bilmemiz gerekiyor
ortalama Moho derinliği. Bu, ızgara adına eklenebilecek ölçek/ofset
eline gelir. Bunu daha önce yapmamıza gerek olmadığına dikkat edin, çünkü ortalama su derinliği
doğrudan verilerden hesaplanır (ayrıca,
o z yukarı pozitif):

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

Şimdi, MBA anomalisini elde etmek için onu deniz dibi anomalisinden çıkarın. Yani:

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

Te = 7 km ile elastik bir levha bat.grd'nin Moho yerçekimi etkisini hesaplamak için, yoğunluğu
2700, yoğunluğu 3300 olan bir manto üzerinde, ortalama 9 km derinlikte

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

Şimdi deniz dibini ve Moho'nun etkilerini eklerseniz, tam yerçekimi tepkisini alırsınız.
izostatik modelinizin Burada sadece Parker genişlemesindeki ilk terimi kullanacağız.

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 EKLE = model.grd

Aynı sonuç, bir sonraki komutla doğrudan elde edilebilir. Ancak, DİKKAT
devamındaki. Bir hatadan mı yoksa bir sınırlamadan mı kaynaklandığını henüz bilmiyorum, ama
gerçek şu ki, aşağıdaki ve önceki komutlar yalnızca şu durumlarda aynı sonucu verir: -E1
kullanıldı. Parker genişletmesinde daha yüksek batimetri güçleri için yalnızca yukarıdaki örnek
doğru sonucu vermek için dikişler.

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

Ve a bölgesinin altında, 50 km derinlikte bir yükün ürettiği jeoid anomalisi ne olurdu?
bat.grd tarafından verilen batimetri, 9 km derinlikte bir Moho ve aynı yoğunlukta
önce?

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

topo.grd batimetrisi ile faa.grd serbest hava anomalisi arasındaki kabulü hesaplamak için
6 km etkili 10 km ortalama kalınlığa sahip bir kabuğun elastik plaka modelini kullanan ızgara
3 km'lik bir bölgedeki elastik kalınlık, ortalama su derinliği:

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

topo.grd batimetrisi ve geoid.grd geoid ızgarası arasındaki kabulü hesaplamak için
40 km'de üst ve yüzey altı yükün aynısı olan "alttan yükleme" (LFB) modeli,
ama şimdi ızgaraların coğrafi olduğunu ve frekans yerine dalga boylarını istediğimizi varsayarsak:

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

kullanarak 2000 km uzunluğunda bir profil boyunca bir LFB'nin yerçekimi teorik kabulünü hesaplamak için
yukarıdakiyle aynı parametreler

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

REFERANSLAR

Luis, JF ve MC Neves. 2006, Azor Platosu'nun izostatik telafisi: bir 3D
Kabul ve tutarlılık analizi. J. Jeotermal Volc. Araş. Cilt 156, Sayılar 1-2, Sayfalar
10-22, http://dx.doi.org/10.1016/j.jvolgeores.2006.03.010 Parker, RL, 1972, Hızlı
potansiyel anomalilerin hesaplanması, Geophys. J., 31, 447-455. Wessel. P., 2001, Küresel
ızgaralı Geosat/ERS-1 altimetrisinden çıkarılan deniz miktarlarının dağılımı, J. Geophys. Res.,
106(B9), 19,431-19,441, http://dx.doi.org/10.1029/2000JB000083

onworks.net hizmetlerini kullanarak gravfftgmt'yi çevrimiçi kullanın