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 SVN: ------- v. ( ) $Id$ WARNING: find_namelist: namelist "_pars" is missing!! bcx1,bcx2= a a2 s a2 s : a a2 s a2 s bcy1,bcy2= p p p p p : p p p p p bcz1,bcz2= p p p p p : p p p p p lperi= F T T SVN: start.in,v v. 1.1 (2017/04/26 17:31:11) wlyra WARNING: find_namelist: namelist "_pars" is missing!! SVN: run.in,v v. 1.1 (2017/04/26 17:31:11) wlyra bcx1,bcx2= a a2 s a2 s : a a2 s a2 s bcy1,bcy2= p p p p p : p p p p p bcz1,bcz2= p p p p p : p p p p p initialize_mpicomm: enabled MPI WARNING: run: quad precision suppressed! 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$ 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$ SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ The verbose level is ip= 10 (ldebug= F ) This is a 2-D run nxgrid, nygrid, nzgrid= 16 48 1 Lx, Ly, Lz= 2.1000000000000001 6.2831853071795862 2.0000000000000000 Vbox= 19.132299260361840 rsnap: read snapshot var.dat in 1.2542100000001444E-004 seconds rprint_list: nname= 7 rprint_list: ix,iy,iz,iz2= -1 -1 -1 -1 rprint_list: nnamev= 0 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= -1 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_shock: run through parse list rprint_viscosity: run through parse list units_general: unit_velocity= 1.0000000000000000 units_general: unit_density= 1.0000000000000000 units_general: unit_length= 1.0000000000000000 units_general: unit_magnetic= 3.5449077018110318 units_eos: unit_temperature= 1.2027220936797545E-008 units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 -5.9914645471079817 5.0000000000000003E-002 2.5000000000000005E-003 1.0000000000000000 alpha_fine = 7.2973525643000000E-003 sigma_Thomson_cgs = 6.6524587321600003E-025 initialize_grid, r_int,r_ext= 0.40000000000000002 2.5000000000000000 timestep based on CFL cond; cdt= 0.40000000000000002 dVol_rel1= 9.5785440613026837E-004 dA_xy_rel1= 9.5785440613026837E-004 dA_xz_rel1= 6.2500000000000000E-002 dA_yz_rel1= 2.0833333333333332E-002 dL_y_rel1= 2.0833333333333332E-002 initialize_gravity: non-smoothed newtonian gravity xmask_den= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 zmask_den= 1.0000000000000000 diffusion: mesh hyperdiffusion initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. xmask_hyd= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 zmask_hyd= 1.0000000000000000 WARNING: initialize_energy: llocal_iso=T. Make sure you have the appropriate INITIAL_CONDITION in Makefile.local! select_eos_variable: Using rho and cs2 viscous force: nu_hyper3_mesh/pi^5 *(Deltav)^6/Deltaq viscous force: nu_shock*(XXXXXXXXXXX) xmask_vis= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 choose_pencils: finding out which pencils are needed for the pencil case. write_pencil_info: pencil information written to the file pencils.list. Lz= 2.0000000000000000 z= 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 pde: ENTER SVN: ------- v. ( ) $Id$ calc_pencils_hydro: call gij_etc calc_pencils_hydro: max(advec_uu) = 1.2722218725854067E-014 calc_pencils_energy: max(advec_cs2) = 2.5986041829730091 duu_dt: SOLVE Bcs for ux, x: < a>, y: < p>, z: < p> Bcs for uy, x: < a2>, y: < p>, z: < p> Bcs for uz, x: < s>, y: < p>, z: < p> dlnrho_dt: SOLVE Bcs for lnrho, x: < a2>, y: < p>, z: < p> dlnrho_dt: diffrho_hyper3_mesh= 20.000000000000000 dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000 ---it-------t----------dt----------umax-----rhomax-rhomin--nshift-- 0 0.000 1.340E-01 1.868E+00 1.0 1.0 3.017 5 0.670 1.340E-01 1.868E+00 1.0 1.0 3.017 10 1.340 1.339E-01 1.868E+00 1.1 1.0 3.016 15 2.009 1.339E-01 1.868E+00 1.2 0.9 3.014 20 2.678 1.338E-01 1.868E+00 1.2 0.9 3.013 25 3.347 1.338E-01 1.869E+00 1.2 0.9 3.013 30 4.016 1.337E-01 1.868E+00 1.2 0.9 3.012 35 4.685 1.337E-01 1.868E+00 1.3 0.8 3.010 40 5.353 1.337E-01 1.868E+00 1.3 0.8 3.010 45 6.021 1.337E-01 1.868E+00 1.2 0.8 3.010 Simulation finished after 50 time-steps Writing final snapshot at time t = 6.6897194356860323 Wall clock time [hours] = 4.144E-05 (+/- 2.7778E-13) Wall clock time/timestep/meshpoint [microsec] = 3.884830 Wall clock time/timestep/local meshpoint [microsec] = 15.53932 Rhs wall clock time/timestep/local meshpoint [microsec] = 6.346752 After substep wall clock time/timestep/local meshpoint [microsec] = 5.269751 Before+After boundary wall clock time/timestep/local meshpoint [microsec] = 0.8995443 Diagnostics wall clock time/timestep/local meshpoint [microsec] = 0.000000 Timestep wall clock time/timestep/local meshpoint [microsec] = 12.97399 Timestep (without waiting) wall clock time/timestep/local meshpoint [microsec] = 12.97399 Timestep wall clock time/timestep/meshpoint [microsec] = 3.243499 Timestep (without waiting) wall clock time/timestep/meshpoint [microsec] = 3.243499 Diagnostics wall clock time/timestep/meshpoint [microsec] = 0.000000 Rhs wall clock time/timestep/meshpoint [microsec] = 1.586688 After substep wall clock time/timestep/meshpoint [microsec] = 1.317438 Before+After boundary wall clock time/timestep/meshpoint [microsec] = 0.2248861 Maximum used memory per cpu [MBytes] = 35.129 Maximum used memory [MBytes] = 138.684 real 0.48 user 0.68 sys 0.12