=========================================================== 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 Department of Earth and Space Sciences University of California NeoKinema: Version 3.0 of 13 December 2012 ----------------------------------------------------------- Run began on 2021. 5. 4 at 18:35: 5 ----------------------------------------------------------- Begin reading parameters for this run from file p_NSHM-WUS_2013002.nki _NSHM-WUS_2013002 [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) 2.00E+08 A0 = area of continuum whose stiffness & isotropy get unit weight (in m**2) 45 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 0.00 12.00 minimum and maximum locking depths of intraplate faults, in km 4.00 22.00 minimum and maximum locking depths of subduction zones, in km F that faults are treated as stress direction data f_UCERF3_FM3p2_and_NSHM2014_NoCA_Slippery20121126.nki preceding line = filename of fault offset rates f_UCERF3_FM3p2dsc_and_NSHM2014_NoCA_Haller20121126.dig preceding line = filename of digitised fault traces sUCERF3.nki preceding line = filename of principal stress directions wus5.omeS_reduction04DForFaults.gps preceding line = filename of geodetic velocities none 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) NSHM-WUS_grid_DForFaults.feg preceding line = filename of finite element grid b_v2_NSHM-WUS_grid_DForFaults.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 traces from f_UCERF3_FM3p2dsc_and_NSHM2014_NoCA_Haller20121126.dig Allocated trace = 0.075 MB, total 0.075 MB Allocated trace_is = 0.100 MB, total 0.175 MB Allocated trace_loc = 0.115 MB, total 0.289 MB Allocated f_dig_degrees = 0.029 MB, total 0.318 MB 6544 fault-trace points were read - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Reading fault offset-rate data from f_UCERF3_FM3p2_and_NSHM2014_NoCA_Slippery20121126.nki Allocated which_trace = 0.004 MB, total 0.322 MB Allocated f_sense = 0.001 MB, total 0.323 MB Allocated fault_name = 0.052 MB, total 0.375 MB Allocated f_dat_dip_degrees = 0.004 MB, total 0.379 MB Allocated f_offset_rate = 0.004 MB, total 0.383 MB Allocated f_offset_rate_bracketed = 0.001 MB, total 0.384 MB Allocated f_old_model_offset_rate = 0.004 MB, total 0.389 MB Allocated f_model_offset_rate = 0.004 MB, total 0.393 MB Allocated f_offset_rate_sigma_ = 0.004 MB, total 0.397 MB Allocated f_offset_rate_floor = 0.004 MB, total 0.401 MB Allocated f_offset_rate_ceiling = 0.004 MB, total 0.405 MB Allocated f_divide = 0.008 MB, total 0.413 MB Allocated f_dat_shadow = 0.001 MB, total 0.415 MB Allocated f_creeping = 0.001 MB, total 0.416 MB Allocated trace_has_dipslip_rate = 0.007 MB, total 0.423 MB Allocated trace_has_strikeslip_rate = 0.007 MB, total 0.430 MB Allocated f_locking_depth_m_max = 0.004 MB, total 0.434 MB Allocated f_locking_depth_m_min = 0.004 MB, total 0.438 MB 1045 offset-rate data (& shadow strike-slip components) were registered - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Reading stress-direction data from sUCERF3.nki Allocated s_ref = 0.018 MB, total 0.456 MB Allocated s_loc = 0.018 MB, total 0.475 MB Allocated s_site = 0.007 MB, total 0.482 MB Allocated s_azim = 0.002 MB, total 0.484 MB Allocated s_sigma_ = 0.002 MB, total 0.487 MB 638 stress-direction sites were read - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Reading geodetic-velocity data from wus5.omeS_reduction04DForFaults.gps Allocated external_benchmark_index = 0.023 MB, total 0.509 MB Allocated internal_benchmark_index = 0.023 MB, total 0.532 MB Allocated benchmark_theta = 0.023 MB, total 0.555 MB Allocated benchmark_phi = 0.023 MB, total 0.577 MB Allocated benchmark_uvec = 0.068 MB, total 0.645 MB Allocated benchmark_is = 0.090 MB, total 0.735 MB Allocated benchmark_G = 0.271 MB, total 1.006 MB Allocated benchmark_vw = 0.045 MB, total 1.051 MB Allocated benchmark_unlocked_vw = 0.045 MB, total 1.096 MB Allocated benchmark_reframed_vw = 0.045 MB, total 1.141 MB Allocated benchmark_model_vw = 0.045 MB, total 1.186 MB Allocated benchmark_covariance = 0.090 MB, total 1.277 MB Allocated benchmark_normal = 0.090 MB, total 1.367 MB Allocated benchmark_name = 0.451 MB, total 1.818 MB 5915 geodetic-velocity data were read - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Reading finite element grid file NSHM-WUS_grid_DForFaults.feg Allocated xyz_nod = 0.125 MB, total 1.943 MB Allocated mu_nod = 0.042 MB, total 1.985 MB Allocated vw = 0.083 MB, total 2.069 MB Allocated vw_interseismic = 0.083 MB, total 2.152 MB Allocated u-flag = 0.083 MB, total 2.235 MB Checking that no node lies on any fault trace Allocated node = 0.246 MB, total 2.482 MB Allocated a_ = 0.082 MB, total 2.564 MB Allocated crack_index = 0.164 MB, total 2.728 MB Allocated ele_strainrate = 0.246 MB, total 2.975 MB Allocated ele_azim = 0.082 MB, total 3.057 MB Allocated ele_sigma = 0.082 MB, total 3.139 MB Allocated ele_stressed = 0.021 MB, total 3.160 MB Allocated boxed = 0.021 MB, total 3.180 MB 10940 nodes and 21534 elements were read - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Reading velocity boundary conditions from file b_v2_NSHM-WUS_grid_DForFaults.nki Allocated boundary_node = 0.000 MB, total 3.181 MB Allocated condition = 0.001 MB, total 3.181 MB 100 velocity boundary conditions were read Table of Euler poles used to compute type-4 velocity boundary conditions: (Note that all are relative to current reference plate NA.) PLATE N_latitude E_longitude Degrees/Ma Note ----- ---------- ----------- ---------- ------------ NA-NA 0.000 0.000 0.0000 stationary PA-NA -49.890 102.989 0.7660 built-in ----- ---------- ----------- ---------- ------------ - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Successfully read all input datasets =============================================================================== Finding all data locations in grid coordinates Allocated neighbor = 0.246 MB, total 3.428 MB Allocated center = 0.246 MB, total 3.674 MB Finding geodetic benchmarks Counting fault segments Allocated seg_def = 0.038 MB, total 3.712 MB Allocated seg_end = 0.114 MB, total 3.825 MB Allocated seg_end_is = 0.151 MB, total 3.977 MB Allocated seg_eta_ = 0.019 MB, total 3.996 MB Allocated seg_kappa_ = 0.019 MB, total 4.015 MB Allocated seg_u_ = 0.019 MB, total 4.034 MB Recording fault segments Found all data locations in grid coordinates Count of geodetic benchmarks inside the .feg area: 5648 Counting active cracks Allocated local_crack = 0.286 MB, total 4.319 MB Recording active cracks Delta_node = 172 Allocated ABCDEF matrix = 29.419 MB, total 33.738 MB Interpolating stress directions Solving for velocities: Active Unfaulted& e1h_azim_error Refinement dV/V RMS(V) brackets Stressed Boxed (N1; in sigmas) 0 1.00000 6.7442E-10 0 60.84% 0.00% 0.394 1 0.02718 6.7540E-10 0 60.84% 0.00% 2.723 2 0.00823 6.7578E-10 0 60.84% 0.00% 2.562 3 0.00518 6.7592E-10 0 60.84% 0.00% 2.434 4 0.00356 6.7599E-10 0 60.84% 0.00% 2.358 5 0.00256 6.7602E-10 0 60.84% 0.00% 2.319 6 0.00190 6.7603E-10 0 60.84% 0.00% 2.294 7 0.00145 6.7602E-10 0 60.84% 0.00% 2.282 8 0.00113 6.7600E-10 0 60.84% 0.00% 2.273 9 0.00090 6.7598E-10 0 60.84% 0.00% 2.270 10 0.00073 6.7596E-10 0 60.84% 0.00% 2.270 11 0.00060 6.7594E-10 0 60.84% 0.00% 2.272 12 0.00050 6.7592E-10 0 60.84% 0.00% 2.272 13 0.00043 6.7590E-10 0 60.84% 0.00% 2.272 14 0.00037 6.7589E-10 0 60.84% 0.00% 2.272 15 0.00142 6.7590E-10 100 60.84% 0.00% 2.272 16 0.00182 6.7594E-10 85 60.84% 0.00% 2.299 17 0.00234 6.7605E-10 71 60.84% 0.00% 2.337 18 0.00295 6.7623E-10 56 60.84% 0.00% 2.361 19 0.00255 6.7638E-10 35 60.84% 0.00% 2.406 20 0.00346 6.7655E-10 26 60.84% 0.00% 2.422 21 0.00456 6.7681E-10 23 60.84% 0.00% 2.433 22 0.00506 6.7710E-10 21 60.84% 0.00% 2.447 23 0.00534 6.7744E-10 15 60.84% 0.00% 2.475 24 0.00442 6.7773E-10 12 60.84% 0.00% 2.497 25 0.00247 6.7791E-10 6 60.84% 0.00% 2.520 26 0.00212 6.7802E-10 6 60.84% 0.00% 2.533 27 0.00182 6.7810E-10 7 60.84% 0.00% 2.541 28 0.00100 6.7817E-10 7 60.84% 0.00% 2.551 29 0.00103 6.7820E-10 6 60.84% 0.00% 2.556 30 0.01235 6.7787E-10 8 60.84% 18.15% 2.558 31 0.00341 6.7762E-10 26 60.84% 21.12% 1.183 32 0.00169 6.7746E-10 17 60.84% 21.53% 0.956 33 0.00151 6.7727E-10 11 60.84% 21.62% 0.916 34 0.00085 6.7719E-10 8 60.84% 21.69% 0.913 35 0.00090 6.7716E-10 7 60.84% 21.71% 0.910 36 0.00039 6.7715E-10 2 60.84% 21.72% 0.910 37 0.00026 6.7714E-10 2 60.84% 21.74% 0.910 38 0.00019 6.7713E-10 1 60.84% 21.75% 0.909 39 0.00016 6.7712E-10 1 60.84% 21.76% 0.909 40 0.00012 6.7712E-10 0 60.84% 21.76% 0.909 41 0.00010 6.7711E-10 0 60.84% 21.77% 0.909 42 0.00008 6.7711E-10 0 60.84% 21.77% 0.909 43 0.00007 6.7711E-10 0 60.84% 21.77% 0.908 44 0.00006 6.7711E-10 0 60.84% 21.77% 0.909 45 0.00006 6.7711E-10 0 60.84% 21.78% 0.909 Largest change is 9.37E-13 = 0.030 mm/a at node 3662 If convergence is not satisfactory, even after 30~45 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_NSHM-WUS_2013002.nko ================================================================ Continuum errors: N0 = 0.043, N1 = 0.467, N2 = 1.391 Stress errors: N0 = 0.103, N1 = 0.909, N2 = 2.464 Offset-rate errors: N0 = 0.044, N1 = 0.438, N2 = 0.906 Potency-rate errors: N0 = 0.129, N1 = 0.967, N2 = 1.582 Geodetic errors: N0 = 0.378, N1 = 2.195, N2 = 2.304 ================================================================ Mean (absolute) value of continuum strain rate = 2.33E-16 Root-mean-square value of continuum strain rate = 6.96E-16 Present mu_ of 5.00E-16 falls in this range. Writing f_NSHM-WUS_2013002.nko Writing h_NSHM-WUS_2013002.nko Writing e_NSHM-WUS_2013002.nko Writing s_NSHM-WUS_2013002.nko Writing v_NSHM-WUS_2013002.out Writing v_interseismic_NSHM-WUS_2013002.out Run ended on 2021. 5. 4 at 19:43:26 =============================================================