Initializing AMReX (26.09-123-g5a08d94fc6ab)... MPI initialized with 4 MPI processes MPI initialized with thread support level 0 AMReX (26.09-123-g5a08d94fc6ab) initialized Successfully read inputs file ... ERF git hash: b0d27c25 AMReX git hash: 26.09-123-g5a08d94fc Problem name (from inputs file) is: "WPS" Reading NC init file names at level 0 and index 0 : wrfinput_chisholmview_d01 Reading NC bdy file name wrfbdy_chisholmview_d01_width1 Level 0 is land Coriolis frequency, f = 0.0001458423212 1/s Thermal stratification based on gradient of potential temperature Reading Times wrfinput datetime 0 : 2013-11-21_06:00:00 1385013600 Have read start_time string at level 0 is 2013-11-21_06:00:00 Have read start_time number at level 0 is 1385013600 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: Kessler Terrain Type: StaticFittedMesh Buildings Type: None Mesh Type: VariableDz ABL Driver Type: Geostrophic Wind (-0,0.002187634818,9.562455688e-11) 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.004999999888 isotropic mixing : 1 Pr_t : 0.3333333433 Sc_t : 1 WRF lateral boundary density forcing: enabled WRF boundary density nudge factor: 10 WRF lateral boundary cloud-water/cloud-ice forcing: disabled 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) (99,99,175) (0,0,0)) ((100,0,0) (199,99,175) (0,0,0)) ((0,100,0) (99,199,175) (0,0,0)) ((100,100,0) (199,199,175) (0,0,0)) ) Reading T00 from wrfinput Reading P00 from wrfinput Reading TLP from wrfinput Reading TISO from wrfinput Reading TLP_STRAT from wrfinput Reading P_STRAT from wrfinput WRF base state parameters (T00, P00, TLP, TISO, TLP_STRAT, P_STRAT) are: (290, 100000, 50, 200, -11, 0) Base state layer interfaces: z_iso = 12901.83398 m Loading initial data from NetCDF file at level 0 Reading from file wrfinput_chisholmview_d01 Will read density from ALB + AL variables Have read (nx,ny,nz) = 200 200 176 Have read (parent_ilo,parent_jlo) = 0 0 Have read refinement ratio = 1 Box in file ((0,0,0) (199,199,175) (0,0,0)) Box to fill ((0,0,0) (199,199,175) (0,0,0)) Checking for ALB ... Reading ALB Checking for AL ... Reading AL Checking for U ... Reading U Checking for V ... Reading V Checking for W ... Reading W Checking for THM ... Reading THM Checking for PH ... Reading PH Checking for PHB ... Reading PHB Checking for PB ... Reading PB Checking for P ... Reading P Checking for PSFC ... Reading PSFC Checking for MUB ... Reading MUB Checking for MAPFAC_U ... Reading MAPFAC_U Checking for MAPFAC_V ... Reading MAPFAC_V Checking for MAPFAC_M ... Reading MAPFAC_M Checking for SST ... Reading SST Checking for TSK ... Reading TSK Checking for LANDMASK ... Reading LANDMASK Checking for C1H ... Reading C1H Checking for C2H ... Reading C2H Checking for RDNW ... Reading RDNW Checking for XLAT_V ... Reading XLAT_V Checking for XLONG_U ... Reading XLONG_U Checking for QVAPOR ... Reading QVAPOR Checking for QCLOUD ... Reading QCLOUD Checking for QRAIN ... Reading QRAIN Checking for IVGTYP ... Reading IVGTYP Checking for ISLTYP ... Reading ISLTYP Terrain has min value = -50.12749481 and max value = 3509.155518 Top of mesh has min value = 15684.13086 and max value = 16560.73047 Max of second-to-highest row = 16382.98828 which is greater than average of top row so defaulting to alternate approach Warning: ProbHi(2) will be ignored; we are setting top of domain to 16471.85938 Constructing nodal heights (z_phys_nd) Nodal reconstruction at k = 0: 50 CG iterations, 0 refinements, regularization 0.0002500000119 nodal heights in [-135.6538086, 3694.466797] m vs WRF z-faces in [-70.00720215, 3612.088867] m max |avg4(nodal) - WRF z-faces| = 35.23657227 m (direct interpolation gives 366.8664551 m) max deviation from direct interpolation = 629.3427124 m (cap 733.7329102 m) Nodal reconstruction at k = 1: 50 CG iterations, 0 refinements, regularization 0.0002500000119 nodal heights in [-116.8156128, 3710.351807] m vs WRF z-faces in [-51.53600693, 3628.003418] m max |avg4(nodal) - WRF z-faces| = 35.2121582 m (direct interpolation gives 366.6989746 m) max deviation from direct interpolation = 629.1079712 m (cap 733.3979492 m) Verifying reconstructed nodal heights Max |avg4(nodal z) - WRF z-face| over the reconstructed levels: 35.23657227 m Min nodal thickness of the first layer: 15.11303711 m