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: "RICO" Level 0 is land Coriolis frequency, f = 3.76e-05 1/s Thermal stratification based on gradient of potential temperature Maximum simulation length based on stop_time: 3600 s (elapsed) Adding stop time 3600 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 : 1 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 : Upwind_3rd moistscal_vert_adv_type : Upwind_3rd 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 : 0.3333333333 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,24) (0,0,0)) ((16,0,0) (31,15,24) (0,0,0)) ((0,16,0) (15,31,24) (0,0,0)) ((16,16,0) (31,31,24) (0,0,0)) ) max of zphys_nd 4000 max of zlevels 4000 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.15248101072114 kg/m^3 WARNING: HSE Newton iterations did not converge to tolerance! HSE Newton tol: 1.455192e-11 1.000000e-12 WARNING: HSE Newton iterations did not converge to tolerance! HSE Newton tol: 1.455192e-11 1.000000e-12 WARNING: HSE Newton iterations did not converge to tolerance! HSE Newton tol: 1.455192e-11 1.000000e-12 WARNING: HSE Newton iterations did not converge to tolerance! HSE Newton tol: 1.455192e-11 1.000000e-12 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, 12800] x [0, 12800] x [0, 4000] 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 12800 0 12800 0 4000 ) Multiplicity type: sampled Particles per cell: 1 Vapour/Condensate Species material: H2O (distribution: mass_constant, value=4.1887902e-42) Aerosols material: NH4HSO4 (distribution: mass_constant, value=2.77260785e-17) Number of physical particles per cell: 1 Number of super droplets per cell: 1 Adding particles in box with volume 6.5536e+11. Particle container size (level 0): 25600 Total number of superdroplets per cell: 1 Initialized 25600 super-droplets at level 0 in super-droplets moisture model. Based on cfl of one Compressible dt at level 0 would be: 1.129064034 Anelastic dt at level 0 would be: 41.06430732 Fixed dt at level 0 is: 2.5 smallest even ratio is: 4 Using specified ustar, tstar, qstar for MOST = 0.001229 0.001094 0.001133 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.77260785e-17, 2.77260785e-17, 2.77260785e-17 radius [m]: 1.877586467e-07, 1.877586467e-07, 1.877586467e-07 multiplicity: 0.7961281402, 1.145500115, 0.9679862913 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]: NH4HSO4: 2.77260785e-17, 2.77260785e-17, 2.77260785e-17 Writing super_droplets_moisture_g_lnR_00000.txt Max value of |p_hse - p_eos| is 2226.564233 with max in cell (15,15,0)