grdflexuregmt - אונליין בענן

תָכְנִית:

שֵׁם

grdflexure - חישוב דפורמציה כיפופית של משטחים תלת-ממדיים עבור ריאולוגיות שונות

תַקצִיר

גרדפלקסור טופוגרד rm/rl[/ri]/rw Te[u] להתעקש [ Nx/Ny/Nxy ] [ pדגים ] [ yצעיר ] [
נו_א[/ח_א/נו_ם] ] [ רשימה ] [ [f|q|s|nx/ny][+a|d|h|l][+e|n|m][+tרוחב][+w[סִיוֹמֶת]][+z[p]] [
בטא ] [ -Tt0[u][/t1[u]/dt[u]|n][+l] ] [ [רָמָה] ] [ wd] [ zm] [ -fg ]

הערה: אין רווח בין דגל האפשרות לארגומנטים המשויכים.

תיאור

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

נדרש טיעונים

טופוגרד
קובץ רשת בינארי דו-ממדי עם הטופוגרפיה של העומס (במטרים); ראה קובץ רשת
פורמטים למטה. אם -T משמש, טופוגרד עשוי להיות תבנית שם קובץ עם סמל צף
פורמט נקודה (תחביר C) ושם קובץ טעינה שונה יוגדר וייטען עבור
כל שלב זמן. זמני הטעינה תואמים אפוא את הזמנים שניתנו באמצעות -T (אבל לא
בכל עת צריך קובץ מתאים). לחלופין, תן טופוגרד as
=*flist*, איפה פליסט היא טבלת ASCII עם אחת טופוגרד שם הקובץ וזמן הטעינה לכל
רשומה. זמני טעינה אלה יכולים להיות שונים מזמני ההערכה שניתנו באמצעות -T.
לפורמט זמן טעינה, ראה -T.

-Drm/rl[/ri]/rw
קובע צפיפות עבור מעטפת, עומס, מילוי (אופציונלי, אחרת מניחים שהיא שווה ל-
צפיפות העומס), ומים או אוויר. אם ri שונה מ rl לאחר מכן בקירוב
יימצא פתרון. אם ri לא ניתן אז ברירת המחדל היא rl.

-ETe קובע את עובי הפלטה האלסטית (במטרים); הוסף k עבור ק"מ. אם הגומי
עובי עולה על 1e10 הוא יתפרש כקשיחות כיפוף D (כברירת מחדל)
D מחושב מ Te, מודול יאנג, ויחס פואסון; ראה -C לשנות את אלה
ערכים).

-Gאוסף
If -T מוגדר אז grdfile חייב להיות תבנית שם קובץ המכילה קוד צף
פורמט נקודה (תחביר C). אם תבנית שם הקובץ מכילה גם %s (עבור
שם יחידה) או %c (עבור אות היחידה) אז נשתמש בזמן המתאים (ביחידות
צוין ב -T) כדי ליצור את שמות הקבצים הבודדים, אחרת אנו משתמשים בזמן ב
שנים ללא יחידה.

אופציונאלי טיעונים

-ANx/Ny/Nxy
ציין כוחות דחיסה או התארכות במישור בכיווני x ו-y, כ
כמו גם כל כוח גזירה [אין כוחות במישור]. דחיסה מסומנת על ידי כוח שלילי
ערכים, בעוד שכוחות מתיחה מוגדרים באמצעות ערכים חיוביים.

-Cpדגים
שנה את הערך הנוכחי של יחס פואסון [0.25].

-סייצעיר
שנה את הערך הנוכחי של מודול יאנג [7.0e10 N/m^2].

-Fנו_א[/ח_א/נו_ם]
ציין מודל פירמוביסקוסי בשילוב עם עובי פלטה אלסטי
שצוין באמצעות -Eתן רק צמיגות אחת (נו_א) עבור פלטה אלסטית מעל א
חצי חלל צמיג, או גם להוסיף את עובי האסתנוספירה (ח_א) ו
צמיגות מעטפת תחתונה (נו_ם), כאשר הצמיגות הראשונה כעת היא זו של ה-
אסתנוספרה. תן צמיגות ב-Pa*s. אם נעשה שימוש, תן את עובי
אסתנוספרה במטר; הוסף k עבור ק"מ.

-N[f|q|s|ניו יורק/ניו יורק][+a|[+d|h|l][+e|n|m][+tרוחב][+w[סִיוֹמֶת]][+z[p]]
בחרו או שאלו לגבי מידות רשת מתאימות עבור FFT והגדירו אופציונליות
פרמטרים. שליטה על ממד ה-FFT:
-Nf יכריח את ה-FFT להשתמש במידות האמיתיות של הנתונים.

-Nq אשאל לגבי מידות מתאימות יותר, אדווח עליהן ואז אמשיך.

-Ns יציג רשימה של ממדים אופציונליים, ולאחר מכן ייצא.

-Nניו יורק/ניו יורק יבצע FFT על גודל המערך ניו יורק/ניו יורק (חייב להיות גדול מגודל קובץ הרשת). ברירת מחדל
בוחר ממדים >= נתונים אשר ממטבים את המהירות והדיוק של FFT. אם FFT
ממדים > ממדי קובץ רשת, הנתונים מורחבים ומתחדדים לאפס.

שליטה על ניתוק מגמות של נתונים: הוספת שינויים להסרת מגמה ליניארית:
+d: ביטול מגמה של נתוני שינוי, כלומר הסרת מגמה ליניארית מתאימה ביותר [ברירת מחדל].

+a: הסר רק את הערך הממוצע.

+h: הסר רק את הערך האמצעי, כלומר 0.5 * (מקסימום + מינימום).

+l: עזבו את הנתונים בשקט.

שליטה בהרחבה ובצמצום נתונים: השתמשו במתווים כדי לשלוט באופן שבו ההרחבה
ויש לבצע צמצום:
+e מרחיב את הרשת על ידי הטלת סימטריה של נקודת קצה [ברירת מחדל],

+m מרחיב את הרשת על ידי הטלת סימטריית מראה קצה

+n מכבה את הרחבת הנתונים.

ההתאמה בין הנתונים לקצה רשת ה-FFT מתבצעת [100%]. שינוי
אחוז זה באמצעות +tרוחב. כאשר +n בתוקף, ההצטמצמות מוחלת
במקום זאת לשולי הנתונים מכיוון שאין הרחבה זמינה [0%].

כתיבת בקרה של תוצאות זמניות: לחקירה מפורטת ניתן לכתוב את
רשת ביניים מועברת ל-FFT הקדמי; סביר להניח שזה היה
מנותק במגמה, מוארך על ידי סימטריה נקודתית לאורך כל הקצוות, ומתחדד. הוסף
+w[סִיוֹמֶת] שממנו ייווצרו שמות קבצי הפלט (למשל, קידומת_ingrid.ext)
[מחודד], היכן שלוחה היא סיומת הקובץ שלך. לבסוף, תוכל לשמור את הרשת המורכבת
נוצר על ידי ה-FFT הקדמי על ידי הוספה +zכברירת מחדל אנו כותבים את ה- ממשי ו-
רכיבים דמיוניים ל אינגריד_רִיאָל.שלוחה ו אינגריד_תמונה.שלוחה. לְצַרֵף p כדי לחסוך
במקום זאת, הצורה הקוטבית של גודל ופאזה לקבצים אינגריד_mag.שלוחה ו
אינגריד_שָׁלָב.שלוחה.

-Lרשימה רשמו את השמות וזמני ההערכה של כל הרשתות שנוצרו בטקסט
פילה רשימה. דורש -T.

-Mtm ציין מודל ויסקואלסטי בשילוב עם עובי פלטה אלסטית
שצוין באמצעות -Eהוסף את זמן מקסוול tm עבור המודל הוויסקואלסטי (ב-).

-Sבטא ציין חלק חפיר מורעב בטווח 0-1, כאשר 1 פירושו שהחפיר מורעב לחלוטין
ממולא בחומר בעל צפיפות ri בעוד ש-0 פירושו שהוא מלא רק בחומר של
צפיפות rw (כלומר, רק מים) [1].

-Tt0[u][/t1[u]/dt[u]|n][+ל]
לפרט t0, t1, ותוספת זמן (dt) עבור רצף חישובים [ברירת המחדל היא
צעד אחד, ללא תלות בזמן]. עבור זמן ספציפי יחיד, פשוט התחל
זמן t0היחידה היא שנים; הוסף k עבור קיר ו M עבור מיר. עבור זמן לוגריתמי
קנה מידה, הוסף +l ולציין n צעדים במקום dtלחלופין, תן קובץ
עם הזמנים הרצויים בעמודה הראשונה (ייתכן שלזמנים אלה יהיו יחידות נפרדות
(מצורף, אחרת אנו מניחים שנה). לאחר מכן אנו כותבים קובץ רשת מודל נפרד עבור
כל שלב זמן נתון.

-Wwd הגדר עומק ייחוס למשטח הלא-מעוות במטרים [0]. הוסף k ל
לציין ק"מ.

-Zzm ציין עומק ייחוס למשטח מכופף (למשל, Moho) במטרים; הוסף k עבור ק"מ. חובה
היה חיובי. [0].

-V[רָמָה] (יותר ...)
בחר רמת מילוליות [c].

-fg רשתות גיאוגרפיות (ממדים של קו אורך, קו רוחב) יומרו למטרים
באמצעות קירוב של "כדור הארץ השטוח" תוך שימוש בפרמטרי האליפסואידים הנוכחיים.

-^ or רק -
הדפס הודעה קצרה על התחביר של הפקודה, ואז צא (הערה: ב-Windows
להשתמש רק -).

-+ or רק +
הדפס הודעת שימוש נרחבת (עזרה), כולל הסבר על כל
אפשרות ספציפית למודול (אך לא האפשרויות הנפוצות של GMT), ואז יוצאת.

-? or לא טיעונים
לאחר מכן הדפס הודעת שימוש מלאה (עזרה), כולל הסבר על האפשרויות
יציאות.

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

--show-datadir
הדפס את הנתיב המלא לספריית השיתוף של GMT וצא.

GRID קובץ פורמטים

כברירת מחדל, GMT כותב את הרשת כצף דיוק יחיד ב-NETCDF של תלונת COARDS
פורמט קובץ. עם זאת, GMT מסוגל לייצר קבצי רשת בהרבה רשתות נפוצות אחרות
פורמטים של קבצים וגם מקל על מה שנקרא "אריזה" של רשתות, כתיבת נקודה צפה
נתונים כמספרים שלמים של 1 או 2 בתים. כדי לציין את הדיוק, קנה המידה וההיסט, המשתמש צריך
להוסיף את הסיומת =id[/סולם/לקזז[/נאן]], איפה id הוא מזהה בן שתי אותיות של הרשת
סוג ודיוק, ו סולם ו לקזז הם גורם קנה מידה אופציונלי וקיזוז להיות
מיושם על כל ערכי הרשת, ו נאן הוא הערך המשמש לציון נתונים חסרים. במקרה
שתי הדמויות id לא מסופק, כמו ב =/סולם מ id=nf מונחת. מתי
קריאת רשתות, בדרך כלל הפורמט מזוהה אוטומטית. אם לא, אותה סיומת
ניתן להוסיף לשמות קבצי רשת קלט. לִרְאוֹת grdconvert ו-Section grid-fil-פורמט של
עיון טכני של GMT וספר בישול למידע נוסף.

בעת קריאת קובץ netCDF המכיל רשתות מרובות, GMT יקרא, כברירת מחדל, את
רשת דו-ממדית ראשונה שיכולה למצוא בקובץ הזה. לשדל את GMT לקרוא אחר
משתנה רב ממדי בקובץ הגריד, צרף ?varname לשם הקובץ, איפה
varname הוא שם המשתנה. שימו לב שייתכן שתצטרכו לברוח מהמשמעות המיוחדת
of ? בתוכנית המעטפת שלך על ידי הצבת קו נטוי אחורי לפניה, או על ידי הצבת
שם קובץ וסיומת בין מרכאות או מרכאות כפולות. ה ?varname ניתן להשתמש גם בסיומת
עבור רשתות פלט כדי לציין שם משתנה שונה מברירת המחדל: "z". לִרְאוֹת
grdconvert ו-Sections modifiers-for-CF ו-grid-fil-format של GMT Technical
עיון וספר בישול למידע נוסף, במיוחד כיצד לקרוא חיבורים של 3-,
רשתות 4 או 5 מימדיות.

GRID מרחק יחידות

אם לרשת אין מטר כיחידה אופקית, הוסף +uיחידה לקובץ הקלט
שם להמרה מהיחידה שצוינה למטר. אם הרשת שלך היא גיאוגרפית, המר
מרחקים למטרים על ידי אספקה -fg במקום.

שיקולים

רשתות netCDF COARDS יזוהו אוטומטית כגיאוגרפיות. עבור רשתות אחרות
רשתות גיאוגרפיות שבהן ברצונך להמיר מעלות למטרים, בחר -fgאם הנתונים
קרובים לאחד הקטבים, כדאי לשקול להקרין את קובץ הרשת על גבי מלבני
מערכת קואורדינטות באמצעות פרויקט grd.

צלחת כֶּפֶף אורים

פתרון ה-FFT לכיפוף הפלטה דורש שצפיפות המילוי תהיה שווה לצפיפות העומס.
זה נכון בדרך כלל רק ישירות מתחת לעומס; מעבר לעומס המילוי נוטה
להיות משקעים בעלי צפיפות נמוכה יותר או אפילו מים (או אוויר). וסל [2001] הציע
קירוב המאפשר מפרט של צפיפות מילוי שונה מה-
צפיפות עומס תוך מתן אפשרות לפתרון FFT. בעיקרון, כיפוף הפלטה הוא
נפתר באמצעות צפיפות המילוי כצפיפות העומס האפקטיבית אך האמפליטודות הן
מותאם לפי הגורם A = sqrt ((rm - ri)/(rm - rl)), שהוא ההפרש התאורטי
באמפליטודה עקב עומס נקודתי באמצעות שתי צפיפויות עומס שונות. ה-
הקירוב טוב מאוד אבל מתקלקל לעומסים גדולים על לוחות חלשים, פיה
מצב לא שכיח.

דוגמאות

כדי לחשב את כיפוף הפלטה האלסטית מהעומס topo.nc, עבור לוח בעובי 10 ק"מ עם
צפיפויות אופייניות, נסה

gmt grdflexure topo.nc -Gflex.nc -E10k -D2700/3300/1035

כדי לחשב את התגובה הפירמוביסקוסית לסדרה של טעינות מצטברות הנתונות לפי שם הקובץ
וזמן טעינה בטבלה l.lis בזמן יחיד 1 Ma באמצעות הריאולוגיה שצוינה
ערכים, נסו

כיפוף גרדיום gmt -T1M =l.lis -D3300/2800/2800/1000 -E5k -Gflx/smt_fv_%03.1f_%s.nc -F2e20 -Nf+a

ביבליוגרפיה

קתלס, LM, 1975, השמיים צְמִיגוּת of מה היא של כדור הארץ מעטפת, הוצאת אוניברסיטת פרינסטון.

וסל. פ., 2001, תפוצה עולמית של הרי ים כפי שנגזרה מ-Geosat/ERS-1 ברשת.
גובה-מטריה, J. Geophys. Res., 106(B9), 19,431-19,441,
http://dx.doi.org/10.1029/2000JB000083

השתמש ב- grdflexuregmt באופן מקוון באמצעות שירותי onworks.net