Initializing AMReX (26.06-26-ge95d50b2b424)... MPI initialized with 4 MPI processes MPI initialized with thread support level 0 AMReX (26.06-26-ge95d50b2b424) initialized Successfully read inputs file ... Fixing the random seed ERF git hash: 38449378-dirty AMReX git hash: 26.06-26-ge95d50b2b4 Problem name (from inputs file) is: "SDM_Congestus3D" Level 0 is land Thermal stratification based on gradient of potential temperature Maximum simulation length based on stop_time: 8500 s (elapsed) Adding stop time 8500 to start_time 0 Null land surface model! SOLVER CHOICE: At level 0 : compressible with implicit substepping vert_implicit_fac at level : 0 1 1 0 (theta = 1, moisture = 1, tke = 1, momenta = 1) use_coriolis : 0 use_gravity : 1 Moisture Model: None Terrain Type: None Buildings Type: None Mesh Type: ConstantDz ABL Driver Type: None Radiation Model: None Gradp_type : 0 Buoyancy_type at level 0 : 1 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 : 0 alpha_C : 0 dynamic_viscosity : 0 Not using any molecular diffusivity, i.e. using the modeled turbulent diffusivity Sponge choices: Rayleigh damping : None w damping : 0 (coeff = -1 * dz/dt**2) Turbulence Settings at level 0 Using Smagorinsky LES model at level 0 Smagorinsky uses Richardson number correction with Ri_crit = 0.25 Cs : 0.17 isotropic mixing : 1 Pr_t : 0.3333333333 Sc_t : 1 Super-Droplets moisture model! 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) (15,15,49) (0,0,0)) ((16,0,0) (31,15,49) (0,0,0)) ((0,16,0) (15,31,49) (0,0,0)) ((16,16,0) (31,31,49) (0,0,0)) ) max of zphys_nd 10000 max of zlevels 10000 input_sounding file location : input_sounding Successfully opened the input_sounding file. Now reading... Successfully read the 0th input_sounding file... ideal sounding init: surface dry air density = 1.14096930957446 kg/m^3 Note: using stable Redistribute() for particle container super_droplets_moisture. SuperDropletPC(super_droplets_moisture): added 1 int-type attribute(s). SuperDropletPC(super_droplets_moisture): added 6 real-type attribute(s). Using fixed seed for SuperDropletPC random engine. SuperDropletsMoist: diagnostics_interval: 1 cloud/rain radius: 4e-05 [m] kinematic mode: false dimensionality: three_d include phase change: true include particle advection: true include coalescence: true Recycle particles: false number of substeps (phase change): 1 initial phase change relaxation: false SuperDropletPC(super_droplets_moisture): Density scaling: true Nucleate particles: false Inactive particles threshold: 0.01 Recycling bounding box: [0, 6400] x [0, 6400] x [0, 10000] Advect with flow: true Advect with gravity: true Prescribed advection: false Random initial placement: false Coalescence bin size: (1,1,1) Include Brownian coaslescence: false Coalescence kernel: Halls Mass change time integrator: backward_euler (cfl = 1000) Terminal velocity model: CloudRainShima SuperDropletPC(super_droplets_moisture) InitializeParticles at level 0: SuperDropletPC(super_droplets_moisture) Initialization: Initial number density: -1 Initial super-droplets number density: -1 Particle box: (RealBox 0 6400 0 6400 0 10000 ) Multiplicity type: sampled Particles per cell: 1 Vapour/Condensate Species material: H2O (distribution: mass_constant, value=4.1887902e-42) Aerosols material: NH42SO4 (distribution: mass_constant, value=2.649201195e-17) Number of physical particles per cell: 1 Number of super droplets per cell: 1 Adding particles in box with volume 4.096e+11. Particle container size (level 0): 51200 Total number of superdroplets per cell: 1 Initialized 51200 super-droplets at level 0 in super-droplets moisture model. Based on cfl of one Compressible dt at level 0 would be: 0.5775405593 Anelastic dt at level 0 would be: 2009.677145 Fixed dt at level 0 is: 0.5 Fixed fast dt at level 0 is: 0.1 smallest even ratio is: 6 Using specified ustar, tstar, qstar for MOST = 0.28 0 0 If there is water, Charnock relation with C_a=0.0185 will be used SuperDropletPC(super_droplets_moisture) attributes (min, max, avg): mass [kg]: 2.649201195e-17, 2.649201195e-17, 2.649201195e-17 radius [m]: 1.849305995e-07, 1.849305995e-07, 1.849305995e-07 multiplicity: 0.4229960445, 1.133256556, 0.7376726678 velocity components [m/s]: x: 0, 0, 0 y: 0, 0, 0 z: 0, 0, 0 terminal velocity [m/s]: 0, 0, 0 species masses [kg]: H2O: 4.1887902e-42, 4.1887902e-42, 4.1887902e-42 aerosol masses [kg]: NH42SO4: 2.649201195e-17, 2.649201195e-17, 2.649201195e-17 Writing super_droplets_moisture_g_lnR_00000.txt Max value of |p_hse - p_eos| is 2367.317186 with max in cell (15,15,0)