Initializing AMReX (26.06-15-g07311ca5127b)... MPI initialized with 4 MPI processes MPI initialized with thread support level 0 AMReX (26.06-15-g07311ca5127b) initialized Successfully read inputs file ... ERF git hash: 339cea60-dirty AMReX git hash: 26.06-15-g07311ca512 Problem name (from inputs file) is: "Density Current" Level 0 is land Level 1 is land Using constant kinematic diffusion coefficients momentum : 75 m^2/s temperature : 75 m^2/s scalar : 0 m^2/s Thermal stratification based on gradient of potential temperature Thermal stratification based on gradient of potential temperature Turning off native vertical implicit solve from vert_implicit flag. Maximum simulation length based on stop_time: 900 s (elapsed) Adding stop time 900 to start_time 0 Null land surface model! SOLVER CHOICE: At level 0 : compressible with implicit substepping At level 1 : anelastic with no substepping and fixed density vert_implicit_fac : 0 0 0 use_coriolis : 0 use_gravity : 1 Moisture Model: None Terrain Type: StaticFittedMesh Buildings Type: None Mesh Type: StretchedDz ABL Driver Type: None Coupling Type: One-way Radiation Model: None Gradp_type : 0 Buoyancy_type at level 0 : 2 Buoyancy_type at level 1 : 3 Advection Choices: dycore_horiz_adv_type : Upwind_3rd dycore_vert_adv_type : Upwind_3rd dryscal_horiz_adv_type : Upwind_3rd dryscal_vert_adv_type : Upwind_3rd moistscal_horiz_adv_type : WENO3 moistscal_vert_adv_type : WENO3 Diffusion choices: rho0_trans : 1 alpha_T : 75 alpha_C : 0 dynamic_viscosity : 75 Sponge choices: Rayleigh damping : None w damping : 0 (coeff = -1 * dz/dt**2) Turbulence Settings at level 0 Using DNS model at level 0 Turbulence Settings at level 1 Using DNS model at level 1 Refinement ratio at level 1 set to be 2 1 2 Creating new distribution map on level: 0 BA FROM SCRATCH AT LEVEL 0 (BoxArray maxbox(4) m_ref->m_hash_sig(0) ((0,0,0) (63,0,63) (0,0,0)) ((64,0,0) (127,0,63) (0,0,0)) ((128,0,0) (191,0,63) (0,0,0)) ((192,0,0) (255,0,63) (0,0,0)) ) max of zphys_nd 6400 max of zlevels 6400 Problem name "density current" does not add any velocity perturbations. Creating new distribution map on level: 1 BA FROM SCRATCH AT LEVEL 1 (BoxArray maxbox(5) m_ref->m_hash_sig(0) ((0,0,0) (15,0,111) (0,0,0)) ((16,0,0) (31,0,111) (0,0,0)) ((32,0,0) (47,0,111) (0,0,0)) ((48,0,0) (63,0,111) (0,0,0)) ((64,0,0) (79,0,111) (0,0,0)) ) Building the hybrid FFT solver to be used on domain ((0,0,0) (79,0,111) (0,0,0)) Problem name "density current" does not add any velocity perturbations. Problem name "density current" does not add any velocity perturbations. Projecting initial velocity field at level 1 Max/L2 norm of divergence before solve in subdomain 0 at level 1 : 0 0 and volume-weighted sum 0 Subtracting 0 from rhs in subdomain 0 Sum after subtraction 0 in subdomain 0 Based on cfl of one Compressible dt at level 0 would be: 0.2882620166 Anelastic dt at level 0 would be undefined Fixed dt at level 0 is: 1 Fixed fast dt at level 0 is: 0.25 smallest even ratio is: 4 Based on cfl of one Compressible dt at level 1 would be: 0.1440724779 Anelastic dt at level 1 would be undefined Fixed dt at level 1 is: 1 smallest grid spacing at level 1 = 50 dt at level 1 = 1 Viscous CFL is 0.015 Min/max value of dp0/dx are zero Min/max value of dp0/dy are zero Min/max value of dp0/dz + rho0*|g| are less than 1e-10 Min/max value of dp0/dx are zero Min/max value of dp0/dy are zero Min/max value of dp0/dz + rho0*|g| are -6.888949067 1.456612608e-13 with min at face (79,0,112) with max at face (79,0,59)