Pencil1: iapp, nprocs, ncpus= 0 4 4 SVN: ------- v. ( ) $Id$ SVN: start.in,v v. 1.1 (2014/02/13 02:24:37) wlyra SVN: run.in,v v. 1.2 (2014/04/15 10:45:53) wlyra bcx1,bcx2= a s s a2 a2 : a s s a2 a2 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 get_snapdir: dir=data 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: baroclinic_run.f90 v. 19193 (2012-06-30 12:55:46) wdobler SVN: ------- v. ( ) $Id$ SVN: ------- v. ( ) $Id$ The verbose level is ip= 9 (ldebug= F ) This is a 2-D run nxgrid, nygrid, nzgrid= 64 64 1 Lx, Ly, Lz= 9.5999999999999996 6.2831840000000003 0.52000000000000002 Vbox= 313.65654528000005 read persistent ID INITIAL_BLOCK_ID rsnap: read snapshot var.dat in 3.5748481750488281E-003 seconds rprint_list: nname= 18 rprint_list: ix,iy,iz,iz2= -1 -1 -1 -1 rprint_list: nnamev= 0 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_viscosity: run through parse list setup_slices: slice_position = p setup_slices: ix_loc,iy_loc,iz_loc, (video files) = 4 4 4 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 -5.0751738152338266 5.0000000000000003E-002 1.7857142857142861E-003 0.28571428571428570 initialize_grid, r_int,r_ext= 0.40000000000000002 10.000000000000000 timestep based on CFL cond; cdt= 0.40000000000000002 WARNING: initialize_diagnostics: volume averages are calculated wrongly for nonequidistant grids unless lproper_averages=T.! dVol_rel1= 4.7695360195360195E-005 dA_xy_rel1= 4.7695360195360195E-005 dA_xz_rel1= 1.5625000000000000E-002 dA_yz_rel1= 1.5625000000000000E-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 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 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. Lz= 0.52000000000000002 xmask_hyd= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 zmask_hyd= 1.0000000000000000 select_eos_variable: Using rho and ss initialize_energy: nheatc_max,iheatcond= 4 hyper3-mesh nothing nothing nothing heat conduction: hyperdiffusivity of ss viscous force: nu_hyper3_mesh/pi^5 *(Deltav)^6/Deltaq xmask_vis= 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 1.0000000000000000 Lz= 0.52000000000000002 z= 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 Lz= 0.52000000000000002 z= 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 z= 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 Lz= 0.52000000000000002 z= 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 choose_pencils: finding out which pencils are needed for the pencil case. write_pencil_info: pencil information written to the file pencils.list. pde: ENTER SVN: ------- v. ( ) $Id$ it= 1 t= 0.0000000000000000E+000 This is a simulation that is restarting from a snapshot. The list of modes will be read from file, to regenerate the mode structure. calc_pencils_energy: max(advec_cs2) = 18.968230746103675 WARNING: calc_pencils_viscosity: viscous heating not implemented for lvisc_hyper3_mesh! duu_dt: SOLVE Bcs for ux, x: < a>, y: < p>, z: < p> Bcs for uy, x: < s>, y: < p>, z: < p> Bcs for uz, x: < s>, y: < p>, z: < p> duu_dt: max(advec_uu) = 40.209417734726827 dlnrho_dt: SOLVE Bcs for lnrho, x: < a2>, y: < p>, z: < p> 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: < a2>, y: < p>, z: < p> denergy_dt: lnTT,cs2,cp1= -4.1588830833596715 6.2500000000000012E-003 0.0000000000000000E+000 calc_heatcond_hyper3_mesh: chi_hyper3= 0.0000000000000000E+000 calc_heatcond_hyper3: added thdiff calc_heat_cool: lgravz, lgravr, lgravx, lspherical_coords= F T F F dspecial_dt: SOLVE dspecial_dt ---it-------t----------dt----------umax--------rhomin--------rhom--------rhomax--------rho2m--------oz2m----------ozm---------uxuym--------urms---------orms-------potturb2m---potturbmax---potturbmin---gpotturbx2m--gpotturby2m- 0 0.0000 8.672E-03 1.579E+00 3.162E-01 3.024E-01 1.581E+00 1.923E-01 5.207E-02 7.631E-02 0.000E+00 4.379E-01 2.282E-01 2.396E-08 2.284E-03 -2.284E-03 5.881E-06 1.652E-07 10 0.0853 8.399E-03 1.579E+00 3.162E-01 3.024E-01 1.579E+00 1.923E-01 5.209E-02 7.632E-02 1.455E-05 4.379E-01 2.282E-01 2.395E-08 2.287E-03 -2.287E-03 5.895E-06 1.652E-07 20 0.1686 8.270E-03 1.578E+00 3.162E-01 3.024E-01 1.564E+00 1.923E-01 5.210E-02 7.633E-02 4.851E-05 4.379E-01 2.283E-01 2.394E-08 2.288E-03 -2.288E-03 5.908E-06 1.652E-07 30 0.2510 8.202E-03 1.578E+00 3.162E-01 3.024E-01 1.537E+00 1.923E-01 5.211E-02 7.634E-02 9.135E-05 4.379E-01 2.283E-01 2.394E-08 2.288E-03 -2.288E-03 5.921E-06 1.652E-07 40 0.3328 8.165E-03 1.578E+00 3.162E-01 3.024E-01 1.542E+00 1.922E-01 5.212E-02 7.634E-02 1.362E-04 4.379E-01 2.283E-01 2.393E-08 2.288E-03 -2.288E-03 5.932E-06 1.652E-07 50 0.4144 8.145E-03 1.577E+00 3.162E-01 3.024E-01 1.555E+00 1.922E-01 5.211E-02 7.634E-02 1.790E-04 4.379E-01 2.283E-01 2.393E-08 2.287E-03 -2.287E-03 5.943E-06 1.653E-07 60 0.4958 8.132E-03 1.577E+00 3.162E-01 3.024E-01 1.562E+00 1.921E-01 5.206E-02 7.634E-02 2.177E-04 4.378E-01 2.282E-01 2.392E-08 2.286E-03 -2.286E-03 5.953E-06 1.653E-07 70 0.5771 8.124E-03 1.577E+00 3.162E-01 3.024E-01 1.563E+00 1.921E-01 5.200E-02 7.633E-02 2.514E-04 4.378E-01 2.280E-01 2.392E-08 2.284E-03 -2.284E-03 5.962E-06 1.653E-07 80 0.6583 8.118E-03 1.577E+00 3.162E-01 3.023E-01 1.563E+00 1.920E-01 5.191E-02 7.631E-02 2.802E-04 4.378E-01 2.278E-01 2.391E-08 2.282E-03 -2.282E-03 5.969E-06 1.654E-07 90 0.7394 8.113E-03 1.577E+00 3.162E-01 3.023E-01 1.559E+00 1.920E-01 5.180E-02 7.629E-02 3.043E-04 4.377E-01 2.276E-01 2.391E-08 2.280E-03 -2.280E-03 5.976E-06 1.654E-07 Simulation finished after 91 time-steps Writing final snapshot at time t = 0.74753400719229035 Wall clock time [hours] = 2.863E-03 (+/- 2.7778E-10) Wall clock time/timestep/meshpoint [microsec] = 27.65039 Maximum used memory per cpu [MBytes] = 10.512 Maximum used memory [MBytes] = 39.496 real 12.77 user 14.90 sys 2.36