זוהי הפקודה talwani3dgmt שניתן להריץ בספק האירוח החינמי של OnWorks באמצעות אחת מתחנות העבודה המקוונות המרובות שלנו, כגון Ubuntu Online, Fedora Online, אמולטור מקוון של Windows או אמולטור מקוון של MAC OS
תָכְנִית:
שֵׁם
talwani3d - חישוב חריגות שיפוע כבידה באוויר חופשי, גיאואיד או אנכי על פני דו-ממד
גופים
תַקצִיר
talwani3d [ שולחן מודל ] [ ] [ rho ] ] [ f|n|v ] [ grdfile ] [ תוֹסֶפֶת ] [ [h][v] ] [
קובץ מסלול ] [ באזור ] [ רָמָה|obsgrid ] [ [רָמָה] ] [ -דוּ] [ -fg ] [ -i]
[ -o] [ -r ] [ -x[[-]n]]
הערה: אין רווח בין דגל האפשרות לארגומנטים המשויכים.
תיאור
talwani3d יקרא את הרב-קטע שולחן מודל מקובץ או קלט רגיל. הקובץ הזה
מכיל קווי מתאר של גוף תלת מימדי בשונה z-מפלסים, עם קו מתאר אחד לכל קטע. ה
כותרת הקטע חייבת להכיל את הפרמטרים zlevel rho, הקובע את z רמת קווי המתאר
והצפיפות של פרוסה זו (ייתכן שצפיפות פרוסה בודדת תידחה על ידי קבוע
ניגודיות צפיפות הניתנת באמצעות -D). אנו יכולים לחשב חריגות על רשת במרחק שווה (על ידי
ציון רשת חדשה עם -R ו -I או לספק רשת תצפית עם גבהים) או ב
נקודות פלט שרירותיות שצוינו באמצעות -N. בחר מתוך חריגות באוויר חופשי, כוח משיכה אנכי
חריגות שיפוע, או חריגות גיאואידיות. זמינות אפשרויות לבקרת יחידות צירים ו
כיוון.
נדרש טיעונים
שולחן מודל
הקובץ המתאר את קווי המתאר האופקיים של הגופים. קווי מתאר יהיו
נסגר אוטומטית אם עדיין לא סגור, וקודקודים חוזרים יהיו
בוטלו.
-Iשינק[יחידה][=|+][/yinc[יחידה[+]]
x_inc [ובאופציה y_inc] הוא מרווח הרשת. לחלופין, הוסף סיומת
שינוי. גיאוגרפי (מעלות) קואורדינטות: צרף m כדי לציין דקות קשת או s
כדי לציין שניות קשת. אם אחת מהיחידות e, f, k, M, n or u מצורף
במקום זאת, ההנחה היא שהתוספת ניתנת במטר, רגל, ק"מ, מייל, ימי
מייל או רגל סקר בארה"ב, בהתאמה, ויומרו למקבילה
מעלות קו אורך בקו הרוחב האמצעי של האזור (ההמרה תלויה ב
PROJ_ELLIPSOID). אם /y_inc נתון אך מוגדר ל-0 הוא יאופס שווה ל x_inc;
אחרת הוא יומר לקו רוחב מעלות. תעשיות קואורדינטות: אם = is
מצורף ואז המקסימום המתאים x (מזרח) או y (צפון) עשוי להיות מותאם מעט
כדי להתאים בדיוק לתוספת הנתונה [כברירת מחדל, התוספת עשויה להיות מותאמת
מעט כדי להתאים לתחום הנתון]. לבסוף, במקום לתת תוספת אפשר
ציין את מספר of צמתים רצוי על ידי הוספה + למספר השלם שסופק
טַעֲנָה; לאחר מכן, התוספת מחושבת מחדש ממספר הצמתים וה-
תְחוּם. ערך התוספת המתקבל תלוי אם בחרת ב-a
רשת רשומה קו רשת או פיקסל; ראה פורמטים של קבצי אפליקציה לפרטים.
הערה: אם -Rgrdfile נמצא בשימוש אז מרווח הרשת כבר אותחל; להשתמש
-I לעקוף את הערכים.
-ר[יחידה]xmin/xmax/ymin/ymax[ר] (יותר ...)
ציין את אזור העניין.
אופציונאלי טיעונים
-A השמיים z-ציר צריך להיות חיובי כלפי מעלה [ברירת המחדל היא למטה].
-Dיחידה מגדיר ניגודיות צפיפות קבועה שעוקפת כל הגדרה בקובץ הדגם, ב-kg/m^3.
-Ff|n|v
ציין את רכיב שדה הכבידה הרצוי. לבחור בין f (אנומליה באוויר חופשי)
[בְּרִירַת מֶחדָל], n (גיאואיד) או v (שיפוע כבידה אנכי).
-Ggrdfile
ציין את השם של קובץ רשת הפלט; ראה פורמטים של קבצי GRID להלן). נדרש
כאשר רשת במרחק שווה משתמע לתפוקה. אם -N משמש אז הפלט הוא
נכתב ל-stdout.
-מ[ח][v]
מגדיר יחידות בשימוש. לְצַרֵף h כדי לציין מרחקים אופקיים הם בק"מ [מ'], וכן
לצרף z כדי לציין שהמרחקים האנכיים הם בק"מ [מ'].
-Nקובץ מסלול
מציין מיקומים שבהם נרצה לחשב את הערך החזוי. כאשר אפשרות זו
בשימוש, אין רשתות ורשומות נתוני הפלט נכתבות ל-stdout.
-V[רָמָה] (יותר ...)
בחר רמת מילוליות [c].
-Zרָמָה|obsgrid
הגדר את רמת התצפית כקבוע או תן שם של רשת עם
רמות תצפית. אם משתמשים באחרון, הרשת קובעת את רשת הפלט
אזור [0].
-דוּ[ncols][t] (יותר ...)
בחר קלט בינארי מקורי. [ברירת המחדל היא 2 עמודות קלט].
-fg רשתות גיאוגרפיות (מידות קו אורך, קו רוחב) יומרו לק"מ באמצעות a
קירוב "אדמה שטוחה" באמצעות הפרמטרים האליפסואידים הנוכחיים.
-h[i|o][n][+c][+d][+rהֶעָרָה][+rכותרת] (יותר ...)
דלג או הפק רשומות כותרות. לא בשימוש עם נתונים בינאריים.
-iצווארונים[ל][ססולם][oלקזז][,...] (יותר ...)
בחר עמודות קלט (0 הוא העמודה הראשונה).
-oצווארונים[,...] (יותר ...)
בחר עמודות פלט (0 הוא העמודה הראשונה).
-r (יותר ...)
הגדר רישום צומת פיקסל [קו רשת].
-איקס[[-]n] (יותר ...)
הגבל את מספר הליבות בשימוש באלגוריתמים מרובי הליכי (OpenMP נדרש).
-:[i|o] (יותר ...)
החלף את העמודה הראשונה והשנייה בקלט ו/או פלט.
-^ or רק -
הדפס הודעה קצרה על התחביר של הפקודה, ואז צא (הערה: ב-Windows
להשתמש רק -).
-+ or רק +
הדפס הודעת שימוש נרחבת (עזרה), כולל הסבר על כל
אפשרות ספציפית למודול (אך לא האפשרויות הנפוצות של GMT), ואז יוצאת.
-? or לא טיעונים
לאחר מכן הדפס הודעת שימוש מלאה (עזרה), כולל הסבר על האפשרויות
יציאות.
--גִרְסָה
הדפס גרסת GMT וצא.
--show-datadir
הדפס את הנתיב המלא לספריית השיתוף של GMT וצא.
יחידות
ליחידת מרחק מפה, צרף יחידה d עבור דרגת קשת, m עבור דקת קשת, ו s עבור קשת
שני, או e עבור מד [ברירת מחדל], f לרגל, k עבור ק"מ, M עבור מייל חוקי, n עבור ימיים
מייל, ו u עבור רגל סקר בארה"ב. כברירת מחדל אנו מחשבים מרחקים כאלה באמצעות כדור
קירוב עם עיגולים גדולים. הכנה מראש - למרחק (או שהיחידה אינה מרחק הוא
נתון) כדי לבצע חישובי "אדמה שטוחה" (מהירים יותר אך פחות מדויקים) או להקדים + ל
לבצע חישובים גיאודזיים מדויקים (איטיים יותר אך מדויקים יותר).
דוגמאות
כדי לחשב את חריגות האוויר החופשי על רשת על גוף דו-ממדי שעבר קווי מתאר
נשמר ב-body.txt, תוך שימוש ב-1.7 g/cm^3 כניגודיות הצפיפות, נסה
gmt talwani3d -R-200/200/-200/200 -I2 -G3dgrav.nc body.txt -D1700 -Fg
כדי להשיג את האנומליה של שיפוע הכבידה האנכי לאורך המסילה ב-crossing.txt עבור
אותו דגם, נסה
gmt talwani3d -Ncrossing.txt body.txt -D1700 -Fv > vgg_crossing.txt
לבסוף, האנומליה הגיאואידית לאורך אותו מסילה ב-crossing.txt עבור אותו דגם היא
חזר על ידי
gmt talwani3d -Ncrossing.txt body.txt -D1700 -Fn > n_crossing.txt
ביבליוגרפיה
Talwani, M., and M. Ewing (1960), חישוב מהיר של משיכה כבידה של
גופים תלת מימדיים בעלי צורה שרירותית, גיאופיזיקה, 25(203-225).
השתמש ב-talwani3dgmt באינטרנט באמצעות שירותי onworks.net