Pencil1: iapp, nprocs, ncpus= 0 4 4 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 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 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 diffusion: mesh hyperdiffusion diffusion: shock diffusion 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) 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 ( 7) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil grho ( 10) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil ee ( 48) is requested, but does not appear to be required! pencil_consistency_check: possible overcalculation... pencil cv1 ( 66) 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.0000000000000000E+000 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.0000000000000000E+000 dlnrho_dt: max(diffus_diffrho3) = 0.0000000000000000E+000 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.0000000000000000E+000 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.0000000000000000E+000 ---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.20E-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.033 0.0212 20 0.384 1.919E-02 1.20E-05 -1.43E-06 9.983E-01 9.986E-01 0.75 1.180E+00 7.144E-03 1.786E-02 1.002E-01 2.510E-02 0.032 0.0214 30 0.576 1.920E-02 1.70E-05 -2.49E-06 9.983E-01 9.984E-01 0.79 1.171E+00 7.144E-03 1.786E-02 1.002E-01 2.509E-02 0.031 0.0215 40 0.768 1.922E-02 2.12E-05 -3.89E-06 9.983E-01 9.982E-01 0.82 1.161E+00 7.144E-03 1.786E-02 1.002E-01 2.509E-02 0.030 0.0217 50 0.960 1.920E-02 2.43E-05 -5.61E-06 9.983E-01 9.980E-01 0.83 1.148E+00 7.144E-03 1.786E-02 1.002E-01 2.508E-02 0.030 0.0218 60 1.152 1.920E-02 2.74E-05 -7.62E-06 9.983E-01 9.979E-01 0.83 1.139E+00 7.144E-03 1.786E-02 1.002E-01 2.508E-02 0.030 0.0220 70 1.344 1.925E-02 2.93E-05 -9.98E-06 9.983E-01 9.978E-01 0.83 1.127E+00 7.144E-03 1.786E-02 1.002E-01 2.508E-02 0.030 0.0222 80 1.537 1.926E-02 3.04E-05 -1.25E-05 9.983E-01 9.977E-01 0.84 1.120E+00 7.143E-03 1.786E-02 1.001E-01 2.507E-02 0.030 0.0222 90 1.729 1.925E-02 3.08E-05 -1.51E-05 9.983E-01 9.976E-01 0.85 1.123E+00 7.143E-03 1.786E-02 1.001E-01 2.507E-02 0.029 0.0223 Simulation finished after 91 time-steps Writing final snapshot at time t = 1.7484486295741770 Wall clock time [hours] = 7.977E-03 (+/- 2.7778E-10) Wall clock time/timestep/meshpoint [microsec] = 19.26138 Maximum used memory per cpu [MBytes] = 15.207 Maximum used memory [MBytes] = 60.328 real 79.29 user 106.19 sys 17.62