כוח המשיכה - מקוון בענן

תָכְנִית:

שֵׁם

כוח המשיכה - מחשב את שדה הכבידה של כדור הארץ

תַקצִיר

כוח משיכה [ -n שם ] [ -d dir ] [ -G | -D | -A | -H ] [ -c lat h ] [ -w ] [ -p הקודם ] [ -v
] [ --מפריד תגובה commentdelim ] [ --גִרְסָה | -h | - עזרה ] [ --קובץ קלט בקובץ |
--מחרוזת קלט סרגל ] [ --קו-מפריד קו ] [ --קובץ פלט אוסף ]

תיאור

כוח משיכה קורא במיקומים על קלט סטנדרטי ומדפיס את שדה הכבידה על
פלט סטנדרטי.

שורת הקלט היא בצורה lat לון h. lat ו לון הם קו הרוחב והאורך
מבוטא במעלות עשרוניות או במעלות, דקות ושניות; לפרטים על המותר
פורמטים עבור קווי רוחב וקו אורך, עיין בסעיף "קואורדינטות גיאוגרפיות" של
GeoConvert(1). h הוא הגובה מעל אליפסואיד במטרים; כמות זו היא אופציונלית
וברירת המחדל היא 0. לחלופין, ניתן לחשב את שדה הכבידה בנקודות שונות
מעגל של קו רוחב (קבוע lat ו h) דרך ה- -c אוֹפְּצִיָה; במקרה זה רק את
יש לציין קו אורך בקווי הקלט. הכמויות המודפסות כפופות ל
מה היא -G (ברירת מחדל), -D, -A, או -H אפשרויות.

כל דגמי הכבידה הנתמכים, מלבד grs80, משתמשים ב-WGS84 כאליפסואיד הייחוס a
= 6378137 מ', f = 1/298.257223563, אוֹמֶגָה = 7292115e-11 rad/s, ו GM = 3986004.418e8
m^3/s^2.

אפשרויות

-n השתמש במודל שדה הכבידה שם במקום ברירת המחדל "egm96". ראה "מודלים".

-d לקרוא מודלים של כוח הכבידה מ dir במקום ברירת המחדל. ראה "מודלים".

-G מחשב את התאוצה עקב כוח הכבידה (כולל התאוצה הצנטריפוגלית עקב
סיבוב כדור הארץ) g. הפלט מורכב מ gx gy gz (הכל ב-m/s^2), איפה
מה היא x, y, ו z הרכיבים נמצאים בכיוון מזרח, צפון ומעלה,
בהתאמה. בְּדֶרֶך כְּלַל gz הוא שלילי.

-D לחשב את הפרעת הכבידה דלתא = g - גמא, שם גמא הוא כוח המשיכה "הרגיל".
עקב אליפסואיד ההתייחסות. הפלט מורכב מ deltax deltay דלתאז (הכל ב
mGal, 1 mGal = 10^-5 m/s^2), כאשר ה- x, y, ו z הרכיבים נמצאים במזרח,
כיווני צפון ולמעלה, בהתאמה. ציין זאת deltax = gx, כי גאמקס =
0.

-A מחשב את אנומליית הכבידה. הפלט מורכב מ-3 פריטים Dg01 xi ETA, שם
Dg01 הוא ב-mGal (1 mGal = 10^-5 m/s^2) ו xi ו ETA נמצאים בשניות קשת. ה
אנומליה כבידה משווה את שדה הכבידה g at P עם כוח המשיכה הרגיל
גמא at Q שם P נמצא אנכית למעלה Q ופוטנציאל הכבידה ב P
שווה לפוטנציאל הרגיל ב Q. Dg01 נותן את ההבדל בגדלים של
שני הווקטורים הללו ו xi ו ETA תן את ההבדל בכיוונים שלהם (כמו צפון
ורכיבים מזרחיים). החישוב משתמש בקירוב כדורי כדי להתאים את
תוצאות של תוכניות הסינתזה של ה-NGA.

-H חשב את גובה הגיאואיד מעל אליפסואיד הייחוס (במטרים). בזה
מקרה, h צריך להיות אפס. התוצאות תואמות במדויק את התוצאות של ה-NGA
תוכניות סינתזה. GeoidEval(1) יכול לחשב גבהים גיאואידים הרבה יותר מהר על ידי
אינטרפולציה על רשת של תוצאות מחושבות מראש; אולם התוצאות מ GeoidEval(1)
מדויקים רק לכמה מילימטרים.

-c להעריך את השדה במעגל של קו רוחב שניתן על ידי lat ו h במקום לקרוא את אלה
כמויות מקווי הקלט. במקרה הזה, כוח משיכה יכול לחשב את השדה
הרבה יותר מהר. אם מחושבים גבהים גיאואידים (ה -H אפשרות), ואז
h חייב להיות אפס.

-w בקלט ובפלט, קו האורך קודם לקו הרוחב (חוץ מזה שבקלט זה יכול להיות
נדחק על ידי מציין חצי כדור, N, S, E, W).

-p הגדר את דיוק הפלט ל הקודם. כברירת מחדל הקודם הוא 5 להאצה עקב
כוח הכבידה, 3 עבור הפרעת הכבידה ואנומליה, ו-4 עבור גובה הגיאואיד.

-v הדפס מידע על מודל הכבידה על שגיאת תקן לפני עיבוד
קֶלֶט.

--מפריד תגובה
הגדר את מפריד ההערות ל commentdelim (למשל, "#" או "//"). אם מוגדר, הקלט
קווים ייסרקו עבור המפריד הזה, ואם יימצא, את המפריד ושאר
השורה תוסר לפני העיבוד ויצורף לאחר מכן לפלט
קו (מופרד ברווח).

--גִרְסָה
גרסת הדפסה ויציאה.

-h שימוש בהדפסה, נתיב הכבידה והשם המוגדרים כברירת מחדל, ויציאה.

- עזרה
להדפיס את התיעוד המלא ולצאת.

--קובץ קלט
קרא קלט מהקובץ בקובץ במקום מקלט סטנדרטי; שם קובץ של "-"
מייצג קלט סטנדרטי.

--מחרוזת קלט
קרא קלט מהמחרוזת סרגל במקום מקלט רגיל. כל ההתרחשויות
של תו מפריד השורות (ברירת המחדל היא נקודה-פסיק) ב סרגל מומרים ל
שורות חדשות לפני תחילת הקריאה.

--קו-מפריד
הגדר את תו מפריד השורות ל קו. כברירת מחדל זהו נקודה-פסיק.

--קובץ פלט
לכתוב פלט לקובץ אוסף במקום לפלט סטנדרטי; שם קובץ של "-"
מייצג פלט סטנדרטי.

מודלים

כוח משיכה מחשב את שדה הכבידה באמצעות אחד מהדגמים הבאים

egm84, דגם כבידה של כדור הארץ 1984. ראה
http://earth-info.nga.mil/GandG/wgs84/gravitymod/wgs84_180/wgs84_180.html
egm96, דגם כבידה של כדור הארץ 1996. ראה
http://earth-info.nga.mil/GandG/wgs84/gravitymod/egm96/egm96.html
egm2008, דגם כבידה של כדור הארץ 2008. ראה
http://earth-info.nga.mil/GandG/wgs84/gravitymod/egm2008
wgs84, World geodetic system 1984. זה מחזיר את הנורמלי
כוח הכבידה עבור אליפסואיד WGS84.
grs80, מערכת התייחסות גיאודטית 1980. זה מחזיר את הנורמלי
כוח הכבידה עבור אליפסואיד GRS80.

מודלים אלה מקרובים את שדה הכבידה מעל פני כדור הארץ. על ידי
ברירת מחדל, נעשה שימוש במודל הכבידה "egm96". זה עשוי להשתנות על ידי הגדרת הסביבה
המשתנה "GEOGRAPHICLIB_GRAVITY_NAME" או עם -n אוֹפְּצִיָה.

דגמי הכבידה ייטענו מתוך ספרייה שצוינה בזמן הקומפילציה. זה אולי
השתנה על ידי הגדרת משתני הסביבה "GEOGRAPHICLIB_GRAVITY_PATH" או
"GEOGRAPHICLIB_DATA", או עם ה- -d אוֹפְּצִיָה. ה -h אפשרות מדפיסה את ברירת המחדל של כוח המשיכה
נתיב ושם. להשתמש ב -v אפשרות לברר את שם הנתיב המלא של קובץ הנתונים.

הוראות להורדה והתקנה של דגמי כבידה זמינות בכתובת
<http://geographiclib.sf.net/html/gravity.html#כוח המשיכה>.

הסביבה

GEOGRAPHICLIB_GRAVITY_NAME
עוקף את שם הכבידה המוגדר כברירת מחדל בזמן ההידור של "egm96". ה -h אפשרות מדווחת על
ערך של GEOGRAPHICLIB_GRAVITY_NAME, אם מוגדר, אחרת הוא מדווח על זמן ההידור
ערך. אם ה -n שם אם כן, נעשה שימוש באפשרות שם מקבל עדיפות.

GEOGRAPHICLIB_GRAVITY_PATH
עוקף את נתיב הכבידה המוגדר כברירת מחדל בזמן הקומפילציה. זה בדרך כלל
"/usr/local/share/GeographicLib/gravity" במערכות דמויות Unix ו
"C:/ProgramData/GeographicLib/כבידה" במערכות Windows. ה -h אפשרות מדווחת על
ערך של GEOGRAPHICLIB_GRAVITY_PATH, אם מוגדר, אחרת הוא מדווח על זמן ההידור
ערך. אם ה -d dir אם כן, נעשה שימוש באפשרות dir מקבל עדיפות.

GEOGRAPHICLIB_DATA
דרך נוספת לעקוף את נתיב הכבידה המוגדר כברירת מחדל בזמן הקומפילציה. אם זה מוגדר (ואם
GEOGRAPHICLIB_GRAVITY_PATH לא מוגדר), ואז $GEOGRAPHICLIB_DATAנעשה שימוש בכוח הכבידה.

שגיאות

שורת קלט לא חוקית תדפיס הודעת שגיאה לפלט סטנדרטי שמתחיל ב
"ERROR:" וגורמים כוח משיכה כדי להחזיר קוד יציאה של 1. עם זאת, שגיאה אינה גורמת
כוח משיכה לחסל; השורות הבאות יומרו.

דוגמאות

שדה הכבידה מ-EGM2008 בראש הר האוורסט

echo 27:59:17N 86:55:32E 8820 | Gravity -n egm2008
=> -0.00001 0.00103 -9.76782

השתמש ב-Gravity באינטרנט באמצעות שירותי onworks.net