=========================================================== 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.1 of 14 August 2007 ----------------------------------------------------------- Run began on 2007. 8.14 at 21:46:13 ----------------------------------------------------------- Begin reading parameters for this run from file p_ptb2007008.nki _ptb2007008 [name token for use in building output file names] 6.40E+04 L0 = length of fault trace whose offset rate gets unit weight (in m) 1.46E+09 A0 = area of continuum whose stiffness & isotropy get unit weight (in m**2) 20 number of refinements of each velocity solution 3.70E-16 mu_ = scalar measure of typical anelastic strain rates in continuum (/s) 3.00E-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 11.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 ..\common\f_ptb_1.nki preceding line = filename of fault offset rates ..\common\f_ptb_1.dig preceding line = filename of digitised fault traces ..\common\s_ptb.nki preceding line = filename of principal stress directions ..\common\ptb11.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 F conservative_geodetic adjustment? (using geologic slip rates) ..\common\PTB3nbd.feg preceding line = filename of finite element grid ..\common\b_PTB3nbd.nki preceding line = filename of boundary conditions EU 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 ..\common\f_ptb_1.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.025 MB, total 0.027 MB Allocated f_offset_rate = 0.002 MB, total 0.029 MB Allocated f_old_model_offset_rate = 0.002 MB, total 0.031 MB Allocated f_model_offset_rate = 0.002 MB, total 0.033 MB Allocated f_offset_rate_sigma_ = 0.002 MB, total 0.035 MB Allocated f_divide = 0.004 MB, total 0.039 MB Allocated f_dat_shadow = 0.000 MB, total 0.039 MB Allocated f_creeping = 0.000 MB, total 0.040 MB Allocated trace_has_strikeslip_rate = 0.000 MB, total 0.040 MB 514 offset-rate data (& shadow strike-slip components) were registered - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Reading fault traces from ..\common\f_ptb_1.dig Allocated trace = 0.015 MB, total 0.055 MB Allocated trace_is = 0.020 MB, total 0.075 MB Allocated trace_loc = 0.006 MB, total 0.081 MB 1322 fault-trace points were read - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Reading stress-direction data from ..\common\s_ptb.nki Allocated s_ref = 0.025 MB, total 0.106 MB Allocated s_loc = 0.025 MB, total 0.131 MB Allocated s_site = 0.010 MB, total 0.141 MB Allocated s_azim = 0.003 MB, total 0.145 MB Allocated s_sigma_ = 0.003 MB, total 0.148 MB 876 stress-direction sites were read - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Reading geodetic-velocity data from ..\common\ptb11.gps Allocated external_benchmark_index = 0.006 MB, total 0.154 MB Allocated internal_benchmark_index = 0.006 MB, total 0.160 MB Allocated benchmark_theta = 0.006 MB, total 0.166 MB Allocated benchmark_phi = 0.006 MB, total 0.171 MB Allocated benchmark_uvec = 0.017 MB, total 0.189 MB Allocated benchmark_is = 0.023 MB, total 0.212 MB Allocated benchmark_G = 0.070 MB, total 0.282 MB Allocated benchmark_vw = 0.012 MB, total 0.293 MB Allocated benchmark_unlocked_vw = 0.012 MB, total 0.305 MB Allocated benchmark_reframed_vw = 0.012 MB, total 0.316 MB Allocated benchmark_model_vw = 0.012 MB, total 0.328 MB Allocated benchmark_covariance = 0.023 MB, total 0.351 MB Allocated benchmark_name = 0.116 MB, total 0.467 MB 1521 geodetic-velocity data were read - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Reading finite element grid file ..\common\PTB3nbd.feg Allocated xyz_nod = 0.018 MB, total 0.485 MB Allocated mu_nod = 0.006 MB, total 0.492 MB Allocated vw = 0.012 MB, total 0.504 MB Allocated u-flag = 0.012 MB, total 0.516 MB Checking that no node lies on any fault trace Allocated node = 0.034 MB, total 0.550 MB Allocated a_ = 0.011 MB, total 0.561 MB Allocated crack_index = 0.023 MB, total 0.584 MB Allocated ele_strainrate = 0.034 MB, total 0.618 MB Allocated ele_azim = 0.011 MB, total 0.629 MB Allocated ele_sigma = 0.011 MB, total 0.640 MB Allocated ele_stressed = 0.003 MB, total 0.643 MB Allocated boxed = 0.003 MB, total 0.646 MB 1594 nodes and 2964 elements were read - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Reading velocity boundary conditions from file ..\common\b_PTB3nbd.nki Allocated boundary_node = 0.001 MB, total 0.647 MB Allocated condition = 0.002 MB, total 0.648 MB 222 velocity boundary conditions were read - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Successfully read all input datasets =============================================================================== Finding all data locations in grid coordinates Allocated neighbor = 0.034 MB, total 0.682 MB Allocated center = 0.034 MB, total 0.716 MB Finding geodetic benchmarks Counting fault segments Allocated seg_def = 0.013 MB, total 0.729 MB Allocated seg_end = 0.038 MB, total 0.767 MB Allocated seg_end_is = 0.050 MB, total 0.817 MB Allocated seg_eta_ = 0.006 MB, total 0.823 MB Allocated seg_kappa_ = 0.006 MB, total 0.830 MB Allocated seg_u_ = 0.006 MB, total 0.836 MB Recording fault segments Found all data locations in grid coordinates Count of geodetic benchmarks inside the .feg area: 1504 Allocated covariance_mps2 matrix = 69.031 MB, total 69.867 MB Adding reference-frame-loosening rotations of 10.00 degree/Ma Inverting covariance matrix to get normal matrix Allocated normal matrix = 69.031 MB, total 138.898 MB Counting active cracks Allocated local_crack = 0.069 MB, total 138.967 MB Recording active cracks Delta_node = 1581 Allocated ABCDEF matrix = 76.679 MB, total 215.646 MB Interpolating stress directions Solving for velocities Unfaulted& e1h_azim_error Refinement dV/V RMS(V) Stressed Boxed (L1; in sigmas) 0 1.00000 4.9622E-10 45.11% 0.00% 0.702 1 0.13191 4.8343E-10 45.11% 0.00% 1.872 2 0.01208 4.8230E-10 45.11% 0.00% 1.899 3 0.00536 4.8189E-10 45.11% 0.00% 1.926 4 0.00359 4.8175E-10 45.11% 0.00% 1.931 5 0.00222 4.8165E-10 45.11% 0.00% 1.935 6 0.00138 4.8156E-10 45.11% 0.00% 1.936 7 0.00075 4.8150E-10 45.11% 0.00% 1.937 8 0.00035 4.8147E-10 45.11% 0.00% 1.938 9 0.00020 4.8145E-10 45.11% 0.00% 1.938 10 0.07168 4.7582E-10 45.11% 15.28% 1.939 11 0.00792 4.7541E-10 45.11% 16.53% 0.674 12 0.00253 4.7527E-10 45.11% 16.63% 0.493 13 0.00137 4.7524E-10 45.11% 16.67% 0.475 14 0.00064 4.7524E-10 45.11% 16.67% 0.472 15 0.00047 4.7525E-10 45.11% 16.67% 0.472 16 0.00041 4.7526E-10 45.11% 16.67% 0.472 17 0.00030 4.7527E-10 45.11% 16.67% 0.472 18 0.00025 4.7527E-10 45.11% 16.67% 0.471 19 0.00027 4.7526E-10 45.11% 16.67% 0.471 20 0.00027 4.7525E-10 45.11% 16.67% 0.471 Largest change is 2.85E-12 = 0.090 mm/a at node 1103 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. Correcting geodetic velocities, to best-fit the current model, by adding rotation of 0.0399 deg/Ma about ( 55.676N, 108.073E) Writing g_ptb2007008.nko ============================================================= Continuum errors: L0 = 0.033, L1 = 0.593, L2 = 1.071 Stress errors: L0 = 0.023, L1 = 0.471, L2 = 0.722 Fault errors: L0 = 0.070, L1 = 0.683, L2 = 1.077 Geodetic errors: L0 = 0.396, L1 = 1.981, L2 = 1.683 Global errors: L0 = 0.077, L1 = 0.741, L2 = 1.102 ============================================================= Mean (absolute) value of continuum strain rate = 2.19E-16 Root-mean-square value of continuum strain rate = 3.96E-16 Present mu_ of 3.70E-16 falls in this range. Writing f_ptb2007008.nko Writing h_ptb2007008.nko Writing e_ptb2007008.nko Writing s_ptb2007008.nko Writing v_ptb2007008.out Run ended on 2007. 8.14 at 22:25:32 =============================================================