Initializing AMReX (26.10-76-g347a3a2d2675)... MPI initialized with 4 MPI processes MPI initialized with thread support level 0 OMP initialized with 2 OMP threads AMReX (26.10-76-g347a3a2d2675) initialized Successfully read inputs file ... ERF git hash: ec7123d9b AMReX git hash: 26.10-76-g347a3a2d2 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.000145842318 1/s PBL thermal stratification variable: thetav 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 Flat fitted terrain: false Buildings Type: None Mesh Type: VariableDz ABL Driver Type: Geostrophic Wind (-0,0.00218763477,-1.339539959e-19) 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.005 isotropic mixing : 1 Pr_t : 0.3333333333 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)) ((0,100,0) (99,199,175) (0,0,0)) ((100,0,0) (199,99,175) (0,0,0)) ((100,100,0) (199,199,175) (0,0,0)) ) Found MOST at face (2) : Constructing primitive vars for MOST.. 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.83486 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 ... Reading XLAT Checking for XLONG ... Reading XLONG 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.12749566 and max value = 3509.155697 Top of mesh has min value = 15684.13091 and max value = 16560.73164 Max of second-to-highest row = 16382.9893 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.86047 Constructing nodal heights (z_phys_nd) Nodal reconstruction at k = 0: exact de-averaging nodal heights in [-432.7008646, 3723.991744] m vs WRF z-faces in [-70.00720601, 3612.089067] m max |avg4(nodal) - WRF z-faces| = 9.094947018e-13 m (direct interpolation gives 366.8664961 m) max deviation from direct interpolation = 797.3277235 m Verifying reconstructed nodal surface heights Max |avg4(nodal z) - WRF z-face| on the surface: 9.094947018e-13 m Building an ERF grid from the mean wrfinput layer profile: dz0 17.78478816 m (file's own 17.78478816 m), dz_top 210.0059197 m (file's own 190.0667248 m) the 823.8773063 m by which z_top exceeds the file's mean column depth is absorbed aloft: layers are stretched by 0% at the surface and 10.49062899% at the top layer-to-layer stretch ratio above the surface layer is in [1.004572843, 1.031214525] Verifying grid integrity Min nodal thickness of the first layer: 13.76396705 m Cells whose height lies outside the wrfinput column, and are therefore extrapolated rather than interpolated: below the source column, bottom plane : 29393 of 40000 (73.4825% of that plane) below the source column, interior : 0 above the source column : 61425 total : 90818 of 7040000 (these are extrapolated from the two nearest source levels; theta linearly, every other primitive held) The state read from WRF is being rebalanced! Loading boundary data from NetCDF file Reading Times wrfbdy datetime 0 : 2013-11-21_06:00:00 1385013600 wrfbdy datetime 1 : 2013-11-21_09:00:00 1385024400 wrfbdy datetime 2 : 2013-11-21_12:00:00 1385035200 wrfbdy datetime 3 : 2013-11-21_15:00:00 1385046000 wrfbdy datetime 4 : 2013-11-21_18:00:00 1385056800 wrfbdy datetime 5 : 2013-11-21_21:00:00 1385067600 wrfbdy datetime 6 : 2013-11-22_00:00:00 1385078400 wrfbdy datetime 7 : 2013-11-22_03:00:00 1385089200 start_bdy_time 1385013600 final_bdy_time 1385100000 WRF input and boundary hydrostatic rebalance: enabled Reading time slice 0 from wrfbdy_chisholmview_d01_width1 Reading time slice 1 from wrfbdy_chisholmview_d01_width1 Reading time slice 2 from wrfbdy_chisholmview_d01_width1 start_bdy_time is 1385013600 from wrfbdy but note that time variable in simulation is elapsed time Problem name "wps" does not add any velocity perturbations. Problem name "wps" does not add any velocity perturbations. Problem name "wps" does not add any velocity perturbations. Problem name "wps" does not add any velocity perturbations. Problem name "wps" does not add any velocity perturbations. Problem name "wps" does not add any velocity perturbations. NOTE: terrain-fitted mesh with an LES closure (horizontal eddy viscosity K_h) and an explicit vertical-diffusion stage (erf.vert_implicit_fac has a zero): the terrain-metric term K_h h^2 d2/dz2 stays explicit there and limits the time step on steep slopes. To relax it, use erf.vert_implicit_fac = 1 1 1 with erf.implicit_terrain_metric = true, or erf.smag2d_slope_limiter = true for Smagorinsky2D; erf.diffusive_dt_check reports the limit. Based on cfl of one Compressible dt at level 0 would be: 8.781986051 (terrain aware) 121.0229903 (terrain unaware) Anelastic dt at level 0 would be: 32.13141142 (terrain aware) 346.1086493 (terrain unaware) Fixed dt at level 0 is: 0.5 smallest even ratio is: 4 The smooth-flow limit will be included for variable roughness models over sea: 1 If there is water, Charnock relation with C_a=0.0185 will be used MOST reference height at level 0: 6.942394594 to 8.930187725 Variable sea roughness (type charnock) Using MOST with specified surface temperature (land: TSK, sea: SST) Max value of |p_hse - p_eos| is 7.275957614e-11, less than the tolerance 1.007220393e-08 Min/max value of dp0/dx are -0.003206702234 0.0008144358082 with min at face (0,109,0) with max at face (1,110,1) Min/max value of dp0/dy are -0.001466169527 0.0002670507391 with min at face (101,0,1) with max at face (3,200,1) Min/max value of dp0/dz + rho0*|g| are less than 1e-08 Min/Max value of x-gradient of full (moist) pressure are -0.1346856626 0.3885345443 with min at face (0,118,0) with max at face (0,109,0) Min/Max value of y-gradient of full (moist) pressure are -0.2816734687 0.1403411517 with min at face (101,0,0) with max at face (60,0,0) Min/max value of moist dp/dz + rho_m*|g| are -1.274003125e-11 2.015276834e-11 with min at face (133,12,8) with max at face (170,167,2) Min/max value of -dp'/dz + buoyancy are -2.10416129e-11 1.336009081e-11 with min at face (170,167,2) with max at face (126,28,3) Writing native 3D plotfile WPS_Test-OMP_plt00000 3DPlotfile write time = 0.59531653 seconds. TIME= 0 MASS = 1.18689828e+17 RHO THETA = 3.755908636e+19 RHO SCALAR = 0 RHO QTOTAL = 2.54058725e+14 2013-11-21 06:00:00.000000 (0 s elapsed) Coarse STEP 1 starts ... Based on cfl of one Compressible dt at level 0 would be: 8.781986051 (terrain aware) 121.0229903 (terrain unaware) Anelastic dt at level 0 would be: 32.13141142 (terrain aware) 346.1086493 (terrain unaware) Fixed dt at level 0 is: 0.5 smallest even ratio is: 4 [Level 0 step 1] ADVANCE from elapsed time = 0 to 0.5 with dt = 0.5 Making slow rhs at time 0 for fast variables advancing from 0 to 0.1666666666667 Fast time integration at level 0 from 0 to 0.1666666666667 with dt = 0.1666666666667 Time integration of scalars at level 0 from 0 to 0.1666666666667 with dt = 0.1666666666667 using RHS created at 0 Making slow rhs at time 0.1666666666667 for fast variables advancing from 0 to 0.25 Fast time integration at level 0 from 0 to 0.25 with dt = 0.25 Time integration of scalars at level 0 from 0 to 0.25 with dt = 0.25 using RHS created at 0.1666666666667 Making slow rhs at time 0.25 for fast variables advancing from 0 to 0.5 Fast time integration at level 0 from 0 to 0.25 with dt = 0.25 Fast time integration at level 0 from 0.25 to 0.5 with dt = 0.25 Time integration of scalars at level 0 from 0 to 0.5 with dt = 0.5 using RHS created at 0.25 Done with advance_dycore at level 0 Testing new state for NaNs after advance_microphysics [Level 0 step 1] Advanced 7040000 cells Cloud fraction 0 0.1997 Coarse STEP 1 ends. TIME = 0.5 DT = 0.5 Testing new state and vels for NaNs at end of timestep Timestep time = 13.67556701 seconds. TIME= 0.5 MASS = 1.186898345e+17 RHO THETA = 3.75590962e+19 RHO SCALAR = 0 RHO QTOTAL = 2.540583229e+14 Writing native 3D plotfile WPS_Test-OMP_plt00001 3DPlotfile write time = 0.580116143 seconds. 2013-11-21 06:00:00.500000 (0.5 s elapsed) Coarse STEP 2 starts ... Based on cfl of one Compressible dt at level 0 would be: 8.782352267 (terrain aware) 121.0226648 (terrain unaware) Anelastic dt at level 0 would be: 42.83107213 (terrain aware) 288.3729509 (terrain unaware) Fixed dt at level 0 is: 0.5 smallest even ratio is: 4 [Level 0 step 2] ADVANCE from elapsed time = 0.5 to 1 with dt = 0.5 Making slow rhs at time 0.5 for fast variables advancing from 0.5 to 0.6666666666667 Fast time integration at level 0 from 0.5 to 0.6666666666667 with dt = 0.1666666666667 Time integration of scalars at level 0 from 0.5 to 0.6666666666667 with dt = 0.1666666666667 using RHS created at 0.5 Making slow rhs at time 0.6666666666667 for fast variables advancing from 0.5 to 0.75 Fast time integration at level 0 from 0.5 to 0.75 with dt = 0.25 Time integration of scalars at level 0 from 0.5 to 0.75 with dt = 0.25 using RHS created at 0.6666666666667 Making slow rhs at time 0.75 for fast variables advancing from 0.5 to 1 Fast time integration at level 0 from 0.5 to 0.75 with dt = 0.25 Fast time integration at level 0 from 0.75 to 1 with dt = 0.25 Time integration of scalars at level 0 from 0.5 to 1 with dt = 0.5 using RHS created at 0.75 Done with advance_dycore at level 0 Testing new state for NaNs after advance_microphysics [Level 0 step 2] Advanced 7040000 cells Cloud fraction 0.5 0.195925 Coarse STEP 2 ends. TIME = 1 DT = 0.5 Testing new state and vels for NaNs at end of timestep Timestep time = 12.09037355 seconds. TIME= 1 MASS = 1.18689841e+17 RHO THETA = 3.755909611e+19 RHO SCALAR = 0 RHO QTOTAL = 2.540579231e+14 Writing native 3D plotfile WPS_Test-OMP_plt00002 3DPlotfile write time = 0.563109966 seconds. 2013-11-21 06:00:01.000000 (1 s elapsed) Coarse STEP 3 starts ... Based on cfl of one Compressible dt at level 0 would be: 8.782339408 (terrain aware) 121.0223694 (terrain unaware) Anelastic dt at level 0 would be: 53.69394632 (terrain aware) 292.3379659 (terrain unaware) Fixed dt at level 0 is: 0.5 smallest even ratio is: 4 [Level 0 step 3] ADVANCE from elapsed time = 1 to 1.5 with dt = 0.5 Making slow rhs at time 1 for fast variables advancing from 1 to 1.166666666667 Fast time integration at level 0 from 1 to 1.166666666667 with dt = 0.1666666666667 Time integration of scalars at level 0 from 1 to 1.166666666667 with dt = 0.1666666666667 using RHS created at 1 Making slow rhs at time 1.166666666667 for fast variables advancing from 1 to 1.25 Fast time integration at level 0 from 1 to 1.25 with dt = 0.25 Time integration of scalars at level 0 from 1 to 1.25 with dt = 0.25 using RHS created at 1.166666666667 Making slow rhs at time 1.25 for fast variables advancing from 1 to 1.5 Fast time integration at level 0 from 1 to 1.25 with dt = 0.25 Fast time integration at level 0 from 1.25 to 1.5 with dt = 0.25 Time integration of scalars at level 0 from 1 to 1.5 with dt = 0.5 using RHS created at 1.25 Done with advance_dycore at level 0 Testing new state for NaNs after advance_microphysics [Level 0 step 3] Advanced 7040000 cells Cloud fraction 1 0.193875 Coarse STEP 3 ends. TIME = 1.5 DT = 0.5 Testing new state and vels for NaNs at end of timestep Timestep time = 12.13776954 seconds. TIME= 1.5 MASS = 1.186898474e+17 RHO THETA = 3.755909706e+19 RHO SCALAR = 0 RHO QTOTAL = 2.540575256e+14 Writing native 3D plotfile WPS_Test-OMP_plt00003 3DPlotfile write time = 0.592021908 seconds. Total Time: 52.68809144 Unused ParmParse Variables: [TOP]::amr.checkpoint_files_output(nvals = 1) :: [0] AMReX (26.10-76-g347a3a2d2675) finalized