gravfftgmt - آنلاین در ابر

برنامه:

نام

gravfft - محاسبه جاذبه گرانشی سطوح سه بعدی در عدد موج (یا
دامنه فرکانس

خلاصه

گرافت اینگرید [ 2 ] مجموعه [ n/طول موج/mean_depth/tbw ] [ چگالی|روگرید ] [
n_terms ] [ [f[+]|g|v|n|e] ] [ w|b|c|t |k ] [ [f|q|s|nx/ny][+a|d|h
|l][+e|n|m][+tعرض][+w[پسوند]][+z[p]] [ ] te/rl/rm/rw[+m] ] [سطح] ] [ wd] [
zm[zl] ] [ -fg ]

توجه داشته باشید: هیچ فاصله ای بین پرچم گزینه و آرگومان های مرتبط مجاز نیست.

شرح

گرافت می توان در سه حالت اصلی استفاده کرد. حالت 1: به سادگی پتانسیل جغرافیایی ناشی از آن را محاسبه کنید
سطح داده شده در فایل topo.grd. به کنتراست چگالی نیاز دارد (-D) و احتمالاً الف
سطوح مختلف مشاهده (-W). FFT 2 بعدی شبکه را به جلو می برد و استفاده می کند
روش پارکر کامل تا شرایط انتخاب شده. حالت 2: ژئوپتانسیل را محاسبه کنید
پاسخ به دلیل خم شدن فایل توپوگرافی این FFT 2 بعدی رو به جلو را می گیرد
شبکه و از روش کامل پارکر اعمال شده برای مدل ایزواستاتیک انتخاب شده استفاده کنید. در
مدل‌های موجود عبارتند از «بارگیری از بالا» یا مدل صفحه الاستیک و «بارگیری از
زیر" که پاسخ صفحه به یک بار زیرسطحی (مناسب برای گرما) را نشان می دهد
مدل سازی نقطه ای - اگر آنها را باور دارید). در هر دو مورد، پارامترهای مدل با تنظیم می شوند -T
و -Z گزینه ها. حالت 3: پذیرش یا انسجام بین دو شبکه را محاسبه کنید. خروجی
میانگین در جهت شعاعی است. به صورت اختیاری، پذیرش مدل نیز ممکن است باشد
محاسبه شد. ابعاد افقی grdfiles بر حسب متر در نظر گرفته شده است.
شبکه های جغرافیایی ممکن است با مشخص کردن استفاده شود -fg گزینه ای که مقیاس درجه به متر است.
اگر شبکه‌هایی با ابعاد کیلومتر دارید، می‌توانید با استفاده از آن به متر تغییر دهید گریدیت or
مقیاس خروجی با گرم ریاضی. با توجه به تعداد گزینه هایی که این برنامه ارائه می دهد، می باشد
بیان گزینه ها و استدلال های مورد نیاز دشوار است. بستگی به چی داره
داری انجام میدی؛ برای راهنمایی بیشتر به مثال ها مراجعه کنید.

مورد نیاز ادله

اینگرید فایل شبکه دو بعدی دوبعدی که باید روی آن کار کرد. (فرمت های فایل GRID را در زیر ببینید). برای
عملیات متقاطع، همچنین فایل شبکه دوم را می دهد ingrd2.

-Gمجموعه
نام فایل شبکه خروجی یا جدول طیف 1 بعدی را مشخص کنید (نگاه کنید به -E) (دیدن
فرمت های فایل GRID در زیر).

اختیاری ادله

-Cn/طول موج/mean_depth/tbw
فقط منحنی های پذیرش نظری مدل انتخاب شده را محاسبه کنید و از آن خارج شوید. n
و طول موج برای محاسبه (n * طول موج) طول کل پروفیل در استفاده می شود
متر است. میانگین_عمق میانگین عمق آب است. ضمیمه dataflag (یک یا دو) از tbw in
هر سفارش t = از مدل "از بالا" استفاده کنید، b = از مدل "از پایین" استفاده کنید. به صورت اختیاری مشخص کنید
w برای نوشتن طول موج به جای فرکانس.

-Dچگالی|روگرید
کنتراست چگالی را در سراسر سطح تنظیم می کند. به عنوان مثال، برای محاسبه گرانش استفاده می شود
جذب لایه آب که بعداً می تواند با ناهنجاری هوای آزاد ترکیب شود
برای دریافت ناهنجاری Bouguer. در این مورد استفاده نکنید -T. همچنین به طور ضمنی تنظیم می کند
-N+h. از طرف دیگر، اگر a
کنتراست چگالی متغیر مورد نیاز است.

-En_terms
تعداد اصطلاحات استفاده شده در بسط پارکر (حد 10 است، در غیر این صورت اصطلاحات بسته به
n برنامه را منفجر می کند) [پیش فرض = 3]

-F[f[+]|g|v|n|e]
میدان ژئوپتانسیل مورد نظر را مشخص کنید: ژئوئید را به جای گرانش محاسبه کنید
f = ناهنجاری های هوای آزاد (mGal) [پیش فرض]. ضمیمه + برای اضافه کردن در دال ضمنی
هنگام حذف مقدار میانگین از توپوگرافی. این نیاز به توپوگرافی صفر دارد
به معنای عدم وجود ناهنجاری توده ای است.

g = ناهنجاری های ژئوئیدی (m).

v = گرادیان گرانش عمودی (VGG؛ 1 Eotvos = 0.1 mGal/km).

e = انحرافات شرقی قائم (میکرو رادیان).

n = انحرافات شمالی قائم (میکرو رادیان).

-Iw|b|c|t |k
استفاده کنید ingrd2 و ingrd1 (شبکه ای با توپوگرافی / عمق سنجی) برای برآورد
admittance|coherence و آن را در stdout بنویسید (-G در صورت تنظیم نادیده گرفته می شود). این شبکه باید
شامل گرانش یا ژئوئید برای همان منطقه از ingrd1. پیش فرض محاسبه می کند
پذیرش خروجی شامل 3 یا 4 ستون است. فرکانس (طول موج)، پذیرش
(انسجام) یک نوار خطای سیگما و به صورت اختیاری، پذیرش نظری. ضمیمه
پرچم های داده (یک تا سه) از w|b|c|t. w به جای عدد موج، طول موج می نویسد،
k کیلومتر را برای واحد طول موج [m] انتخاب می کند، c انسجام را به جای پذیرش محاسبه می کند، b
ستون چهارم را با پذیرش نظری "بارگیری از پایین" می نویسد و t
ستون چهارم را با پذیرش نظری "صفحه الاستیک" می نویسد.

-N[f|q|s|nx/ny][+a|[+d|h|l][+e|n|m][+tعرض][+w[پسوند]][+z[p]]
ابعاد شبکه مناسب برای FFT را انتخاب کنید یا پرس و جو کنید و اختیاری تنظیم کنید
مولفه های. بعد FFT را کنترل کنید:
-Nf FFT را مجبور می کند که از ابعاد واقعی داده ها استفاده کند.

-نق ابعاد مناسب تر را جویا می شود، آن ها را گزارش می دهد، سپس ادامه می دهد.

-Ns لیستی از ابعاد اختیاری را ارائه می دهد، سپس خارج می شود.

-Nnx/ny FFT را در اندازه آرایه انجام خواهد داد nx/ny (باید >= اندازه فایل شبکه ای باشد). پیش فرض
ابعاد >= داده را انتخاب می کند که سرعت و دقت FFT را بهینه می کند. اگر FFT
dimensions > ابعاد فایل شبکه، داده ها گسترش یافته و به صفر کاهش می یابد.

کاهش روند داده ها را کنترل کنید: اصلاح کننده هایی را برای حذف روند خطی اضافه کنید:
+d: داده‌های Detrend، یعنی حذف روند خطی مناسب [پیش‌فرض].

+a: فقط مقدار میانگین را حذف کنید.

+h: فقط مقدار متوسط ​​یعنی 0.5 * (حداکثر + دقیقه) را حذف کنید.

+l: داده ها را به حال خود رها کنید.

کنترل گسترش و کاهش داده ها: از اصلاح کننده ها برای کنترل نحوه برنامه افزودنی استفاده کنید
و Tapering باید انجام شود:
+e شبکه را با اعمال تقارن نقطه لبه [پیش‌فرض] گسترش می‌دهد،

+m شبکه را با تحمیل تقارن آینه لبه گسترش می دهد

+n پسوند داده را خاموش می کند.

Tapering از لبه داده به لبه شبکه FFT [100٪] انجام می شود. تغییر دادن
این درصد از طریق +tعرضاست. وقتی که +n در واقع، مخروطی اعمال می شود
به جای آن به حاشیه داده ها، زیرا هیچ برنامه افزودنی در دسترس نیست [0٪].

نوشتن کنترل نتایج موقت: برای بررسی دقیق می توانید بنویسید
شبکه میانی به FFT رو به جلو منتقل می شود. این احتمالا بوده است
بدون روند، با تقارن نقطه ای در تمام لبه ها گسترش یافته و مخروطی شده است. ضمیمه
+w[پسوند] که نام(های) فایل خروجی از آن ایجاد می شود (یعنی، ingrid_prefix.ext)
[مخروطی]، که در آن EXT پسوند فایل شماست در نهایت، ممکن است شبکه پیچیده را ذخیره کنید
تولید شده توسط FFT فوروارد با الحاق +z. به طور پیش فرض واقعی و را می نویسیم
اجزای خیالی به اینگرید_واقعیEXT و اینگرید_تصویر.EXT. ضمیمه p برای صرفه جویی در
در عوض شکل قطبی قدر و فاز به فایل ها اینگرید_mag.EXT و
اینگرید_فاز.EXT.

-Q یک شبکه با توپوگرافی خمشی (با z مثبت بالا) می نویسد که میانگین آن
تنظیم شده توسط -Zzm و پارامترهای مدل توسط -T (و خروجی توسط -G). آن است
"موهو وزن سنجی". -Q به طور ضمنی تنظیم می کند -N+h

-S گرانش یا شبکه ژئوئیدی پیش بینی شده را به دلیل بار زیر صفحه ای که توسط آن تولید می شود را محاسبه می کند
آب سنجی فعلی و مدل نظری پارامترهای لازم تنظیم شده است
در داخل -T و -Z گزینه ها. تعداد قدرت ها در بسط پارکر محدود به
1. نمونه ای را در ادامه مطلب ببینید.

-Tte/rl/rm/rw[+m]
جبران ایزواستاتیک را از بار توپوگرافی (فایل شبکه ورودی) روی یک محاسبه کنید
صفحه الاستیک با ضخامت te. همچنین چگالی برای بار، گوشته و آب در داخل اضافه کنید
واحدهای SI عمق متوسط ​​گوشته را از طریق -Z. اگر ضخامت الاستیک > 1e10 باشد
به عنوان صلبیت خمشی تفسیر می شود (به طور پیش فرض از محاسبه می شود te و
مدول جوان). به صورت اختیاری، ضمیمه کنید +m برای نوشتن یک شبکه با پتانسیل جغرافیایی Moho
اثر (نگاه کنید به -F) از مدل انتخاب شده توسط -T. اگر te = 0 سپس پاسخ Airy است
برگشت -T+m به طور ضمنی تنظیم می کند -N+h

-Wwd عمق آب (یا ارتفاع مشاهده) را نسبت به توپوگرافی [0] تنظیم کنید. ضمیمه k به
کیلومتر را نشان می دهد.

-Zzm[zl]
Moho [و متورم] عمق جبران متوسط. برای مدل "بار از بالا" شما
فقط باید ارائه کند zm، اما برای "بارگیری از پایین" فراموش نکنید zl.

-V[سطح] (بیشتر ...)
سطح پرحرفی [c] را انتخاب کنید.

-fg شبکه های جغرافیایی (ابعاد طول و عرض جغرافیایی) به متر تبدیل خواهند شد
از طریق یک تقریب "زمین مسطح" با استفاده از پارامترهای بیضی فعلی.

-^ or تنها -
یک پیام کوتاه در مورد نحو دستور چاپ کنید، سپس از آن خارج شوید (توجه: در ویندوز
فقط استفاده کنید -).

-+ or تنها +
یک پیام کاربردی (راهنما) گسترده، از جمله توضیح هر یک را چاپ کنید
گزینه مخصوص ماژول (اما نه گزینه های رایج GMT)، سپس خارج می شود.

-? or نه استدلال
سپس یک پیام استفاده کامل (راهنما)، از جمله توضیح گزینه ها را چاپ کنید
خارج می شود.

- نسخه
نسخه GMT را چاپ کنید و خارج شوید.

--show-datadir
مسیر کامل را به دایرکتوری اشتراک گذاری GMT ​​چاپ کنید و از آن خارج شوید.

GRID فایل فرمها

به‌طور پیش‌فرض GMT شبکه را به‌عنوان شناور دقیق تکی در یک netCDF شکایت COARDS می‌نویسد.
فرمت فایل. با این حال، GMT می‌تواند فایل‌های شبکه‌ای را در بسیاری از شبکه‌های رایج دیگر تولید کند
فرمت های فایل و همچنین به اصطلاح "بسته بندی" شبکه ها، نوشتن نقطه شناور را تسهیل می کند
داده ها به صورت اعداد صحیح 1 یا 2 بایتی. برای تعیین دقت، مقیاس و افست، کاربر باید
پسوند را اضافه کنید =id[/مقیاس/چاپ افست[/نان]]، جایی که id یک شناسه دو حرفی از شبکه است
نوع و دقت و مقیاس و چاپ افست ضریب مقیاس اختیاری و افست هستند
به همه مقادیر شبکه اعمال می شود، و نان مقداری است که برای نشان دادن داده های از دست رفته استفاده می شود. در صورت
دو شخصیت id ارائه نشده است، همانطور که در =/مقیاس از یک id=nf مفروض است. چه زمانی
شبکه های خواندن، قالب به طور کلی به طور خودکار تشخیص داده می شود. اگر نه، همان پسوند
را می توان به نام فایل های شبکه ورودی اضافه کرد. دیدن grdconvert و بخش grid-file-format of the
مرجع فنی GMT ​​و کتاب آشپزی برای اطلاعات بیشتر.

هنگام خواندن یک فایل netCDF که حاوی چندین شبکه است، GMT به طور پیش‌فرض، آن را می‌خواند
اولین شبکه 2 بعدی که می توانید در آن فایل پیدا کنید. برای ترغیب GMT به خواندن دیگری
متغیر چند بعدی در فایل شبکه، پیوست ?نام متغیر به نام فایل، جایی که
نام متغیر نام متغیر است. توجه داشته باشید که ممکن است لازم باشد از معنای خاص فرار کنید
of ? در برنامه پوسته خود با قرار دادن یک بک اسلش در مقابل آن یا با قرار دادن علامت
نام فایل و پسوند بین نقل قول یا دو نقل قول. در ?نام متغیر پسوند نیز می تواند استفاده شود
برای شبکه های خروجی برای تعیین نام متغیر متفاوت از پیش فرض: "z". دیدن
grdconvert و بخش‌های اصلاح‌کننده-برای-CF و قالب-فایل شبکه ای GMT ​​فنی
مرجع و کتاب آشپزی برای اطلاعات بیشتر، به ویژه در مورد نحوه خواندن پیوندهای 3-،
شبکه های 4 یا 5 بعدی.

GRID فاصله یاتاقان ها

اگر شبکه دارای متر به عنوان واحد افقی نیست، ضمیمه کنید +uواحد به فایل ورودی
نام برای تبدیل از واحد مشخص شده به متر. اگر شبکه شما جغرافیایی است، تبدیل کنید
فواصل متر با تامین -fg به جای آن.

بحث در مورد

شبکه های netCDF COARDS به طور خودکار به عنوان جغرافیایی شناسایی می شوند. برای شبکه های دیگر
شبکه های جغرافیایی که می خواهید درجه ها را به متر تبدیل کنید، انتخاب کنید -fg. اگر داده ها
نزدیک به هر یک از قطب ها هستند، باید در نظر بگیرید که فایل شبکه را روی یک مستطیل شکل دهید
سیستم مختصات با استفاده از پروژه grd.

صفحه FLEXURE

راه حل FFT برای خمش صفحه الاستیک نیاز به چگالی پر کردن برای برابر بار دارد
تراکم این معمولاً فقط مستقیماً در زیر بار صادق است. فراتر از بار
پر کردن رسوبات با چگالی کمتر یا حتی آب (یا هوا) است. وسل [2001] پیشنهاد کرد
تقریبی که اجازه می دهد تا مشخصات چگالی پرکننده متفاوت از
چگالی بار در حالی که هنوز برای یک راه حل FFT اجازه می دهد. اساساً خمش صفحه است
برای استفاده از چگالی پر شدن به عنوان چگالی بار موثر حل شده است، اما دامنه ها هستند
با یک عامل تنظیم می شود A = sqrt ((rm - ri)/(rm - rl))، که تفاوت نظری است
در دامنه ناشی از بار نقطه ای با استفاده از دو چگالی بار مختلف. در
تقریب بسیار خوب است، اما برای بارهای بزرگ در صفحات ضعیف، یک پری خراب می شود
وضعیت غیر معمول

مثال ها

برای محاسبه اثر لایه آب بالای عمق سنجی bat.grd با استفاده از 2700 و 1035
برای چگالی پوسته و آب و نوشتن نتیجه در water_g.grd (محاسبه کردن
به چهارمین قدرت آب سنجی در بسط پارکر):

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

اکنون آن را به ناهنجاری هوای آزاد خود کم کنید faa.grd و ناهنجاری Bouguer را دریافت خواهید کرد. شما
ممکن است تعجب کنیم که چرا ما کم می کنیم و اضافه نمی کنیم. پس از همه ناهنجاری بوگر وانمود می کند
برای تصحیح کمبود جرم ارائه شده توسط لایه آب، پس باید اضافه کنیم زیرا
چگالی آب کمتر از سنگ های زیر است. پاسخ به نحوه تأثیرات جاذبه بستگی دارد
با روش پارکر و جنبه های عملی استفاده از FFT محاسبه شده است.

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

آیا یک ناهنجاری MBA می خواهید؟ خوب سهم گوشته پوسته را محاسبه کرده و به آن اضافه کنید
ناهنجاری کف دریا با فرض یک پوسته با ضخامت 6 کیلومتر با چگالی 2700 و یک گوشته با 3300
چگالی می‌توانیم با استفاده از 600، دستور مورد استفاده برای محاسبه ناهنجاری لایه آب را تکرار کنیم
(3300 - 2700) به عنوان کنتراست چگالی. اما ما اکنون مشکل داریم زیرا باید بدانیم
میانگین عمق موهو در آن زمان است که مقیاس / افست قابل اضافه شدن به نام شبکه است
در دست می آید توجه داشته باشید که قبلاً نیازی به انجام این کار نداشتیم زیرا عمق آب متوسط ​​بود
به طور مستقیم از داده ها محاسبه می شود (به علامت منفی افست نیز توجه کنید
که z مثبت است):

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

حال، آن را به ناهنجاری کف دریا تفریق کنید تا ناهنجاری MBA به دست آید. به این معنا که:

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

برای محاسبه اثر گرانش موهو یک صفحه الاستیک bat.grd با Te = 7 کیلومتر، چگالی
2700، بیش از یک گوشته با چگالی 3300، در عمق متوسط ​​9 کیلومتر

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

اگر اکنون افکت های sea-bottom و Moho را اضافه کنید، پاسخ جاذبه کامل را دریافت خواهید کرد
مدل ایزواستاتیک شما ما در اینجا فقط از اولین عبارت در بسط پارکر استفاده خواهیم کرد.

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

همین نتیجه را می توان مستقیماً با دستور بعدی به دست آورد. با این حال، توجه به
به شرح زیر. من هنوز نمی دانم که آیا این به دلیل یک اشکال است یا به دلیل برخی محدودیت ها، اما
واقعیت این است که دستورات زیر و قبلی فقط نتیجه یکسانی را می دهند اگر -E1
استفاده می شود. برای قدرت های بالاتر عمق سنجی در بسط پارکر، فقط مثال بالا
درزها برای دادن نتیجه صحیح

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

و ناهنجاری ژئوئیدی که توسط بار در عمق 50 کیلومتری، زیر ناحیه a ایجاد می شود، چه خواهد بود؟
که عمق سنجی آن توسط bat.grd، یک موهو در عمق 9 کیلومتری و همان چگالی ارائه شده است.
قبل از؟

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

برای محاسبه میزان پذیرش بین عمق سنجی topo.grd و ناهنجاری هوای آزاد faa.grd
شبکه با استفاده از مدل صفحه الاستیک از پوسته 6 کیلومتر ضخامت متوسط ​​با 10 کیلومتر موثر
ضخامت الاستیک در منطقه 3 کیلومتری به طور میانگین عمق آب:

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

برای محاسبه میزان پذیرش بین عمق سنجی topo.grd و geoid.grd geoid grid با
مدل "بارگیری از پایین" (LFB) با بار مشابه بالا و بار زیرسطحی در 40 کیلومتر،
اما با فرض اینکه اکنون شبکه ها در جغرافیا هستند و ما به جای فرکانس، طول موج می خواهیم:

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

برای محاسبه پذیرش نظری گرانش یک LFB در طول یک پروفیل به طول 2000 کیلومتر با استفاده از
همان پارامترهای بالا

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

مراجع

لوئیس، جی اف و ام سی نئوس. 2006، جبران ایزواستاتیک فلات آزور: یک 3 بعدی
پذیرش و تحلیل انسجام J. Geothermal Volc. Res. جلد 156، شماره 1-2، صفحات
10 22، http://dx.doi.org/10.1016/j.jvolgeores.2006.03.010 Parker, RL, 1972, The Rapid
محاسبه ناهنجاری های بالقوه، ژئوفیز. ج.، 31، 447-455. وسل پ.، 2001، جهانی
توزیع seamounts استنباط از شبکه بندی Geosat/ERS-1 ارتفاع سنجی، J. Geophys. پاسخ،
106 (B9)، 19,431،19,441-XNUMX،XNUMX، http://dx.doi.org/10.1029/2000JB000083

با استفاده از خدمات onworks.net از gravfftgmt به صورت آنلاین استفاده کنید