Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Adding NWM namelists and parameter tables #794

Merged
merged 1 commit into from
Nov 25, 2024
Merged

Adding NWM namelists and parameter tables #794

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
226 changes: 226 additions & 0 deletions src/template/examples/nwm/namelists/v2.0/analysis_assim/hydro.namelist
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,226 @@
&HYDRO_nlist
!!!! ---------------------- SYSTEM COUPLING ----------------------- !!!!

! Specify what is being coupled: 1=HRLDAS (offline Noah-LSM), 2=WRF, 3=NASA/LIS, 4=CLM
sys_cpl = 1

!!!! ------------------- MODEL INPUT DATA FILES ------------------- !!!!

! Specify land surface model gridded input data file (e.g.: 'geo_em.d01.nc')
GEO_STATIC_FLNM = './DOMAIN/geo_em.d01_1km.nc'

! Specify the high-resolution routing terrain input data file (e.g.: 'Fulldom_hires.nc')
GEO_FINEGRID_FLNM = './DOMAIN/Fulldom_hires_netcdf_250m.nc'

! Specify the spatial hydro parameters file (e.g.: 'hydro2dtbl.nc')
! If you specify a filename and the file does not exist, it will be created for you.
HYDROTBL_F = './DOMAIN/HYDRO_TBL_2D.nc'

! Specify spatial metadata file for land surface grid. (e.g.: 'GEOGRID_LDASOUT_Spatial_Metadata.nc')
LAND_SPATIAL_META_FLNM = './DOMAIN/WRF_Hydro_NWM_geospatial_data_template_land_GIS.nc'

! Specify the name of the restart file if starting from restart...comment out with '!' if not...
RESTART_FILE = './RESTART/HYDRO_RST.2018-06-01_00:00_DOMAIN1'

!!!! --------------------- MODEL SETUP OPTIONS -------------------- !!!!

! Specify the domain or nest number identifier...(integer)
IGRID = 1

! Specify the restart file write frequency...(minutes)
! A value of -99999 will output restarts on the first day of the month only.
rst_dt = 60

! Reset the LSM soil states from the high-res routing restart file (1=overwrite, 0=no overwrite)
! NOTE: Only turn this option on if overland or subsurface rotuing is active!
rst_typ = 1

! Restart file format control
rst_bi_in = 0 !0: use netcdf input restart file (default)
!1: use parallel io for reading multiple restart files, 1 per core
rst_bi_out = 0 !0: use netcdf output restart file (default)
!1: use parallel io for outputting multiple restart files, 1 per core

! Restart switch to set restart accumulation variables to 0 (0=no reset, 1=yes reset to 0.0)
RSTRT_SWC = 1

! Specify baseflow/bucket model initialization...(0=cold start from table, 1=restart file)
GW_RESTART = 1

!!!! -------------------- MODEL OUTPUT CONTROL -------------------- !!!!

! Specify the output file write frequency...(minutes)
out_dt = 60

! Specify the number of output times to be contained within each output history file...(integer)
! SET = 1 WHEN RUNNING CHANNEL ROUTING ONLY/CALIBRATION SIMS!!!
! SET = 1 WHEN RUNNING COUPLED TO WRF!!!
SPLIT_OUTPUT_COUNT = 1

! Specify the minimum stream order to output to netcdf point file...(integer)
! Note: lower value of stream order produces more output.
order_to_write = 1

! Flag to turn on/off new I/O routines: 0 = deprecated output routines (use when running with Noah LSM),
! 1 = with scale/offset/compression, ! 2 = with scale/offset/NO compression,
! 3 = compression only, 4 = no scale/offset/compression (default)
io_form_outputs = 2

! Realtime run configuration option:
! 0=all (default), 1=analysis, 2=short-range, 3=medium-range, 4=long-range, 5=retrospective,
! 6=diagnostic (includes all of 1-4 outputs combined)
io_config_outputs = 1

! Option to write output files at time 0 (restart cold start time): 0=no, 1=yes (default)
t0OutputFlag = 1

! Options to output channel & bucket influxes. Only active for UDMP_OPT=1.
! Nonzero choice requires that out_dt above matches NOAH_TIMESTEP in namelist.hrldas.
! 0=None (default), 1=channel influxes (qSfcLatRunoff, qBucket)
! 2=channel+bucket fluxes (qSfcLatRunoff, qBucket, qBtmVertRunoff_toBucket)
! 3=channel accumulations (accSfcLatRunoff, accBucket) *** NOT TESTED ***
output_channelBucket_influx = 2

! Output netcdf file control
CHRTOUT_DOMAIN = 1 ! Netcdf point timeseries output at all channel points (1d)
! 0 = no output, 1 = output
CHANOBS_DOMAIN = 0 ! Netcdf point timeseries at forecast points or gage points (defined in Routelink)
! 0 = no output, 1 = output at forecast points or gage points.
CHRTOUT_GRID = 0 ! Netcdf grid of channel streamflow values (2d)
! 0 = no output, 1 = output
! NOTE: Not available with reach-based routing
LSMOUT_DOMAIN = 0 ! Netcdf grid of variables passed between LSM and routing components (2d)
! 0 = no output, 1 = output
! NOTE: No scale_factor/add_offset available
RTOUT_DOMAIN = 1 ! Netcdf grid of terrain routing variables on routing grid (2d)
! 0 = no output, 1 = output
output_gw = 0 ! Netcdf GW output
! 0 = no output, 1 = output
outlake = 1 ! Netcdf grid of lake values (1d)
! 0 = no output, 1 = output
frxst_pts_out = 0 ! ASCII text file of forecast points or gage points (defined in Routelink)
! 0 = no output, 1 = output

!!!! ------------ PHYSICS OPTIONS AND RELATED SETTINGS ------------ !!!!

! Specify the number of soil layers (integer) and the depth of the bottom of each layer... (meters)
! Notes: In Version 1 of WRF-Hydro these must be the same as in the namelist.input file.
! Future versions will permit this to be different.
NSOIL=4
ZSOIL8(1) = -0.10
ZSOIL8(2) = -0.40
ZSOIL8(3) = -1.00
ZSOIL8(4) = -2.00

! Specify the grid spacing of the terrain routing grid...(meters)
DXRT = 250.0

! Specify the integer multiple between the land model grid and the terrain routing grid...(integer)
AGGFACTRT = 4

! Specify the channel routing model timestep...(seconds)
DTRT_CH = 300

! Specify the terrain routing model timestep...(seconds)
DTRT_TER = 10

! Switch to activate subsurface routing...(0=no, 1=yes)
SUBRTSWCRT = 1

! Switch to activate surface overland flow routing...(0=no, 1=yes)
OVRTSWCRT = 1

! Specify overland flow routing option: 1=Seepest Descent (D8) 2=CASC2D (not active)
! NOTE: Currently subsurface flow is only steepest descent
rt_option = 1

! Switch to activate channel routing...(0=no, 1=yes)
CHANRTSWCRT = 1

! Specify channel routing option: 1=Muskingam-reach, 2=Musk.-Cunge-reach, 3=Diff.Wave-gridded
channel_option = 2

! Specify the reach file for reach-based routing options (e.g.: 'Route_Link.nc')
route_link_f = './DOMAIN/RouteLink_NHDPLUS.nc'

! If using channel_option=2, activate the compound channel formulation? (Default=.FALSE.)
compound_channel = .TRUE.

! Specify the lake parameter file (e.g.: 'LAKEPARM.nc').
! Note REQUIRED if lakes are on.
route_lake_f = './DOMAIN/LAKEPARM_NHDPLUS.nc'

! Switch to activate baseflow bucket model...(0=none, 1=exp. bucket, 2=pass-through)
GWBASESWCRT = 1

! Groundwater/baseflow 2d mask specified on land surface model grid (e.g.: 'GWBASINS.nc')
!Note: Only required if baseflow model is active (1 or 2) and UDMP_OPT=0.
!gwbasmskfil = './DOMAIN/GWBASINS.nc'

! Groundwater bucket parameter file (e.g.: 'GWBUCKPARM.nc')
GWBUCKPARM_file = './DOMAIN/GWBUCKPARM_CONUS.nc'

! User defined mapping, such NHDPlus: 0=no (default), 1=yes
UDMP_OPT = 1

! If on, specify the user-defined mapping file (e.g.: 'spatialweights.nc')
udmap_file = './DOMAIN/spatialweights_250m_all_basins.nc'

/

&NUDGING_nlist

! Path to the 'timeslice' observation files.
timeSlicePath = './nudgingTimeSliceObs/'

nudgingParamFile = './DOMAIN/nudgingParams.nc'

! Nudging restart file = 'nudgingLastObsFile'
! nudgingLastObsFile defaults to '', which will look for nudgingLastObs.YYYY-mm-dd_HH:MM:SS.nc
! **AT THE INITALIZATION TIME OF THE RUN**. Set to a missing file to use no restart.
nudgingLastObsFile = './RESTART/nudgingLastObs.2018-06-01_00:00:00.nc'

!! Parallel input of nudging timeslice observation files?
readTimesliceParallel = .TRUE.

! temporalPersistence defaults to true, only runs if necessary params present.
temporalPersistence = .TRUE.

! The total number of last (obs, modeled) pairs to save in nudgingLastObs for
! removal of bias. This is the maximum array length. (This option is active when persistBias=FALSE)
! (Default=960=10days @15min obs resolution, if all the obs are present and longer if not.)
nLastObs = 480

! If using temporalPersistence the last observation persists by default.
! This option instead persists the bias after the last observation.
persistBias = .TRUE.

! AnA (FALSE) vs Forecast (TRUE) bias persistence.
! If persistBias: Does the window for calculating the bias end at
! model init time (=t0)?
! FALSE = window ends at model time (moving),
! TRUE = window ends at init=t0(fcst) time.
! (If commented out, Default=FALSE)
! Note: Perfect restart tests require this option to be .FALSE.
biasWindowBeforeT0 = .FALSE.

! If persistBias: Only use this many last (obs, modeled) pairs. (If Commented out, Default=-1*nLastObs)
! > 0: apply an age-based filter, units=hours.
! = 0: apply no additional filter, use all available/usable obs.
! < 0: apply an count-based filter, units=count
maxAgePairsBiasPersist = 3

! If persistBias: The minimum number of last (obs, modeled) pairs, with age less than
! maxAgePairsBiasPersist, required to apply a bias correction. (default=8)
minNumPairsBiasPersist = 1

! If persistBias: give more weight to observations closer in time? (default=FALSE)
invDistTimeWeightBias = .TRUE.

! If persistBias: 'No constructive interference in bias correction?', Reduce the bias adjustment
! when the model and the bias adjustment have the same sign relative to the modeled flow at t0?
! (default=FALSE)
! Note: Perfect restart tests require this option to be .FALSE.
noConstInterfBias = .TRUE.

/
72 changes: 72 additions & 0 deletions src/template/examples/nwm/namelists/v2.0/analysis_assim/namelist.hrldas
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,72 @@
&NOAHLSM_OFFLINE

HRLDAS_SETUP_FILE = './DOMAIN/wrfinput_d01_1km.nc'
INDIR = './forcing'
SPATIAL_FILENAME = './DOMAIN/soil_veg_properties_ASM.nc'
OUTDIR = './'

START_YEAR = 2018
START_MONTH = 06
START_DAY = 01
START_HOUR = 00
START_MIN = 00

RESTART_FILENAME_REQUESTED = './RESTART/RESTART.2018060100_DOMAIN1'

! Specification of simulation length in days OR hours
!KDAY = 1
KHOUR = 3

! Physics options (see the documentation for details)
DYNAMIC_VEG_OPTION = 4
CANOPY_STOMATAL_RESISTANCE_OPTION = 1
BTR_OPTION = 1
RUNOFF_OPTION = 3
SURFACE_DRAG_OPTION = 1
FROZEN_SOIL_OPTION = 1
SUPERCOOLED_WATER_OPTION = 1
RADIATIVE_TRANSFER_OPTION = 3
SNOW_ALBEDO_OPTION = 1
PCP_PARTITION_OPTION = 1
TBOT_OPTION = 2
TEMP_TIME_SCHEME_OPTION = 3
GLACIER_OPTION = 2
SURFACE_RESISTANCE_OPTION = 4

! Timesteps in units of seconds
FORCING_TIMESTEP = 3600
NOAH_TIMESTEP = 3600
OUTPUT_TIMESTEP = 3600

! Land surface model restart file write frequency
RESTART_FREQUENCY_HOURS = 1

! Split output after split_output_count output times.
SPLIT_OUTPUT_COUNT = 1

! Soil layer specification
NSOIL=4
soil_thick_input(1) = 0.10
soil_thick_input(2) = 0.30
soil_thick_input(3) = 0.60
soil_thick_input(4) = 1.00

! Forcing data measurement height for winds, temp, humidity
ZLVL = 10.0

! Restart file format options
rst_bi_in = 0 !0: use netcdf input restart file
!1: use parallel io for reading multiple restart files (1 per core)
rst_bi_out = 0 !0: use netcdf output restart file
!1: use parallel io for outputting multiple restart files (1 per core)
/

&WRF_HYDRO_OFFLINE

! Specification of forcing data: 1=HRLDAS-hr format, 2=HRLDAS-min format, 3=WRF,
! 4=Idealized, 5=Ideal w/ spec. precip,
! 6=HRLDAS-hr format w/ spec. precip,
! 7=WRF w/ spec. precip
FORC_TYP = 2

/
Loading