gravfftgmt - ออนไลน์ในคลาวด์

โครงการ:

ชื่อ

gravfft - คำนวณแรงดึงดูดของพื้นผิวสามมิติในเวฟนัมเบอร์ (หรือ
ความถี่) โดเมน

เรื่องย่อ

Grafft อิงกริด [ อิงกริด2 ] ออกจากไฟล์ [ n/ความยาวคลื่น/mean_ความลึก/tbw ] [ |โรกริด ] [
n_terms -f[+]|g|v|n|e- w|b|c|t |k -f|q|s|nx/นิวยอร์ก][+a|d|h
|l][+e|n|m][+tความกว้าง][+w[วิภัตติ-+z[p- เท/rl/rm/rw[+m-ระดับ- wd] [
zm[zl- -ฉ ]

หมายเหตุ ไม่อนุญาตให้มีช่องว่างระหว่างแฟล็กตัวเลือกและอาร์กิวเมนต์ที่เกี่ยวข้อง

DESCRIPTION

Grafft ใช้งานได้สามโหมดหลัก โหมด 1: เพียงคำนวณศักยภาพทางภูมิศาสตร์เนื่องจาก
พื้นผิวที่ระบุในไฟล์ topo.grd ต้องการคอนทราสต์ความหนาแน่น (-D) และอาจจะเป็น
ระดับการสังเกตที่แตกต่างกัน (-W). มันจะนำ FFT ไปข้างหน้า 2 มิติของกริดและใช้
วิธีการของ Parker แบบเต็มจนถึงเงื่อนไขที่เลือก โหมด 2: คำนวณศักยภาพทางภูมิศาสตร์
การตอบสนองเนื่องจากการงอของไฟล์ภูมิประเทศ มันจะนำ FFT ไปข้างหน้า 2 มิติของ
กริดและใช้วิธีการของ Parker แบบเต็มที่ใช้กับแบบจำลอง isostatic ที่เลือก NS
รุ่นที่มีจำหน่ายคือ "โหลดจากด้านบน" หรือรุ่นแผ่นยางยืดและ "โหลดจาก
ด้านล่าง" ซึ่งอธิบายการตอบสนองของเพลตต่อโหลดใต้พื้นผิว (เหมาะสมสำหรับความร้อน
การสร้างแบบจำลองเฉพาะจุด - ถ้าคุณเชื่อ) ในทั้งสองกรณี พารามิเตอร์โมเดลถูกกำหนดด้วย -T
และ -Z ตัวเลือก. โหมด 3: คำนวณการรับเข้าหรือการเชื่อมโยงกันระหว่างสองตาราง ผลลัพธ์
คือค่าเฉลี่ยในแนวรัศมี ทางเลือก การรับเข้าแบบจำลองอาจเป็น
คำนวณแล้ว ขนาดแนวนอนของ grdfiles จะถือว่าอยู่ในหน่วยเมตร
กริดทางภูมิศาสตร์อาจใช้โดยการระบุ -ฉ ตัวเลือกที่ปรับขนาดองศาเป็นเมตร
หากคุณมีตารางที่มีขนาดเป็นกม. คุณสามารถเปลี่ยนค่านี้เป็นเมตรได้โดยใช้ กริดดิต or
ปรับขนาดเอาต์พุตด้วย กริด. จากจำนวนตัวเลือกที่โปรแกรมนี้เสนอให้ คือ
ยากที่จะระบุว่าตัวเลือกคืออะไรและมีข้อโต้แย้งใดที่จำเป็น ขึ้นอยู่กับว่า
คุณกำลังทำ; ดูตัวอย่างสำหรับคำแนะนำเพิ่มเติม

ที่จำเป็น อาร์กิวเมนต์

อิงกริด ไฟล์กริดไบนารี 2 มิติที่จะเปิดใช้งาน (ดูรูปแบบไฟล์ GRID ด้านล่าง) สำหรับ
การดำเนินการข้ามสเปกตรัม ให้ไฟล์กริดที่สองด้วย ingrd2.

-Gออกจากไฟล์
ระบุชื่อไฟล์กริดเอาต์พุตหรือตารางสเปกตรัม 1 มิติ (ดู -E) (ดู
รูปแบบไฟล์ GRID ด้านล่าง)

ตัวเลือก อาร์กิวเมนต์

-Cn/ความยาวคลื่น/mean_ความลึก/tbw
คำนวณเฉพาะเส้นโค้งการรับเข้าตามทฤษฎีของแบบจำลองที่เลือกและทางออก n
และ ความยาวคลื่น ใช้ในการคำนวณ (n * ความยาวคลื่น) ความยาวโปรไฟล์ทั้งหมดเป็น
เมตร หมายถึง_ความลึก คือ ความลึกของน้ำเฉลี่ย ผนวก dataflags (หนึ่งหรือสอง) ของ TBW in
คำสั่งใดๆ t = ใช้โมเดล "จากด้านบน" b = ใช้รูปแบบ "จากด้านล่าง" เลือกระบุ
w เพื่อเขียนความยาวคลื่นแทนความถี่

-D |โรกริด
ตั้งค่าคอนทราสต์ของความหนาแน่นทั่วทั้งพื้นผิว ใช้สำหรับคำนวณแรงโน้มถ่วง
แรงดึงดูดของชั้นน้ำที่สามารถรวมเข้ากับความผิดปกติของอากาศอิสระได้ในภายหลัง
เพื่อรับความผิดปกติของ Bouguer ในกรณีนี้อย่าใช้ -T. นอกจากนี้ยังกำหนดโดยปริยาย
-N+ชม. อีกทางหนึ่ง ระบุกริดที่ลงทะเบียนร่วมด้วยคอนทราสต์ของความหนาแน่นหาก a
ต้องการคอนทราสต์ของความหนาแน่นแบบแปรผัน

-En_terms
จำนวนคำศัพท์ที่ใช้ในส่วนขยาย Parker (จำกัดคือ 10 มิฉะนั้นข้อกำหนดขึ้นอยู่กับ
n จะระเบิดโปรแกรมออก) [ค่าเริ่มต้น = 3]

-F[f[+]|g|v|n|e]
ระบุเขตข้อมูล geopotential ที่ต้องการ: คำนวณ geoid มากกว่าแรงโน้มถ่วง
f = ความผิดปกติของอากาศอิสระ (mGal) [ค่าเริ่มต้น] ผนวก + เพิ่มในแผ่นโดยนัย
เมื่อเอาค่าเฉลี่ยออกจากภูมิประเทศ สิ่งนี้ต้องการภูมิประเทศเป็นศูนย์
หมายถึงไม่มีความผิดปกติของมวล

g = ความผิดปกติของจีออยด์ (ม.)

v = การไล่ระดับความโน้มถ่วงในแนวตั้ง (VGG; 1 Eotvos = 0.1 mGal/km)

e = การโก่งตัวทางทิศตะวันออกของแนวตั้ง (ไมโครเรเดียน)

n = การโก่งตัวทางทิศเหนือของแนวตั้ง (ไมโครเรเดียน)

-Iw|ข|ค|t |k
ใช้ ingrd2 และ ingrd1 (ตารางที่มีภูมิประเทศ/บาธมาตรศาสตร์) ที่จะประมาณการ
การรับเข้าเรียน | การเชื่อมโยงกันและเขียนไปยัง stdout (-G ละเว้นหากตั้งค่าไว้) ตารางนี้ควร
มีแรงโน้มถ่วงหรือ geoid สำหรับภูมิภาคเดียวกันของ ingrd1. คอมพิวเตอร์เริ่มต้น
การรับเข้า ผลลัพธ์ประกอบด้วย 3 หรือ 4 คอลัมน์ ความถี่ (ความยาวคลื่น) การรับเข้า
(การเชื่อมโยงกัน) แถบข้อผิดพลาดซิกมาหนึ่งแถบและทางเลือกคือการรับเข้าตามทฤษฎี ผนวก
dataflags (หนึ่งถึงสาม) จาก w|b|c|t. w เขียนความยาวคลื่นแทนจำนวนคลื่น
k เลือกกิโลเมตรสำหรับหน่วยความยาวคลื่น [m], c คำนวณการเชื่อมโยงกันแทนการรับเข้า b
เขียนคอลัมน์ที่สี่ด้วยการรับเข้าตามทฤษฎี "โหลดจากด้านล่าง" และ t
เขียนคอลัมน์ที่สี่ด้วยการรับเข้าตามทฤษฎี "แผ่นยางยืด"

-N[ฉ|q|s|nx/นิวยอร์ก][+a|[+d|h|l][+e|n|m][+tความกว้าง][+w[วิภัตติ]][+z[p]]
เลือกหรือสอบถามเกี่ยวกับขนาดกริดที่เหมาะสมสำหรับ FFT และตั้งค่าตัวเลือกเพิ่มเติม
พารามิเตอร์ ควบคุมมิติ FFT:
-Nf จะบังคับให้ FFT ใช้มิติข้อมูลจริง

-Nq จะสอบถามเกี่ยวกับมิติข้อมูลที่เหมาะสมกว่า รายงานเหล่านั้น จากนั้นดำเนินการต่อ

-NS จะแสดงรายการมิติข้อมูลเสริม จากนั้นออก

-Nnx/นิวยอร์ก จะทำ FFT กับขนาดอาร์เรย์ nx/นิวยอร์ก (ต้องเป็น >= ขนาดไฟล์กริด) ค่าเริ่มต้น
เลือกขนาด >= ข้อมูลที่ปรับความเร็วและความแม่นยำของ FFT ให้เหมาะสมที่สุด ถ้าFFT
มิติ > ขนาดไฟล์กริด ข้อมูลถูกขยายและลดลงเหลือศูนย์

ควบคุมการปฏิเสธข้อมูล: ผนวกตัวดัดแปลงเพื่อลบแนวโน้มเชิงเส้น:
+d: ลบข้อมูลแนวโน้ม กล่าวคือ ลบแนวโน้มเชิงเส้นที่เหมาะสมที่สุด [ค่าเริ่มต้น]

+a: ลบเฉพาะค่ากลางเท่านั้น

+h: ลบเฉพาะค่ากลาง เช่น 0.5 * (สูงสุด + นาที)

+l: ปล่อยให้ข้อมูลอยู่คนเดียว

ควบคุมการขยายและการเรียวของข้อมูล: ใช้ตัวปรับแต่งเพื่อควบคุมวิธีการขยาย
และจะต้องทำการเทเปอร์:
+e ขยายตารางโดยกำหนดความสมมาตรของจุดขอบ [ค่าเริ่มต้น]

+m ขยายตารางโดยกำหนดความสมมาตรของกระจกขอบ

+n ปิดส่วนขยายข้อมูล

ทำการเรียวจากขอบข้อมูลไปยังขอบกริด FFT [100%] เปลี่ยน
เปอร์เซ็นต์นี้ผ่าน +tความกว้าง. เมื่อ +n มีผลใช้การเรียว
แทนขอบข้อมูลเนื่องจากไม่มีส่วนขยาย [0%]

ควบคุมการเขียนผลลัพธ์ชั่วคราว: สำหรับการตรวจสอบโดยละเอียด คุณสามารถเขียน
กริดกลางถูกส่งไปยัง FFT ไปข้างหน้า นี้น่าจะเป็น
detrended ขยายโดยสมมาตรจุดตามขอบทั้งหมดและเรียว ผนวก
+w[วิภัตติ] ซึ่งชื่อไฟล์เอาต์พุตจะถูกสร้างขึ้น (เช่น ingrid_prefix.ext)
[เรียว] โดยที่ ต่อ เป็นนามสกุลไฟล์ของคุณ สุดท้าย คุณอาจบันทึกตารางที่ซับซ้อนได้
ผลิตโดยไปข้างหน้า FFT โดยผนวก +z. โดยค่าเริ่มต้นเราเขียนของจริงและ
องค์ประกอบจินตภาพถึง อิงกริด_จริง.ต่อ และ อิงกริด_ภาพ.ต่อ. ผนวก p เพื่อประหยัด
แทนรูปแบบขั้วของขนาดและเฟสเป็นไฟล์ อิงกริด_แม็ก.ต่อ และ
อิงกริด_เฟส.ต่อ.

-Q เขียนตารางที่มีภูมิประเทศดัดโค้ง (โดยมีค่า z เป็นบวก) ซึ่งมีค่าเฉลี่ย
ถูกกำหนดโดย -Zzm และพารามิเตอร์โมเดลโดย -T (และส่งออกโดย -G). นั่นคือ
"กราวิเมตริก โมโฮ" -Q กำหนดโดยปริยาย -N+ชม

-S คำนวณความโน้มถ่วงที่คาดการณ์ไว้หรือตาราง geoid เนื่องจากโหลดแผ่นย่อยที่ผลิตโดย
การวัดปริมาณน้ำในปัจจุบันและแบบจำลองทางทฤษฎี มีการตั้งค่าพารามิเตอร์ที่จำเป็น
ภายใน -T และ -Z ตัวเลือก. จำนวนอำนาจในการขยาย Parker ถูก จำกัด ไว้ที่
1. ดูตัวอย่างเพิ่มเติมด้านล่าง

-Tเท/rl/rm/rw[+ม.]
คำนวณการชดเชย isostatic จากภาระภูมิประเทศ (ไฟล์กริดอินพุต) บน an
แผ่นยางยืดหนา te. ผนวกความหนาแน่นของน้ำหนักบรรทุก เสื้อคลุม และน้ำเข้าด้วย
หน่วย SI ให้ความลึกของเสื้อคลุมเฉลี่ยผ่าน -Z. ถ้าความหนายืดหยุ่น > 1e10 ก็
จะถูกตีความว่าเป็นความแข็งแกร่งของการดัด (โดยค่าเริ่มต้นจะคำนวณจาก te และ
โมดูลัสหนุ่ม) หรือไม่ก็ได้ ต่อท้าย +m เพื่อเขียนตารางด้วยศักยภาพทางภูมิศาสตร์ของ Moho
ผล (ดู -F) จากรุ่นที่เลือกโดย -T. ถ้า te = 0 ดังนั้นการตอบสนองที่โปร่งสบายคือ
กลับ -T+ม กำหนดโดยปริยาย -N+ชม

-Wwd ตั้งค่าความลึกของน้ำ (หรือความสูงของการสังเกต) ที่สัมพันธ์กับภูมิประเทศ [0] ผนวก k ไปยัง
ระบุกม.

-Zzm[zl]
Moho [และบวม] ความลึกของการชดเชยโดยเฉลี่ย สำหรับรุ่น "โหลดจากด้านบน" คุณ
ต้องจัดให้เท่านั้น zmแต่สำหรับการ "โหลดจากด้านล่าง" อย่าลืม zl.

-วี[ระดับ] (มากกว่า ... )
เลือกระดับการใช้คำฟุ่มเฟือย [c]

-ฉ กริดทางภูมิศาสตร์ (ขนาดลองจิจูด ละติจูด) จะถูกแปลงเป็นเมตร
ผ่านการประมาณ "Flat Earth" โดยใช้พารามิเตอร์ทรงรีปัจจุบัน

-^ or เพียงแค่ -
พิมพ์ข้อความสั้น ๆ เกี่ยวกับไวยากรณ์ของคำสั่ง จากนั้นออก (หมายเหตุ: บน Windows
ใช้เพียง -).

-+ or เพียงแค่ +
พิมพ์ข้อความการใช้งาน (ช่วยเหลือ) อย่างกว้างขวาง รวมถึงคำอธิบายของทุก ๆ
ตัวเลือกเฉพาะโมดูล (แต่ไม่ใช่ตัวเลือกทั่วไปของ GMT) จากนั้นออก

-? or ไม่ ข้อโต้แย้ง
พิมพ์ข้อความการใช้งาน (วิธีใช้) ที่สมบูรณ์ รวมทั้งคำอธิบายตัวเลือก จากนั้น
ทางออก

--รุ่น
พิมพ์เวอร์ชัน GMT และออก

--show-datadir
พิมพ์เส้นทางแบบเต็มไปยังไดเรกทอรีแชร์ GMT และออก

GRID ไฟล์ รูปแบบ

โดยค่าเริ่มต้น GMT จะเขียนกริดออกเป็นความแม่นยำเดียวใน netCDF . ร้องเรียน COARDS
รูปแบบไฟล์. อย่างไรก็ตาม GMT สามารถสร้างไฟล์กริดในกริดอื่น ๆ ที่ใช้กันทั่วไปได้
รูปแบบไฟล์และยังอำนวยความสะดวกที่เรียกว่า "การบรรจุ" ของกริดการเขียนจุดลอยตัว
ข้อมูลเป็นจำนวนเต็ม 1 หรือ 2 ไบต์ ในการระบุความแม่นยำ มาตราส่วน และออฟเซ็ต ผู้ใช้ควร
เพิ่มคำต่อท้าย =id[/ขนาด/ชดเชย[/น่าน]], ที่ไหน id เป็นตัวระบุสองตัวอักษรของกริด
ชนิดและความแม่นยำ และ ขนาด และ ชดเชย เป็นสเกลแฟกเตอร์เสริมและออฟเซ็ตเป็น
นำไปใช้กับค่ากริดทั้งหมดและ น่าน คือค่าที่ใช้ระบุข้อมูลที่ขาดหายไป เผื่อ
ตัวละครทั้งสอง id ไม่ได้ระบุไว้เช่นใน =/ขนาด กว่า a id=nf จะถือว่า เมื่อไหร่
อ่านตาราง รูปแบบโดยทั่วไปจะรู้จักโดยอัตโนมัติ ถ้าไม่ใช่ คำต่อท้ายเดียวกัน
สามารถเพิ่มลงในชื่อไฟล์กริดอินพุต ดู แปลงกริด และรูปแบบไฟล์กริดมาตราของ
ข้อมูลอ้างอิงทางเทคนิค GMT และตำราอาหารสำหรับข้อมูลเพิ่มเติม

เมื่ออ่านไฟล์ netCDF ที่มีหลายกริด GMT จะอ่านตามค่าเริ่มต้น
ตาราง 2 มิติแรกที่หาได้ในไฟล์นั้น เพื่อเกลี้ยกล่อม GMT ให้อ่านอีก
ตัวแปรหลายมิติในไฟล์กริด ผนวก ?นามสกุล ไปที่ชื่อไฟล์โดยที่
นามสกุล เป็นชื่อของตัวแปร โปรดทราบว่าคุณอาจต้องหนีความหมายพิเศษ
of ? ในโปรแกรมเชลล์ของคุณโดยใส่แบ็กสแลชไว้ข้างหน้า หรือโดยการวาง
ชื่อไฟล์และคำต่อท้ายระหว่างเครื่องหมายคำพูดหรือเครื่องหมายคำพูดคู่ NS ?นามสกุล ต่อท้ายก็ใช้ได้นะ
สำหรับกริดเอาต์พุตเพื่อระบุชื่อตัวแปรที่แตกต่างจากค่าเริ่มต้น: "z" ดู
แปลงกริด และส่วนดัดแปลงสำหรับ CF และรูปแบบไฟล์กริดของ GMT Technical
อ้างอิงและตำราอาหารสำหรับข้อมูลเพิ่มเติม โดยเฉพาะอย่างยิ่งในการอ่าน splices ของ 3-,
กริด 4 หรือ 5 มิติ

GRID ระยะทาง UNITS

หากกริดไม่มีมิเตอร์เป็นหน่วยแนวนอน ให้ผนวก +uหน่วย ไปยังไฟล์อินพุต
ชื่อที่จะแปลงจากหน่วยที่ระบุเป็นเมตร หากกริดของคุณเป็นพื้นที่ทางภูมิศาสตร์ ให้แปลง
ระยะทางเป็นเมตรโดยการจัดหา -ฉ แทน.

ที่ต้องคำนึงถึง

กริด netCDF COARDS จะรับรู้โดยอัตโนมัติตามพื้นที่ทางภูมิศาสตร์ สำหรับกริดอื่นๆ
กริดทางภูมิศาสตร์ที่คุณต้องการแปลงองศาเป็นเมตร ให้เลือก -ฉ. ถ้าข้อมูล
อยู่ใกล้กับขั้วใดขั้วหนึ่ง คุณควรพิจารณาฉายไฟล์กริดลงบนสี่เหลี่ยม
ระบบพิกัดโดยใช้ โครงการ grd.

PLATE เฟล็กเชอร์

สารละลาย FFT สำหรับการดัดงอของแผ่นยางยืดต้องการความหนาแน่นของวัสดุเติมเพื่อให้เท่ากับน้ำหนักบรรทุก
ความหนาแน่น. โดยทั่วไปแล้วสิ่งนี้เป็นจริงภายใต้การโหลดเท่านั้น เกินภาระ
infill มีแนวโน้มที่จะเป็นตะกอนที่มีความหนาแน่นต่ำกว่าหรือแม้กระทั่งน้ำ (หรืออากาศ) Wessel [2001] เสนอ
การประมาณที่ช่วยให้สามารถระบุความหนาแน่นของ infill ที่แตกต่างจาก
ความหนาแน่นของโหลดในขณะที่ยังคงอนุญาตให้ใช้โซลูชัน FFT โดยทั่วไป ความโค้งของแผ่นคือ
แก้ไขโดยใช้ความหนาแน่นของ infill เป็นความหนาแน่นของโหลดที่มีประสิทธิภาพ แต่แอมพลิจูดคือ
ปรับด้วยปัจจัย A = sqrt ((rm - ri)/(rm - rl)) ซึ่งเป็นข้อแตกต่างทางทฤษฎี
ในแอมพลิจูดเนื่องจากการโหลดแบบจุดโดยใช้ความหนาแน่นโหลดสองแบบที่แตกต่างกัน NS
การประมาณนั้นดีมาก แต่พังทลายลงสำหรับการบรรทุกขนาดใหญ่บนจานที่อ่อนแอ นางฟ้า
สถานการณ์ที่ไม่ปกติ

ตัวอย่าง

ในการคำนวณผลกระทบของชั้นน้ำเหนือ bat.grd bathymetry โดยใช้ 2700 และ 1035
สำหรับความหนาแน่นของเปลือกโลกและน้ำ และเขียนผลลัพธ์บน water_g.grd (คำนวณค่า
สู่พลังที่สี่ของการวัดปริมาณน้ำในการขยายตัวของ Parker):

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

ตอนนี้ลบไปที่ faa.grd ความผิดปกติในอากาศฟรีของคุณและคุณจะได้รับความผิดปกติของ Bouguer คุณ
อาจสงสัยว่าทำไมเราลบแล้วไม่บวก หลังจากที่ทุกความผิดปกติของ Bouguer แสร้งทำเป็น
เพื่อแก้ไขการขาดมวลที่นำเสนอโดยชั้นน้ำดังนั้นเราควรเพิ่มเพราะ
น้ำมีความหนาแน่นน้อยกว่าหินด้านล่าง คำตอบจะขึ้นอยู่กับผลกระทบของแรงโน้มถ่วง
คำนวณโดยวิธีของ Parker และแง่มุมเชิงปฏิบัติของการใช้ FFT

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

ต้องการความผิดปกติ MBA หรือไม่? คำนวณการมีส่วนร่วมของเปลือกโลกแล้วเพิ่มลงใน
ความผิดปกติของก้นทะเล สมมติว่าเปลือกโลกหนา 6 กม. มีความหนาแน่น 2700 และเสื้อคลุมด้วย3300
ความหนาแน่น เราสามารถทำซ้ำคำสั่งที่ใช้ในการคำนวณความผิดปกติของชั้นน้ำโดยใช้600
(3300 - 2700) ตามความหนาแน่นที่ตัดกัน แต่ตอนนี้เรามีปัญหาเพราะเราต้องรู้
ความลึกเฉลี่ยของ Moho นั่นคือเมื่อสเกล/ออฟเซ็ตที่สามารถต่อท้ายชื่อกริดได้
มาอยู่ในมือ สังเกตว่าเราไม่จำเป็นต้องทำอย่างนั้นมาก่อนเพราะความลึกของน้ำเฉลี่ยคือ
คำนวณโดยตรงจากข้อมูล (โปรดสังเกตเครื่องหมายลบของการชดเชยเนื่องจากข้อเท็จจริง
ที่ z เป็นบวกขึ้น):

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

ทีนี้ ลบมันด้วยความผิดปกติก้นทะเลเพื่อให้ได้ความผิดปกติ MBA นั่นคือ:

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

ในการคำนวณเอฟเฟกต์แรงโน้มถ่วงของ Moho ของแผ่นยืดหยุ่น bat.grd ด้วย Te = 7 กม. ความหนาแน่นของ
2700 บนเสื้อคลุมที่มีความหนาแน่น 3300 ที่ความลึกเฉลี่ย 9 กม.

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

หากคุณเพิ่มเอฟเฟกต์ของก้นทะเลและ Moho ในตอนนี้ คุณจะได้รับการตอบสนองต่อแรงโน้มถ่วงเต็มรูปแบบ
ของแบบจำลอง isostatic ของคุณ เราจะใช้เฉพาะเทอมแรกในส่วนขยาย Parker

gmt Grafft bat.grd -D1665 -Gwater_g.grd -E1
gmt Grafft bat.grd -Gelastic.grd -T7000/2700/3300/1035+m -Z9000 -E1
gmt grdmath water_g.grd elastic.grd เพิ่ม = model.grd

ผลลัพธ์เดียวกันสามารถรับได้โดยตรงจากคำสั่งถัดไป อย่างไรก็ตามให้ความสนใจกับ
ต่อไปนี้ ฉันยังไม่รู้ว่าเป็นเพราะข้อบกพร่องหรือเพราะข้อจำกัดบางอย่าง แต่
ความจริงก็คือคำสั่งต่อไปนี้และคำสั่งก่อนหน้าจะให้ผลลัพธ์เดียวกัน if . เท่านั้น -E1
ถูกนำมาใช้. สำหรับพลังการอาบน้ำที่สูงขึ้นในการขยายตัวของ Parker เฉพาะตัวอย่างข้างต้นเท่านั้น
ตะเข็บเพื่อให้ได้ผลลัพธ์ที่ถูกต้อง

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

และสิ่งที่จะเป็นความผิดปกติ geoid ที่เกิดจากโหลดที่ความลึก 50 กม. ใต้พื้นที่
โดย bat.grd ซึ่งเป็น Moho ที่ความลึก 9 กม. และมีความหนาแน่นเท่ากับ
ก่อน?

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

เพื่อคำนวณการอนุญาติให้เข้าระหว่าง topo.grd bathymetry และ faa.grd free-air anomaly
ตารางโดยใช้แบบจำลองแผ่นยางยืดของเปลือกโลกที่มีความหนาเฉลี่ย 6 กม. ที่มีประสิทธิภาพ 10 กม
ความหนายืดหยุ่นในพื้นที่ 3 กม. ความลึกของน้ำเฉลี่ย:

gmt gravfft topo.grd faa.grd -มัน -T10000/2700/3300/1035 -Z9000

ในการคำนวณการอนุญาติระหว่าง topo.grd bathymetry และ geoid.grd geoid grid ด้วย
รุ่น "โหลดจากด้านล่าง" (LFB) ที่เหมือนกันกับด้านบนและโหลดใต้พื้นผิวที่ 40 กม.
แต่สมมติว่าตอนนี้กริดอยู่ในพื้นที่ และเราต้องการความยาวคลื่นแทนความถี่:

GMT Grafft 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 การชดเชยไอโซสแตติกของที่ราบสูงอะซอเรส: a 3D
การรับเข้าและการวิเคราะห์การเชื่อมโยงกัน เจ. ความร้อนใต้พิภพ Volc. ความละเอียด เล่มที่ 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, Global
การกระจายของภูเขาทะเลที่อนุมานจากการวัดระยะสูง Geosat/ERS-1, J. Geophys ความละเอียด
106(B9), 19,431-19,441, http://dx.doi.org/10.1029/2000JB000083

ใช้ gravfftgmt ออนไลน์โดยใช้บริการ onworks.net