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
rsnap: read snapshot var.dat in 0.0000000000000000 seconds
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.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000
alpha_fine = 7.2973525643000000E-003
sigma_Thomson_cgs = 6.6524587321600003E-025
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: minmax(ad)= 0.62035049089940009 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
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
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: 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.0000000000000000
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.0000000000000000
---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.053E-04 (+/- 2.7778E-07)
Wall clock time/timestep/meshpoint [microsec] = 18.62536
Wall clock time/timestep/local meshpoint [microsec] = 18.62536
Rhs wall clock time/timestep/local meshpoint [microsec] = 58.34318
Maximum used memory per cpu [MBytes] = 9.012
Maximum used memory [MBytes] = 9.012
real 11.53
user 10.88
sys 0.12