gravfftgmt - ອອນລາຍໃນຄລາວ

ໂຄງການ:

NAME

gravfft - ຄິດໄລ່ແຮງໂນ້ມຖ່ວງຂອງພື້ນຜິວ 3-D ໃນ wavenumber (ຫຼື
ຄວາມຖີ່) ໂດເມນ

ສະຫຼຸບສັງລວມ

gravfft ingrid [ ເສັ້ນ​ທາງ 2 ] outfile [ n/wavelength/mean_depth/tbw ] [ ຄວາມ​ຫນາ​ແຫນ້ນ|rhogrid ] [
n_terms ][f[+]|g|v|n|e] ] [ w|b|c|t |k ][f|q|s|nx/ny][+a|d|h
|l][+e|n|m][+twidth][+w[បច្ច័យ]][+z[p]] [ ] [ te/rl/rm/rw[+m] ] [ລະດັບ] ] [ wd] [
zm[zl] ] [ -fg ]

ຫມາຍ​ເຫດ​: ບໍ່ອະນຸຍາດໃຫ້ມີຊ່ອງຫວ່າງລະຫວ່າງທຸງທາງເລືອກ ແລະອາກິວເມັນທີ່ກ່ຽວຂ້ອງ.

ລາຍລະອຽດ

gravfft ສາ​ມາດ​ນໍາ​ໃຊ້​ເປັນ​ສາມ​ທິ​ບາຍ​ຮູບ​ແບບ​ຕົ້ນ​ຕໍ​. ໂຫມດ 1: ພຽງແຕ່ຄິດໄລ່ geopotential ເນື່ອງຈາກ
ພື້ນຜິວທີ່ໃຫ້ໄວ້ໃນໄຟລ໌ topo.grd. ຕ້ອງການຄວາມຄົມຊັດຄວາມຫນາແຫນ້ນ (-D) ແລະ​ອາດ​ຈະ​ເປັນ​
ລະ​ດັບ​ການ​ສັງ​ເກດ​ທີ່​ແຕກ​ຕ່າງ​ກັນ (-W). ມັນຈະໃຊ້ເວລາ 2-D forward FFT ຂອງຕາຂ່າຍໄຟຟ້າແລະນໍາໃຊ້
ວິທີການຂອງ Parker ເຕັມໄປເຖິງເງື່ອນໄຂທີ່ເລືອກ. ຮູບແບບທີ 2: ຄິດໄລ່ທ່າແຮງທາງພູມສາດ
ການຕອບສະໜອງເນື່ອງຈາກການຢືດຕົວຂອງໄຟລ໌ພູມສັນຖານ. ມັນຈະໃຊ້ເວລາ 2-D forward FFT ຂອງ
ຕາຂ່າຍໄຟຟ້າແລະນໍາໃຊ້ວິທີການ Parker ເຕັມທີ່ນໍາໃຊ້ກັບຮູບແບບ isostatic ທີ່ເລືອກ. ໄດ້
ຮູບແບບທີ່ມີຢູ່ແມ່ນ "ການໂຫຼດຈາກເທິງ", ຫຼືຮູບແບບແຜ່ນ elastic, ແລະ "ການໂຫຼດຈາກ
ຂ້າງລຸ່ມນີ້" ເຊິ່ງກວມເອົາການຕອບສະຫນອງຂອງແຜ່ນຕໍ່ການໂຫຼດຂອງພື້ນຜິວຍ່ອຍ (ເຫມາະສົມສໍາລັບຄວາມຮ້ອນ
ການສ້າງແບບຈໍາລອງຈຸດ - ຖ້າທ່ານເຊື່ອວ່າພວກເຂົາ). ໃນທັງສອງກໍລະນີ, ຕົວກໍານົດການຂອງຕົວແບບແມ່ນຖືກກໍານົດດ້ວຍ -T
ແລະ -Z ທາງເລືອກ. ຮູບແບບທີ 3: ຄິດໄລ່ການຮັບ ຫຼື ຄວາມສອດຄ່ອງລະຫວ່າງສອງຕາຂ່າຍ. ຜົນຜະລິດ
ແມ່ນສະເລ່ຍໃນທິດທາງ radial. ທາງເລືອກ, ການຍອມຮັບແບບຈໍາລອງອາດຈະເປັນ
ຄິດໄລ່. ຂະຫນາດຕາມລວງນອນຂອງ grdfiles ແມ່ນສົມມຸດວ່າເປັນແມັດ.
ຕາຂ່າຍທາງພູມສາດອາດຈະຖືກນໍາໃຊ້ໂດຍການລະບຸ -fg ທາງ​ເລືອກ​ທີ່​ຂະ​ຫຍາຍ​ອົງ​ສາ​ເປັນ​ແມັດ​.
ຖ້າທ່ານມີຕາຂ່າຍໄຟຟ້າທີ່ມີຂະຫນາດເປັນກິໂລແມັດ, ທ່ານສາມາດປ່ຽນເປັນແມັດໄດ້ໂດຍໃຊ້ grdedit or
ຂະຫນາດຜົນຜະລິດທີ່ມີ grdmath. ໃຫ້ຈໍານວນຂອງທາງເລືອກທີ່ໂຄງການນີ້ສະເຫນີ, ແມ່ນ
ຍາກທີ່ຈະລະບຸສິ່ງທີ່ເປັນທາງເລືອກແລະການໂຕ້ຖຽງທີ່ຕ້ອງການ. ມັນຂຶ້ນກັບສິ່ງທີ່
ເຈົ້າກໍາລັງເຮັດ; ເບິ່ງຕົວຢ່າງສໍາລັບຄໍາແນະນໍາເພີ່ມເຕີມ.

ຕ້ອງການ ການໂຕ້ຖຽງ

ingrid ໄຟລ໌ຕາຕະລາງສອງມິຕິ 2-D ທີ່ຈະດໍາເນີນການ. (ເບິ່ງຮູບແບບໄຟລ໌ຕາຕະລາງຂ້າງລຸ່ມນີ້). ສໍາລັບ
ການດໍາເນີນງານຂ້າມ spectral, ຍັງໃຫ້ໄຟລ໌ຕາຂ່າຍໄຟຟ້າທີສອງ ສ່ວນປະກອບ2.

-Goutfile
ລະບຸຊື່ຂອງໄຟລ໌ຕາຕະລາງຜົນຜະລິດ ຫຼືຕາຕະລາງ 1-D spectrum (ເບິ່ງ -E) (ເບິ່ງ
ຮູບແບບໄຟລ໌ຕາຕະລາງຂ້າງລຸ່ມນີ້).

ທາງເລືອກ ການໂຕ້ຖຽງ

-Cn/wavelength/mean_depth/tbw
ຄິດໄລ່ພຽງແຕ່ເສັ້ນໂຄ້ງການຍອມຮັບທາງທິດສະດີຂອງຕົວແບບທີ່ເລືອກແລະອອກ. n
ແລະ wavelength ຖືກນໍາໃຊ້ເພື່ອຄິດໄລ່ (n * wavelength) ຄວາມຍາວຂອງໂປຣໄຟລ໌ທັງຫມົດໃນ
ແມັດ. mean_depth ແມ່ນຄວາມເລິກຂອງນ້ໍາສະເລ່ຍ. ຕື່ມຂໍ້ມູນໃສ່ທຸງ (ໜຶ່ງ ຫຼື ສອງ) ຂອງ tbw in
ຄໍາສັ່ງໃດໆ. t = ໃຊ້ຕົວແບບ "ຈາກເທິງ", b = ໃຊ້ຕົວແບບ "ຈາກລຸ່ມນີ້". ກໍານົດທາງເລືອກ
w ຂຽນ wavelength ແທນຄວາມຖີ່.

-Dຄວາມ​ຫນາ​ແຫນ້ນ|rhogrid
ຕັ້ງຄ່າຄວາມຄົມຊັດຂອງຄວາມຫນາແຫນ້ນໃນທົ່ວພື້ນຜິວ. ຕົວຢ່າງ, ໃຊ້ເພື່ອຄິດໄລ່ແຮງໂນ້ມຖ່ວງ
ຄວາມດຶ່ງດູດຂອງຊັ້ນນ້ໍາທີ່ຕໍ່ມາສາມາດຖືກລວມເຂົ້າກັບຄວາມຜິດປົກກະຕິຂອງອາກາດທີ່ບໍ່ເສຍຄ່າ
ເພື່ອໃຫ້ໄດ້ຮັບຄວາມຜິດປົກກະຕິ Bouguer. ໃນກໍລະນີນີ້, ຢ່າໃຊ້ -T. ມັນຍັງກໍານົດ implicitly
-N+ຊ. ອີກທາງເລືອກ, ລະບຸຕາໜ່າງທີ່ລົງທະບຽນຮ່ວມກັບຄວາມໜາແໜ້ນຂອງຄວາມຄົມຊັດ ຖ້າ a
ຄວາມຫນາແຫນ້ນຂອງຕົວປ່ຽນແປງແມ່ນຕ້ອງການ.

-En_terms
ຈໍາ​ນວນ​ຂອງ​ຄໍາ​ສັບ​ຕ່າງໆ​ທີ່​ນໍາ​ໃຊ້​ໃນ​ການ​ຂະ​ຫຍາຍ Parker (ຂໍ້​ຈໍາ​ກັດ​ແມ່ນ 10​, ຖ້າ​ບໍ່​ດັ່ງ​ນັ້ນ​ຂໍ້​ກໍາ​ນົດ​ຂຶ້ນ​ກັບ​
n ຈະລະເບີດອອກໂຄງການ) [Default = 3]

-F[f[+]|g|v|n|e]
ລະບຸຊ່ອງຂໍ້ມູນພູມສາດທີ່ຕ້ອງການ: ຄິດໄລ່ geoid ຫຼາຍກວ່າແຮງໂນ້ມຖ່ວງ
f = ຜິດປົກກະຕິທາງອາກາດ (mGal) [ຄ່າເລີ່ມຕົ້ນ]. ຕື່ມໃສ່ + ຕື່ມໃນຝາອັດປາກຂຸມ implied
ເມື່ອເອົາຄ່າສະເລ່ຍອອກຈາກພູມສັນຖານ. ອັນນີ້ຕ້ອງການພູມສັນຖານສູນ
ຫມາຍຄວາມວ່າບໍ່ມີຄວາມຜິດປົກກະຕິຂອງມະຫາຊົນ.

g = ຜິດປົກກະຕິ Geoid (m).

v = Gravity Gradient (VGG; 1 Eotvos = 0.1 mGal/km).

e = ທິດຕາເວັນອອກ deflections ຂອງແນວຕັ້ງ (micro-radian).

n = ການເໜັງຕີງທາງເໜືອຂອງແນວຕັ້ງ (micro-radian).

-Iw|b|c|t |k
ການນໍາໃຊ້ ສ່ວນປະກອບ2 ແລະ ສ່ວນປະກອບ1 (ຕາຂ່າຍໄຟຟ້າທີ່ມີພູມສັນຖານ/bathymetry) ເພື່ອຄາດຄະເນ
admittance|ຄວາມສອດຄ່ອງກັນ ແລະຂຽນໃສ່ stdout (-G ບໍ່ສົນໃຈຖ້າຕັ້ງ). ຕາຂ່າຍໄຟຟ້ານີ້ຄວນ
ມີແຮງໂນ້ມຖ່ວງ ຫຼື geoid ສໍາລັບພາກພື້ນດຽວກັນຂອງ ສ່ວນປະກອບ1. ຄອມພິວເຕີເລີ່ມຕົ້ນ
ເສັງເຂົ້າ. ຜົນຜະລິດມີ 3 ຫຼື 4 ຖັນ. ຄວາມຖີ່ (ຄື້ນຄວາມຖີ່), ການຮັບ
(ຄວາມສອດຄ່ອງກັນ) ແຖບຄວາມຜິດພາດຫນຶ່ງ sigma ແລະ, ທາງເລືອກ, ການຍອມຮັບທາງທິດສະດີ. ຕື່ມໃສ່
dataflags (ຫນຶ່ງຫາສາມ) ຈາກ w|b|c|t. w ຂຽນ wavelength ແທນ wavenumber,
k ເລືອກກິໂລແມັດສຳລັບຫົວໜ່ວຍຄວາມຍາວຄື້ນ [m], c ຄິດ​ໄລ່​ຄວາມ​ສອດ​ຄ່ອງ​ແທນ​ທີ່​ຈະ​ເປັນ​ການ​ຍອມ​ຮັບ​, b
ຂຽນຖັນທີສີ່ດ້ວຍ "ການໂຫຼດຈາກຂ້າງລຸ່ມ" ການຍອມຮັບທາງທິດສະດີ, ແລະ t
ຂຽນຖັນທີສີ່ດ້ວຍ "ແຜ່ນ elastic" ການຍອມຮັບທາງທິດສະດີ.

-N[f|q|s|nx/ny][+a|[+d|h|l][+e|n|m][+twidth][+w[បច្ច័យ]][+z[p]]
ເລືອກຫຼືສອບຖາມກ່ຽວກັບຂະຫນາດຕາຂ່າຍໄຟຟ້າທີ່ເຫມາະສົມສໍາລັບ FFT ແລະກໍານົດທາງເລືອກ
ຕົວກໍານົດການ. ຄວບ​ຄຸມ​ຂະ​ຫນາດ FFT​:
-Nf ຈະບັງຄັບໃຫ້ FFT ໃຊ້ຂະຫນາດຕົວຈິງຂອງຂໍ້ມູນ.

- ນຄ ຈະສອບຖາມກ່ຽວກັບຂະຫນາດທີ່ເຫມາະສົມກວ່າ, ລາຍງານເຫຼົ່ານັ້ນ, ຫຼັງຈາກນັ້ນສືບຕໍ່.

- ນ ຈະນໍາສະເຫນີບັນຊີລາຍຊື່ຂອງຂະຫນາດທາງເລືອກ, ຫຼັງຈາກນັ້ນອອກ.

-Nnx/ny ຈະເຮັດ FFT ໃນຂະຫນາດ array nx/ny (ຕ້ອງເປັນ >= ຂະໜາດໄຟລ໌ຕາຂ່າຍ). ຄ່າເລີ່ມຕົ້ນ
ເລືອກຂະໜາດ >= ຂໍ້ມູນທີ່ເພີ່ມປະສິດທິພາບຄວາມໄວ ແລະຄວາມຖືກຕ້ອງຂອງ FFT. ຖ້າ FFT
ຂະ​ຫນາດ​> ຂະ​ຫນາດ​ໄຟລ​໌​ຕາ​ຂ່າຍ​ໄຟ​ຟ້າ​, ຂໍ້​ມູນ​ແມ່ນ​ຂະ​ຫຍາຍ​ແລະ tapered ເປັນ​ສູນ​.

ຄວບ​ຄຸມ​ການ​ຫຼຸດ​ລົງ​ຂອງ​ຂໍ້​ມູນ​: ເພີ່ມ​ເຕີມ​ຕົວ​ປັບ​ປຸງ​ສໍາ​ລັບ​ການ​ລົບ​ແນວ​ໂນ້ມ​ເສັ້ນ​:
+d: Detrend data, ie ເອົາ best-fitting linear trend [Default].

+a: ເອົາພຽງແຕ່ຄ່າສະເລ່ຍ.

+h: ພຽງແຕ່ເອົາຄ່າກາງ, ie 0.5 * (max + min).

+l: ປ່ອຍຂໍ້ມູນໄວ້ຢ່າງດຽວ.

ຄວບ​ຄຸມ​ການ​ຂະ​ຫຍາຍ​ແລະ tapering ຂອງ​ຂໍ້​ມູນ​: ໃຊ້​ຕົວ​ປັບ​ປຸງ​ເພື່ອ​ຄວບ​ຄຸມ​ວິ​ທີ​ການ​ຂະ​ຫຍາຍ​
ແລະ tapering ແມ່ນ​ຈະ​ໄດ້​ຮັບ​ການ​ປະ​ຕິ​ບັດ​:
+e ຂະຫຍາຍຕາຂ່າຍໂດຍການວາງສົມມາຕຣິກຈຸດຂອບ [ຄ່າເລີ່ມຕົ້ນ],

+m ຂະຫຍາຍຕາໜ່າງໂດຍການວາງຂອບບ່ອນແລກປ່ຽນຄວາມສົມມາ

+n ປິດການຂະຫຍາຍຂໍ້ມູນ.

Tapering ແມ່ນປະຕິບັດຈາກຂອບຂໍ້ມູນໄປຫາຂອບຕາຂ່າຍໄຟຟ້າ FFT [100%]. ປ່ຽນແປງ
ອັດຕາສ່ວນນີ້ຜ່ານ +twidthທີ່ຢູ່ ເມື່ອ​ໃດ​ +n ມີຜົນບັງຄັບໃຊ້, tapering ຖືກນໍາໃຊ້
ແທນທີ່ຂອບຂໍ້ມູນເນື່ອງຈາກບໍ່ມີການຂະຫຍາຍ [0%].

ຄວບຄຸມການຂຽນຜົນໄດ້ຮັບຊົ່ວຄາວ: ສໍາລັບການສືບສວນລາຍລະອຽດທ່ານສາມາດຂຽນໄດ້
ຕາຂ່າຍໄຟຟ້າລະດັບປານກາງຖືກສົ່ງກັບ FFT ໄປຂ້າງຫນ້າ; ນີ້ເປັນໄປໄດ້
detreded, ຂະຫຍາຍໂດຍຈຸດ, symmetry ຕາມແຄມທັງຫມົດ, ແລະ tapered. ຕື່ມໃສ່
+w[បច្ច័យ] ຈາກ​ທີ່​ຊື່​ໄຟລ​໌​ຜົນ​ຜະ​ລິດ​ຈະ​ໄດ້​ຮັບ​ການ​ສ້າງ​ຕັ້ງ (ເຊັ່ນ​: ingrid_prefix.ext)
[tapered], ບ່ອນທີ່ ext ແມ່ນ​ການ​ຂະ​ຫຍາຍ​ໄຟລ​໌​ຂອງ​ທ່ານ​. ສຸດທ້າຍ, ທ່ານອາດຈະຊ່ວຍປະຢັດຕາຂ່າຍໄຟຟ້າສະລັບສັບຊ້ອນ
ຜະລິດໂດຍ FFT ຂ້າງຫນ້າໂດຍການເພີ່ມເຕີມ +z. ໂດຍຄ່າເລີ່ມຕົ້ນພວກເຮົາຂຽນທີ່ແທ້ຈິງແລະ
ອົງ​ປະ​ກອບ​ຂອງ​ຈິນ​ຕະ​ນາ​ການ​ ingrid_ຈິງ.ext ແລະ ingrid_ຮູບ.ext. ຕື່ມໃສ່ p ບັນ​ທຶກ
ແທນທີ່ຈະເປັນຮູບແບບຂົ້ວໂລກຂອງຂະຫນາດແລະໄລຍະກັບໄຟລ໌ ingrid_ແມັກ.ext ແລະ
ingrid_ໄລຍະ.ext.

-Q ຂຽນຕາຕະລາງທີ່ມີພູມສັນຖານ flexural (ກັບ z ບວກຂຶ້ນ) ເຊິ່ງສະເລ່ຍ
ຖືກກໍານົດໂດຍ -Zzm ແລະຕົວກໍານົດການແບບຈໍາລອງໂດຍ -T (ແລະຜົນຜະລິດໂດຍ -G). ນັ້ນແມ່ນ
"gravimetric Moho". -Q ກໍານົດ implicitly -N+ຊ

-S ຄິດໄລ່ແຮງໂນ້ມຖ່ວງທີ່ຄາດຄະເນຫຼືຕາຂ່າຍໄຟຟ້າ geoid ເນື່ອງຈາກການໂຫຼດ subplate ທີ່ຜະລິດໂດຍ
bathymmetry ໃນ​ປະ​ຈຸ​ບັນ​ແລະ​ຕົວ​ແບບ​ທາງ​ທິດ​ສະ​ດີ​. ຕົວກໍານົດການທີ່ຈໍາເປັນແມ່ນຖືກກໍານົດ
ພາຍໃນ -T ແລະ -Z ທາງເລືອກ. ຈໍານວນຂອງອໍານາດໃນການຂະຫຍາຍຕົວ Parker ແມ່ນຈໍາກັດ
1. ເບິ່ງຕົວຢ່າງເພີ່ມເຕີມລົງລຸ່ມ.

-Tte/rl/rm/rw[+m]
ຄິດ​ໄລ່​ການ​ຊົດ​ເຊີຍ isostatic ຈາກ​ການ​ໂຫຼດ topography (ໄຟລ​໌​ຕາ​ຂ່າຍ​ໄຟ​ຟ້າ​ການ​ປ້ອນ​ຂໍ້​ມູນ​) ໃນ​
ແຜ່ນ elastic ຂອງຄວາມຫນາ te. ນອກຈາກນີ້ຍັງເພີ່ມຄວາມຫນາແຫນ້ນສໍາລັບການໂຫຼດ, mantle, ແລະນ້ໍາໃນ
ໜ່ວຍ SI. ໃຫ້ຄວາມເລິກສະເລ່ຍຂອງ mantle ຜ່ານ -Z. ຖ້າຄວາມຫນາ elastic ແມ່ນ> 1e10 ມັນ
ຈະ​ໄດ້​ຮັບ​ການ​ຕີ​ລາ​ຄາ​ເປັນ flexural rigidity (ໂດຍ​ໃນ​ຕອນ​ຕົ້ນ​ມັນ​ແມ່ນ​ຄໍາ​ນວນ​ຈາກ​ te ແລະ
ໂມດູລອ່ອນ). ທາງເລືອກ, ຕື່ມໃສ່ +m ຂຽນຕາຂ່າຍໄຟຟ້າທີ່ມີທ່າແຮງທາງດ້ານພູມສາດຂອງ Moho
ຜົນກະທົບ (ເບິ່ງ -F) ຈາກຕົວແບບທີ່ເລືອກໂດຍ -Tທີ່ຢູ່ ຖ້າຫາກວ່າ te = 0 ຫຼັງຈາກນັ້ນການຕອບສະຫນອງ Airy ແມ່ນ
ກັບຄືນມາ. -T+ມ ກໍານົດ implicitly -N+ຊ

-Wwd ກຳນົດຄວາມເລິກຂອງນ້ຳ (ຫຼືຄວາມສູງການສັງເກດ) ທຽບກັບພູມສັນຖານ [0]. ຕື່ມໃສ່ k to
ຊີ້ km.

-Zzm[zl]
Moho [ແລະບວມ] ຄວາມເລິກການຊົດເຊີຍສະເລ່ຍ. ສໍາລັບຮູບແບບ "ໂຫຼດຈາກເທິງ" ທ່ານ
ພຽງແຕ່ມີໃຫ້ zm, ແຕ່ສໍາລັບ "ໂຫລດຈາກຂ້າງລຸ່ມນີ້" ຢ່າລືມ zl.

-V[ລະດັບ] (ຫຼາຍ ... )
ເລືອກລະດັບ verbosity [c].

-fg ຕາຂ່າຍທາງພູມສາດ (ຂະໜາດຂອງເສັ້ນແວງ, ເສັ້ນຂະໜານ) ຈະຖືກປ່ຽນເປັນແມັດ
ໂດຍຜ່ານການປະມານ "ໂລກຮາບພຽງ" ໂດຍໃຊ້ຕົວກໍານົດການ ellipsoid ໃນປັດຈຸບັນ.

-^ or ພຽງແຕ່ -
ພິມຂໍ້ຄວາມສັ້ນກ່ຽວກັບ syntax ຂອງຄໍາສັ່ງ, ຫຼັງຈາກນັ້ນອອກ (ຫມາຍເຫດ: ໃນ Windows
ໃຊ້ພຽງແຕ່ -).

-+ or ພຽງແຕ່ +
ພິມຂໍ້ຄວາມການນໍາໃຊ້ຢ່າງກວ້າງຂວາງ (ການຊ່ວຍເຫຼືອ), ລວມທັງຄໍາອະທິບາຍຂອງໃດໆ
ທາງເລືອກໂມດູນສະເພາະ (ແຕ່ບໍ່ແມ່ນຕົວເລືອກທົ່ວໄປ GMT), ຫຼັງຈາກນັ້ນອອກ.

-? or no ກະທູ້ທີ່
ພິມຂໍ້ຄວາມການນໍາໃຊ້ທີ່ສົມບູນ, ລວມທັງຄໍາອະທິບາຍທາງເລືອກ, ຫຼັງຈາກນັ້ນ
ອອກ.

- ການປ່ຽນແປງ
ພິມສະບັບ GMT ແລະອອກ.

--show-datadir
ພິມເສັ້ນທາງເຕັມໄປຫາ GMT share directory ແລະອອກ.

GRID ເອກະສານ ຮູບແບບ

ຕາມຄ່າເລີ່ມຕົ້ນ GMT ຂຽນອອກຕາຂ່າຍເປັນຄວາມແມ່ນຍໍາອັນດຽວລອຍຢູ່ໃນ netCDF-ການຮ້ອງຮຽນ COARDS
ຮູບແບບໄຟລ໌. ຢ່າງໃດກໍຕາມ, GMT ສາມາດຜະລິດໄຟລ໌ຕາຂ່າຍໄຟຟ້າຢູ່ໃນຕາຂ່າຍໄຟຟ້າອື່ນໆທີ່ໃຊ້ທົ່ວໄປຫຼາຍ
ຮູບ​ແບບ​ໄຟລ​໌​ແລະ​ຍັງ​ສ້າງ​ຄວາມ​ສະ​ດວກ​ທີ່​ເອີ້ນ​ວ່າ "ການ​ຫຸ້ມ​ຫໍ່​" ຕາ​ຂ່າຍ​ໄຟ​ຟ້າ​, ການ​ຂຽນ​ອອກ​ຈຸດ​ທີ່​ເລື່ອນ​ໄດ້​
ຂໍ້ມູນເປັນຈໍານວນ 1- ຫຼື 2-byte. ເພື່ອກໍານົດຄວາມແມ່ນຍໍາ, ຂະຫນາດແລະການຊົດເຊີຍ, ຜູ້ໃຊ້ຄວນ
ເພີ່ມຄໍາຕໍ່ທ້າຍ =id[/ຂະຫນາດ/ຊົດເຊີຍ[/nan]], ຢູ່ໃສ id ແມ່ນຕົວລະບຸສອງຕົວອັກສອນຂອງຕາຂ່າຍໄຟຟ້າ
ປະເພດແລະຄວາມແມ່ນຍໍາ, ແລະ ຂະຫນາດ ແລະ ຊົດເຊີຍ ແມ່ນປັດໄຈຂະຫນາດທາງເລືອກແລະຊົດເຊີຍທີ່ຈະເປັນ
ນຳໃຊ້ກັບຄ່າຕາຂ່າຍທັງໝົດ, ແລະ nan ແມ່ນຄ່າທີ່ໃຊ້ເພື່ອຊີ້ບອກຂໍ້ມູນທີ່ຂາດຫາຍໄປ. ໃນ​ກໍ​ລະ​ນີ
ສອງລັກສະນະ id ບໍ່ໄດ້ສະຫນອງໃຫ້, ຄືກັບໃນ =/ຂະຫນາດ ກວ່າຫນຶ່ງ id=nf ແມ່ນສົມມຸດ. ເມື່ອ​ໃດ​
ຕາ​ຂ່າຍ​ໄຟ​ຟ້າ​ການ​ອ່ານ​, ຮູບ​ແບບ​ໂດຍ​ທົ່ວ​ໄປ​ໄດ້​ຮັບ​ຮູ້​ອັດ​ຕະ​ໂນ​ມັດ​. ຖ້າບໍ່ແມ່ນ, ຄຳຕໍ່ທ້າຍອັນດຽວກັນ
ສາມາດຖືກເພີ່ມໃສ່ຊື່ໄຟລ໌ຕາຕະລາງການປ້ອນຂໍ້ມູນ. ເບິ່ງ grdconvert ແລະພາກສ່ວນຕາຕະລາງໄຟລ໌ຮູບແບບຂອງ
GMT ເອກະສານອ້າງອີງດ້ານວິຊາການແລະ Cookbook ສໍາລັບຂໍ້ມູນເພີ່ມເຕີມ.

ເມື່ອອ່ານໄຟລ໌ netCDF ທີ່ປະກອບດ້ວຍຫຼາຍຕາຂ່າຍ, GMT ຈະອ່ານ, ໂດຍຄ່າເລີ່ມຕົ້ນ, the
ຕາໜ່າງ 2 ມິຕິລະດັບທຳອິດທີ່ສາມາດຊອກຫາໄດ້ໃນໄຟລ໌ນັ້ນ. ເພື່ອ coax GMT ໃນການອ່ານອື່ນ
ຕົວແປຫຼາຍມິຕິໃນໄຟລ໌ຕາຂ່າຍ, ຕໍ່ທ້າຍ ?varname ກັບຊື່ໄຟລ໌, ບ່ອນທີ່
varname ແມ່ນຊື່ຂອງຕົວແປ. ໃຫ້ສັງເກດວ່າທ່ານອາດຈະຈໍາເປັນຕ້ອງຫນີຄວາມຫມາຍພິເສດ
of ? ໃນໂຄງການ shell ຂອງທ່ານໂດຍການວາງ backslash ຢູ່ທາງຫນ້າຂອງມັນ, ຫຼືໂດຍການວາງ
ຊື່ໄຟລ໌ ແລະຄຳຕໍ່ທ້າຍລະຫວ່າງວົງຢືມ ຫຼືວົງຢືມຄູ່. ໄດ້ ?varname suffix ຍັງສາມາດຖືກນໍາໃຊ້
ສໍາລັບຕາຕະລາງຜົນຜະລິດເພື່ອລະບຸຊື່ຕົວແປທີ່ແຕກຕ່າງຈາກຄ່າເລີ່ມຕົ້ນ: "z". ເບິ່ງ
grdconvert ແລະຕົວແກ້ໄຂພາກສ່ວນ-for-CF ແລະຮູບແບບໄຟລ໌ຕາຕະລາງຂອງ GMT Technical
ເອກະສານອ້າງອີງແລະ Cookbook ສໍາລັບຂໍ້ມູນເພີ່ມເຕີມ, ໂດຍສະເພາະກ່ຽວກັບວິທີການອ່ານ splices ຂອງ 3-,
ຕາຂ່າຍ 4-, ຫຼື 5-ມິຕິ.

GRID DISTANCE UNITS

ຖ້າຕາຂ່າຍໄຟຟ້າບໍ່ມີແມັດເປັນຫົວໜ່ວຍລວງນອນ, ໃຫ້ຕື່ມໃສ່ +uຫນ່ວຍບໍລິການ ກັບໄຟລ໌ປ້ອນຂໍ້ມູນ
ຊື່ທີ່ຈະປ່ຽນຈາກຫົວໜ່ວຍທີ່ລະບຸໄປເປັນແມັດ. ຖ້າຕາຂ່າຍໄຟຟ້າຂອງເຈົ້າເປັນພູມສັນຖານ, ແປງ
ໄລຍະຫ່າງກັບແມັດໂດຍການສະຫນອງ -fg ແທນທີ່ຈະເປັນ.

ການພິຈາລະນາ

ຕາຂ່າຍ netCDF COARDS ຈະຖືກຮັບຮູ້ໂດຍອັດຕະໂນມັດເປັນພູມສັນຖານ. ສໍາລັບຕາຂ່າຍໄຟຟ້າອື່ນໆ
ຕາຂ່າຍທາງພູມສາດແມ່ນທ່ານຕ້ອງການປ່ຽນອົງສາເປັນແມັດ, ເລືອກ -fg. ຖ້າ​ຫາກ​ວ່າ​ຂໍ້​ມູນ​
ຢູ່ໃກ້ກັບເສົາທັງສອງ, ທ່ານຄວນພິຈາລະນາການວາງໄຟລ໌ຕາຂ່າຍໄຟຟ້າເປັນສີ່ຫລ່ຽມ
ລະບົບປະສານງານໂດຍໃຊ້ ໂຄງການ grd.

ແຜນການ flexURE

ການແກ້ໄຂ FFT ສໍາລັບການຍືດຫຍຸ່ນແຜ່ນ elastic ຮຽກຮ້ອງໃຫ້ມີຄວາມຫນາແຫນ້ນຂອງ infill ເພື່ອເທົ່າກັບການໂຫຼດ
ຄວາມ​ຫນາ​ແຫນ້ນ. ນີ້ແມ່ນປົກກະຕິແລ້ວພຽງແຕ່ເປັນຄວາມຈິງໂດຍກົງພາຍໃຕ້ການໂຫຼດ; ເກີນການໂຫຼດໄດ້
infill ມັກຈະເປັນຕະກອນທີ່ມີຄວາມຫນາແຫນ້ນຕ່ໍາຫຼືແມ້ກະທັ້ງນ້ໍາ (ຫຼືອາກາດ). Wessel [2001] ສະເໜີ
ການປະມານທີ່ອະນຸຍາດໃຫ້ສໍາລັບຄວາມຫນາແຫນ້ນຂອງ infill ທີ່ແຕກຕ່າງຈາກ
ຄວາມຫນາແຫນ້ນຂອງການໂຫຼດໃນຂະນະທີ່ຍັງອະນຸຍາດໃຫ້ແກ້ໄຂ FFT. ໂດຍພື້ນຖານແລ້ວ, ຄວາມຍືດຫຍຸ່ນຂອງແຜ່ນແມ່ນ
ການແກ້ໄຂສໍາລັບການນໍາໃຊ້ຄວາມຫນາແຫນ້ນຂອງ infill ເປັນຄວາມຫນາແຫນ້ນຂອງການໂຫຼດປະສິດທິພາບແຕ່ຄວາມກວ້າງຂວາງແມ່ນ.
ປັບ​ໂດຍ​ປັດ​ໄຈ​ຫນຶ່ງ​ A = sqrt ((rm - ri)/(rm - rl)), ເຊິ່ງເປັນຄວາມແຕກຕ່າງທາງທິດສະດີ
ໃນຄວາມກວ້າງໃຫຍ່ໄພສານເນື່ອງຈາກການໂຫຼດຈຸດໂດຍໃຊ້ສອງຄວາມຫນາແຫນ້ນຂອງການໂຫຼດທີ່ແຕກຕ່າງກັນ. ໄດ້
approximation ແມ່ນດີຫຼາຍແຕ່ breaks ລົງສໍາລັບການໂຫຼດຂະຫນາດໃຫຍ່ຢູ່ໃນແຜ່ນທີ່ອ່ອນແອ, fairy
ສະຖານະການທີ່ບໍ່ທໍາມະດາ.

ຕົວຢ່າງ

ເພື່ອຄິດໄລ່ຜົນກະທົບຂອງຊັ້ນນ້ໍາຂ້າງເທິງ bathymetry bat.grd ໂດຍໃຊ້ 2700 ແລະ 1035.
ສໍາລັບຄວາມຫນາແຫນ້ນຂອງ crust ແລະນ້ໍາແລະຂຽນຜົນໄດ້ຮັບໃນ water_g.grd (computing up
ກັບພະລັງງານທີ່ສີ່ຂອງ bathymetry ໃນການຂະຫຍາຍຕົວ Parker):

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

ຕອນນີ້ຫັກມັນອອກເປັນ faa.grd ຜິດປົກກະຕິຂອງເຈົ້າ ແລະເຈົ້າຈະໄດ້ຮັບຄວາມຜິດກະຕິ Bouguer. ເຈົ້າ
ອາດຈະສົງໄສວ່າເປັນຫຍັງພວກເຮົາລົບແລະບໍ່ເພີ່ມ. ຫຼັງຈາກທີ່ທັງຫມົດ Bouguer ຜິດລັກທໍາທ່າ
ເພື່ອແກ້ໄຂການຂາດມະຫາຊົນທີ່ນໍາສະເຫນີໂດຍຊັ້ນນ້ໍາ, ດັ່ງນັ້ນພວກເຮົາຄວນເພີ່ມເພາະວ່າ
ນ້ໍາມີຄວາມຫນາແຫນ້ນຫນ້ອຍກວ່າກ້ອນຫີນຂ້າງລຸ່ມນີ້. ຄໍາຕອບແມ່ນຂຶ້ນກັບວິທີການຜົນກະທົບຂອງແຮງໂນ້ມຖ່ວງ
ຄິດໄລ່ໂດຍວິທີການຂອງ Parker ແລະພາກປະຕິບັດຂອງການນໍາໃຊ້ FFT.

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

ຕ້ອງການ MBA ຜິດປົກກະຕິ? ດີຄິດໄລ່ການປະກອບສ່ວນຂອງ crust mantle ແລະເພີ່ມມັນໃສ່
ຄວາມຜິດກະຕິຂອງທະເລລຸ່ມ. ສົມມຸດວ່າເປືອກເປືອກຫນາ 6 ກິໂລແມັດຂອງຄວາມຫນາແຫນ້ນ 2700 ແລະ mantle ກັບ 3300.
ຄວາມຫນາແຫນ້ນພວກເຮົາສາມາດເຮັດຊ້ໍາຄໍາສັ່ງທີ່ໃຊ້ໃນການຄິດໄລ່ຄວາມຜິດປົກກະຕິຂອງຊັ້ນນ້ໍາ, ໂດຍໃຊ້ 600
(3300 - 2700) ເປັນຄວາມຫນາແຫນ້ນຂອງກົງກັນຂ້າມ. ແຕ່ດຽວນີ້ພວກເຮົາມີບັນຫາເພາະວ່າພວກເຮົາຕ້ອງການຮູ້
ຄວາມເລິກຂອງ Moho ສະເລ່ຍ. ນັ້ນແມ່ນເວລາທີ່ຂະຫນາດ / ຊົດເຊີຍທີ່ສາມາດຖືກຕໍ່ທ້າຍກັບຊື່ຂອງຕາຂ່າຍໄຟຟ້າ
ມາຢູ່ໃນມື. ໃຫ້ສັງເກດວ່າພວກເຮົາບໍ່ຈໍາເປັນຕ້ອງເຮັດແນວນັ້ນກ່ອນເພາະວ່າຄວາມເລິກຂອງນ້ໍາແມ່ນ
ຄິດ​ໄລ່​ໂດຍ​ກົງ​ຈາກ​ຂໍ້​ມູນ (ສັງ​ເກດ​ເຫັນ​ຍັງ​ສັນ​ຍານ​ທາງ​ລົບ​ຂອງ​ການ​ຊົດ​ເຊີຍ​ເນື່ອງ​ຈາກ​ຄວາມ​ຈິງ​
ທີ່ z ເປັນ​ທາງ​ບວກ​)​:

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

ດຽວນີ້, ເອົາມັນອອກເປັນຄວາມຜິດກະຕິຂອງທະເລລຸ່ມເພື່ອໃຫ້ໄດ້ MBA ຜິດປົກກະຕິ. ນັ້ນ​ແມ່ນ:

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

ເພື່ອຄິດໄລ່ຜົນກະທົບແຮງໂນ້ມຖ່ວງຂອງ Moho ຂອງແຜ່ນ elastic bat.grd ດ້ວຍ Te = 7 ກິໂລແມັດ, ຄວາມຫນາແຫນ້ນຂອງ
2700, ຫຼາຍກວ່າ mantle ຂອງຄວາມຫນາແຫນ້ນ 3300, ໃນຄວາມເລິກສະເລ່ຍຂອງ 9 ກິໂລແມັດ.

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

ຖ້າທ່ານເພີ່ມຜົນກະທົບຂອງທະເລລຸ່ມແລະ Moho, ທ່ານຈະໄດ້ຮັບຄໍາຕອບຕໍ່ແຮງໂນ້ມຖ່ວງຢ່າງເຕັມທີ່
ຂອງຕົວແບບ isostatic ຂອງທ່ານ. ພວກເຮົາຈະໃຊ້ທີ່ນີ້ພຽງແຕ່ໄລຍະທໍາອິດໃນການຂະຫຍາຍຕົວ 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 elastic.grd ADD = model.grd

ຜົນໄດ້ຮັບດຽວກັນສາມາດໄດ້ຮັບໂດຍກົງໂດຍຄໍາສັ່ງຕໍ່ໄປ. ຢ່າງໃດກໍຕາມ, ຈ່າຍເອົາໃຈໃສ່
ຕໍ່ໄປນີ້. ຂ້າ​ພະ​ເຈົ້າ​ຍັງ​ບໍ່​ທັນ​ຮູ້​ວ່າ​ມັນ​ເປັນ​ຍ້ອນ​ການ bug ຫຼື​ເນື່ອງ​ຈາກ​ການ​ຈໍາ​ກັດ​ບາງ​ຢ່າງ​, ແຕ່​ວ່າ​
ຄວາມຈິງແມ່ນວ່າຄໍາສັ່ງຕໍ່ໄປນີ້ແລະກ່ອນຫນ້າພຽງແຕ່ໃຫ້ຜົນໄດ້ຮັບດຽວກັນຖ້າ -E1
ຖືກນໍາໃຊ້. ສໍາລັບອໍານາດທີ່ສູງຂຶ້ນຂອງ bathymetry ໃນການຂະຫຍາຍຕົວ Parker, ພຽງແຕ່ຕົວຢ່າງຂ້າງເທິງ
seams ເພື່ອໃຫ້ຜົນໄດ້ຮັບທີ່ຖືກຕ້ອງ.

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

ແລະສິ່ງທີ່ຈະເປັນຜິດປົກກະຕິ geoid ຜະລິດໂດຍການໂຫຼດຢູ່ທີ່ຄວາມເລິກ 50 ກິໂລແມັດ, ຂ້າງລຸ່ມນີ້ພາກພື້ນ
ການ​ອາບ​ນ້ຳ​ແມ່ນ​ໃຫ້​ໂດຍ bat.grd, Moho ຢູ່​ທີ່​ຄວາມ​ເລິກ 9 ກິ​ໂລ​ແມັດ ແລະ​ຄວາມ​ໜາ​ແໜ້ນ​ຄື​ກັນ.
ກ່ອນ?

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

ເພື່ອຄິດໄລ່ການເຂົ້າລະຫວ່າງຫ້ອງອາບນໍ້າ topo.grd ແລະ faa.grd ຜິດປົກກະຕິ
ຕາຂ່າຍໄຟຟ້າໃຊ້ຮູບແບບແຜ່ນ elastic ຂອງ crust ຄວາມຫນາຂອງ 6 ກິໂລແມັດຫມາຍຄວາມວ່າປະສິດທິພາບ 10 ກິໂລແມັດ
ຄວາມຫນາຂອງ elastic ໃນພາກພື້ນຂອງ 3 ກິໂລແມັດຫມາຍຄວາມວ່າຄວາມເລິກນ້ໍາ:

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

ເພື່ອຄິດໄລ່ການເຂົ້າລະຫວ່າງ topo.grd bathymetry ແລະ geoid.grd ຕາຂ່າຍໄຟຟ້າ geoid ກັບ
ຮູບແບບ "Loading from below" (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

ຂໍ້ມູນອ້າງອິງ

Luis, JF ແລະ MC Neves. 2006, ການຊົດເຊີຍ isostatic ຂອງພູພຽງ Azores: ເປັນ 3D
ການ​ເຂົ້າ​ມາ​ແລະ​ການ​ວິ​ເຄາະ​ຄວາມ​ສອດ​ຄ່ອງ​. J. Volc ຄວາມຮ້ອນໃຕ້ດິນ Res. ເຫຼັ້ມທີ 156, ສະບັບທີ 1-2, ໜ້າ
10-22, http://dx.doi.org/10.1016/j.jvolgeores.2006.03.010 Parker, RL, 1972, ໄວ
ການຄິດໄລ່ຄວາມຜິດປົກກະຕິທີ່ອາດຈະເກີດຂື້ນ, Geophys. J., 31, 447-455. ເວສເຊວ. P., 2001, ທົ່ວໂລກ
ການແຜ່ກະຈາຍຂອງພູເຂົາທະເລທີ່ສົມມຸດຕິຖານມາຈາກ Geosat/ERS-1 altimetry, J. Geophys. ສະບັບ.,
106(B9), 19,431-19,441, http://dx.doi.org/10.1029/2000JB000083

ໃຊ້ gravfftgmt ອອນໄລນ໌ໂດຍໃຊ້ບໍລິການ onworks.net