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 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 "special_init_pars" is missing!! SVN: start.in,v v. 1.1 (2009-09-12 16:05:39) wlyra SVN: run.in,v v. 1.2 (2009-09-24 08:59:37) wlyra 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: ------- 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$ SVN: ------- v. ( ) $Id$ The verbose level is ip= 14 (ldebug= F ) This is a 3-D run nxgrid, nygrid, nzgrid= 32 32 16 Lx, Ly, Lz= 0.40000000000000002 1.6000000000000001 0.20000000000000001 Vbox= 0.12800000000000003 rsnap: read snapshot var.dat in 4.8361100000000000E-004 seconds setup_slices: slice_position = p setup_slices: ix_loc,iy_loc,iz_loc,iz2_loc, (video files) = 4 4 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.5449077018110318 units_eos: unit_temperature= 3.4363488390850127E-009 units_eos: cp, lnTT0, cs0, pp0, Rgas= 1.0000000000000000 -3.6888794541139358 0.10000000000000001 7.1428571428571461E-003 0.28571428571428570 alpha_fine = 7.2973525643000000E-003 sigma_Thomson_cgs = 6.6524587321600003E-025 diffusion: mesh hyperdiffusion diffusion: shock diffusion WARNING: initialize_density: For diffusion energy/momentum correction should use lmassdiff_fix=T! initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. select_eos_variable: Using rho and ss initialize_energy: nheatc_max,iheatcond= 4 hyper3-mesh shock nothing nothing heat conduction: hyperdiffusivity of ss heat conduction: shock viscous force: nu_hyper3_mesh/pi^5 *(Deltav)^6/Deltaq viscous force: nu_shock*(XXXXXXXXXXX) sh: 1: Syntax error: Bad fd number pencil_consistency_check: checking pencil case pencil_consistency_check: checking requested pencils pencil_consistency_check: performing small pencil check pencil_consistency_check: the small pencil check reported no problems pencil_consistency_check: performing full pencil check (takes a while) pencil_consistency_check: possible overcalculation... pencil rho ( 13) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil grho ( 16) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil ee ( 55) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cv1 ( 73) is requested, but does not appear to be required! pencil_consistency_check: checking dependence on pencil initialization pencil_consistency_check: results are independent of pencil initialization pencil_consistency_check: checking diagnostics pencils pencil_consistency_check: performing small pencil check pencil_consistency_check: the small pencil check reported no problems pencil_consistency_check: performing full pencil check (takes a while) pencil_consistency_check: set lpencil_check_diagnos_opti=T in run.in to report pencil_consistency_check: optimisation possibilities for diagnostics pencils pencil_consistency_check: checking dependence on pencil initialization pencil_consistency_check: diagnostics are independent of pencil initialization pencil_consistency_check: all tests passed pde: ENTER SVN: ------- v. ( ) $Id$ calc_pencils_energy: max(advec_cs2) = 139.18980252797598 WARNING: calc_pencils_viscosity: viscous heating not implemented for lvisc_hyper3_mesh! duu_dt: SOLVE Bcs for ux, x: < she>, y: < p>, z: < p> Bcs for uy, x: < she>, y: < p>, z: < p> Bcs for uz, x: < she>, y: < p>, z: < p> coriolis_cartesian: add Coriolis force; Omega= 1.0000000000000000 duu_dt: max(advec_uu) = 0.0000000000000000 dlnrho_dt: SOLVE Bcs for lnrho, x: < she>, y: < p>, z: < p> dlnrho_dt: diffrho_shock= 1.0000000000000000 dlnrho_dt: diffrho_hyper3_mesh= 5.0000000000000000 dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000 denergy_dt: SOLVE denergy_dt Bcs for ss, x: < she>, y: < p>, z: < p> denergy_dt: lnTT,cs2,cp1= -3.6946269386926391 9.9426900061287231E-003 1.1728390000000000 calc_heatcond_shock: use shock diffusion calc_heatcond_shock: chi_shock= 1.0000000000000000 calc_heatcond_shock: added thdiff calc_heatcond_hyper3_mesh: chi_hyper3= 0.0000000000000000 calc_heatcond_hyper3: added thdiff calc_heat_cool: lgravz, lgravr, lgravx, lspherical_coords= F F F F dspecial_dt: SOLVE dSPECIAL_dt shearing: Sshear,Sshear1= -1.5000000000000000 -1.5000000000000000 shearing: qshear,qshear0= 1.5000000000000000 0.0000000000000000 ---it-------t--------dt-------uxm-------uym-------rhom-------rho2m---rhomin---rhomax-------ppm-------ethm--------csm--------TTm-----TTmax--TTmin-- 0 0.000 1.921E-02 0.00E+00 0.00E+00 9.983E-01 9.991E-01 0.73 1.189E+00 7.143E-03 1.786E-02 1.002E-01 2.511E-02 0.034 0.0210 10 0.192 1.920E-02 8.15E-06 -5.18E-07 9.983E-01 9.989E-01 0.73 1.186E+00 7.143E-03 1.786E-02 1.002E-01 2.510E-02 0.034 0.0212 20 0.384 1.919E-02 1.20E-05 -1.42E-06 9.983E-01 9.986E-01 0.75 1.181E+00 7.144E-03 1.786E-02 1.002E-01 2.510E-02 0.033 0.0213 30 0.576 1.920E-02 1.70E-05 -2.49E-06 9.983E-01 9.984E-01 0.78 1.174E+00 7.144E-03 1.786E-02 1.002E-01 2.509E-02 0.032 0.0215 40 0.768 1.921E-02 2.13E-05 -3.87E-06 9.983E-01 9.982E-01 0.81 1.164E+00 7.144E-03 1.786E-02 1.002E-01 2.509E-02 0.030 0.0216 50 0.960 1.920E-02 2.43E-05 -5.59E-06 9.983E-01 9.981E-01 0.82 1.151E+00 7.144E-03 1.786E-02 1.002E-01 2.509E-02 0.031 0.0217 60 1.152 1.920E-02 2.75E-05 -7.61E-06 9.983E-01 9.980E-01 0.82 1.142E+00 7.144E-03 1.786E-02 1.002E-01 2.508E-02 0.031 0.0219 70 1.344 1.924E-02 2.94E-05 -9.98E-06 9.983E-01 9.978E-01 0.82 1.130E+00 7.144E-03 1.786E-02 1.002E-01 2.508E-02 0.030 0.0221 80 1.537 1.926E-02 3.05E-05 -1.25E-05 9.983E-01 9.977E-01 0.83 1.123E+00 7.144E-03 1.786E-02 1.002E-01 2.508E-02 0.030 0.0222 90 1.729 1.925E-02 3.10E-05 -1.51E-05 9.983E-01 9.977E-01 0.84 1.126E+00 7.144E-03 1.786E-02 1.001E-01 2.507E-02 0.029 0.0222 Simulation finished after 91 time-steps Writing final snapshot at time t = 1.7483417219707393 Wall clock time [hours] = 1.147E-03 (+/- 2.7778E-13) Wall clock time/timestep/meshpoint [microsec] = 2.769996 Wall clock time/timestep/local meshpoint [microsec] = 11.07998 Rhs wall clock time/timestep/local meshpoint [microsec] = 22.31000 Maximum used memory per cpu [MBytes] = 31.621 Maximum used memory [MBytes] = 124.531 real 12.57 user 47.58 sys 1.25