=========================================================== 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 2013. 1. 2 at 12:31:38 ----------------------------------------------------------- Begin reading parameters for this run from file p_NSHM-WUS_2013001.nki _NSHM-WUS_2013001 [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_FM3p1_and_NSHM2014_NoCA_Slippery20121126.nki preceding line = filename of fault offset rates f_UCERF3_FM3p1dsc_and_NSHM2014_NoCA_Haller20121126.dig preceding line = filename of digitised fault traces sUCERF3.nki preceding line = filename of principal stress directions wus5.omeS_reduction06CForFaults.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_CForFaults.feg preceding line = filename of finite element grid b_v2_NSHM-WUS_grid_CForFaults.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_FM3p1dsc_and_NSHM2014_NoCA_Haller20121126.dig Allocated trace = 0.075 MB, total 0.075 MB Allocated trace_is = 0.100 MB, total 0.174 MB Allocated trace_loc = 0.115 MB, total 0.289 MB Allocated f_dig_degrees = 0.029 MB, total 0.318 MB 6531 fault-trace points were read - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Reading fault offset-rate data from f_UCERF3_FM3p1_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.374 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.388 MB Allocated f_model_offset_rate = 0.004 MB, total 0.392 MB Allocated f_offset_rate_sigma_ = 0.004 MB, total 0.396 MB Allocated f_offset_rate_floor = 0.004 MB, total 0.400 MB Allocated f_offset_rate_ceiling = 0.004 MB, total 0.404 MB Allocated f_divide = 0.008 MB, total 0.413 MB Allocated f_dat_shadow = 0.001 MB, total 0.414 MB Allocated f_creeping = 0.001 MB, total 0.415 MB Allocated trace_has_dipslip_rate = 0.007 MB, total 0.422 MB Allocated trace_has_strikeslip_rate = 0.007 MB, total 0.429 MB Allocated f_locking_depth_m_max = 0.004 MB, total 0.433 MB Allocated f_locking_depth_m_min = 0.004 MB, total 0.437 MB 1038 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.474 MB Allocated s_site = 0.007 MB, total 0.481 MB Allocated s_azim = 0.002 MB, total 0.483 MB Allocated s_sigma_ = 0.002 MB, total 0.486 MB 638 stress-direction sites were read - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Reading geodetic-velocity data from wus5.omeS_reduction06CForFaults.gps Allocated external_benchmark_index = 0.022 MB, total 0.508 MB Allocated internal_benchmark_index = 0.022 MB, total 0.530 MB Allocated benchmark_theta = 0.022 MB, total 0.552 MB Allocated benchmark_phi = 0.022 MB, total 0.573 MB Allocated benchmark_uvec = 0.066 MB, total 0.639 MB Allocated benchmark_is = 0.088 MB, total 0.727 MB Allocated benchmark_G = 0.263 MB, total 0.989 MB Allocated benchmark_vw = 0.044 MB, total 1.033 MB Allocated benchmark_unlocked_vw = 0.044 MB, total 1.077 MB Allocated benchmark_reframed_vw = 0.044 MB, total 1.120 MB Allocated benchmark_model_vw = 0.044 MB, total 1.164 MB Allocated benchmark_covariance = 0.088 MB, total 1.252 MB Allocated benchmark_normal = 0.088 MB, total 1.339 MB Allocated benchmark_name = 0.438 MB, total 1.777 MB 5735 geodetic-velocity data were read - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Reading finite element grid file NSHM-WUS_grid_CForFaults.feg Allocated xyz_nod = 0.124 MB, total 1.901 MB Allocated mu_nod = 0.041 MB, total 1.942 MB Allocated vw = 0.083 MB, total 2.025 MB Allocated vw_interseismic = 0.083 MB, total 2.108 MB Allocated u-flag = 0.083 MB, total 2.190 MB Checking that no node lies on any fault trace Allocated node = 0.244 MB, total 2.435 MB Allocated a_ = 0.081 MB, total 2.516 MB Allocated crack_index = 0.163 MB, total 2.679 MB Allocated ele_strainrate = 0.244 MB, total 2.923 MB Allocated ele_azim = 0.081 MB, total 3.005 MB Allocated ele_sigma = 0.081 MB, total 3.086 MB Allocated ele_stressed = 0.020 MB, total 3.106 MB Allocated boxed = 0.020 MB, total 3.127 MB 10846 nodes and 21346 elements were read - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Reading velocity boundary conditions from file b_v2_NSHM-WUS_grid_CForFaults.nki Allocated boundary_node = 0.000 MB, total 3.127 MB Allocated condition = 0.001 MB, total 3.128 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.244 MB, total 3.372 MB Allocated center = 0.244 MB, total 3.617 MB Finding geodetic benchmarks Counting fault segments Allocated seg_def = 0.038 MB, total 3.654 MB Allocated seg_end = 0.113 MB, total 3.767 MB Allocated seg_end_is = 0.151 MB, total 3.918 MB Allocated seg_eta_ = 0.019 MB, total 3.936 MB Allocated seg_kappa_ = 0.019 MB, total 3.955 MB Allocated seg_u_ = 0.019 MB, total 3.974 MB Recording fault segments Found all data locations in grid coordinates Count of geodetic benchmarks inside the .feg area: 5468 Counting active cracks Allocated local_crack = 0.284 MB, total 4.259 MB Recording active cracks Delta_node = 169 Allocated ABCDEF matrix = 28.658 MB, total 32.917 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.6926E-10 0 60.55% 0.00% 0.393 1 0.02647 6.7043E-10 0 60.55% 0.00% 2.711 2 0.00862 6.7098E-10 0 60.55% 0.00% 2.602 3 0.00545 6.7128E-10 0 60.55% 0.00% 2.510 4 0.00382 6.7144E-10 0 60.55% 0.00% 2.464 5 0.00274 6.7154E-10 0 60.55% 0.00% 2.433 6 0.00202 6.7159E-10 0 60.55% 0.00% 2.421 7 0.00150 6.7162E-10 0 60.55% 0.00% 2.418 8 0.00114 6.7162E-10 0 60.55% 0.00% 2.415 9 0.00089 6.7162E-10 0 60.55% 0.00% 2.415 10 0.00070 6.7161E-10 0 60.55% 0.00% 2.417 11 0.00056 6.7160E-10 0 60.55% 0.00% 2.418 12 0.00046 6.7158E-10 0 60.55% 0.00% 2.421 13 0.00039 6.7157E-10 0 60.55% 0.00% 2.422 14 0.00031 6.7155E-10 0 60.55% 0.00% 2.422 15 0.00150 6.7159E-10 97 60.55% 0.00% 2.422 16 0.00197 6.7163E-10 82 60.55% 0.00% 2.451 17 0.00246 6.7173E-10 70 60.55% 0.00% 2.485 18 0.00291 6.7191E-10 55 60.55% 0.00% 2.518 19 0.00255 6.7202E-10 37 60.55% 0.00% 2.563 20 0.00338 6.7218E-10 29 60.55% 0.00% 2.567 21 0.00437 6.7240E-10 25 60.55% 0.00% 2.578 22 0.00484 6.7268E-10 20 60.55% 0.00% 2.593 23 0.00526 6.7300E-10 20 60.55% 0.00% 2.622 24 0.00450 6.7324E-10 12 60.55% 0.00% 2.647 25 0.00250 6.7337E-10 8 60.55% 0.00% 2.670 26 0.00214 6.7348E-10 8 60.55% 0.00% 2.676 27 0.00181 6.7357E-10 9 60.55% 0.00% 2.687 28 0.00111 6.7364E-10 8 60.55% 0.00% 2.697 29 0.00126 6.7369E-10 7 60.55% 0.00% 2.702 30 0.01291 6.7329E-10 6 60.55% 18.86% 2.705 31 0.00346 6.7302E-10 25 60.55% 22.00% 1.252 32 0.00160 6.7287E-10 12 60.55% 22.34% 0.986 33 0.00100 6.7275E-10 7 60.55% 22.42% 0.953 34 0.00069 6.7267E-10 5 60.55% 22.47% 0.949 35 0.00048 6.7263E-10 3 60.55% 22.51% 0.943 36 0.00033 6.7260E-10 0 60.55% 22.51% 0.938 37 0.00023 6.7258E-10 1 60.55% 22.51% 0.937 38 0.00076 6.7257E-10 2 60.55% 22.53% 0.937 39 0.00023 6.7257E-10 0 60.55% 22.53% 0.939 40 0.00015 6.7257E-10 0 60.55% 22.53% 0.939 41 0.00010 6.7257E-10 0 60.55% 22.54% 0.939 42 0.00008 6.7257E-10 0 60.55% 22.54% 0.939 43 0.00007 6.7257E-10 0 60.55% 22.54% 0.938 44 0.00007 6.7257E-10 0 60.55% 22.55% 0.938 45 0.00007 6.7256E-10 0 60.55% 22.55% 0.938 Largest change is 1.18E-12 = 0.037 mm/a at node 8542 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_2013001.nko ================================================================ Continuum errors: N0 = 0.044, N1 = 0.459, N2 = 1.356 Stress errors: N0 = 0.105, N1 = 0.938, N2 = 2.564 Offset-rate errors: N0 = 0.045, N1 = 0.454, N2 = 0.932 Potency-rate errors: N0 = 0.132, N1 = 1.002, N2 = 1.621 Geodetic errors: N0 = 0.375, N1 = 2.182, N2 = 2.295 ================================================================ Mean (absolute) value of continuum strain rate = 2.30E-16 Root-mean-square value of continuum strain rate = 6.78E-16 Present mu_ of 5.00E-16 falls in this range. Writing f_NSHM-WUS_2013001.nko Writing h_NSHM-WUS_2013001.nko Writing e_NSHM-WUS_2013001.nko Writing s_NSHM-WUS_2013001.nko Writing v_NSHM-WUS_2013001.out Writing v_interseismic_NSHM-WUS_2013001.out Run ended on 2013. 1. 2 at 15:38:26 =============================================================