Il s'agit de la commande grdgravmag3dgmt qui peut être exécutée dans le fournisseur d'hébergement gratuit OnWorks en utilisant l'un de nos multiples postes de travail en ligne gratuits tels que Ubuntu Online, Fedora Online, l'émulateur en ligne Windows ou l'émulateur en ligne MAC OS
PROGRAMME:
Nom
grdgravmag3d - Calculer l'effet de gravité d'une grille par la méthode d'Okabe
SYNOPSIS
grdgravmag3d fichier_grd_top [fichier_grd_bot] [ densité ] [ épaisse ] [ fichier_xy ] [ grille ] [
<...> ] [ incrément ] [ z_obs ] [ [nn_pad]|[pad_dist]|[ ] ] [ région ] [ radius ]
[[niveau] ] [ niveau[b|t] ] [ -fg ] [ -x+une|n|-n ]
Remarque: Aucun espace n'est autorisé entre l'indicateur d'option et les arguments associés.
DESCRIPTION
grdgravmag3d calculera l'anomalie de gravité d'un corps décrit par un ou (éventuellement)
deux grilles La sortie peut être soit le long d'un ensemble donné d'emplacements xy, soit sur une grille. Cette
méthode n'est pas particulièrement rapide mais permet de calculer l'anomalie de complexité arbitraire
formes.
REQUIS ARGUMENTS
fichier_grd_top [fichier_grd_bot]
Fichier de grille dont l'effet de gravité va être calculé. Si deux grilles sont fournies
puis l'effet gravitationnel/magnétique du volume entre eux est calculé.
-Cdensité
Définit la densité corporelle en SI. Cette option est mutuellement exclusive avec -H
-Ffichier_xy
Indiquez les emplacements où l'anomalie sera calculée. Notez que cette option est mutuellement
exclusif avec -G.
-Ggrille
Sortir l'anomalie de gravité aux nœuds de ce fichier de grille.
EN OPTION ARGUMENTS
-Eépaisseur
Pour fournir l'épaisseur de couche en m [Par défaut = 500 m].
-Hf_dec/f_dip/m_int/m_dec/m_dip -H+m -Hx|y|z|h|t -H+je|+g|+r|+f|+n
Définit les paramètres pour le calcul de l'anomalie magnétique (peut être utilisé plusieurs fois).
f_dec/f_dip -> déclinaison/inclinaison géomagnétique
m_int/m_dec/m_dip -> intensité/déclinaison/inclinaison magnétique du corps
OU pour un mode grille
+m où 'magfile' est le nom du fichier d'intensité magnétique.
Pour calculer un composant, spécifiez l'un des éléments suivants :
x|X|e|E pour calculer la composante EW.
y|Y|n|N pour calculer la composante NS.
z|Z pour calculer la composante verticale.
h|H pour calculer la composante horizontale.
t|T|f|F pour calculer le champ total.
Pour une inclinaison et une déclinaison variables, utilisez IGRF. Définissez l'un des -H+je|+g|+r|+f|+n
pour faire ça
-Ixinc[unité][=|+][/yinc[unité][=|+]]
x_inc [et éventuellement y_inc] est l'espacement de la grille. Ajoutez éventuellement un suffixe
modificateur. Géographique (degrés) les coordonnées: Ajouter m pour indiquer les minutes d'arc ou s
pour indiquer les secondes d'arc. Si l'une des unités e, f, k, M, n or u est annexé
au lieu de cela, l'incrément est supposé être donné en mètre, pied, km, mille, nautique
mile ou pied d'arpentage américain, respectivement, et sera converti en l'équivalent
degrés de longitude à la latitude moyenne de la région (la conversion dépend de
PROJ_ELLIPSOID). Si /y_inc est donné mais mis à 0 il sera remis à zéro égal à x_inc;
sinon, il sera converti en degrés de latitude. Tous les coordonnées: Si = is
puis le max correspondant x (est) ou y (nord) peut être légèrement ajusté
pour s'adapter exactement à l'incrément donné [par défaut, l'incrément peut être ajusté
légèrement pour s'adapter au domaine donné]. Enfin, au lieu de donner un incrément, vous pouvez
spécifie le nombre of nœuds souhaité en ajoutant + à l'entier fourni
argument; l'incrément est alors recalculé à partir du nombre de nœuds et le
domaine. La valeur d'incrément résultante dépend de la sélection ou non d'un
grille enregistrée sur le quadrillage ou sur les pixels ; voir App-file-formats pour plus de détails.
Remarque : si -Rfichier grd est utilisé alors l'espacement de la grille a déjà été initialisé ; utilisation
-I pour remplacer les valeurs.
-Lz_obs
Définit le niveau d'observation [Par défaut = 0]. C'est la hauteur (z) à laquelle les anomalies
sont calculés.
-Q[nn_pad]|[pad_dist]|[ ]
étendre le domaine of calcul avec respect à sortie -R région.
-Qnn_pad étend artificiellement la largeur du bord extérieur des alvéoles pour avoir un
fausse largeur de n_pad *dx[/dy].
-Qpad_dist étendre la région par west-pad, east+pad, etc.
-Qrégion Même syntaxe que -R.
-R[unité]x min/xmax/ymin/ymax[r] (plus ...)
Spécifiez la région d'intérêt. Remarque : cela remplace la région de grille source
(Par défaut : utiliser la même région que l'entrée)
-Sradius
Définir le rayon de recherche en km (valable uniquement dans le mode deux grilles OU lorsque -E) [Par défaut = 30
km]. Cette option sert à accélérer le calcul en ne calculant pas l'effet de
prismes plus éloignés que radius du nœud courant.
-V[niveau] (plus ...)
Sélectionnez le niveau de verbosité [c].
-Zniveau[b|t]
niveau du plan de référence [Par défaut = 0]. Utilisez cette option lorsque les triangles décrivent
une surface non fermée et le volume est défini à partir de chaque triangle et ce
Niveau de référence. Un exemple sera la profondeur d'eau pour calculer une anomalie de Bouguer.
Utilisez le -Zb or Zt pour fermer le corps en son bas (par exemple, pour calculer l'effet
d'un dôme) ou à son sommet (pour calculer l'effet d'un cuillère).
-fg Les grilles géographiques (dimensions de longitude, latitude) seront converties en mètres
via une approximation « Terre plate » en utilisant les paramètres actuels de l'ellipsoïde.
-x+une|n|-n
Choisissez le nombre de processeurs utilisés en multi-threading (Uniquement disponible avec
versions multi-threads).
+a Utilisez tous les processeurs disponibles.
n Utilisez n processeurs (pas plus que le maximum disponible bien sûr).
-n Utilisez (tous - n) des processeurs.
-^ or juste -
Imprime un court message sur la syntaxe de la commande, puis quitte (REMARQUE : sous Windows
utiliser juste -).
-+ or juste +
Imprimez un message d'utilisation détaillé (aide), y compris l'explication de tout
option spécifique au module (mais pas les options communes GMT), puis se ferme.
-? or aucune arguments
Imprimez un message d'utilisation (aide) complet, y compris l'explication des options, puis
sorties.
--version
Imprimer la version GMT et quitter.
--show-datadir
Affichez le chemin complet vers le répertoire de partage GMT et quittez.
GRID DISTANCE UNITÉS
Si la grille n'a pas de mètre comme unité horizontale, ajoutez +uunité au fichier d'entrée
nom à convertir de l'unité spécifiée en mètre. Si votre grille est géographique, convertissez
distances aux mètres en fournissant -fg à la place.
EXEMPLES
Supposons que vous vouliez calculer l'effet de gravité du fantôme "Sandy Island" avec
ce n'est pas un mont sous-marin fantôme
gmt grdgravmag3d Sandy_bat.grd -C1700 -Z-4300 -fg -I1m -Gsandy_okb.grd -V
Pour calculer la composante verticale due à une aimantation stockée dans mag.grd sur une zone
défini par la surface bat.grd, en utilisant une déclinaison et une inclinaison variables à condition que le
l'IGRF et à l'aide de 4 processeurs, faites :
gmt grdgravmag3d bat.grd -E10000 -Gcomp_Z.grd -Hz -H+n -H+mmag.grd -x4 -V -S50
Utilisez grdgravmag3dgmt en ligne en utilisant les services onworks.net
