gravfftgmt - 클라우드 온라인

프로그램:

이름

gravfft - 파수에서 3차원 표면의 중력 인력을 계산합니다(또는
주파수) 도메인

개요

중력 잉그리드 [ 잉그리드 ] 아웃파일 [ n/파장/평균_깊이/tbw ] [ density|로그리드 ] [
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] ] [ -fg ]

참고 : 옵션 플래그와 관련 인수 사이에는 공백이 허용되지 않습니다.

기술

중력 세 가지 주요 모드로 사용할 수 있습니다. 모드 1: 단순히 다음으로 인한 지리 전위를 계산합니다.
topo.grd 파일에 제공된 표면. 밀도 대비가 필요합니다(-D) 그리고 아마도
다른 관찰 수준(-W). 그리드의 2-D 순방향 FFT를 사용하고
선택한 용어까지 전체 Parker의 방법. 모드 2: 지리 전위 계산
지형 파일의 굴곡으로 인한 응답. 그것은 2-D 순방향 FFT를 취할 것입니다
그리드를 선택하고 선택한 등압 모델에 적용된 전체 Parker 방법을 사용합니다. 그만큼
사용 가능한 모델은 "상단에서 하중" 또는 탄성판 모델과 "하향에서 하중"입니다.
표면 아래 하중에 대한 플레이트의 응답을 설명하는 "아래"(뜨거운
스폿 모델링 - 당신이 그들을 믿는다면). 두 경우 모두 모델 매개변수는 다음과 같이 설정됩니다. -T
과 -Z 옵션. 모드 3: 두 그리드 사이의 어드미턴스 또는 일관성을 계산합니다. 출력
방사형 방향의 평균입니다. 선택적으로 모델 어드미턴스는 다음과 같을 수도 있습니다.
계획된. grdfile의 수평 치수는 미터 단위로 가정합니다.
지리적 그리드는 다음을 지정하여 사용할 수 있습니다. -fg 각도를 미터로 조정하는 옵션입니다.
km 단위의 그리드가 있는 경우 다음을 사용하여 이를 미터로 변경할 수 있습니다. 감사 or
출력 크기 조정 수학. 이 프로그램이 제공하는 선택의 수를 고려할 때,
옵션이 무엇이고 필수 인수가 무엇인지 명시하기 어렵습니다. 그것은 무엇에 달려 있습니다
당신은 일을; 추가 지침은 예제를 참조하십시오.

필요한 인수

잉그리드 작업할 2D 이진 그리드 파일입니다. (아래 그리드 파일 형식 참조). 을 위한
교차 스펙트럼 작업, 두 번째 그리드 파일 제공 잉그라드2.

-G아웃파일
출력 그리드 파일 또는 1D 스펙트럼 테이블의 이름을 지정합니다(참조: -E). (만나다
아래의 그리드 파일 형식).

선택 사항 인수

-Cn/파장/평균_깊이/tbw
선택한 모델의 이론적 어드미턴스 곡선만 계산하고 종료합니다. n
과 파장 (n * 파장)에서 총 프로파일 길이를 계산하는 데 사용됩니다.
미터. 평균 깊이 평균 수심입니다. 다음의 데이터 플래그(하나 또는 둘)를 추가합니다. TBW in
어떤 주문. t = "위에서" 모델 사용, b = "아래에서" 모델을 사용합니다. 선택적으로 지정
w 주파수 대신 파장을 씁니다.

-Ddensity|로그리드
표면 전체의 밀도 대비를 설정합니다. 예를 들어 중력을 계산하는 데 사용됩니다.
나중에 자유 공기 이상과 결합될 수 있는 수층의 인력
Bouguer 변칙을 얻기 위해. 이 경우 사용하지 마십시오. -T. 또한 암시적으로 설정합니다.
-N+시. 또는 밀도 대비가 있는 공동 등록 그리드를 지정하십시오.
가변 밀도 콘트라스트가 필요합니다.

-En_terms
Parker 확장에 사용되는 용어 수(한도는 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
잉그라드2 과 잉그라드1 (지형/수심 측량이 포함된 그리드) 추정
어드미턴스|일관성을 확인하고 stdout(-G 설정된 경우 무시됨). 이 그리드는
동일한 지역에 대한 중력 또는 지오이드를 포함 잉그라드1. 기본 컴퓨팅
입장. 출력에는 3개 또는 4개의 열이 포함됩니다. 주파수(파장), 어드미턴스
(일관성) XNUMX 시그마 오차 막대 및 선택적으로 이론적 어드미턴스. 추가
데이터 플래그(XNUMX에서 XNUMX까지) w|b|c|t. w 파수 대신 파장을 씁니다.
k 파장 단위 [m]에 대해 km를 선택하고, c 어드미턴스 대신 일관성을 계산합니다. b
"아래에서 로드" 이론적인 어드미턴스로 네 번째 열을 작성하고 t
"탄성판"이론적 어드미턴스로 네 번째 칼럼을 씁니다.

-N[f|q|s|nx/뉴욕][+a|[+d|h|l][+e|n|m][+t폭][+우[접미사]][+지[p]]
FFT에 적합한 그리드 치수를 선택하거나 문의하고 옵션으로 설정
매개변수. FFT 차원 제어:
-Nf FFT가 데이터의 실제 차원을 사용하도록 합니다.

-Nq 더 적합한 치수에 대해 문의하고 보고한 다음 계속합니다.

-Ns 선택적 측정기준 목록을 표시한 다음 종료합니다.

-Nnx/뉴욕 배열 크기에 대해 FFT를 수행합니다. nx/뉴욕 (>= 그리드 파일 크기여야 함). 기본
FFT의 속도와 정확도를 최적화하는 차원 >= 데이터를 선택합니다. FFT인 경우
차원 > 그리드 파일 차원, 데이터가 확장되고 XNUMX으로 테이퍼됩니다.

데이터의 추세 제거 제어: 선형 추세를 제거하기 위한 수정자 추가:
+d: 데이터 추세 제거, 즉 가장 적합한 선형 추세 제거 [기본값].

+a: 평균값만 제거합니다.

+h: 중간 값, 즉 0.5 *(최대 + 최소)만 제거합니다.

+l: 데이터를 그대로 두십시오.

데이터 확장 및 테이퍼링 제어: 수정자를 사용하여 확장
테이퍼링을 수행해야 합니다.
+e 에지 포인트 대칭을 부과하여 그리드 확장 [기본값],

+m 에지 미러 대칭을 적용하여 그리드 확장

+n 데이터 확장을 끕니다.

데이터 에지에서 FFT 그리드 에지로 테이퍼링이 수행됩니다[100%]. 변화
이 비율을 통해 +t폭. 언제 +n 테이퍼링이 적용됩니다.
확장자를 사용할 수 없으므로 데이터 여백에 대신 적용됩니다[0%].

임시 결과 쓰기 제어: 자세한 조사를 위해 다음을 작성할 수 있습니다.
중간 그리드가 순방향 FFT로 전달되고 있습니다. 이것은 아마도
모든 모서리를 따라 점 대칭으로 확장되고 테이퍼됩니다. 추가
+w[접미사] 출력 파일 이름이 생성될 것입니다(즉, ingrid_prefix.ext)
[테이퍼드], 어디에 내선 파일 확장자입니다. 마지막으로 복잡한 그리드를 저장할 수 있습니다.
다음을 추가하여 순방향 FFT에 의해 생성 +z. 기본적으로 우리는 실제와
허수 성분 잉그리드_진짜.내선 과 잉그리드_이미지.내선. 추가 p 저장
대신 파일에 대한 크기 및 위상의 극성 형태 잉그리드_잡지.내선 과
잉그리드_단계.내선.

-Q 굴곡 지형(z 양수가 위로)이 있는 그리드를 작성합니다.
에 의해 설정되었습니다 -Zzm 및 모델 매개변수 -T (및 출력 -G). 그것은
"중량 모호". -Q 암시적으로 설정 -N+시

-S 에 의해 생성된 서브플레이트 하중으로 인해 예상 중력 또는 지오이드 그리드를 계산합니다.
현재 수심 측량 및 이론적 모델. 필요한 매개변수가 설정됩니다.
이내 -T 과 -Z 옵션. Parker 확장의 권한 수는 다음으로 제한됩니다.
1. 더 아래에 있는 예를 참조하십시오.

-T테/rl/rm/rw[+엠]
지형 하중(입력 그리드 파일)에서 등압 보상을 계산합니다.
두께의 탄성 판 te. 또한 하중, 맨틀 및 물에 대한 밀도를 추가합니다.
SI 단위. 평균 맨틀 깊이 제공 -Z. 탄성 두께가 > 1e10이면
굴곡 강성으로 해석됩니다(기본적으로 다음에서 계산됨). te 과
젊은 계수). 선택적으로 추가 +m Moho의 지오포텐셜로 그리드 작성하기
효과(참조 -F) 선택한 모델에서 -T. 면 te = 0이면 Airy 응답은
반환되었습니다. -티+엠 암시적으로 설정 -N+시

-Wwd 지형 [0]을 기준으로 수심(또는 관찰 높이)을 설정합니다. 추가 k 에
km를 표시합니다.

-Zzm[zl]
Moho [and swell] 평균 보상 깊이. "위에서 로드" 모델의 경우
제공해야만 zm, 그러나 "아래에서 로드"를 위해 잊지 마세요 zl.

-V[수평] (더 ~)
상세 수준 [c]를 선택합니다.

-fg 지리적 그리드(경도, 위도 치수)가 미터로 변환됩니다.
현재 타원체 매개변수를 사용하여 "평평한 지구" 근사치를 통해.

-^ or 다만 -
명령 구문에 대한 짧은 메시지를 인쇄한 다음 종료합니다(참고: Windows
그냥 사용 -).

-+ or 다만 +
설명을 포함하여 광범위한 사용(도움말) 메시지를 인쇄하십시오.
모듈별 옵션(GMT 공통 옵션 아님)을 선택한 다음 종료됩니다.

-? or 아니 인수
옵션 설명을 포함하여 전체 사용법(도움말) 메시지를 인쇄한 다음
출구.

--번역
GMT 버전을 인쇄하고 종료합니다.

--show-datadir
GMT 공유 디렉토리의 전체 경로를 인쇄하고 종료합니다.

GRID FILE 형식

기본적으로 GMT는 COARDS 호환 netCDF에서 단정밀도 부동 소수점으로 그리드를 작성합니다.
파일 형식. 그러나 GMT는 일반적으로 사용되는 다른 많은 그리드에서 그리드 파일을 생성할 수 있습니다.
파일 형식 및 또한 부동 소수점 작성, 그리드의 "패킹"을 용이하게 합니다.
데이터를 1바이트 또는 2바이트 정수로 표시합니다. 정밀도, 스케일 및 오프셋을 지정하려면 사용자가
접미사를 추가 =id[/규모/오프셋[/할머니]], 어디 id 그리드의 두 글자 식별자입니다.
유형 및 정밀도, 규모 과 오프셋 선택적 스케일 팩터 및 오프셋은 다음과 같습니다.
모든 그리드 값에 적용 할머니 누락된 데이터를 나타내는 데 사용되는 값입니다. 경우에
두 캐릭터 id 와 같이 제공되지 않습니다. =/규모 ~보다 id=nf 가정됩니다. 언제
그리드를 읽을 때 형식은 일반적으로 자동으로 인식됩니다. 그렇지 않은 경우 동일한 접미사
입력 그리드 파일 이름에 추가할 수 있습니다. 보다 grd변환 및 섹션 그리드 파일 형식의
자세한 내용은 GMT 기술 참조 및 요리 책을 참조하십시오.

여러 그리드가 포함된 netCDF 파일을 읽을 때 GMT는 기본적으로 다음을 읽습니다.
해당 파일에서 찾을 수 있는 첫 번째 2차원 그리드입니다. GMT가 다른 것을 읽도록 유도하기 위해
그리드 파일의 다차원 변수, 추가 ?이름 파일 이름에, 여기서
이름 변수의 이름입니다. 특별한 의미를 피해야 할 수도 있습니다.
of ? 쉘 프로그램에서 백슬래시를 앞에 놓거나
따옴표 또는 큰따옴표 사이에 파일 이름과 접미사를 입력합니다. NS ?이름 접미사도 사용할 수 있습니다
출력 그리드가 기본값과 다른 변수 이름을 지정하려면: "z". 보다
grd변환 및 GMT 기술의 CF 및 그리드 파일 형식용 섹션 수정자
특히 3-의 접합을 읽는 방법에 대한 자세한 내용은 참조 및 요리책,
4차원 또는 5차원 격자.

GRID 거리 단위

그리드에 수평 단위로 미터가 없으면 추가 +u단위 입력 파일에
지정된 단위에서 미터로 변환할 이름입니다. 그리드가 지리적인 경우 변환
공급하여 미터까지의 거리 -fg 대신.

고려 사항

netCDF COARDS 그리드는 자동으로 지리적으로 인식됩니다. 다른 그리드의 경우
지리적 그리드는 도를 미터로 변환하려는 경우 선택합니다. -fg. 만약 데이터가
두 극에 가깝다면 그리드 파일을 직사각형에 투영하는 것을 고려해야 합니다.
좌표계 grd 프로젝트.

플레이트 굴곡

탄성 플레이트 굴곡에 대한 FFT 솔루션은 충전 밀도가 하중과 같아야 합니다.
밀도. 이는 일반적으로 부하 바로 아래에서만 해당됩니다. 부하를 넘어
충전물은 밀도가 낮은 퇴적물 또는 심지어 물(또는 공기)인 경향이 있습니다. Wessel [2001] 제안
다른 채우기 밀도의 사양을 허용하는 근사
FFT 솔루션을 허용하면서 로드 밀도. 기본적으로 플레이트 굴곡은
채우기 밀도를 유효 하중 밀도로 사용하기 위해 풀었지만 진폭은 다음과 같습니다.
요인으로 조정 A = sqrt((rm - ri)/(rm - rl)), 이는 이론적 차이입니다.
두 가지 다른 하중 밀도를 사용하는 점 하중으로 인한 진폭. NS
근사치는 매우 좋지만 약한 판에 큰 하중을 가하면 분해됩니다.
드문 상황.

사용 예

2700 및 1035를 사용하여 bat.grd 수심 측량법 위의 수층 효과를 계산하려면
지각과 물의 밀도에 대해 계산하고 그 결과를 water_g.grd에 기록합니다(계산
Parker 확장에서 수심 측량법의 XNUMX제곱으로):

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

이제 그것을 free-air anomaly faa.grd에서 빼면 Bouguer anomaly가 됩니다. 너
왜 더하지 않고 빼는지 궁금할 수 있습니다. 결국 Bouguer 변칙은 척
수층에 의해 나타나는 질량 결핍을 수정하기 위해 추가해야 합니다.
물은 아래의 암석보다 밀도가 낮습니다. 답은 중력 효과가 나타나는 방식에 달려 있습니다.
Parker의 방법과 FFT 사용의 실용적인 측면으로 계산됩니다.

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

MBA 이상을 원하십니까? 지각 맨틀 기여도를 잘 계산하고
해저 이상 현상. 밀도가 6인 2700km 두께의 지각과 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

Te = 7km인 탄성판 bat.grd의 Moho 중력 효과를 계산하려면 밀도
2700, 맨틀 밀도 3300, 평균 깊이 9km

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 추가 = model.grd

다음 명령으로 동일한 결과를 직접 얻을 수 있습니다. 그러나 다음 사항에 주의하십시오.
다음과 같은. 버그 때문인지 어떤 제약 때문인지는 아직 모르겠지만
사실 다음과 이전 명령은 다음과 같은 경우에만 동일한 결과를 제공합니다. -E1
사용. Parker 확장에서 더 높은 수심 측량의 경우 위의 예만
솔기가 올바른 결과를 제공합니다.

GMT Gravfft bat.grd -Gmodel.grd -T7000/2700/3300/1035 -Z9000 -E1

그리고 a 지역 아래 50km 깊이에서 하중에 의해 생성된 지오이드 변칙은 무엇입니까?
그의 수심 측량은 bat.grd, 9km 깊이의 Moho 및 동일한 밀도로 제공됩니다.
전?

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

topo.grd 수심 측량과 faa.grd 자유 대기 이상 사이의 어드미턴스를 계산하려면
유효 6km의 평균 두께 10km 지각의 탄성판 모델을 사용하는 그리드
평균 수심 3km 영역의 탄성 두께:

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

topo.grd 수심 측량과 geoid.grd 지오이드 그리드 사이의 어드미턴스를 계산하려면
위와 동일한 "하부하중"(LFB) 모델 및 40km에서의 지표하중,
그러나 이제 그리드가 지리적이고 주파수 대신 파장을 원한다고 가정합니다.

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

다음을 사용하여 2000km 길이의 프로파일을 따라 LFB의 중력 이론적 어드미턴스를 계산하려면
위와 같은 매개변수

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

참조

루이스, JF 및 MC Neves. 2006, 아조레스 고원의 등압 보상: 3D
어드미턴스 및 일관성 분석. J. 지열 화산. 해상도 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 고도계, J. Geophys에서 추론한 해산 분포. 해상도,
106(B9), 19,431-19,441, http://dx.doi.org/10.1029/2000JB000083

onworks.net 서비스를 사용하여 온라인에서 gravfftgmt 사용