=========================================================== 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 (with Guadalupe pole for NA-PA) ----------------------------------------------------------- Run began on 2008.11. 7 at 16:23: 7 ----------------------------------------------------------- Begin reading parameters for this run from file p_GCN2008088.nki _GCN2008088 [name token for use in building output file names] 2.00E+04 L0 = length of fault trace whose offset rate gets unit weight (in m) 8.00E+08 A0 = area of continuum whose stiffness & isotropy get unit weight (in m**2) 20 number of refinements of each velocity solution 5.00E-16 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 fGCN_merged_WGCEPFM2p2_200810.nki preceding line = filename of fault offset rates fGCN_merged_WGCEPFM2p2_200810.dig preceding line = filename of digitised fault traces s_Gorda-Cal-Nev.nki preceding line = filename of principal stress directions GPS2006_NA_over2kmMergedFM2p2.gps preceding line = filename of geodetic velocities GPS2006_NA_over2kmMergedFM2p2.gp2 preceding line = filename of velocity-covariance matrix F that velocity reference frame of geodetic data is allowed to float F conservative_geodetic adjustment? (using geologic slip rates) GCN8p9.feg preceding line = filename of finite element grid b_GCN8p9.nki preceding line = filename of boundary conditions NA 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 fGCN_merged_WGCEPFM2p2_200810.nki Allocated which_trace = 0.010 MB, total 0.010 MB Allocated f_sense = 0.002 MB, total 0.012 MB Allocated fault_name = 0.119 MB, total 0.131 MB Allocated f_offset_rate = 0.010 MB, total 0.141 MB Allocated f_old_model_offset_rate = 0.010 MB, total 0.151 MB Allocated f_model_offset_rate = 0.010 MB, total 0.160 MB Allocated f_offset_rate_sigma_ = 0.010 MB, total 0.170 MB Allocated f_divide = 0.019 MB, total 0.189 MB Allocated f_dat_shadow = 0.002 MB, total 0.191 MB Allocated f_creeping = 0.002 MB, total 0.194 MB Allocated trace_has_dipslip_rate = 0.004 MB, total 0.198 MB Allocated trace_has_strikeslip_rate = 0.004 MB, total 0.202 MB Allocated f_locking_depth_m_max = 0.010 MB, total 0.211 MB Allocated f_locking_depth_m_min = 0.010 MB, total 0.221 MB 2444 offset-rate data (& shadow strike-slip components) were registered - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Reading fault traces from fGCN_merged_WGCEPFM2p2_200810.dig Allocated trace = 0.147 MB, total 0.368 MB Allocated trace_is = 0.196 MB, total 0.565 MB Allocated trace_loc = 0.067 MB, total 0.631 MB 12868 fault-trace points were read - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Reading stress-direction data from s_Gorda-Cal-Nev.nki Allocated s_ref = 0.059 MB, total 0.691 MB Allocated s_loc = 0.059 MB, total 0.750 MB Allocated s_site = 0.024 MB, total 0.773 MB Allocated s_azim = 0.008 MB, total 0.781 MB Allocated s_sigma_ = 0.008 MB, total 0.789 MB 2068 stress-direction sites were read - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Reading geodetic-velocity data from GPS2006_NA_over2kmMergedFM2p2.gps Allocated external_benchmark_index = 0.005 MB, total 0.794 MB Allocated internal_benchmark_index = 0.005 MB, total 0.799 MB Allocated benchmark_theta = 0.005 MB, total 0.804 MB Allocated benchmark_phi = 0.005 MB, total 0.809 MB Allocated benchmark_uvec = 0.015 MB, total 0.825 MB Allocated benchmark_is = 0.020 MB, total 0.845 MB Allocated benchmark_G = 0.061 MB, total 0.905 MB Allocated benchmark_vw = 0.010 MB, total 0.916 MB Allocated benchmark_unlocked_vw = 0.010 MB, total 0.926 MB Allocated benchmark_reframed_vw = 0.010 MB, total 0.936 MB Allocated benchmark_model_vw = 0.010 MB, total 0.946 MB Allocated benchmark_covariance = 0.020 MB, total 0.966 MB Allocated benchmark_name = 0.101 MB, total 1.067 MB 1324 geodetic-velocity data were read - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Reading finite element grid file GCN8p9.feg Allocated xyz_nod = 0.074 MB, total 1.141 MB Allocated mu_nod = 0.025 MB, total 1.165 MB Allocated vw = 0.049 MB, total 1.215 MB Allocated u-flag = 0.049 MB, total 1.264 MB Checking that no node lies on any fault trace Allocated node = 0.145 MB, total 1.408 MB Allocated a_ = 0.048 MB, total 1.457 MB Allocated crack_index = 0.096 MB, total 1.553 MB Allocated ele_strainrate = 0.145 MB, total 1.697 MB Allocated ele_azim = 0.048 MB, total 1.746 MB Allocated ele_sigma = 0.048 MB, total 1.794 MB Allocated ele_stressed = 0.012 MB, total 1.806 MB Allocated boxed = 0.012 MB, total 1.818 MB 6452 nodes and 12627 elements were read - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Reading velocity boundary conditions from file b_GCN8p9.nki Allocated boundary_node = 0.001 MB, total 1.819 MB Allocated condition = 0.002 MB, total 1.821 MB 259 velocity boundary conditions were read - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Successfully read all input datasets =============================================================================== Finding all data locations in grid coordinates Allocated neighbor = 0.145 MB, total 1.965 MB Allocated center = 0.145 MB, total 2.110 MB Finding geodetic benchmarks Counting fault segments Allocated seg_def = 0.052 MB, total 2.162 MB Allocated seg_end = 0.157 MB, total 2.319 MB Allocated seg_end_is = 0.209 MB, total 2.529 MB Allocated seg_eta_ = 0.026 MB, total 2.555 MB Allocated seg_kappa_ = 0.026 MB, total 2.581 MB Allocated seg_u_ = 0.026 MB, total 2.607 MB Recording fault segments Found all data locations in grid coordinates Count of geodetic benchmarks inside the .feg area: 1210 Allocated covariance_mps2 matrix = 44.681 MB, total 47.288 MB Reading geodetic covariance matrix from GPS2006_NA_over2kmMergedFM2p2.gp2 2069595 entries were processed from this file. Inverting covariance matrix to get normal matrix Allocated normal matrix = 44.681 MB, total 91.969 MB Counting active cracks Allocated local_crack = 0.316 MB, total 92.285 MB Recording active cracks Delta_node = 6158 Allocated ABCDEF matrix = 1185.217 MB, total 1277.502 MB Interpolating stress directions Solving for velocities Unfaulted& e1h_azim_error Refinement dV/V RMS(V) Stressed Boxed (L1; in sigmas) 0 1.00000 7.6823E-10 63.29% 0.00% 0.309 1 0.14544 7.7248E-10 63.29% 0.00% 2.815 2 0.01092 7.7332E-10 63.29% 0.00% 2.371 3 0.00668 7.7337E-10 63.29% 0.00% 2.278 4 0.00462 7.7336E-10 63.29% 0.00% 2.239 5 0.00331 7.7340E-10 63.29% 0.00% 2.230 6 0.00247 7.7347E-10 63.29% 0.00% 2.228 7 0.00187 7.7354E-10 63.29% 0.00% 2.229 8 0.00145 7.7356E-10 63.29% 0.00% 2.231 9 0.00117 7.7359E-10 63.29% 0.00% 2.237 10 0.01527 7.7132E-10 63.29% 11.36% 2.240 11 0.00403 7.7144E-10 63.29% 13.00% 1.094 12 0.00164 7.7140E-10 63.29% 13.22% 0.877 13 0.00131 7.7150E-10 63.29% 13.27% 0.860 14 0.00103 7.7145E-10 63.29% 13.29% 0.857 15 0.00072 7.7147E-10 63.29% 13.29% 0.857 16 0.00056 7.7144E-10 63.29% 13.30% 0.858 17 0.00048 7.7149E-10 63.29% 13.32% 0.859 18 0.00037 7.7144E-10 63.29% 13.33% 0.860 19 0.00031 7.7147E-10 63.29% 13.33% 0.857 20 0.00027 7.7141E-10 63.29% 13.34% 0.858 Largest change is 2.97E-12 = 0.094 mm/a at node 3924 If convergence is not satisfactory, even after 20-40 iterations, then run program Analyze_Velocity_Evolution to get insight. Also consider setting input parameter conservative_geodetic_adjustment to TRUE. Finished computing velocities at nodes. Writing g_GCN2008088.nko ================================================================ Continuum errors: N0 = 0.055, N1 = 0.585, N2 = 1.401 Stress errors: N0 = 0.109, N1 = 0.858, N2 = 1.786 Offset-rate errors: N0 = 0.018, N1 = 0.210, N2 = 0.643 Potency-rate errors: N0 = 0.103, N1 = 0.844, N2 = 1.711 Geodetic errors: N0 = 0.281, N1 = 1.786, N2 = 1.638 ================================================================ Mean (absolute) value of continuum strain rate = 2.92E-16 Root-mean-square value of continuum strain rate = 7.01E-16 Present mu_ of 5.00E-16 falls in this range. Writing f_GCN2008088.nko Writing h_GCN2008088.nko Writing e_GCN2008088.nko Writing s_GCN2008088.nko Writing v_GCN2008088.out Run ended on 2008.11. 9 at 3:58:34 =============================================================