Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified Authorization required, but no authorization protocol specified ===================================================================== WARNING: using SINGLE PRECISION, you'd better know what you're doing! ===================================================================== SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ initialize_mpicomm: enabled MPI SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: gravity_r.f90,v v. 1.1 (2018/08/24 15:48:10) wlyra SVN: ------- v. ( ) $Id SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ bcx1,bcx2= a s s a2 a2 s s a : a s s a2 a2 s s a bcy1,bcy2= s a s a2 a2 s s a : s a s a2 a2 s s a bcz1,bcz2= s s a s s s s s : s s a a2 a2 s s a SVN: ------- v. ( ) $Id The verbose level is ip= 13 (ldebug= F ) This is a 3-D run nxgrid, nygrid, nzgrid= 64 64 32 Lx, Ly, Lz= 2.20000005 2.20000005 2.20000005 Vbox= 10.6480007 rsnap: read snapshot var.dat in 3.7175444592834506E-002 seconds rprint_list: nname= 18 rprint_list: ix,iy,iz,iz2= -1 -1 -1 -1 rprint_list: nnamev= 5 rprint_list: nnamez= -1 rprint_list: nnamez= -1 rprint_list: nnamey= -1 rprint_list: nnamex= -1 rprint_list: nnamer= 0 rprint_list: nnamexz= -1 rprint_list: nnamexy= -1 rprint_list: nnamerz= 7 rprint_general: run through parse list rprint_hydro: run through parse list rprint_hydro: run through parse list rprint_density: run through parse list rprint_viscosity: run through parse list setup_slices: slice_position = e setup_slices: ix_loc,iz_loc,iz2_loc, (video files) = 35 4 19 units_general: unit_velocity= 1.0000000000000000 units_general: unit_density= 1.0000000000000000 units_general: unit_length= 1.0000000000000000 units_general: unit_magnetic= 3.5449078083038330 units_eos: unit_temperature= 4.8108880235694808E-009 units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.00000000 1.19209221E-06 0.816497028 0.400000453 0.399999976 alpha_fine = 7.2973525643000000E-003 sigma_Thomson_cgs = 6.6524587321600003E-025 timestep based on CFL cond; cdt= 0.400000006 initialize_gravity: smoothed 1/r potential in spherical shell xmask_den= 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 zmask_den= 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 diffusion: nothing (i.e. no mass diffusion) initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. lT0mu_total is True so setting T00_total and T0i_total to True lT0mu_total is True so setting T00_total and T0i_total to True lT0mu_total is True so setting T00_total and T0i_total to True initialize_hydro: set outer radius rdampext=r_ext xmask_hyd= 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 zmask_hyd= 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 WARNING: initialize_hydro: need to set orms in print.in to get othresh. Disabled othresh! lT0mu_total is True so setting T00_total and T0i_total to True select_eos_variable: Using lnrho and ss initialize_energy: set boundary temperatures for spherical shell problem initialize_energy: g0,mpoly,beta1 1.20000005 1.00000000 1.50000012 initialize_energy: TT_int, TT_ext= 2.50000000 1.00000000 initialize_energy: cs2_ext, cs2_int= 0.666666627 1.66666651 initialize_energy: nheatc_max,iheatcond= 4 K-const nothing nothing nothing heat conduction: K=cte xmask_mag= 0 0.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 0.00000000 ymask_mag= 0 0.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 zmask_mag= 0 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 resistivity: shell viscous force: nu*(del2u+graddivu/3+2S.glnrho) xmask_vis= 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 1.00000000 choose_pencils: finding out which pencils are needed for the pencil case. write_pencil_info: pencil information written to the file pencils.list. acCompile: No changes in input; running make in case DSL code has changed; logging to /home/pencil/GPU/samples/geodynamo/ac_compilation.log Thu Nov 12 13:22:59 2026 : Done setupConfig && acCompile Thu Nov 12 13:22:59 2026 : acGridInit: n[xyz]grid: (64,64,32) n[xyz]local (64,32,16) Thu Nov 12 13:22:59 2026 : acGridInit: Creating device Thu Nov 12 13:23:03 2026 : acGridInit: decomp: (1,2,2) Thu Nov 12 13:23:03 2026 : acGridInit: Stencil order: 6 Thu Nov 12 13:23:12 2026 : acGridInit: Done Started autotuning: watch progress in ac_autotuning.log Thu Nov 12 13:23:45 2026 : DONE initializeGPU WARNING: run: Pencil consistency check not supported on GPUs. You can consider running it with a CPU-only build! pde: ENTER SVN: ------- v. ( ) $Id$ ----it--------t--------dt-------urms------umax------brms------bmax------epsK------epsM-----rhom------ssm------ethm-----dtc-----dtu-----dtb-----dtnu---dtchi---dteta-- 0 0.000 7.08E-03 0.01547 0.02717 4.064E-05 1.331E-04 4.402E-05 5.488E-09 1.2077 0.02953 0.9739 0.9802 0.0198 0.0001 0.0027 0.0037 0.0192 2 0.014 7.09E-03 0.01172 0.04736 3.978E-05 1.329E-04 1.981E-05 5.061E-09 1.2076 0.02954 0.9738 0.9851 0.0355 0.0001 0.0027 0.0037 0.0192 4 0.028 7.09E-03 0.01257 0.06671 3.897E-05 1.327E-04 1.563E-05 4.672E-09 1.2075 0.02958 0.9736 0.9865 0.0450 0.0001 0.0027 0.0038 0.0192 6 0.043 7.10E-03 0.01394 0.07307 3.819E-05 1.325E-04 1.715E-05 4.317E-09 1.2073 0.02966 0.9734 0.9863 0.0487 0.0000 0.0028 0.0038 0.0193 8 0.057 7.11E-03 0.01505 0.07590 3.746E-05 1.323E-04 1.844E-05 3.992E-09 1.2072 0.02977 0.9732 0.9890 0.0568 0.0000 0.0028 0.0038 0.0193 10 0.071 7.12E-03 0.01573 0.07496 3.676E-05 1.320E-04 1.738E-05 3.696E-09 1.2070 0.02990 0.9729 0.9905 0.0561 0.0000 0.0028 0.0038 0.0193 12 0.085 7.13E-03 0.01616 0.07487 3.610E-05 1.317E-04 1.615E-05 3.424E-09 1.2069 0.03005 0.9727 0.9917 0.0556 0.0000 0.0028 0.0038 0.0194 14 0.099 7.14E-03 0.01655 0.07907 3.547E-05 1.314E-04 1.585E-05 3.176E-09 1.2067 0.03021 0.9725 0.9930 0.0583 0.0000 0.0028 0.0039 0.0194 16 0.114 7.15E-03 0.01693 0.07939 3.487E-05 1.311E-04 1.588E-05 2.949E-09 1.2066 0.03037 0.9723 0.9942 0.0590 0.0000 0.0028 0.0039 0.0194 Simulation finished after 17 time-steps Writing final snapshot at time t = 0.12093318637926131 Wall clock time [hours] = 4.006E-02 (+/- 2.7778E-13) Wall clock time/timestep/meshpoint [microsec] = 64.71691 Wall clock time/timestep/local meshpoint [microsec] = 258.8676 Rhs wall clock time/timestep/local meshpoint [microsec] = 5.554025 After substep wall clock time/timestep/local meshpoint [microsec] = 2.0790407E-02 Before+After boundary wall clock time/timestep/local meshpoint [microsec] = 1.387534 Diagnostics wall clock time/timestep/local meshpoint [microsec] = 9.357957 Timestep wall clock time/timestep/local meshpoint [microsec] = 12.39201 Timestep (without waiting) wall clock time/timestep/local meshpoint [microsec] = 6.964512 Timestep wall clock time/timestep/meshpoint [microsec] = 3.098003 Timestep (without waiting) wall clock time/timestep/meshpoint [microsec] = 1.741128 Diagnostics wall clock time/timestep/meshpoint [microsec] = 2.339489 Rhs wall clock time/timestep/meshpoint [microsec] = 1.388506 After substep wall clock time/timestep/meshpoint [microsec] = 5.1976019E-03 Before+After boundary wall clock time/timestep/meshpoint [microsec] = 0.3468834 Maximum used memory per cpu [MBytes] = 1134.211 Maximum used memory [GBytes] = 4.429 real 145.16 user 333.45 sys 40.91