Gravfftgmt - عبر الإنترنت في السحابة

برنامج:

اسم

Gravfft - حساب جاذبية الأسطح ثلاثية الأبعاد في العدد الموجي (أو
التردد) المجال

موجز

com.gravfft إنغريد [ ingrid2 ] ملف [ ن/الطول الموجي/متوسط_العمق/تبو ] [ كثافة|rhogrid ] [
n_terms ] [[f[+]|g|v|n|e]] [ w|b|c|t |k ] [[f|q|s|نكس / نيويورك][+a|d|h
|l][+e|n|m][+tعرض][+w[لاحقة]] [+z[p]] [] [ ر.ه./rl/rm/rw[+m]] [[مستوى]] [ wd] [
zm[zl]] [ -fg ]

ملحوظة: لا يسمح بمسافة بين علامة الخيار والوسيطات المرتبطة بها.

الوصف

com.gravfft يمكن استخدامها في ثلاثة أوضاع رئيسية. الوضع 1: ما عليك سوى حساب الإمكانات الجغرافية الناتجة عن ذلك
السطح الوارد في ملف topo.grd. يتطلب تباين الكثافة (-D) وربما أ
مستوى ملاحظة مختلف (-W). وسوف يستغرق 2-D إلى الأمام FFT من الشبكة واستخدامها
طريقة باركر الكاملة حتى الشروط المختارة. الوضع 2: حساب الإمكانات الجغرافية
الاستجابة بسبب انثناء الملف الطبوغرافي. سوف يستغرق الأمر 2-D للأمام FFT لـ
الشبكة واستخدم طريقة باركر الكاملة المطبقة على النموذج المتوازن المختار. ال
النماذج المتوفرة هي نموذج "التحميل من الأعلى"، أو نموذج اللوحة المرنة، ونموذج "التحميل من الأعلى".
أدناه" والذي يمثل استجابة اللوحة لحمل تحت السطح (مناسب للسخونة
النمذجة الفورية - إذا كنت تصدقهم). وفي كلتا الحالتين، يتم تعيين معلمات النموذج باستخدام -T
و -Z خيارات. الوضع 3: حساب القبول أو التماسك بين شبكتين. الإخراج
هو المتوسط ​​في الاتجاه الشعاعي. بشكل اختياري، قد يكون القبول النموذجي أيضًا
محسوب. من المفترض أن تكون الأبعاد الأفقية لملفات grd بالمتر.
يمكن استخدام الشبكات الجغرافية عن طريق تحديد -fg الخيار الذي يقيس الدرجات إلى الأمتار.
إذا كان لديك شبكات بأبعاد بالكيلومترات، فيمكنك تغيير ذلك إلى أمتار باستخدام grdedit or
قم بتوسيع حجم الإخراج مع com.grdmath. ونظرًا لعدد الخيارات التي يقدمها هذا البرنامج، فهو
من الصعب تحديد ما هي الخيارات وما هي الحجج المطلوبة. يعتمد على ماذا
انت تفعل؛ انظر الأمثلة لمزيد من التوجيه.

مطلوبة الحجج

إنغريد ملف شبكي ثنائي الأبعاد ليتم التشغيل عليه. (راجع تنسيقات ملفات GRID أدناه). ل
العمليات عبر الطيفية، كما تعطي ملف الشبكة الثاني ingrd2.

-Gملف
حدد اسم ملف شبكة الإخراج أو جدول الطيف ثنائي الأبعاد (انظر -E). (نرى
تنسيقات ملفات GRID أدناه).

اختياري الحجج

-Cن/الطول الموجي/متوسط_العمق/تبو
حساب فقط منحنيات القبول النظرية للنموذج المحدد والخروج. n
و الطول الموجي تُستخدم لحساب (n * الطول الموجي) الطول الإجمالي للملف الجانبي
متر. متوسط_العمق هو متوسط ​​عمق الماء. إلحاق علامات البيانات (واحدة أو اثنتين) من TBW in
اي طلب. t = استخدم النموذج "من الأعلى"، b = استخدم نموذج "من الأسفل". تحديد اختياريا
w لكتابة الطول الموجي بدلا من التردد.

-Dكثافة|rhogrid
يضبط تباين الكثافة عبر السطح. تستخدم، على سبيل المثال، لحساب الجاذبية
جاذبية طبقة الماء التي يمكن دمجها لاحقًا مع شذوذ الهواء الحر
للحصول على شذوذ Bouguer. في هذه الحالة لا تستخدم -T. كما أنه يحدد ضمنا
-ن+ح. بدلاً من ذلك، حدد شبكة مسجلة بشكل مشترك مع تباينات الكثافة إذا كان أ
مطلوب تباين الكثافة المتغيرة.

-En_terms
عدد المصطلحات المستخدمة في توسيع Parker (الحد الأقصى هو 10، وإلا فإن المصطلحات تعتمد على ذلك
n سوف يفجر البرنامج) [افتراضي = 3]

-F[f[+]|g|v|n|e]
حدد مجال الجهد الأرضي المطلوب: قم بحساب الجيوديوم بدلاً من الجاذبية
f = شذوذات الهواء الحر (mGal) [افتراضي]. ألحق + لإضافة في بلاطة ضمنا
عند إزالة القيمة المتوسطة من التضاريس. وهذا يتطلب صفر التضاريس
ليعني عدم وجود شذوذ جماعي.

g = الشذوذ الجيود (م).

v = تدرج الجاذبية العمودية (VGG؛ 1 Eotvos = 0.1 ملي جالون/كم).

e = الانحرافات الشرقية للعمودي (الراديان الصغير).

n = الانحرافات الشمالية للعمودي (الراديان الصغير).

-إو|ب|ج|ر |k
استعمل ingrd2 و ingrd1 (شبكة مع التضاريس / قياس الأعماق) للتقدير
القبول | التماسك وكتابته إلى stdout (-G تم تجاهله إذا تم ضبطه). ينبغي لهذه الشبكة
تحتوي على الجاذبية أو الجيود لنفس المنطقة ingrd1. الحسابات الافتراضية
القبول. يحتوي الإخراج على 3 أو 4 أعمدة. التردد (الطول الموجي)، القبول
(التماسك) شريط خطأ سيجما واحد، واختياريًا، الاعتراف النظري. ألحق
علامات البيانات (من واحد إلى ثلاثة) من w|b|c|t. w يكتب الطول الموجي بدلا من wavenumber،
k يختار كم لوحدة الطول الموجي [م]، c يحسب التماسك بدلا من القبول، b
يكتب عمودًا رابعًا مع القبول النظري "التحميل من الأسفل"، و t
يكتب العمود الرابع مع القبول النظري "لوحة مرنة".

-N [f | q | s |نكس / نيويورك] [+ أ | [+ د | ح | ل] [+ ه | ن | م] [+ رعرض] [+ ث [لاحقة]] [+ z [p]]
اختر أو استفسر عن أبعاد الشبكة المناسبة لـ FFT وقم بتعيينها اختياريًا
حدود. التحكم في البعد FFT:
-نف سيجبر FFT على استخدام الأبعاد الفعلية للبيانات.

-نق سوف يستفسر عن أبعاد أكثر ملاءمة ، ويبلغ عنها ، ثم يستمر.

-ن سيقدم قائمة بالأبعاد الاختيارية ، ثم الخروج.

-Nنكس / نيويورك سوف تفعل FFT على حجم الصفيف نكس / نيويورك (يجب أن يكون> ​​= حجم ملف الشبكة). تقصير
يختار الأبعاد> = البيانات التي تعمل على تحسين سرعة ودقة FFT. إذا كان FFT
الأبعاد> أبعاد ملف الشبكة ، يتم تمديد البيانات وتدقيقها إلى الصفر.

التحكم في تعطيل البيانات: إلحاق مُعدِّلات لإزالة اتجاه خطي:
+d: Detrend data ، أي إزالة الاتجاه الخطي الأنسب [افتراضي].

+a: قم بإزالة متوسط ​​القيمة فقط.

+h: قم بإزالة القيمة المتوسطة فقط ، أي 0.5 * (حد أقصى + حد أدنى).

+l: اترك البيانات وشأنها.

التحكم في امتداد البيانات واستدقاقها: استخدم المعدِّلات للتحكم في كيفية الامتداد
والاستدقاق يجب أن يتم:
+e يوسع الشبكة بفرض تناظر نقطة الحافة [افتراضي] ،

+m يوسع الشبكة بفرض تناظر مرآة الحواف

+n يوقف تمديد البيانات.

يتم إجراء الاستدقاق من حافة البيانات إلى حافة شبكة FFT [100٪]. يتغير
هذه النسبة عبر +tعرض. عندما +n في الواقع ، يتم تطبيق التناقص
بدلاً من ذلك إلى هوامش البيانات حيث لا يوجد امتداد متاح [0٪].

التحكم في كتابة النتائج المؤقتة: للحصول على تحقيق مفصل ، يمكنك كتابة
تمرير الشبكة الوسيطة إلى FFT الأمامي ؛ من المحتمل أن يكون هذا
منحرف ، ممتد بالتناظر النقطي على طول جميع الحواف ، ومدبب. ألحق
+w[لاحقة] الذي سيتم إنشاء اسم (أسماء) ملف الإخراج منه (على سبيل المثال ، ingrid_prefix.ext)
[مدبب] ، أين تحويلة هو امتداد الملف الخاص بك. أخيرًا ، يمكنك حفظ الشبكة المعقدة
التي تنتجها FFT إلى الأمام عن طريق إلحاق +z. بشكل افتراضي نكتب الحقيقي و
مكونات وهمية إنغريد_حقيقي.تحويلة و إنغريد_imag.تحويلة. ألحق p لإنقاذ
بدلا من الشكل القطبي من حيث الحجم والمرحلة إلى الملفات إنغريد_ماج.تحويلة و
إنغريد_مرحلة.تحويلة.

-Q يكتب شبكة ذات تضاريس الانحناء (مع z موجب لأعلى) متوسطها
تم تعيينه بواسطة -Zzm ومعلمات النموذج بواسطة -T (والإخراج بواسطة -G). هذا هو
"قياس الجاذبية موهو". -Q مجموعات ضمنية -ن+ح

-S يحسب الجاذبية المتوقعة أو الشبكة الجغرافية بسبب حمل اللوحة الفرعية الناتج عن
قياس الأعماق الحالي والنموذج النظري. يتم تعيين المعلمات الضرورية
في غضون -T و -Z خيارات. يقتصر عدد القوى في توسعة باركر على
1. انظر مثالاً في الأسفل.

-Tر.ه./rl/rm/rw[+ م]
حساب التعويض المتوازن من تحميل التضاريس (ملف شبكة الإدخال) على
لوحة مرنة من سمك te. قم أيضًا بإلحاق كثافات الحمل والوشاح والماء فيها
وحدات SI. إعطاء متوسط ​​عمق الوشاح عبر -Z. إذا كان سمك المرونة> 1e10
سيتم تفسيرها على أنها صلابة الانثناء (افتراضيًا يتم حسابها من te و
عارض شاب). اختياريا، إلحاق +m لكتابة شبكة مع الإمكانات الجغرافية لموهو
تأثير (انظر -F) من النموذج المحدد بواسطة -T. إذا te = 0 فإن الاستجابة Airy هي
عاد. -ت+م مجموعات ضمنية -ن+ح

-Wwd اضبط عمق المياه (أو ارتفاع المراقبة) بالنسبة للتضاريس [0]. ألحق k إلى
تشير كم.

-Zzm[zl]
موهو [وتنتفخ] متوسط ​​أعماق التعويض. لنموذج "التحميل من الأعلى" لك
يجب فقط أن توفر zmلكن بالنسبة إلى "التحميل من الأسفل" فلا تنسَ zl.

-الخامس[مستوى] (أكثر ...)
حدد مستوى الإسهاب [ج].

-fg سيتم تحويل الشبكات الجغرافية (أبعاد خطوط الطول والعرض) إلى أمتار
عن طريق تقريب "الأرض المسطحة" باستخدام المعلمات الإهليلجية الحالية.

-^ or م -
اطبع رسالة قصيرة حول صيغة الأمر ، ثم قم بالخروج (ملاحظة: في Windows
استخدم فقط -).

-+ or م +
اطبع رسالة (تعليمات) استخدام شاملة ، بما في ذلك شرح أي منها
خيار خاص بالوحدة النمطية (ولكن ليس خيارات GMT الشائعة) ، ثم يخرج.

-? or لا الحجج
اطبع رسالة استخدام (مساعدة) كاملة ، بما في ذلك شرح الخيارات ، بعد ذلك
المخارج.

--الإصدار
اطبع نسخة GMT واخرج.

--إظهار- datadir
طباعة المسار الكامل إلى دليل مشاركة GMT والخروج.

غريد FILE صيغ

بشكل افتراضي يكتب GMT خارج الشبكة عندما يطفو الدقة المفردة في netCDF شكوى COARDS
تنسيق الملف. ومع ذلك ، فإن GMT قادرة على إنتاج ملفات الشبكة في العديد من الشبكات الأخرى الشائعة الاستخدام
وتنسيقات الملفات ويسهل أيضًا ما يسمى بـ "حزم" الشبكات ، وكتابة النقطة العائمة
البيانات كأعداد صحيحة 1 أو 2 بايت. لتحديد الدقة والمقياس والإزاحة ، يجب على المستخدم
أضف اللاحقة =id[/مقياس/عوض[/نان]]، أين id هو معرف من حرفين للشبكة
النوع والدقة ، و مقياس و عوض عامل مقياس اختياري ويقابله
يتم تطبيقه على جميع قيم الشبكة ، و نان هي القيمة المستخدمة للإشارة إلى البيانات المفقودة. في حالة
الشخصيتين id لم يتم توفيرها ، كما في =/مقياس من أ id=nf يفترض. متي
قراءة الشبكات ، يتم التعرف على التنسيق تلقائيًا بشكل عام. إذا لم يكن كذلك ، نفس اللاحقة
يمكن إضافتها إلى أسماء ملفات الشبكة المدخلة. ارى com.grdconvert وقسم شبكة ملف بتنسيق
GMT المرجع الفني وكتاب الطبخ لمزيد من المعلومات.

عند قراءة ملف netCDF يحتوي على شبكات متعددة ، فإن GMT ستقرأ افتراضيًا ملف
أول شبكة ثنائية الأبعاد يمكن العثور عليها في هذا الملف. لإقناع GMT في قراءة أخرى
متغير متعدد الأبعاد في ملف الشبكة ، إلحاق ?فارمين إلى اسم الملف ، حيث
فارمين هو اسم المتغير. لاحظ أنك قد تحتاج إلى الهروب من المعنى الخاص
of ? في برنامج shell الخاص بك عن طريق وضع شرطة مائلة للخلف أمامه ، أو بوضع الرمز
اسم الملف واللاحقة بين علامات الاقتباس أو علامات الاقتباس المزدوجة. ال ?فارمين يمكن أيضًا استخدام اللاحقة
لشبكات الإخراج لتحديد اسم متغير مختلف عن الافتراضي: "z". ارى
com.grdconvert ومعدلات الأقسام لـ CF وتنسيق ملف الشبكة لـ GMT Technical
المرجع وكتاب الطبخ لمزيد من المعلومات ، لا سيما حول كيفية قراءة وصلات 3 ،
شبكات ذات 4 أو خماسية الأبعاد.

غريد مسافه: بعد وحدات

إذا كانت الشبكة لا تحتوي على عداد كوحدة أفقية ، قم بإلحاقه +uوحدة إلى ملف الإدخال
اسم للتحويل من الوحدة المحددة إلى متر. إذا كانت شبكتك جغرافية ، قم بالتحويل
مسافات إلى أمتار عن طريق التوريد -fg بدلا من ذلك.

الاعتبارات

سيتم التعرف على شبكات netCDF COARDS تلقائيًا على أنها جغرافية. للشبكات الأخرى
الشبكات الجغرافية التي تريد تحويل الدرجات إلى أمتار فيها ، حدد -fg. إذا كانت البيانات
بالقرب من أي من القطبين ، يجب أن تفكر في إسقاط ملف الشبكة على مستطيل
نظام الإحداثيات باستخدام grdproject.

طبق فليكسور

يتطلب حل FFT لثني الصفائح المرنة أن تكون كثافة الحشو مساوية للحمل
كثافة. عادةً ما يكون هذا صحيحًا فقط أسفل الحمل مباشرةً؛ وراء الحمل
يميل الحشو إلى أن يكون رواسب منخفضة الكثافة أو حتى ماء (أو هواء). اقترح ويسيل [2001].
تقريب يسمح بمواصفات كثافة ملء مختلفة عن
كثافة الحمل مع استمرار السماح بحل FFT. في الأساس ، انثناء اللوحة هو
تم حلها لاستخدام كثافة الملء باعتبارها كثافة الحمل الفعالة ولكن السعات
تعديلها بعامل A = الجذر التربيعي ((rm - ri) / (rm - rl)) ، وهو الاختلاف النظري
في السعة بسبب تحميل نقطة باستخدام كثافتي الحمل المختلفتين. ال
التقريب جيد جدًا ولكنه ينهار للأحمال الكبيرة على الأطباق الضعيفة ، إنها خرافية
حالة غير شائعة.

أمثلة

لحساب تأثير الطبقة المائية فوق قياس الأعماق الخفافيش باستخدام 2700 و 1035
لكثافة القشرة والماء وكتابة النتيجة على water_g.grd (حساب
إلى القوة الرابعة لقياس الأعماق في توسعة باركر):

gmt جراففت بات.grd -D1665 -Gwater_g.grd -E4

الآن قم بطرحها إلى شذوذ الهواء الحر faa.grd وستحصل على شذوذ Bouguer. أنت
قد يتساءل لماذا نقوم بالطرح وعدم الإضافة. بعد كل شيء، يتظاهر شذوذ Bouguer
لتصحيح النقص الكتلي الذي تمثله طبقة الماء فيجب أن نضيف لأنه
الماء أقل كثافة من الصخور أدناه. تعتمد الإجابة على طريقة تأثير الجاذبية
محسوبة بطريقة باركر والجوانب العملية لاستخدام FFT.

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

تريد شذوذ ماجستير في إدارة الأعمال؟ حسنًا ، احسب مساهمة القشرة الأرضية وأضفها إلى
شذوذ قاع البحر. بافتراض قشرة سمكها 6 كم كثافتها 2700 ووشاحها 3300
يمكننا تكرار الأمر المستخدم لحساب شذوذ طبقة الماء باستخدام 600
(3300 - 2700) كتباين الكثافة. ولكن لدينا الآن مشكلة لأننا بحاجة إلى أن نعرف
متوسط ​​عمق موهو. وذلك عندما يتم إلحاق المقياس/الإزاحة باسم الشبكة
يأتي في متناول اليد. لاحظ أننا لم نكن بحاجة إلى القيام بذلك من قبل لأن متوسط ​​عمق الماء كان كذلك
محسوبة مباشرة من البيانات (لاحظ أيضًا الإشارة السلبية للإزاحة نظرًا لحقيقة
أن z هو إيجابي يصل):

بتوقيت جرينتش غراففت بات.grd=nf/1/-6000 -D600 -Gmoho_g.grd

الآن، قم بطرحها إلى شذوذ قاع البحر للحصول على شذوذ ماجستير إدارة الأعمال. إنه:

بتوقيت جرينتش 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

إذا أضفت الآن تأثيرات قاع البحر وموهو، فستحصل على استجابة الجاذبية الكاملة
من النموذج المتوازن الخاص بك. سنستخدم هنا الفصل الأول فقط في توسعة باركر.

gmt جراففت بات.grd -D1665 -Gwater_g.grd -E1
gmt Gravfft Bat.grd -Gelastic.grd -T7000/2700/3300/1035+m -Z9000 -E1
بتوقيت جرينتش grdmath water_g.grd elastic.grd ADD = model.grd

يمكن الحصول على نفس النتيجة مباشرة عن طريق الأمر التالي. ومع ذلك، انتبه إلى
الأتى. لا أعرف حتى الآن ما إذا كان ذلك بسبب خطأ أو بسبب بعض القيود، ولكن
والحقيقة هي أن الأوامر التالية والسابقة تعطي نفس النتيجة فقط إذا -E1
يستخدم. للحصول على قوى أعلى لقياس الأعماق في توسعة باركر، فقط المثال أعلاه
طبقات لإعطاء النتيجة الصحيحة.

gmt جراففت بات.grd -Gmodel.grd -T7000/2700/3300/1035 -Z9000 -E1

وما هو الشذوذ الجيودي الناتج عن حمل على عمق 50 كم، أسفل المنطقة
يتم إعطاء قياس الأعماق بواسطة Bat.grd، وهو Moho على عمق 9 كم وبنفس الكثافة
قبل؟

بتوقيت جرينتش جراففت topo.grd -Gswell_geoid.grd -T7000/2700/3300/1035 -Fg -Z9000/50000 -S -E1

لحساب القبول بين قياس الأعماق topo.grd وشذوذ الهواء الحر faa.grd
شبكة باستخدام نموذج اللوحة المرنة لقشرة يبلغ متوسط ​​سمكها 6 كم وفعالية 10 كم
سمك مرن في منطقة 3 كم متوسط ​​عمق الماء:

بتوقيت جرينتش جراففت topo.grd faa.grd -It -T10000/2700/3300/1035 -Z9000

لحساب القبول بين قياس الأعماق topo.grd وشبكة Geoid.grd مع
نموذج "التحميل من الأسفل" (LFB) بنفس الحمولة العلوية وتحت السطحية عند 40 كم،
لكن بافتراض أن الشبكات الآن جغرافية ونريد أطوال موجية بدلاً من التردد:

بتوقيت جرينتش جراففت 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

المراجع

لويس، JF وMC نيفيز. 2006، التعويض التوازني لهضبة جزر الأزور: ثلاثي الأبعاد
تحليل القبول والتماسك. J. الطاقة الحرارية الأرضية Volc. الدقة. المجلد 156، الأعداد 1-2، الصفحات
10-22، http://dx.doi.org/10.1016/j.jvolgeores.2006.03.010 باركر، ر.ل، 1972، السريع
حساب الشذوذات المحتملة، جيوفيس. ج.، 31، 447-455. فيسيل. ص، 2001، العالمية
توزيع الجبال البحرية المستنتجة من قياس الارتفاع الشبكي Geosat/ERS-1، J. Geophys. الدقة.،
106(ب9)، 19,431-19,441، http://dx.doi.org/10.1029/2000JB000083

استخدم gravfftgmt عبر الإنترنت باستخدام خدمات onworks.net