========== STRATEGIC PARAMETERS (DEFINE THE REAL-EARTH PROBLEM) ======== MODEL 90-03: NORTHERN OPTION; LIKE 89-28, WITH IMPROVED I.C. AND B.C.'S CRUST MANTLE \ PARAMETER (LINE ABOVE IS TITLE) 0.85 0.85 COEFFICIENT OF FRICTION 6.60E+03 4.26E+04 PRE-EXPONENTIAL SHEAR STRESS CONSTANT FOR CREEP 2.50E+08 N/A PRE-EXPONENTIAL FOR LOWER CRUST, BELOW CONRAD 0.333333 0.285714 STRAIN-RATE EXPONENT FOR CREEP (1/N) 0.333333 N/A STRAIN-RATE EXPONENT FOR LOWER CRUST, BELOW CONRAD 16000. 18314. B FOR CREEP =(ACTIVATION ENERGY)/R/N (IN K) 9194. N/A B FOR CREEP OF LOWER CRUST, BELOW CONRAD 0.0E-05 17.1E-05 C FOR CREEP = DERIVITIVE OF B WITH RESPECT TO DEPTH 0.0E-05 N/A C FOR CREEP OF LOWER CRUST, BELOW CONRAD 10.E+09 10.E+09 MAXIMUM SHEAR STRESS UNDER ANY CONDITIONS 10.E+09 N/A MAXIMUM SHEAR STRESS FOR LOWER CRUST, BELOW CONRAD 2.50E+05 4.06E+05 THERMAL CONDUCTIVITY (ENERGY/LENGTH/SEC/DEG) 0.012 0.012 THERMAL DIFFUSIVITY (LENGTH**2/SEC) 6.060E-06 1.60E-07 RADIOACTIVE HEAT PRODUCTION (ENERGY/VOL/SEC) 39.1E+05 70.1E+05 THICKNESS OF LAYER IN NORMAL CONTINENT 1223. 1573. CONVECTION TEMPERATURE (TMAX) IN DEGREES KELVIN 2.88 3.30 DENSITY, CORRECTED TO 0 K 24.0E-06 31.0E-06 VOLUMETRIC THERMAL EXPANSION (1/VOL)*(D.VOL/D.T) 7.21E+05 9.32E+05 MEAN P-WAVE VELOCITY (VP) AT 0 K AND HIGH PRESSURE -1.00E-04 -1.00E-04 (1/VP)*(D.VP/D.T): TEMPERATURE SENSITIVITY OF VP 0.0 0.05 (DELTA.VP/VP): STRAIN(EZZ)-INDUCED ANISOTROPY IN VP 1.00 1.00 CHARACTERISTIC STRAIN TO DEVELOP ANISOTROPY 3.157 DENSITY OF ASTHENOSPHERE (ADJUST TO CORRECT ALL ELEVATIONS) 1.03 DENSITY OF GROUNDWATER AND LAKES/OCEANS 0.0 EFFECTIVE PRESSURE (BIOT) COEFFICIENT, 0.0 TO 1.0 980. GRAVITATIONAL ACCELERATION (LENGTH/SEC**2) 6371.E+05 RADIUS OF EARTH (EFFECTIVELY DEFINES YOUR LENGTH UNIT) -105.00 LONGITUDE OF X/Y ORIGIN IN DEGREES (EAST = +, WEST = -) 39.75 LATITUDE OF X/Y ORIGIN IN DEGREES (NORTH = +, SOUTH = -) 45.0 LATITUDE OF BASE-MAP CONIC PROJECTION TANGENT IN DEG.(NORTH=+) 1 BELOW-INDEX; SELECTS PLATE MODEL FOR BASAL BOUNDARY CONDITIONS: 0=NONE,1=N.AMER./NORTH,2=N.AMER./SOUTH,3=S.AMER.,4=ASIA 283. SURFACE TEMPERATURE IN DEGREES KELVIN 653. TEMPERATURE OF SLAB-TOP SHEAR ZONE AT 1000 KM INLAND, IN KELVIN 275.E+06 SHEAR STRESS LIMIT ON BASE OF CONTINENT (MELANGE STRENGTH) 730.E+05 WIDTH OF CORDILLERA, MEASURED TRENCH TO PLAINS (OR 0.0 FOR NONE) 0.E+06 EXTRA SHEAR STRESS APPLIED TO LEFT MARGIN FOREARC ONLY -300.E+06 EXCESS-WEIGHT/UNIT-AREA OF NEW LITHOSPHERE W/ RESP. TO ASTHENOSPHERE 595.E+06 THERMAL EXCESS-WEIGHT/UNIT-AREA @ 100 MA WRT NEW LITHOSPHERE 0.E+06 NON-ISOSTATIC FLEXURAL UPLIFT BY SLAB (IN FOREARC REGION ONLY) ============== TACTICAL PARAMETERS (HOW TO FIND THE SOLUTION) =========== 14 NUMBER OF ROWS OF 2-ELEMENT QUADRILATERALS 10 NUMBER OF COLUMNS OF 2-ELEMENT QUADRILATERALS 2.6824E+15 BEGINNING OF CALCULATION (POSITIVE SECONDS BEFORE PRESENT) 3.16E+13 SIZE OF TIME STEPS (POSITIVE SECONDS); MAY BE REDUCED BY PROG 2.2090E+15 ENDING OF CALCULATION (POSITIVE SECONDS BEFORE PRESENT) 100.E+05 MAX. HORIZONTAL DISPLACEMENT OF ANY NODE IN 1 TIME STEP 20.E+05 MAX. CHANGE IN LAYER THICKNESS BY PURE SHEAR IN 1 TIME STEP 300.E+05 WIDTH OF LINEAR RAMP SMOOTHING OF SLAB WEIGHT 1000 MAX. NUMBER OF CRUSTAL-THICKNESS SMOOTHINGS, EACH TIMESTEP 10 MAXIMUM ITERATIONS WITHIN VELOCITY SOLUTION IN EACH TIMESTEP 0.0200 ACCEPTABLE RMS FRACTIONAL ERROR (STOPS ITERATION EARLY) 6.00E+25 MAXIMUM AVERAGE VISCOSITY ALLOWED IN ANY LAYER (STRESS*SEC) 1.00E+17 MAXIMUM LAYER/LAYER COUPLING ALLOWED (STRESS/VELOCITY-DIFFERENCE) CRUST MANTLE\ LIMITS ON LAYER THICKNESSES: 10.E+05 10.E+05 MINIMUM THICKNESS (TRIGGERS VOLUME ADDITION) 80.E+05 150.E+05 MAXIMUM THICKNESS (TRIGGERS VOLUME REDUCTION) F ALLREP: SHOULD REPORTS BE PRODUCED AT EVERY TIMESTEP ? 2 NUMBER OF INTERMEDIATE REPORTS (IF ALLREP = F) T THAT DETAILED OUTPUT IS PROVIDED ON DEVICE 9 (USUALLY T) ================== INITIALIZATION PARAMETERS (INITIAL CONDITIONS) ======= T RESTRT: IF = T, THEN RESTART FROM OLD REPORT ON DEVICE 8 1 IF (RESTRT): ORDINAL NUMBER OF OLD REPORT ON DEVICE 8 ------ NEXT LINES ARE USED ONLY IF RESTRT = F -------------------------- 40.0E+05 INITIAL DEPTH OF CONRAD DISCONTINUITY IN CRUST OF PLAINS 54. INITIAL STEADY-STATE HEAT-FLOW OF PLAINS (ENERGY/LENGTH**2/SE 3.40E+05 HEIGHT OF ALTIPLANO IN CORDILLERA (OR 0.0 FOR NONE) 3.0E-8 GROSS ESTIMATE OF DECOLLEMENT VELOCITY BETWEEN CRUST AND MANTLE T OLDGRD: SHALL EXISITNG GRID OF NODES LOCATIONS BE READ (ON 8)? -------- FOLLOWING LINES DEFINE AN AUTOMATICALLY GENERATED GRID, AND --- ----------- ARE USED ONLY IF RESTRT = F AND OLDGRD = F ----------------- 220.E+05 "WIDTH OF GRID FROM "LEFT" (TRENCH SIDE) TO "RIGHT" (INLAND) 475.E+05 "HEIGHT" OF GRID FROM "TOP" (NODE ROW 1) TO "BOTTOM" (LAST ROW) 90. ANGLE GRID IS ROTATED FROM ("RIGHT"=+X,"TOP"=+Y),IN DEG. ===== POST-PROCESSING PLOT CONTROL PARAMETERS (NOT USED BY LARAMY) ===== 999 PLOTFILE # ON DEVICE-8 DATASET (OR 999 FOR LAST, OR 0 FOR ALL) F 25.E+06 125.E+06+1 PLOT MANTLE BASAL SHEAR STRESS F 25.E+06 125.E+06+1 PLOT CRUSTAL BASAL SHEAR STRESS F 3.0E-08 15.0E-08+1 PLOT MANTLE VELOCITY VECTORS F 1.6E-08 8.0E-08+1 PLOT CRUSTAL VELOCITY VECTORS F 1.0E-15 5.0E-15+1 PLOT MANTLE GREATEST PRIN. STRAIN RATES F +1 PLOT CRUSTAL GREATEST PRIN. STRAIN RATE F +1 PLOT MANTLE PRIN. STRESS ANOMALY INTEGR F 1.0E+15 5.0E+15+1 PLOT CRUST PRIN. STRESS ANOMALY INTEGR F PLOT MANTLE GRID OF ELEMENTS T PLOT CRUSTAL GRID OF ELEMENTS F -1 PLOT MANTLE LAYER THICKENING RATE F -1 PLOT CRUSTAL LAYER THICKENING RATE F 20.E+05 80.E+05-1 PLOT MANTLE THICKNESS F 5.E+05 35.E+05+1 PLOT CRUSTAL THICKNESS F 1.E+09 5.E+09+1 PLOT MANTLE BASAL TEMPERATURE F 100. 900.+1 PLOT CRUSTAL BASAL TEMPERATURE F 0.0000001+1 PLOT TELESEISMIC P TRAVEL-TIME RESIDUAL F 0.5E+05 0.001E+05+1 PLOT ISOSTATIC ELEVATIONS F 10. 54.+1 PLOT PALEO-HEAT-FLOW F 0.5E+05 0.001E+05+1 PLOT ELEVATIONS AFTER DELAMINATION F 0.03 0.15+1 PLOT LOG10(CRUSTAL NET STRAIN, TO NOW) F 5. 0.0000001+1 PLOT CRUSTAL ROTATION IN DEGREES F 5.E+05 35.E+05+1 PLOT UPPER CRUSTAL THICKNESS F 5.E+05 25.E+05+1 PLOT LOWER CRUSTAL THICKNESS 2.500E+07 SCALE OF VERSATEC PLOTS (INPUT LENGTH UNITS)/(CM ON VERSATEC) 10 APPROXIMATE NUMBER OF CONTOURS IN EACH PLOT T F PLOT STATE OUTLINES ? RETRO-DEFORMED ? 1 PEN WEIGHT (INTEGER) OF CONTOURS & ELEMENT SIDES 3 PEN WEIGHT (INTEGER) OF STATE OUTLINES, IF ANY. 2 PEN WEIGHT (INTEGER) OF TEXT LABELS & NODES 0.500 RMS LENGTH OF PLOTTED VECTORS, IN INCHES T COLOR? (F GIVES BLACK-AND-WHITE) (APPLIES TO VERSATEC ONLY) 0. 0. (XCENTR,YCENTR) OF PLOTS (in m); MAY BE NEEDED FOR CLIPPING?