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 ... Fixing the random seed ERF git hash: 74cc55ec-dirty AMReX git hash: 26.06-15-g07311ca512 Problem name (from inputs file) is: "SinusoidalMassFlux" Level 0 is land 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 : 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 : WENO5 moistscal_vert_adv_type : WENO5 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) (1,1,49) (0,0,0)) ((2,0,0) (3,1,49) (0,0,0)) ((0,2,0) (1,3,49) (0,0,0)) ((2,2,0) (3,3,49) (0,0,0)) ) max of zphys_nd 2975 max of zlevels 2975 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 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: 100 cloud/rain radius: 4e-05 [m] kinematic mode: false dimensionality: one_d_z include phase change: true include particle advection: true include coalescence: false 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, 400] x [0, 400] x [0, 2975] Advect with flow: true Advect with gravity: false Prescribed advection: true Random initial placement: false Coalescence bin size: (1,1,1) Include Brownian coaslescence: false Coalescence kernel: sedimentation 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: 50000000 Initial super-droplets number density: -1 Particle box: (RealBox 0 400 0 400 0 2975 ) Multiplicity type: sampled Particles per cell: 128 Vapour/Condensate Species material: H2O (distribution: mass_constant, value=4.1887902e-42) Aerosols material: NH42SO4 (distribution: radius_log_normal, min=1e-09, max=5e-06, mean=8e-08, std=1.157196188) Number of physical particles per cell: 2.975e+13 Number of super droplets per cell: 128 Adding particles in box with volume 476000000. Particle container size (level 0): 102400 Total number of superdroplets per cell: 128 Initialized 102400 super-droplets at level 0 in super-droplets moisture model. Based on cfl of one Compressible dt at level 0 would be: 0.2900277086 Anelastic dt at level 0 would be undefined Fixed dt at level 0 is: 0.25 Fixed fast dt at level 0 is: 0.25 smallest even ratio is: 2 SuperDropletPC(super_droplets_moisture) attributes (min, max, avg): mass [kg]: 7.420163307e-24, 9.267419255e-13, 4.183945705e-18 radius [m]: 1.209970972e-09, 6.048161795e-06, 9.767178499e-08 multiplicity: 0.879715783, 1.6018481e+13, 2.357893625e+11 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: 7.420163307e-24, 9.267419255e-13, 4.183945705e-18 Writing super_droplets_moisture_g_lnR_00000.txt Min/max value of x-gradient of base state pressure are zero Min/max value of y-gradient of base state pressure are zero Min/max value of dp0/dz + rho0*|g| are -6.892264537e-13 1.02851061e-12 with min at face (1,1,29) with max at face (1,1,22) Min/max value of x-gradient of full (moist) pressure are zero Min/max value of y-gradient of full (moist) pressure are zero Min/max value of moist dp/dz + rho_m*|g| are -6.892264537e-13 1.02851061e-12 with min at face (1,1,29) with max at face (1,1,22) ParticleData: the following additional Eulerian variables are available to plot: super_droplets_moisture_mass_flux_x super_droplets_moisture_mass_flux_y super_droplets_moisture_mass_flux_z super_droplets_moisture_number_density super_droplets_moisture_sd_number_density super_droplets_moisture_mass_density super_droplets_moisture_radius super_droplets_moisture_mass_density_H2O super_droplets_moisture_mass_flux_x_H2O super_droplets_moisture_mass_flux_y_H2O super_droplets_moisture_mass_flux_z_H2O super_droplets_moisture_aerosol_mass_density_NH42SO4 super_droplets_moisture_aerosol_mass_flux_x_NH42SO4 super_droplets_moisture_aerosol_mass_flux_y_NH42SO4 super_droplets_moisture_aerosol_mass_flux_z_NH42SO4 super_droplets_moisture_species_mass_density_H2O super_droplets_moisture_species_mass_flux_x_H2O super_droplets_moisture_species_mass_flux_y_H2O super_droplets_moisture_species_mass_flux_z_H2O SuperDroplets: the following additional variables are available to plot: accum_NH42SO4 Writing native 3D plotfile SDM_SineMassFlux_plt00000 3DPlotfile write time = 0.019588456 seconds. TIME= 0 MASS = 482896614.4 RHO THETA = 1.460647369e+11 RHO SCALAR = 0 RHO QTOTAL = 5104305.629 Coarse STEP 1 starts ... Based on cfl of one Compressible dt at level 0 would be: 0.2900277086 Anelastic dt at level 0 would be undefined Fixed dt at level 0 is: 0.25 Fixed fast dt at level 0 is: 0.25 smallest even ratio is: 2 [Level 0 step 1] ADVANCE from elapsed time = 0 to 0.25 with dt = 0.25 Making slow rhs at time 0 for fast variables advancing from 0 to 0.08333333333333