SVN: ------- v. ( ) $Id$ SVN: start.in,v v. 1.2 (2013-03-22 03:42:03) wlyra SVN: run.in,v v. 1.2 (2013-03-22 03:42:03) wlyra SVN: ------- v. ( ) $Id$ SVN: centrifugal_balanc v. 19193 (2012-06-30 12:55:46) wdobler 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= 16 16 8 Lx, Ly, Lz= 1.0000000000000000 1.0000000000000000 0.59999999999999998 Vbox= 0.59999999999999998 setup_slices: slice_position = p setup_slices: ix_loc,iy_loc,iz_loc,iz2_loc, (video files) = 4 4 4 11 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 0.0000000000000000E+000 1.0000000000000000 1.0000000000000000 1.0000000000000000 initialize_gravity: no x-gravity initialize_gravity: no y-gravity initialize_gravity: linear z-grav, nu= 1.0000000000000000 diffusion: (d^6/dx^6+d^6/dy^6+d^6/dz^6)lnrho initialize_density: no need to read initial stratification for lanti_shockdiffusion=F. select_eos_variable: Using lnrho and cs2 initialize_dustvelocity: dust_chemistry = nothing recalculated: md0= 1.0000000000000000 initialize_dustvelocity: ad= 0.62035049089940009 initialize_dustvelocity: No betad calculation for draglaw=epstein_cst initialize_dustvelocity: dust geometry = sphere Viscous force (dust): nud*del6ud Viscous force (dust): nud*del2ud initialize_dustdensity: ldustcoagulation,ldustcondensation = F F initialize_dustdensity: ldustcoagulation_simplified,ldustcondensation_simplified = F F advec_ddensity: plain vanilla scheme initialize_dustdensity: diffnd_ndustspec= 0.0000000000000000E+000 dust diffusion: (d^6/dx^6+d^6/dy^6+d^6/dz^6)lnnd dust diffusion: [div(DT*grad(nd))] initialize_shock: prenormalised shock_factor sum= 1.0000000000000000 viscous force: nu_hyper*del6v viscous force: nu*del2v 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: 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) = 648.11111111111109 duu_dt: SOLVE Bcs for ux, x: < she>, y: < p>, z: < a> Bcs for uy, x: < she>, y: < p>, z: < a> Bcs for uz, x: < she>, y: < p>, z: < out> coriolis_cartesian: add Coriolis force; Omega= 1.0000000000000000 dss_dt: max(advec_cs2) = 972.16666666666663 duu_dt: max(advec_uu) = 0.15373795466934773 dlnrho_dt: SOLVE Bcs for lnrho, x: < she>, y: < p>, z: < v3> dlnrho_dt: diffrho_hyper3= 2.9999999999999998E-015 dlnrho_dt: max(diffus_diffrho ) = 0.0000000000000000E+000 dlnrho_dt: max(diffus_diffrho3) = 1.0822817511522633E-007 duud_dt: SOLVE duud_dt Bcs for udx, x: < she>, y: < p>, z: < a> Bcs for udy, x: < she>, y: < p>, z: < a> Bcs for udz, x: < she>, y: < p>, z: < set:set> duud_dt: add Coriolis force; Omega= 1.0000000000000000 dndmd_dt: SOLVE dnd_dt, dmd_dt, dmi_dt Bcs for nd, x: < she>, y: < p>, z: < v3> dspecial_dt: SOLVE dspecial_dt shearing: Sshear,Sshear1= -1.5000000000000000 -1.5000000000000000 shearing: qshear,qshear0= 1.5000000000000000 0.0000000000000000E+000 ---it-------t----------dt----------umax--------rhomin--------rhom--------rhomax--------rho2m--------urms----------ndm---------ndmin--------ndmax--------nd2m---------rhodm-------rhodmin------rhodmax-------epsdm-------epsdmin------epsdmax--- 0 0.0000 9.876E-03 4.168E-02 1.549E+00 1.790E+00 1.917E+00 3.225E+00 1.537E-02 4.597E-02 2.195E-05 1.480E-01 5.766E-03 4.597E-02 2.195E-05 1.480E-01 2.871E-02 1.192E-05 9.551E-02 25 0.2466 9.927E-03 4.321E-02 1.479E+00 1.784E+00 1.968E+00 3.207E+00 1.788E-02 4.604E-02 2.132E-05 1.501E-01 5.792E-03 4.604E-02 2.132E-05 1.501E-01 2.915E-02 1.162E-05 1.013E-01 50 0.4951 9.908E-03 5.826E-02 1.464E+00 1.791E+00 2.096E+00 3.239E+00 2.079E-02 4.625E-02 2.149E-05 1.571E-01 5.896E-03 4.625E-02 2.149E-05 1.571E-01 2.961E-02 1.126E-05 1.057E-01 75 0.7434 9.947E-03 6.895E-02 1.398E+00 1.795E+00 2.146E+00 3.270E+00 1.966E-02 4.665E-02 2.295E-05 1.687E-01 6.006E-03 4.665E-02 2.295E-05 1.687E-01 3.082E-02 1.160E-05 1.200E-01 Simulation finished after 76 time-steps Writing final snapshot at time t = 0.75339681323679542 Wall clock time [hours] = 8.294E-04 (+/- 2.7778E-07) Wall clock time/timestep/meshpoint [microsec] = 19.18431 Maximum used memory per cpu [MBytes] = 8.992 Maximum used memory [MBytes] = 8.992 real 12.37 user 10.95 sys 0.08