WRF v4.2.1 : WARNING : Older v3 input data detected

Any issues with the actual running of the WRF.

WRF v4.2.1 : WARNING : Older v3 input data detected

Postby goddarpb » Wed Mar 10, 2021 1:22 pm

Hello,

I ran across an error while running the wrf.exe from an ndown simulation; it is similar to the issue posted here (and within the links provided by Enzo): viewtopic.php?f=7&t=11255

The error occurs at the beginning of the run:
Max map factor in domain 1 = 1.04. Scale the dt in the model accordingly.
D01: Time step = 30.00000 (s)
D01: Grid Distance = 9.000000 (km)
D01: Grid Distance Ratio dt/dx = 3.333333 (s/km)
D01: Ratio Including Maximum Map Factor = 3.466527 (s/km)
D01: NML defined reasonable_time_step_ratio = 6.000000
---- WARNING : Older v3 input data detected
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 684
---- Error : Cannot use moist theta option with old data
-------------------------------------------

Unlike the postings from Enzo, I am using WRF v4.2.1 and I am not attempting a vertical refinement. In both the wrfinput_d01 and wrfbdy_d01, the files state:
TITLE = " OUTPUT FROM NDOWN_EM V4.2.1 PREPROCESSOR" ;

The only "v3" data I use is from the geo-static data used in the WPS step. I don't think this is the error, in the past I have successfully ran Ndown simulations with the geo-static data v3.

Below is my namelist.input.

Thanks for any troubleshooting ideas,
Paul


------------
&time_control ! time
start_year = 2014, 2014, ! start year
start_month = 08, 08, ! start month
start_day = 27, 27, ! start day
start_hour = 00, 00, ! start hour
start_minute = 00, 00, ! start minute
start_second = 00, 00, ! start second
end_year = 2014, 2014, ! end year
end_month = 09, 09, ! end month
end_day = 30, 30, ! end day
end_hour = 00, 00, ! end hour
end_minute = 00, 00, ! end minute
end_second = 00, 00, ! end second
interval_seconds = 10800, ! interval between wrfout files in seconds
input_from_file = .true., .true., ! whether the nested run will have input files for domains other than domain 1
history_interval = 60, 60, ! history output file interval in minutes (integer only)
frames_per_outfile = 24, 24, ! number of output times bulked into each history file
restart = .false., ! whether this run is a restart
restart_interval = 7200, ! restart output file interval in minutes
io_form_history = 2 ! NetCDF
io_form_restart = 2 ! NetCDF
io_form_input = 2 ! NetCDF
io_form_boundary = 2 ! NetCDF
debug_level = 0 ! debugging level
output_diagnostics = 1 ! Output Max, Min, Mean, SD of common Meteo Variables
auxhist3_outname = "wrfxtrm_d<domain>_<date>" ! Output name of extreme diagnostics
auxhist3_interval = 1440, 1440, ! Output interval in minutes, set to daily output
frames_per_auxhist3 = 1, 1, ! One per file
! io_form_auxinput2 = 2, ! Set for Ndown
io_form_auxhist3 = 2, ! outfile as NetCDF
auxinput4_inname = "wrflowinp_d<domain>", ! sst update name (create from real.exe)
auxinput4_interval = 180, ! sst update in minutes
io_form_auxinput4 = 2, ! NetCDF
! iofields_filename = "my_delete_outfields.txt", "myoutfields.txt", "myoutfields.txt", ! an option to request particular variables to appear in output
! io_form_auxinput24 = 2, ! outfields as netCDF
! auxhist24_interval = 60, 60, ! outfields every ## minutes
io_form_auxinput5 = 2, ! anthropogenic as netCDF
! io_form_auxinput6 = 2, ! biogenic as netCDF
! io_form_auxinput7 = 2, ! biomass burning as netCDF
auxinput5_inname = 'wrfchemi_<hr>z_d<domain>', ! anthropogentic emission filename
! auxinput5_interval_m = 60, 60, ! anthropogenic emission time interval, minutes
auxinput5_interval_d = 30, 30, ! anthropogenic emission time interval (EDGAR-HTAP1 is one month), days
! auxinput6_inname = 'wrfbiochemi_d<domain>', ! biogenic emission filename
! auxinput6_interval_d = 600, 600, ! biogenic emission time interval, days - set high because static
! auxinput7_inname = 'wrffirechemi_d<domain>_<date>', ! biomass burning filename
! auxinput7_interval_m = 60, 60, ! biomass burning time interval, minutes - fire is onehour intervals
! frames_per_auxinput7 = 1, 1, ! biomass burning files per time interval
! io_form_auxinput12 = __ISRESTARTVALUE__, ! restart file type
force_use_old_data = T, ! allow WRFChem4 to run with data from WRFChem3
/

&domains ! domains - dimensions, nesting, parameters
time_step = 30, ! timestep for coarsest domain meteorology, in seconds (90s still hand some banding signatures)
time_step_fract_num = 0, ! numerator for fractional time step
time_step_fract_den = 1, ! denominator for fractional time step
max_dom = 1, ! number of domains
nproc_x = 16, ! good for <10 per and 16x9 is 144 which is 3 cores with 48 tasks per
nproc_y = 9, ! good for <10 per and 16x9 is 144 which is 3 cores with 48 tasks per
e_we = 163, 283, ! westeast dimension
e_sn = 118, 187, ! southnorth dimension
e_vert = 57, 57, ! vertical dimension, then setting eta values with more near the surface
eta_levels =
1.000, 0.997, 0.992, 0.985, 0.978, 0.969, 0.960, 0.950,
0.938, 0.925, 0.910, 0.894, 0.876, 0.857, 0.835, 0.812,
0.787, 0.760, 0.731, 0.700, 0.668, 0.635, 0.600, 0.565,
0.530, 0.494, 0.458, 0.423, 0.388, 0.355, 0.323, 0.293,
0.264, 0.237, 0.212, 0.188, 0.167, 0.147, 0.130, 0.114,
0.099, 0.086, 0.074, 0.064, 0.054, 0.046, 0.039, 0.032,
0.027, 0.022, 0.017, 0.013, 0.010, 0.007, 0.004, 0.002,
0.000,

num_metgrid_levels = 38, ! number of vertical levels in WPS output
num_metgrid_soil_levels = 4, ! number of soil levels or layers in WPS output
dx = 9000, 3000, ! westeast resolution, metres
dy = 9000, 3000, ! southnorth resolution, metres
grid_id = 2, 3, ! grid ID
parent_id = 1, 2, ! parent ID
i_parent_start = 20, 40, ! x coordinate of the lower-left corner
j_parent_start = 20, 40, ! y coordinate of the lower-left corner
parent_grid_ratio = 3, 3, ! nesting ratio relative to the domains parent
! parent_time_step_ratio = 3, 3, ! parent-to-nest time step ratio
feedback = 0, ! feedback from nest to its parent domain
/

&afwa
afwa_diag_opt = 1 ! additional meteo variables in package I modified
/

&diags
solar_diagnostics = 1 ! additional variables including LWP and effective radii
/

&physics ! physics
mp_physics = 40, 40, ! microphysics scheme, 18 = Noaa Severe Storms Lab with CCN vars (works with Mozart), 10 = Morrison 2-moment (no ccn vars), 40 = Morrison 2-moment with aerosol variables (does not work with mozART)
progn = 1, 1, ! prognostic number density, switch to use mix-activate scheme, 0 = off
ra_lw_physics = 4, 4, ! longwave radiation scheme, 4 = RRTMG
ra_sw_physics = 4, 4, ! shortwave radiation scheme, 4 = RRTMG
radt = 5, 5, ! minutes between radiation physics calls, recommended 1 minute per km of dx (e.g. 10 for 10 km grid); use the same value for all nests
sf_sfclay_physics = 5, 5, ! surface layer physics option, 5 = MYNN (Ravan's suggestion), 2 = myjsfcscheme
sf_surface_physics = 2, 2, ! land surface physics option, 2 = NOAH
bl_pbl_physics = 5, 5, ! boundary layer physics option, 5 = MYNN 2.5, 1 = YSU scheme
bldt = 0, 0, ! minutes between boundary-layer physics calls (0=call every timestep)
cu_physics = 5, 5, ! cumulus parameterization option, 5 = Grell 3D, turn off for 4 km domain?, 11 = MSKF
cudt = 0, 0, ! minutes between cumulus physics calls; should be set to 0 when using all cu_physics except Kain-Fritsch, set equal to 10 to match radiation calls
cugd_avedx = 1, ! number of grid boxes over which subsidence is spread, set to 3 for 4km run, 1 for 36km
isfflx = 1, ! heat and moisture fluxes from the surface for real-data cases and when a PBL is used
ifsnow = 0, ! snow-cover effects
icloud = 1, ! cloud effect to the optical depth in radiation
icloud_bl = 1, ! option to couple subgrid-scale clouds from PBL scheme (MYNN only) to radiative scheme
bl_mynn_cloudmix = 1, ! mixing cloud water and ice (qnc and qni are mixed when scalar pblmix = 1)
bl_mynn_cloudpdf = 2, ! Chaboureau and Bechtold (2 = default)
bl_mynn_edmf = 1, ! activate mass-flux in mynn
bl_mynn_mixlength = 2, ! from blending (default = 2?)
scalar_pblmix = 1, ! mix scalar fields consistent with PBL option (exch_h)
tracer_pblmix = 0, ! mix tracer fields consistent with PBL option (exch_h)
grav_settling = 2, ! gravitational settling of fog/cloud droplets, 2 = Fogdes at surface, and Dyunkerke in the atmosphere
surface_input_source = 1, ! where landuse and soil category data come from, 1 = WPS
num_soil_layers = 4, ! number of soil levels or layers in WPS output
sf_urban_physics = 0, 0, ! activate urban canopy model, 1 = single layer, some break the inner domain and no idea why, set to zero for online nests
mp_zero_out = 0, ! this keeps moisture variables above a threshold value, set to zero bc a better method is to set the ??_adv_opt equal to 4 in thr &dynamics section
mp_zero_out_thresh = 1.e-8, ! critical value for moisture variable threshold, below which moisture arrays (except for Qv) are set to zero, will be ignored bc of above 0
cu_rad_feedback = .true., .false., ! sub-grid cloud effect to the optical depth in radiation
cu_diag = 1, 1, ! Additional time-averaged diagnostics from cu_physics
aercu_opt = 0, ! Option to control aerosol interaction in MSKF and Morrison microphysics (mp_physics = 40), = 2 aersol interaction with both MSKF and Morrison
slope_rad = 0, 0, ! use slope-dependent radiation
topo_shading = 0, 0, ! applies neighboring-point shadow effects
num_land_cat = 21, ! number of land categories in input data
sst_update = 1, ! Do an SST update
clean_atm_diag = 1, ! Output SW and LW in clean atmosphere (ie., no aerosols)
/

&fdda ! FDDA - options for grid, obs and spectral nudging -- Many set in WRFotron example namelist.wrf.blueprint_202
/

&dynamics ! dynamics - diffusion, damping options, advection options
rk_ord = 3, ! time-integration scheme option, 3 = Runge-Kutta 3rd order
w_damping = 1, ! vertical velocity damping flag, 1 = with damping
diff_opt = 1, 1, ! turbulence and mixing option, 1 = evaluates 2nd order diffusion term on coordinate surfaces
km_opt = 4, 4, ! eddy coefficient option, 4 = horizontal Smagorinsky first order closure
diff_6th_opt = 0, 0, ! 6th-order numerical diffusion, 0 = none
diff_6th_factor = 0.12, 0.12 ! 6th-order numerical diffusion nondimensional rate
base_temp = 290. ! base state temperature (K)
damp_opt = 3, ! upper-level damping flag, 3 = Rayleigh damping
zdamp = 5000., 5000., ! damping depth (m) from model top
dampcoef = 0.2, 0.2, ! damping coefficient
khdif = 0, 0, ! horizontal diffusion constant (m2/s)
kvdif = 0, 0, ! vertical diffusion constant (m2/s)
non_hydrostatic = .true., .false., ! running the model in nonhydrostatic mode
moist_adv_opt = 4, 4, ! advection options for moisture, 4 = WENO 5th order w/ positive definite, 2 = monotonic
chem_adv_opt = 4, 4, ! advection options for chemistry, 4 = WENO 5th order w/ positive definite, 2 = monotonic
scalar_adv_opt = 4, 4, ! advection options for scalars, 4 = WENO 5th order w/ positive definite, 2 = monotonic
tke_adv_opt = 4, 4, ! advection options for TKE, 4 = WENO 5th order w/ positive definite, 2 = monotonic
tracer_adv_opt = 4, 4, ! advection options for tracers, 4 = WENO 5th order w/ positive definite
momentum_adv_opt = 4, 4, ! advection options for momentum, 4 = WENO 5th order w/ positive definite
do_avgflx_em = 0, 0, ! outputs time-averaged masscoupled advective velocities, 1 = on
/

&bdy_control ! Boundary condition control
spec_bdy_width = 5, ! total number of rows for specified boundary value nudging
spec_zone = 1, ! number of points in specified zone
relax_zone = 4, ! number of points in relaxation zone
specified = .true., .true., ! specified boundary condition
nested = .false., .false., ! nested boundary conditions
/

&grib2
/

&namelist_quilt ! options for asynchronized I/O for MPI applications
nio_tasks_per_group = 0, ! # of processors used for IO quilting per IO group
nio_groups = 1, ! number of quilting groups
/

&chem ! chemistry
kemit = 1, ! number of vertical levels in the emissions input data file
chem_opt = 202, 202, ! chemistry option, 201 = MOZART-MOSAIC (4 bins + simplified SOA + no aqeuous chemistry), 202 = MOZART-MOSAIC (4 bins + VBS SOA + aqeuous chemistry).
! bioemdt = 2.5, 2.5, ! timestep, biogenic, minutes !NEED TO CHANGE D2 and D3
photdt = 10., 10., ! timestep, photolysis, minutes
chemdt = 0.5, 0.166, ! timestep, chemistry, minutes !NEED TO CHANGE D2 and D3
io_style_emissions = 1, ! anthropogenic emissions, files, two 12-h emissions data files used
emiss_inpt_opt = 102, 102, ! =102 RADM2 emission speciation adapted after reading data file to follow the RADM2/SORGAM framework (including isoprene)
emiss_opt = 10, 10, ! anthropogenic emissions, setting, 10 = MOZART (MOZART + aerosols) emissions
chem_in_opt = 0, 0, ! initialize chemistry, we use mozbc
phot_opt = 4, 4, ! photolysis option, 1 = Full TUV, 3 = Madronich F-TUV, 4 = New full TUV scheme
gas_drydep_opt = 1, 1, ! dry deposition of gas species, 1 = on
aer_drydep_opt = 1, 1, ! dry deposition of aerosols, 1 = on
! bio_emiss_opt = 3, 3, ! includes MEGAN biogenic emissions online based upon the weather, land use data
gas_bc_opt = 1, 1, ! gas boundary conditions, 1 = default
gas_ic_opt = 1, 1, ! gas initial conditions, 1 = default
aer_bc_opt = 1, 1, ! aerosol boundary conditions, 1 = default
aer_ic_opt = 1, 1, ! aerosol initial conditions, 1 = default
gaschem_onoff = 1, 1, ! gas phase chemistry, 1 = on
aerchem_onoff = 1, 1, ! aerosol chemistry, 1 = on
wetscav_onoff = 1, 1, ! wet scavenging in stratocumulus clouds, 1 = on
cldchem_onoff = 1, 1, ! aqueous chemistry in stratocumulus clouds, 1 = on
vertmix_onoff = 1, 1, ! vertical turbulent mixing, 1 = on
chem_conv_tr = 1, 0, ! subgrid convective transport, 1 = on
conv_tr_wetscav = 1, 0, ! wet scavenging in cumulus clouds, subgrid, 1 = on
conv_tr_aqchem = 1, 0, ! aqueous chemistry in cumulus clouds, subgrid, 1 = on
seas_opt = 2, ! sea salt emissions, 2 = MOSAIC or MADE/SORGAM sea salt emissions
dust_opt = 3, ! dust emissions, 1 = GOCART dust emissions; 3 = GOCART dust emissions with AFWA modifications
dmsemis_opt = 1, ! include GOCART dms emissions from sea surface
! biomass_burn_opt = 2, 2, ! biomass burning emissions, 2 = MOCART
! plumerisefire_frq = 30, 30, ! time interval for calling the biomass burning plume rise subroutine
! scale_fire_emiss = .true., .true., ! must be equal to .true. when running with FINN emissions
aer_ra_feedback = 1, 1, ! feedback from the aerosols to the radiation schemes, 1 = on
ne_area = 600, ! total number of chemical species, in the chemical name list, best to set to a value larger than all chemical species, for MEGAN -- not needed if not using?
opt_pars_out = 1, ! include optical properties in output
have_bcs_chem = .true., .true., ! gets lateral boundary data from wrfbdy (.true.) or idealized profile (.false.)
have_bcs_upper = .false., .false., ! upper boundary boundary condition for chemical species
aer_op_opt = 4, 4, ! aerosol optical properties, 1 = volume, 2 = approximate Maxwell-Garnet, 3 = complex volume-mixing, 4 = complex Maxwell-Garnet (complex/exact is shown to be better than approx), 5 = complex core-shell
! bbinjectscheme = 2, 2, ! 0 = plumerise (biomass_burn_opt), 1 = all ground level, 2 = flaming evenly in BL (recommended), 3 = flaming top BL, 4 = flaming injected at specific height
/
goddarpb
 
Posts: 4
Joined: Mon Feb 22, 2021 3:43 pm

Return to Runtime Problems

Who is online

Users browsing this forum: Google [Bot] and 9 guests