gravfftgmt - Online sa Cloud

PROGRAMA:

NAME

gravfft - Compute gravitational attraction ng 3-D surface sa wavenumber (o
dalas) domain

SINOPSIS

gravfft ingrid [ ingrid2 ] outfile [ n/wavelength/mean_depth/tbw ] [ Density|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][+tlapad][+w[suffix]][+z[p]] [ ] [ te/rl/rm/rw[+m] ] [ [antas] ] [ wd] [
zm[zl] ] [ -fg ]

tandaan: Walang puwang ang pinapayagan sa pagitan ng flag ng opsyon at ng mga nauugnay na argumento.

DESCRIPTION

gravfft maaaring gamitin sa tatlong pangunahing mga mode. Mode 1: I-compute lang ang geopotential dahil sa
ang ibabaw na ibinigay sa topo.grd file. Nangangailangan ng kaibahan ng density (-D) at posibleng a
iba't ibang antas ng pagmamasid (-W). Aabutin nito ang 2-D forward FFT ng grid at gamitin
ang buong pamamaraan ni Parker hanggang sa mga napiling termino. Mode 2: Compute ng geopotential
tugon dahil sa flexure ng topography file. Aabutin ang 2-D forward FFT ng
grid at gamitin ang buong pamamaraan ng Parker na inilapat sa napiling modelong isostatic. Ang
ang mga available na modelo ay ang "loading from top", o elastic plate model, at ang "loading from
sa ibaba" na tumutukoy sa tugon ng plato sa isang sub-surface load (angkop para sa mainit
spot modeling - kung naniniwala ka sa kanila). Sa parehong mga kaso, ang mga parameter ng modelo ay nakatakda sa -T
at -Z mga pagpipilian. Mode 3: kalkulahin ang pagpasok o pagkakaugnay sa pagitan ng dalawang grids. Ang output
ay ang average sa radial na direksyon. Opsyonal, ang modelong pagpasok ay maaari ding
kalkulado. Ang mga pahalang na sukat ng mga grdfile ay ipinapalagay na nasa metro.
Maaaring gamitin ang mga geograpikal na grid sa pamamagitan ng pagtukoy sa -fg opsyon na nagsusukat ng digri hanggang metro.
Kung mayroon kang mga grid na may mga sukat sa km, maaari mo itong baguhin sa metro gamit grdedit or
sukatin ang output gamit ang grdmath. Dahil sa bilang ng mga pagpipilian na inaalok ng program na ito, ay
mahirap sabihin kung ano ang mga opsyon at kung ano ang mga kinakailangang argumento. Depende kung ano
ang iyong ginagawa; tingnan ang mga halimbawa para sa karagdagang gabay.

KAILANGAN MGA PANGANGATWIRANG

ingrid 2-D binary grid file na paganahin. (Tingnan ang GRID FILE FORMATS sa ibaba). Para sa
cross-spectral operations, ibigay din ang pangalawang grid file ingrd2.

-Goutfile
Tukuyin ang pangalan ng output grid file o ang 1-D spectrum table (tingnan ang -E). (Tingnan
MGA FORMAT NG GRID FILE sa ibaba).

OPSYONAL MGA PANGANGATWIRANG

-Cn/wavelength/mean_depth/tbw
I-compute lang ang theoretical admittance curves ng napiling modelo at exit. n
at haba ng daluyong ay ginagamit upang kalkulahin (n * wavelength) ang kabuuang haba ng profile sa
metro. mean_depth ay ang ibig sabihin ng lalim ng tubig. Idagdag ang mga dataflag (isa o dalawa) ng tbw in
anumang order. t = gamitin ang "mula sa itaas" na modelo, b = gamitin ang "mula sa ibaba" na modelo. Opsyonal na tukuyin
w isulat ang wavelength sa halip na frequency.

-DDensity|rhogrid
Nagtatakda ng contrast ng density sa buong ibabaw. Ginamit, halimbawa, upang kalkulahin ang gravity
atraksyon ng layer ng tubig na maaaring pagsamahin sa ibang pagkakataon sa free-air anomalya
para makuha ang Bouguer anomalya. Sa kasong ito, huwag gamitin -T. Implicitly din itong nagtatakda
-N+h. Bilang kahalili, tukuyin ang isang co-registered grid na may density contrasts kung a
kailangan ang variable density contrast.

-En_terms
Bilang ng mga terminong ginamit sa pagpapalawak ng Parker (ang limitasyon ay 10, kung hindi, ang mga termino ay depende sa
n ay sasabog ang programa) [Default = 3]

-F[f[+]|g|v|n|e]
Tukuyin ang gustong geopotential field: compute geoid sa halip na gravity
f = Mga anomalya sa free-air (mGal) [Default]. Idugtong + upang idagdag sa slab na ipinahiwatig
kapag inaalis ang mean value mula sa topograpiya. Nangangailangan ito ng zero topography
ang ibig sabihin ay walang anomalya sa masa.

g = Geoid anomalya (m).

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

e = East deflections ng vertical (micro-radian).

n = North deflections ng vertical (micro-radian).

-Iw|b|c|t |k
paggamit ingrd2 at ingrd1 (isang grid na may topograpiya/batymetry) upang tantiyahin
admittance|coherence at isulat ito sa stdout (-G hindi pinansin kung itinakda). Ang grid na ito ay dapat
naglalaman ng gravity o geoid para sa parehong rehiyon ng ingrd1. Default na pag-compute
pagpasok. Ang output ay naglalaman ng 3 o 4 na column. Dalas (haba ng daluyong), pagpasok
(pagkakaugnay) isang sigma error bar at, opsyonal, isang teoretikal na pagpasok. Idugtong
dataflags (isa hanggang tatlo) mula sa w|b|c|t. w nagsusulat ng wavelength sa halip na wavenumber,
k pumipili ng km para sa wavelength unit [m], c kinukuwenta ang pagkakaugnay-ugnay sa halip na pagpasok, b
nagsusulat ng ikaapat na hanay na may "naglo-load mula sa ibaba" na teoretikal na pagpasok, at t
nagsusulat ng ika-apat na hanay na may "elastic plate" theoretical admittance.

-N[f|q|s|nx/ny][+a|[+d|h|l][+e|n|m][+tlapad][+w[suffix]][+z[p]]
Pumili o magtanong tungkol sa angkop na mga sukat ng grid para sa FFT at itakda ang opsyonal
mga parameter. Kontrolin ang dimensyon ng FFT:
-Nf pipilitin ang FFT na gamitin ang aktwal na mga sukat ng data.

-Nq ay magtatanong tungkol sa mas angkop na mga sukat, iulat ang mga iyon, pagkatapos ay magpapatuloy.

-Ns ay magpapakita ng listahan ng mga opsyonal na dimensyon, pagkatapos ay lalabas.

-Nnx/ny gagawa ng FFT sa laki ng array nx/ny (dapat na >= laki ng grid file). Default
pinipili ang mga dimensyon >= data na nag-o-optimize ng bilis at katumpakan ng FFT. Kung FFT
mga dimensyon > mga sukat ng grid file, ang data ay pinahaba at naka-tape sa zero.

Kontrolin ang pag-detrending ng data: Idagdag ang mga modifier para sa pag-alis ng linear na trend:
+d: Detrend data, ibig sabihin, alisin ang pinakaangkop na linear trend [Default].

+a: Alisin lamang ang ibig sabihin ng halaga.

+h: Alisin lamang ang mid value, ibig sabihin, 0.5 * (max + min).

+l: Iwanan ang data.

Kontrolin ang extension at pag-taping ng data: Gumamit ng mga modifier para kontrolin kung paano ang extension
at ang tapering ay isasagawa:
+e pinapalawak ang grid sa pamamagitan ng pagpapataw ng edge-point symmetry [Default],

+m pinapalawak ang grid sa pamamagitan ng paglalagay ng edge mirror symmetry

+n in-off ang extension ng data.

Ginagawa ang tapering mula sa gilid ng data hanggang sa gilid ng grid ng FFT [100%]. Baguhin
ang porsyentong ito sa pamamagitan ng +tlapad. Kailan +n ay may bisa, ang tapering ay inilapat
sa halip sa mga margin ng data dahil walang available na extension [0%].

Kontrolin ang pagsulat ng mga pansamantalang resulta: Para sa detalyadong pagsisiyasat maaari mong isulat ang
intermediate grid na ipinapasa sa forward FFT; ito ay malamang na naging
detrended, pinalawak ng point-symmetry sa lahat ng mga gilid, at tapered. Idugtong
+w[suffix] kung saan malilikha ang pangalan ng (mga) file ng output (ibig sabihin, ingrid_prefix.ext)
[tapered], saan ext ay ang iyong file extension. Sa wakas, maaari mong i-save ang kumplikadong grid
ginawa ng forward FFT sa pamamagitan ng pagdaragdag +z. Bilang default, isinusulat namin ang tunay at
haka-haka na mga bahagi sa ingrid_totoo.ext at ingrid_imag.ext. Idugtong p upang i-save
sa halip ang polar form ng magnitude at phase sa mga file ingrid_mag.ext at
ingrid_phase.ext.

-Q Nagsusulat ng grid na may flexural topography (na may z positive up) na ang average
ay itinakda ni -Zzm at mga parameter ng modelo sa pamamagitan ng -T (at output ni -G). Iyon ay ang
"gravimetric Moho". -Q implicitly sets -N+h

-S Kinuwenta ang hinulaang gravity o geoid grid dahil sa isang subplate load na ginawa ng
kasalukuyang bathymetry at ang teoretikal na modelo. Ang mga kinakailangang parameter ay nakatakda
sa loob ng -T at -Z mga pagpipilian. Ang bilang ng mga kapangyarihan sa pagpapalawak ng Parker ay limitado sa
1. Tingnan ang isang halimbawa sa ibaba.

-Tte/rl/rm/rw[+m]
Kalkulahin ang isostatic compensation mula sa topography load (input grid file) sa isang
nababanat na plato ng kapal te. Idagdag din ang mga densidad para sa load, mantle, at tubig
Mga yunit ng SI. Ibigay ang average na lalim ng mantle sa pamamagitan ng -Z. Kung ang elastic na kapal ay > 1e10 ito
ay bibigyang-kahulugan bilang ang flexural rigidity (bilang default ay nakalkula ito mula sa te at
Young modulus). Opsyonal, idagdag +m para magsulat ng grid gamit ang geopotential ng Moho
epekto (tingnan -F) mula sa modelong pinili ni -T. Kung te = 0 kung gayon ang Airy na tugon ay
bumalik. -T+m implicitly sets -N+h

-Wwd Itakda ang lalim ng tubig (o taas ng pagmamasid) na may kaugnayan sa topograpiya [0]. Idugtong k sa
ipahiwatig ang km.

-Zzm[zl]
Moho [at swell] ang average na lalim ng compensation. Para sa "load mula sa itaas" na modelo mo
kailangan lang magbigay zm, ngunit para sa "naglo-load mula sa ibaba" huwag kalimutan zl.

-V[antas] (higit pa ...)
Piliin ang antas ng verbosity [c].

-fg Ang mga geographic na grid (mga sukat ng longitude, latitude) ay gagawing metro
sa pamamagitan ng pagtatantya ng "Flat Earth" gamit ang kasalukuyang mga parameter ng ellipsoid.

-^ or m -
Mag-print ng maikling mensahe tungkol sa syntax ng command, pagkatapos ay lumabas (TANDAAN: sa Windows
gamitin lang -).

-+ or m +
Mag-print ng malawak na paggamit (tulong) na mensahe, kasama ang paliwanag ng anuman
opsyong tukoy sa module (ngunit hindi ang mga karaniwang opsyon sa GMT), pagkatapos ay lalabas.

-? or hindi argumento
Mag-print ng kumpletong mensahe ng paggamit (tulong), kasama ang paliwanag ng mga opsyon, pagkatapos
labasan.

--bersyon
I-print ang bersyon ng GMT at lumabas.

--show-datadir
I-print ang buong path sa GMT share directory at lumabas.

Parilya FILE FORMATS

Bilang default, isinusulat ng GMT ang grid habang lumulutang ang solong katumpakan sa isang netCDF na reklamo ng COARDS
format ng file. Gayunpaman, nakakagawa ang GMT ng mga grid file sa maraming iba pang karaniwang ginagamit na grid
mga format ng file at pinapadali din ang tinatawag na "pag-iimpake" ng mga grid, na nagsusulat ng lumulutang na punto
data bilang 1- o 2-byte na integer. Upang tukuyin ang katumpakan, sukat at offset, dapat ang user
idagdag ang panlapi =id[/sukatan/ginalaw[/nan]], saan id ay isang dalawang-titik na identifier ng grid
uri at katumpakan, at sukatan at ginalaw ay opsyonal na scale factor at offset na
inilapat sa lahat ng mga halaga ng grid, at nan ay ang halaga na ginagamit upang ipahiwatig ang nawawalang data. Kung sakali
ang dalawang karakter id ay hindi ibinigay, tulad ng sa =/sukatan kaysa sa isang id=nf ay ipinapalagay. Kailan
pagbabasa ng mga grids, ang format ay karaniwang awtomatikong kinikilala. Kung hindi, ang parehong suffix
maaaring idagdag sa mga pangalan ng file ng grid ng input. Tingnan mo grdconvert at Section grid-file-format ng
GMT Technical Reference at Cookbook para sa higit pang impormasyon.

Kapag nagbabasa ng netCDF file na naglalaman ng maraming grids, babasahin ng GMT, bilang default, ang
unang 2-dimensional na grid na mahahanap sa file na iyon. Upang hikayatin ang GMT na magbasa ng isa pa
multi-dimensional na variable sa grid file, idagdag ?varname sa pangalan ng file, kung saan
varname ay ang pangalan ng variable. Tandaan na maaaring kailanganin mong takasan ang espesyal na kahulugan
of ? sa iyong shell program sa pamamagitan ng paglalagay ng backslash sa harap nito, o sa pamamagitan ng paglalagay ng
filename at suffix sa pagitan ng mga quote o double quote. Ang ?varname maaari ding gamitin ang panlapi
para sa mga output grid upang tumukoy ng variable na pangalan na iba sa default: "z". Tingnan mo
grdconvert at Mga Seksyon modifier-para-CF at grid-file-format ng GMT Technical
Sanggunian at Cookbook para sa higit pang impormasyon, partikular sa kung paano magbasa ng mga splice ng 3-,
4-, o 5-dimensional na grid.

Parilya DISTANCE UNITS

Kung ang grid ay walang metro bilang horizontal unit, dugtungan +uyunit sa input file
pangalan na iko-convert mula sa tinukoy na unit patungo sa metro. Kung geographic ang iyong grid, mag-convert
mga distansya sa metro sa pamamagitan ng pagbibigay -fg sa halip.

KONSIDERASYON

Ang netCDF COARDS grids ay awtomatikong makikilala bilang geographic. Para sa iba pang mga grids
geographical grids kung gusto mong i-convert ang mga degree sa metro, piliin -fg. Kung ang datos
ay malapit sa alinmang poste, dapat mong isaalang-alang ang pag-project ng grid file sa isang hugis-parihaba
coordinate system gamit ang grdproject.

MAG-PLATE PAGBABAGO

Ang solusyon ng FFT sa elastic plate flexure ay nangangailangan ng infill density na katumbas ng load
densidad. Ito ay karaniwang totoo lamang nang direkta sa ilalim ng pagkarga; lampas sa load ang
Ang infill ay malamang na mga sediment na mas mababa ang density o kahit na tubig (o hangin). Iminungkahi ni Wessel [2001].
isang pagtatantya na nagbibigay-daan para sa pagtutukoy ng isang infill density na naiiba sa
densidad ng pag-load habang pinapayagan pa rin ang isang solusyon sa FFT. Talaga, ang plate flexure ay
nalutas para sa paggamit ng infill density bilang ang epektibong density ng pagkarga ngunit ang mga amplitude ay
inaayos ng isang salik A = sqrt ((rm - ri)/(rm - rl)), na siyang teoretikal na pagkakaiba
sa amplitude dahil sa isang point load gamit ang dalawang magkaibang load density. Ang
approximation ay napakahusay ngunit break down para sa malalaking load sa mahina plates, isang engkanto
hindi pangkaraniwang sitwasyon.

HALIMBAWA

Upang kalkulahin ang epekto ng layer ng tubig sa itaas ng paniki.grd bathymetry gamit ang 2700 at 1035
para sa densidad ng crust at tubig at pagsusulat ng resulta sa water_g.grd (computing up
sa ika-apat na kapangyarihan ng bathymetry sa pagpapalawak ng Parker):

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

Ngayon ay ibawas ito sa iyong free-air anomaly faa.grd at makukuha mo ang Bouguer anomaly. Ikaw
maaaring magtaka kung bakit tayo ay nagbabawas at hindi nagdadagdag. Kung tutuusin ay nagpapanggap ang anomalyang Bouguer
upang itama ang mass deficiency na ipinakita ng layer ng tubig, kaya dapat nating idagdag dahil
ang tubig ay hindi gaanong siksik kaysa sa mga bato sa ibaba. Ang sagot ay nakasalalay sa paraan ng mga epekto ng gravity
kinalkula sa pamamagitan ng pamamaraan ng Parker at mga praktikal na aspeto ng paggamit ng FFT.

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

Gusto ng MBA anomalya? Well kalkulahin ang crust mantle kontribusyon at idagdag ito sa
anomalya sa ilalim ng dagat. Ipagpalagay na ang isang 6 na km makapal na crust ng density ay 2700 at isang mantle na may 3300
density maaari naming ulitin ang utos na ginamit upang makalkula ang anomalya ng layer ng tubig, gamit ang 600
(3300 - 2700) bilang contrast ng density. Pero may problema tayo ngayon dahil kailangan nating malaman
ang ibig sabihin ng lalim ng Moho. Iyon ay kapag ang sukat/offset na maaaring idugtong sa pangalan ng grid
dumating sa kamay. Pansinin na hindi namin kailangang gawin iyon noon dahil ang ibig sabihin ay lalim ng tubig
direktang nakalkula mula sa data (pansinin din ang negatibong tanda ng offset dahil sa katotohanan
na z ay positibo):

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

Ngayon, ibawas ito sa sea-bottom anomalya upang makuha ang MBA anomalya. Yan ay:

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

Para kalkulahin ang Moho gravity effect ng isang elastic plate bat.grd na may Te = 7 km, density ng
2700, higit sa isang mantle na may density na 3300, sa average na lalim na 9 km

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

Kung idaragdag mo ngayon ang sea-bottom at ang mga epekto ng Moho, makukuha mo ang buong gravity response
ng iyong isostatic na modelo. Gagamitin namin dito ang unang termino sa pagpapalawak ng 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

Ang parehong resulta ay maaaring makuha nang direkta sa pamamagitan ng susunod na utos. Gayunpaman, MAGBIGAY NG PANSIN sa
ang mga sumusunod. Hindi ko pa alam kung ito ay dahil sa isang bug o dahil sa ilang limitasyon, ngunit
ang katotohanan ay ang mga sumusunod at ang naunang mga utos ay nagbibigay lamang ng parehong resulta kung -E1
Ginagamit. Para sa mas mataas na kapangyarihan ng bathymetry sa pagpapalawak ng Parker, tanging ang halimbawa sa itaas
tahi upang magbigay ng tamang resulta.

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

At ano ang magiging geoid anomaly na ginawa ng isang load sa lalim na 50 km, sa ibaba ng isang rehiyon
na ang bathymetry ay ibinibigay ng bat.grd, isang Moho sa lalim na 9 km at kapareho ng densidad ng
bago?

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

Upang kalkulahin ang pagpasok sa pagitan ng topo.grd bathymetry at faa.grd free-air anomaly
grid gamit ang elastic plate na modelo ng isang crust na 6 km ang ibig sabihin ng kapal na may 10 km epektibo
nababanat na kapal sa isang rehiyon na 3 km ang ibig sabihin ng lalim ng tubig:

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

Upang kalkulahin ang pagpasok sa pagitan ng topo.grd bathymetry at geoid.grd geoid grid na may
ang modelong "loading from below" (LFB) na kapareho ng nasa itaas at sub-surface load sa 40 km,
ngunit ipagpalagay na ngayon ang mga grids ay nasa geographic at gusto namin ng mga wavelength sa halip na dalas:

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

Upang kalkulahin ang gravity theoretical admittance ng isang LFB kasama ang isang 2000 km ang haba na profile gamit ang
ang parehong mga parameter tulad ng sa itaas

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

Mga sanggunian

Luis, JF at MC Neves. 2006, Ang isostatic compensation ng Azores Plateau: isang 3D
pagsusuri ng pagpasok at pagkakaugnay. J. Geothermal Volc. Res. Tomo 156, Mga Isyu 1-2, Mga Pahina
10-22, http://dx.doi.org/10.1016/j.jvolgeores.2006.03.010 Parker, RL, 1972, Ang mabilis
pagkalkula ng mga potensyal na anomalya, Geophys. J., 31, 447-455. Wessel. P., 2001, Pandaigdig
pamamahagi ng mga seamount na hinuha mula sa gridded Geosat/ERS-1 altimetry, J. Geophys. Res.,
106(B9), 19,431-19,441, http://dx.doi.org/10.1029/2000JB000083

Gamitin ang gravfftgmt online gamit ang mga serbisyo ng onworks.net