=========================================================== A record of a run of program NeoKinema: (NEOtectonic, anelastic velocity solution from KINEMAtic data such as plate rotations, geodetic velocities, and fault slip rates with thin-viscous-shell interpolation guided by a smooth interpolated stress field) by Peter Bird Dept. of Earth and Space Sciences University of California Los Angeles, CA 90095-1567 Version 2.2 of 30 January 2008 ----------------------------------------------------------- Run began on 2008.12. 6 at 12:34:52 ----------------------------------------------------------- Begin reading parameters for this run from file Param22432.nki KRPNS22432 [name token for use in building output file names] 5.00E+05 L0 = length of fault trace whose offset rate gets unit weight (in m) 1.00E+10 A0 = area of continuum whose stiffness & isotropy get unit weight (in m**2) 40 number of refinements of each velocity solution 1.00E-15 mu_ = scalar measure of typical anelastic strain rates in continuum (/s) 3.20E-17 xi_ = small strain-rate increment, /s 20.0 sigma, in degrees, of angle between [heave vectors of dip-slip faults] & [trace-normal direction] 6.371E+03 radius of the planet (R), in km 1.00 12.00 minimum and maximum locking depths of intraplate faults, in km 14.00 40.00 minimum and maximum locking depths of subduction zones, in km F that faults are treated as stress direction data faultphil_08.nki preceding line = filename of fault offset rates Phil_faults_05.dig preceding line = filename of digitised fault traces s_philippines.nki preceding line = filename of principal stress directions Phil_All.gps preceding line = filename of geodetic velocities none preceding line = filename of velocity-covariance matrix T that velocity reference frame of geodetic data is allowed to float T conservative_geodetic adjustment? (using geologic slip rates) phil_orogen.feg preceding line = filename of finite element grid b_phil_org_11.nki preceding line = filename of boundary conditions SU plate defining velocity reference frame for type-4 boundary conditions F that velocity solutions in all iterations will be written to v_log.nko End Parameter Section =============================================================================== Begin reading input data files Reading fault offset-rate data from faultphil_08.nki Allocated which_trace = 0.002 MB, total 0.002 MB Allocated f_sense = 0.000 MB, total 0.002 MB Allocated fault_name = 0.024 MB, total 0.027 MB Allocated f_offset_rate = 0.002 MB, total 0.028 MB Allocated f_old_model_offset_rate = 0.002 MB, total 0.030 MB Allocated f_model_offset_rate = 0.002 MB, total 0.032 MB Allocated f_offset_rate_sigma_ = 0.002 MB, total 0.034 MB Allocated f_divide = 0.004 MB, total 0.038 MB Allocated f_dat_shadow = 0.000 MB, total 0.039 MB Allocated f_creeping = 0.000 MB, total 0.039 MB Allocated trace_has_dipslip_rate = 0.000 MB, total 0.039 MB Allocated trace_has_strikeslip_rate = 0.000 MB, total 0.040 MB Allocated f_locking_depth_m_max = 0.002 MB, total 0.042 MB Allocated f_locking_depth_m_min = 0.002 MB, total 0.044 MB 498 offset-rate data (& shadow strike-slip components) were registered - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Reading fault traces from Phil_faults_05.dig Allocated trace = 0.061 MB, total 0.105 MB Allocated trace_is = 0.081 MB, total 0.186 MB Allocated trace_loc = 0.005 MB, total 0.191 MB 5327 fault-trace points were read - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Reading stress-direction data from s_philippines.nki Allocated s_ref = 0.058 MB, total 0.248 MB Allocated s_loc = 0.058 MB, total 0.306 MB Allocated s_site = 0.023 MB, total 0.329 MB Allocated s_azim = 0.008 MB, total 0.336 MB Allocated s_sigma_ = 0.008 MB, total 0.344 MB 2010 stress-direction sites were read - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Reading geodetic-velocity data from Phil_All.gps Allocated external_benchmark_index = 0.001 MB, total 0.344 MB Allocated internal_benchmark_index = 0.001 MB, total 0.345 MB Allocated benchmark_theta = 0.001 MB, total 0.346 MB Allocated benchmark_phi = 0.001 MB, total 0.346 MB Allocated benchmark_uvec = 0.002 MB, total 0.348 MB Allocated benchmark_is = 0.002 MB, total 0.350 MB Allocated benchmark_G = 0.007 MB, total 0.357 MB Allocated benchmark_vw = 0.001 MB, total 0.358 MB Allocated benchmark_unlocked_vw = 0.001 MB, total 0.359 MB Allocated benchmark_reframed_vw = 0.001 MB, total 0.360 MB Allocated benchmark_model_vw = 0.001 MB, total 0.362 MB Allocated benchmark_covariance = 0.002 MB, total 0.364 MB Allocated benchmark_name = 0.011 MB, total 0.375 MB 150 geodetic-velocity data were read - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Reading finite element grid file phil_orogen.feg Allocated xyz_nod = 0.016 MB, total 0.392 MB Allocated mu_nod = 0.005 MB, total 0.397 MB Allocated vw = 0.011 MB, total 0.408 MB Allocated u-flag = 0.011 MB, total 0.419 MB Checking that no node lies on any fault trace Allocated node = 0.032 MB, total 0.451 MB Allocated a_ = 0.011 MB, total 0.462 MB Allocated crack_index = 0.021 MB, total 0.484 MB Allocated ele_strainrate = 0.032 MB, total 0.516 MB Allocated ele_azim = 0.011 MB, total 0.527 MB Allocated ele_sigma = 0.011 MB, total 0.537 MB Allocated ele_stressed = 0.003 MB, total 0.540 MB Allocated boxed = 0.003 MB, total 0.543 MB 1438 nodes and 2813 elements were read - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Reading velocity boundary conditions from file b_phil_org_11.nki Allocated boundary_node = 0.000 MB, total 0.543 MB Allocated condition = 0.000 MB, total 0.543 MB 49 velocity boundary conditions were read - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Successfully read all input datasets =============================================================================== Finding all data locations in grid coordinates Allocated neighbor = 0.032 MB, total 0.575 MB Allocated center = 0.032 MB, total 0.608 MB Finding geodetic benchmarks Counting fault segments Allocated seg_def = 0.013 MB, total 0.620 MB Allocated seg_end = 0.038 MB, total 0.658 MB Allocated seg_end_is = 0.050 MB, total 0.708 MB Allocated seg_eta_ = 0.006 MB, total 0.714 MB Allocated seg_kappa_ = 0.006 MB, total 0.720 MB Allocated seg_u_ = 0.006 MB, total 0.727 MB Recording fault segments Found all data locations in grid coordinates Count of geodetic benchmarks inside the .feg area: 69 Allocated covariance_mps2 matrix = 0.145 MB, total 0.872 MB Adding reference-frame-loosening rotations of 10.00 degree/Ma Inverting covariance matrix to get normal matrix Allocated normal matrix = 0.145 MB, total 1.017 MB Counting active cracks Allocated local_crack = 0.071 MB, total 1.088 MB Recording active cracks Delta_node = 1351 Allocated ABCDEF matrix = 57.568 MB, total 58.657 MB Interpolating stress directions Solving for velocities Unfaulted& e1h_azim_error Refinement dV/V RMS(V) Stressed Boxed (L1; in sigmas) 0 1.00000 1.6377E-09 53.86% 0.00% 1.375 1 0.13019 1.6307E-09 53.86% 0.00% 1.119 2 0.04796 1.6503E-09 53.86% 0.00% 1.559 3 0.02267 1.6605E-09 53.86% 0.00% 1.442 4 0.01119 1.6668E-09 53.86% 0.00% 1.474 5 0.00642 1.6712E-09 53.86% 0.00% 1.452 6 0.00304 1.6733E-09 53.86% 0.00% 1.458 7 0.00141 1.6742E-09 53.86% 0.00% 1.456 8 0.00188 1.6759E-09 53.86% 0.00% 1.463 9 0.00072 1.6758E-09 53.86% 0.00% 1.468 10 0.00047 1.6762E-09 53.86% 0.00% 1.469 11 0.00041 1.6761E-09 53.86% 0.00% 1.469 12 0.00058 1.6758E-09 53.86% 0.00% 1.468 13 0.00016 1.6758E-09 53.86% 0.00% 1.469 14 0.00020 1.6759E-09 53.86% 0.00% 1.468 15 0.00024 1.6758E-09 53.86% 0.00% 1.467 16 0.00057 1.6754E-09 53.86% 0.00% 1.467 17 0.00022 1.6754E-09 53.86% 0.00% 1.467 18 0.00054 1.6756E-09 53.86% 0.00% 1.467 19 0.00032 1.6754E-09 53.86% 0.00% 1.466 20 0.05251 1.6148E-09 53.86% 11.52% 1.466 21 0.00431 1.6132E-09 53.86% 12.55% 0.616 22 0.00263 1.6142E-09 53.86% 12.73% 0.493 23 0.00448 1.6111E-09 53.86% 12.76% 0.463 24 0.00342 1.6146E-09 53.86% 12.76% 0.452 25 0.00191 1.6135E-09 53.86% 12.76% 0.446 26 0.00265 1.6150E-09 53.86% 12.76% 0.446 27 0.00160 1.6163E-09 53.86% 12.76% 0.445 28 0.00133 1.6163E-09 53.86% 12.76% 0.445 29 0.00093 1.6154E-09 53.86% 12.80% 0.444 30 0.00110 1.6149E-09 53.86% 12.80% 0.444 31 0.00139 1.6152E-09 53.86% 12.80% 0.443 32 0.00122 1.6166E-09 53.86% 12.80% 0.443 33 0.00257 1.6147E-09 53.86% 12.80% 0.443 34 0.00224 1.6164E-09 53.86% 12.80% 0.443 35 0.00150 1.6160E-09 53.86% 12.80% 0.442 36 0.00147 1.6151E-09 53.86% 12.80% 0.442 37 0.00177 1.6168E-09 53.86% 12.80% 0.442 38 0.00242 1.6148E-09 53.86% 12.80% 0.441 39 0.00073 1.6152E-09 53.86% 12.80% 0.441 40 0.00334 1.6180E-09 53.86% 12.80% 0.441 Largest change is 8.58E-12 = 0.271 mm/a at node 1298 If convergence is not satisfactory, even after 20-40 iterations, then run program Analyze_Velocity_Evolution to get insight. Finished computing velocities at nodes. Correcting geodetic velocities, to best-fit the current model, by adding rotation of 0.1065 deg/Ma about (-49.561N, -52.487E) Writing gKRPNS22432.nko ================================================================ Continuum errors: N0 = 0.065, N1 = 0.696, N2 = 1.086 Stress errors: N0 = 0.033, N1 = 0.441, N2 = 0.796 Offset-rate errors: N0 = 0.007, N1 = 0.328, N2 = 0.588 Potency-rate errors: N0 = 0.011, N1 = 0.781, N2 = 1.158 Geodetic errors: N0 = 0.507, N1 = 2.377, N2 = 2.012 ================================================================ Mean (absolute) value of continuum strain rate = 6.96E-16 Root-mean-square value of continuum strain rate = 1.09E-15 Present mu_ of 1.00E-15 falls in this range. Writing fKRPNS22432.nko Writing hKRPNS22432.nko Writing eKRPNS22432.nko Writing sKRPNS22432.nko Writing vKRPNS22432.out Run ended on 2008.12. 6 at 12:40:34 =============================================================