Recent Releases of WRF-ARW

WRF-ARW - WRF Version 4.6.1 (Bug-fix Release)

The WRF model has been update to Version v4.6.1 on October 16, 2024.

Acknowledgements: We would like to thank Jeronimo Bande (IDING SAS), Jesus Fernandez, Andrés Simón-Moral and Josipa Milovac (Instituto de Física de Cantabria, CSIC-Universidad de Cantabria, Spain), Parag Joshi (Brookhaven National Lab), Katia Lamer (Brookhaven National Lab), Tanya Spero and Jeff Willison (U.S. EPA), and Ted Mansell (NOAA/NSSL) for their contributions to this release.

This is a bug fix release. Following are fixes associated with this release.

Physics

Correct the shortwave radiation balance at the wall, particularly the reflected direct and diffuse radiation reaching the wall which leads to underestimation of SW radiation at the wall. (#2101) Details
Fix non-initialized 3D moisture field in pressure- and height-level diagnostic interpolation. (#2104) Details
Fixes for the calculations for PSIM in the mos subroutine and the temperature used in surface energy balance for the green roof. The later follows the OSU1DPBL document located at: [https://ftp.emc.ncep.noaa.gov/mmb/gcp/ldas/nceplsm/OSU1DPBL/OSU1DPBL-userguide.pdf] (#2038) Details
Lech's stability functions are corrected following the original formulation in the single layer urban canopy model. (#2032) Details
Fix for rare undefined zol in sfclayrev. (#2105) Details
Bug fixes for parameters associated with the heat transfer coefficient for the upward heat flow over PV panels. (#2049) Details
Fixed the bug of missing radar reflectivity in WSM6 scheme when do_radar_ref is set to 1. This bug was introduced when WSM6 scheme was moved to physics_mmm/ directory, and shared with other MMM models in v4.6. (#2098) Details

Data Assimilation

Bugfix for duplicated allocation and file closing for radar_rhv_opt = 2 (#2078) Details

Chemistry

Expanded global attributes to the full suite for auxiliary history files. Also added two global attributes when the WRF-CMAQ model is invoked (WRF_CMAQ_OPTION and DIRECT_SW_FEEDBACK). (#2044) Details

Software

Compilation Tests for ARW (#2115) Details
CMake compilation flags as categorical properties (#2088) Details
CMake improve configure_reader.py script (#2087) Details
Fix fseek test (#2055) Details
Restructure netCDF find modules to use modern import target (#2054) Details
CMake Prepare for alternate core selections (#2086) Details
Change internal flag organization in CMake build to not be global (#2053) Details
Breakout deallocation calls into simpler smaller files (#2070) Details
CMake build no longer uses generator expressions in defines (#2056) Details
Break out module_dm external subroutines into separate files (#2069) Details
Restructure alloc space calls into simpler smaller files (#2068) Details
Introduce a modularized testing framework that allows testing locally and natively on HPC systems that lives within the WRF repository (#2095) Details
Add explicit casts to mismatched pointer types in rsl_lite (#2062) Details
Replace nonstandard isnan() call with value input checks instead (#2061) Details
Fix line continuation syntax causing problems with certain compilers (#2060) Details
Remove unsupported -auto C flag from Intel oneAPI stanza (#2059) Details
Fixed a compile problem with .F90 source files on case-insensitive file systems. (#2057) Details

Miscellaneous

The documentation of a few namelists is updated in this PR. (#2074) Details
Corrected end date for em_b_wave test case (#2107) Details

NOTE!!!! : If downloading the files from the section below, please choose either the v4.6.1.tar.gz file, or the v4.6.1.zip file. DO NOT choose those named "Source Code." They do not include the mandatory NoahMP submodule - needed for compiling WRF.

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by islas 6 months ago

WRF-ARW - WRF Version 4.6.0

This is the Version 4.6 release of the WRF Model. This release includes all components of the system: WRF-ARW, WRF-Chem, and WRFDA.

Acknowledgements: We would like to thank Kasra Shamsaei (University of Nevada, Reno), Tim Juliano, Domingo Munoz-Esparza, Branko Kosovic (NCAR/RAL), Jun Tang (Amazon), RenChuanhua and Huangxin (Nanjing University), David Robertson (Rutgers University), Kasra Shamsaei (University of Nevada, Reno), Tim Juliano, Branko Kosovic (NCAR/RAL), Cristina Archer (UDEL), Yulong Ma(UDEL;GWA-MWF), Ahmad Vasel-Be-Hagh(Tennessee Technological University), Zhixiong Chen (Fujian Normal University, China),JJ Guerrette (NCAR/MMM, now at tomorrow.io), Deqin Li (Liaoning Meteorological Bureau of CMA), Cenlin He (NCAR), Joseph Olson (NOAA), Laura Fowler (NCAR), Sergey Osipov (KAUST), Daniel Wesloh (Penn State), Do Ngoc Khanh (Tokyo Institute of Technology), Jordan Schnell (NOAA), Soren Rasmussen and Ryan Cabell (NCAR),Parag Joshi and Katia Lamer (Brookhaven National Lab), and Ted Mansell (NOAA) for their contributions to this release.

New features in WRFV4.6

Physics

Add revised MM5 surface layer scheme from the future shared MMM-physics repository. Detail
Add new Tiedtke cumulus scheme from the future shared MMM-physics repository. Detail
Add WSM6 microphysics scheme from the future shared MMM-physics repository. detail
Add YSU PBL scheme from the future shared MMM-physics repository. Detail
Add GWDO (gwd_opt=1) code from the future shared MMM-physics repository. Detail
Add two options (slucm_distributed_drag and distributed_ahe_opt) to WRF SLUCM (sf_urban_physics = 1) so that spatially varying urban morphological parameters (building height, plan area index, frontal area index, roughness length for momentum, and displacement height) can be considered. Detail
Add salinity effect in MM5 and revised MM5 surface layer schemes. The effect is lower the saturation vapor pressure over ocean by about 2%. Detail
Add options for three new wind farm parameterizations and their ensembles to account for the individual and overlapping sub-grid wakes of wind turbines within a wind farm. It is recommended to use the new wind farm parameterizations, particularly for coarse resolution, high turbine density, and wind directions aligned with the turbine columns. Detail
Add 3-moment option to NSSL microphysics scheme (for rain, graupel, and hail) (mp_physics=18). Changes to snow aggregation parameters reduce stratiform reflectivity bias. Also add namelist options to enable/disable bulk CCN prediction and graupel/hail density prediction, and deprecate old mp_physics options (17, 19, 21, 22). Main default parameter changes are for graupel/hail fall speed options (icdx, icdxhl; changed from 3 to 6), and default maximum gr/hail droplet collection efficiencies (ehw0/ehlw0 changed from 0.5/0.75 to 0.9/0.9) See doc/README.NSSLmp for details. A problem is fixed for high resolution vertical domains with nz > 128, where hydrometeor sedimentation did not work at k > 128 (this bug affected most or all previous releases). Detail

Real-Data initialization

Add a new namelist option "use_sh_qv = T" to allow for choosing to use SH or Qv data directly. If water mixing ratio or specific humidity is available in the metgrid output, they can be used directly in vertical interpolation. Previously these fields are first converted to RH, interpolated vertically, and the water vapor mixing ratio is computed at model levels.Detail

Data assimilation

Add the capability for assimilating GOES-ABI radiance data. Detail
Add a lightning data assimilation scheme to assimilate pseudo vertical velocity, divergence fields, or water vapor retrievals from lightning flash rate data. Detail
Add a background-dependent hydrometeor retrieval scheme for indirect radar reflectivity assimilation. Detail

WRF-FIRE

Add Scott and Burgan (2005) 40 fuel models for WRF-Fire. Detail
Add a new Truncated Gaussian scheme to distribute fire heat and smoke in the atmosphere. Detail.

Chemistry

Add WRF-SoilN-chem, a dynamic ammonia emission model capable of calculating NH3 emission rate interactively with time- and spatial-varying meteorological and soil conditions. Detail

Compiling

Introduce a CMake build system for em_real and em_ideal as an optional compilation method. Detail. Users are welcome to try and provide feedback.
Add CMake build option to WRF-Hydro. General bug fixes and code cleanup are also done. MPI routine is more efficient. All hydro file suffixes are now .F90. Detail

Bug fixes and enhancement in WRFV4.6

Physics

Correct net long wave fluxes for application in modeling urban climate using Single Layer Urban Canopy Model (SLUCM). It slightly improved 2-m temperature in urban area. Details
Change the behavior of the mp_zero_out flag so that it affects only the 'moist' array, whereas previously it also caused the scalar/chem/tracer arrays to also be set to zero for values below threshold. Now there is a separate flag (mp_zero_out_all) if one wishes to reproduce the old behavior. Details
Move the call to shalwater_init from module_sf_sfclayrev.F to module_physics_init.F. This is in preparation to make sfclayrev module ccpp-compliant, and potentially make the initialization and physics available to other surface layer modules. Detail
Fix a bug in Noah-MP for uninitialized leaf mass when LCZ is used. Detail
Update urban LCZ parameter table (URBPARM_LCZ.TBL) with more reasonable values. Detail
Update MYNN-EDMF . Detail
Fix bugs in Noah-MP for snow combination, vegetation fraction scaling, and urban ground heat flux sign. Detail
Remove zero negative latent heat flux limit (atmosphere to water) in the revised MM5 and original MM5 surface layer schemes. Detail
Make the Pleim-Xiu LSM compatible with MODIS LCZ. Also, effects of evaporation from transpiration, soil in both vegetated and non-veg parts and wet leaves on ground temperature all are considered. Previously only the effect from transpiration and evaporation from non-veg soil are considered. Detail
Modify the RUC LSM to add options to the snow model, use real-time vegetation fraction as an indicator for irrigation scheme, and change soil diffusional and hydraulic conductivities. Detail
Correct diffusion coefficients for tke in LES. Detail

Chemistry

A building issues with Intel compilers (ifx/icx) for WRF/WRF-chem is addressed. Detail
Add KPP configure option for alternative libfl name, libfl.so, in addition to libfl.a. Detail
Fix n2o5 heterogenous namelist option for MADE schemes. Detail
A bug was introduced after splitting chem_opt 100 and 106. Currently, the chemistry initialization always calls for module_aerosol_soa_vbs routine, leaving the module_aerosol_soa_vbs_HET and corresponding data constants unitialized. As a result, aerosol concentrations are set to 0 after the first time integration (10**-16). This bug is fixed by Differentiate between module_aerosol_soa_vbs and module_aerosol_soa_vbs_het initialization routines. Detail
Update the registry description for settling velocity of bin 5 for GOCART. Detail

Data assimilation

Change an existing HDF5 lib name and add additional HDF5 libs to make sure WRFDA can be compiled in NCAR HPC Derecho when including RTTOV lib with HDF5. Detail
Fix Bufrlib C code to ensure that WRFDA can be built with intel oneAPI compiler. A run-time segfault bug is also fixed. Detail

Compiling

Add error checks in registry to catch code generation failure. Detail
Add the ability for registry code generation to be out of source. Detail
Fix missing or erroneous dependencies in make rules. Detail
Fix a compiling issue on Cygwin. Detail
Some versions of make appear to fail on the USENETCDFPAR logic in the top level Makefile if the variable is undefined. Thais issue is fixed by making sure it is at least set to 0 by the configure script fixes the issue. Detail
Many files are generated in the repo's directories and should be ignored to avoid accidental addition to commits.Detail
Update g2lib png encode/decode to use correct png_voidp define. Detail
Fix bugs in gen_allocs to use correct conditional compilation logic on allocatables vs pointers. Detail
Remove redundant code files wrf_status_codes.h, wrf_io_flags.h, io_int_idx_tags.h. Detail
Fix issue with generalized compiler version output using the wrong flags. Detail

Miscellaneous

Fix erroneous syntax in dyn_em/module_first_rk_step_part1.F. Detail
Loop order is rearranged (permute) to help memory access; it also breaks down larger loops into smaller ones to allow partial vectorization (GCC). [Detail](loop order (permute) to help memory access; it also breaks down larger loops into smaller ones to allow partial vectorization (GCC))
Adjust #defines to accommodate future CMake build. Detail
Add a print for parameters used when running real.exe using auto_levels_opt = 2 option, which is the default. Details

NOTE!!!! : If downloading the files from the section below, please choose either the v4.6.0.tar.gz file, or the v4.6.0.zip file. DO NOT choose those named "Source Code." They do not include the mandatory NoahMP submodule - needed for compiling WRF.

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by islas 12 months ago

WRF-ARW - WRF Version 4.5.2 (Bug-fix Release)

The WRF model has been updated to Version 4.5.2 on December 22, 2023.

Acknowledgements: We would like to thank Mathieu Landreau (LHEEA, Centrale Nantes, France), Marc Imberger (Technical University of Denmark), Alexander Knyazev (Intel), Wei-keng Liao (Northwestern University), Jun Tang (Amazon), Zhixiao Zhang (University of Oxford), Jacob Shpund (Pacific Northwest National Lab), Timothy Brown (AWS), Songyou Hong (PSL/NOAA), Dan Li (Boston University), David Robertson (Rutgers University), Cenlin He (NCAR RAL) for their contributions to this release.

This is a bug fix release. Following are fixes associated with this release.

Physics

Fix two errors in the WDM6 and WDM7 that have been in the code since V4.0: (1) uninitialized cloud number concentration when cloud is present at the start, and (2) The cloud autoconversion rate was off by a factor of 600. The effect of the first error isn't large, but the effect of the wrong autoconversion rate has resulted in doubling the surface rainfall in some tests (warm rain processes). Details
Fix a typo in computing saturated specific humidity in the single-layer urban canopy model. The effect is likely small since the calculation is relevant only in the case of rain and on impervious surfaces when using default evaporation scheme (which is set in URBPARM.TBL). Details
Fix the calculation order in the immediate melting, which was erroneous and based on the hydrometeor state after melting in previous version. In addtion, the calculation method in the supersaturation dynamical tendency is changed to avoid high values of total mass imbalance.Details
Fix a bug in the code related to trigger condition for CAPE calculation. Details
Add missing wrf_dm_bcast_real for stc2 in fitch parameterization scheme, which affects standing thrust coefficients for wind speeds at turbine hub-heights above cut-out. Details
Fix an issue in which time series for a parent domain were re-initialized when a child domain started later than the initial time for a simulation. Details

Compiling

Add build configuration for new Intel oneAPI LLVM ifx/icx compilers, which will be available on NCAR's new computer (Derecho). Details
The fortran library can be either hdf5_hl_fortran or hdf5hl_fortran for HDF5, which is now checked in Config.pl. Details
Improve the detection and linking with shared and static configurations of the NetCDF-4 libraries (split or unified) inside and outside Spack-Stack to avoid issues when NetCDF-C and NetCDF-Fortran reside in different directories Details
Add compile defaults for Arm Compiler (armclang/armflang) for Linux. Details
Resolves building issues with Intel compilers (ifx/icx). Details
Fixed an error in Makefile to check if USENETCDFPAR is set first. Details

Dynamics

Fix loop ranges in advect_w for lat-long map projection. This fix only affects polar boundary condition. Details

NOTE!!!! : If downloading the files from the section below, please choose either the v4.5.2.tar.gz file, or the v4.5.2.zip file. DO NOT choose those named "Source Code." They do not include the mandatory NoahMP submodule - needed for compiling WRF.

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by kkeene44 over 1 year ago

WRF-ARW - WRF Version 4.5.1 (Bug-fix Release)

The WRF model has been updated to Version 4.5.1 on July 25, 2023.

Acknowledgements: We would like to thank Stephen Sachs (AWS), Greg Thompson (NCAR RAL), Michael Laufer (Israel), Andrea Zonato (Koninklijk Nederlands Meteorologisch Instituut), Cenlin He (NCAR RAL), Becky Adams-Selin (Atmospheric and Environmental Research),Samm Elliott (TempoQuest Inc.) for their contribution to this release.

This is a bug fix release. Following are fixes associated with this release.

Physics

Correct the counter-gradient term for temperature equation at the first level in the KEPS scheme (bl_pbl_physics=17), which causes instabilities and large biases at the first model level over water surfaces. Details
Make Minor changes to Thompson-MP to align with changes made in CCPP for use in UFS. Details
Fix a bug related to FVEG scaling of canopy heat storage in Noah-MP canopy energy balance calculations. Details
Fix problems in HAILCAST quilted processes/tiling and some incorrect SI/non-SI unit conversions. Details
Fix a bug in building height specification in the LCZ urban parameter table. Details
Fix an issue in the revised MM5 (sf_sfclay_physics=1) scheme, where the model could potentially encounter an infinite loop. In specific conditions floating point roundoff errors were preventing a convergence condition from ever being met. Details

Compiling

Make configure script more portable. Details
Fix the ADIOS2 init API to support ADIOS2 v2.9.0. Details
Apply type-mismatch for gcc-10 to all gcc targets. Details

Miscellaneous

Replace a few wrf_messages and wrf_error_fatal calls by macros defined at the top of modules for urban schemes so that the code maintains consistency with those in MPAS. Details
Update references and remove obsolete contact information related to lightning options (lightning_option = 1 and 2). Details
Add User-friendly fatal error messages if the optional CESM aerosol data files required for the MSKF cumulus/Morrison microphysics aerosol interaction option (aercu_opt = 1 or 2) are missing. Details
Remove the file "test/em_real/namelist.input_default", which is more confusing than helpful for users. [Details] (https://github.com/wrf-model/WRF/commit/b55293062e400f82d7be0ebcf38f4ef87381ed83)

NOTE!!!! : If downloading the files from the section below, please choose either the v4.5.1.tar.gz file, or the v4.5.1.zip file. DO NOT choose those named "Source Code." They do not include the mandatory NoahMP submodule - needed for compiling WRF.

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by kkeene44 almost 2 years ago

WRF-ARW - WRF Version 4.5

This is the Version 4.5 release of the WRF Model. This release includes all components of the system: WRF-ARW, WRF-Chem, and WRFDA.

Acknowledgements: We would like to thank Cenlin He (NCAR RAL), Greg Thompson (JCSDA, UCAR), Anders Jensen (RAL, NCAR), Dylan Reynolds (SLF), Joseph Olson (NOAA-GSL), Wayne Angevine (NOAA-CSL), Andrea Zonato (Department of Environmental, Civil and Mechanical Engineering, University of Trento, Trento, Italy), Jun Tang (Amazon), Jeff Adie (NVIDIA Corporation), Hugh Morrison (NCAR-MMM), Brian Vanderwende (NCAR), Jordan Schnell (NOAA), Kengo Miyamoto (RIST), Zhang Chunxi, Maria Frediani (NCAR RAL), Michael Laufer (Toga Networks, a Huawei Company), Erick Fredj(The Jerusalem College of Technology, Rutgers University, Toga Networks, a Huawei Company), Joseph Brodie (Rutgers University), Lori Garzio (Rutgers University), Sen YANG, Deqin LI, Liqiang CHEN (Institute of Atmospheric Environment, China Meteorological Administration, Shenyang) for their contributions to this release.

New features in WRFV4.5

Physics

A new package to compute two-moment prognostics for graupel/hail and a predicted density graupel category is added in the Thompson scheme (mp_physics=38). Other updates to the scheme include a change to the Y-intercept relationship for one-moment graupel and replacement of air temperature for wet-bulb temperature in riming and mixed phase processes (Jensen et al., in review). This code requires the datafile "qr_acr_qg_mp38V1.dat" to run, which can be found on the file system for the NCAR community HPC "Cheyenne" under /glade/work/wrfhelp/WRF_files/ and online at http://www2.mmm.ucar.edu/wrf/src/wrf_files/. Details
A scale-aware or grid-distance-dependent option is added to the new Tiedtke scheme. The code is based on the paper by Wang (2022) with two minor changes removed. The option becomes active automatically when grid sizes is less than 15 km. PR1806 and PR1840
A new k-epsilon-theta^2 PBL scheme, including two additional prognostic equations for dissipation rate and temperature variance, has been included in the WRF model (see Zonato et al.,2022:A New K–ε Turbulence Parameterization for Mesoscale Meteorological Models. Monthly Weather Review, DOI: https://doi.org/10.1175/MWR-D-21-0299.1). This scheme is designed for finer grids. It works with sfclay=2(time step = 90), sfclay=1(time step = 60), diff_opt=1 km_opt=4). Details

Compiling

Add support to build WRF on ARM64 (aarch64) architecture systems with the NVIDIA HPC SDK Compilers (nvfortran,nvc). Details
Provided AMD AOCC Compiler support at the time of configuring WRF for AMD Architectures including Zen1, Zen2 and Zen3. Details

Bug fixes and enhancement in WRFV4.5

Physics

Substantial updates are added to the MYNN-EDMF, including a reorganization of the subroutines as part of the “universalization efforts”, which are meant to keep the scheme similar in all dycores. The major physics changes include (1) improved conservation of scalars which improves performance in medium-range global forecasting, (2) revised subgrid clouds which leverages the prognostic q’2, (3) vacate bl_pbl_physics = 6 (now only run with bl_pbl_physics = 5); will instead regulate the closure level with new namelist variable “bl_mynn_closure” set to 2.5, 2.6, or 3.0, (4) important bug fixes for the TKE budget and the bl_mynn_tkebudget flag is not renamed to tke_budget for use in other schemes (although may not be hooked up for other schemes yet), (5) new downdraft option (hidden hard-coded flag in module_bl_mynn.F) but still under development to best fit within the MYNN-EDMF, (6) and many miscellaneous tuning/minor bug fixes.Details
Fix an error with O31D exceeding vertical dimension kte. In WRFV4.4, we introduced ozone output o3rad when o3input = 0, that is to use ozone profile provided by the RRTMG code. But the use of output array accessed memory with the vertical dimension > kte. This can cause segmentation fault on some systems. Details
Update P3 microphysics. Details
Bug fix for the vertical levels of TKE in the EEPS scheme and more features are added in this scheme. Details
Fixed a bug for lightning diagnostic option 3, an option for lightning potential index calculation. The bug prevented nested case to run, and may produce NaNs in the diagnostics. Details
Fix an issue in roof longwave calculation in the bep_bem urban scheme (module_sf_bep_bem.F), in which a misplaced parenthesis is corrected. Details
The hard-coded tunable canopy heat capacity parameter in Noah-MP is added to MPTABLE for the convenience of tuning. Details
Fixed a dimension error for arrays refl, t and z in module_ltng_crmpr92.F.Details
Update hydro directory to match WRF-Hydro v5.3.0. Details
NoahMP submodule is updated to its v4.5 branch.Details

Compiling

Stop intentionally removing OpenMP flags when compiling certain listed source program files, which causes problems with a Fujitsu Fortran compiler. Details
Change configure to test for RPC package instead of checking hardcoded paths, so that non-standard RPC installs can still be found. Details
Fix numerical divergence on x86 and arm64. It is recommended to use armclang for best performance of WRF on arm64. Details
Correct string test in configure to avoid 'Unexpected Operator' error. Details

Data assimilation

Add a regularized version of WSM6 and its TL/AD for 4DVar with ice-phase hydrometeor analysis variables.Details
Enhance AHI radiance DA, including all-sky observation error model, Level-2 AHI product read, and more diagnostic output. Details

Chemistry

Address a bug introduced in WRF v4.3 that affects users using chem_opts = 201 (MOZART_MOSAIC_4BIN_KPP) and 202 (MOZART_MOSAIC_4BIN_AQ_KPP). Users experienced a stalled simulation upon writing history files. This was caused by the pH diagnostics added in v4.3 (#1342) not being activated with the namelist option mozart_ph_diag = 1. Details
Add an option to use GFS total ozone (2D) in the TUV photolysis (phot_opt=4) to scale the simulated column to the GFS total column in the calculation of photolysis rates. Details
chem_opt=100 is separated from chem_opt=108 and 109. This is because species in chem_opt=100 (HCL, CLNO2, MG, K, CA) are referenced by the aerosol module for chem_opt=108, 109 leading to memory mapping errors. Details

WRF-Fire

Update the firebrand spotting parameterization to (1) allow likelihood of spotting on short grass and (2) fix the location and height of ember release. Details

###Miscellaneous

Fix two numeric operators that are consecutive. Details
Change the exponent letters of the constants those are compared to double precision variables in SUBROUTINE FIND_RC0 from ‘E’ to ‘D’. Details
Remove excessive print statements in CLM. Details
Fix small indexing error in ISHMAEL MP scheme. Details
Check the existence of WRF_CMAQ. Details
Update REAME.namelist to include new options and missed pre-existing options.

NOTE!!!! : If downloading the files from the section below, please choose either the v4.5.tar.gz file, or the v4.5.zip file. DO NOT choose those named "Source Code." They do not include the mandatory NoahMP submodule - needed for compiling WRF.

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by kkeene44 about 2 years ago

WRF-ARW - WRF Version v4.4.2 (Bug-fix Release)

The WRF model has been updated to Version 4.4.2 on December 19, 2022.

Acknowledgements: We would like to thank Cenlin He (NCAR/RAL), Florian Sammüller (Universität Bayreuth), Adam Dury(WeatherQuest)) and Zhang Zhixiao (University of Oxford and University of Utah) for their contributions to this release.

This is a bug fix release. Following are fixes associated with this release.

Physics

Fix a bug for mixing ratio calculation in Jarvis stomatal conductance scheme in Noah-MP. This fix affects NoahMP option OPT_CRS=2. Details
Update urban local climate zone (LCZ) numbers in parameter tables to avoid overlap with existing NLCD land types. This update affects VEGPARM.TBL, LANDUSE.TBL, and MPTABLE.TBL. The LCZ numbers in WRF-urban are changed from 31-41 to 51-61. Details
Fix the units and unit conversions of several variables in subroutines CARBON and CARBON_CROP (NoahMP scheme); update urban pixel identification by including LCZ urban types in module_sf_noahmpdrv.F for surface variable initialization. Details
Refine CAPE and CIN calculations to be consistent with AMS definitions; improve LFC identification for adapting to multi-inversion layers profiles. Details
Change the declaration of several variables in subroutine TWOSTREAM from OUT to INOUT, change the soil ice treatment in subroutine COMBINE, correct the units of three accumulated flux variables in Registry, removed the hard-coded 4 soil layers and replace it by the parameter "NSOIL". These changes only affect NoahMP land surface scheme. Details
Change the declaration of QSFC from OUT to INOUT in sfclay and sfclayrev modules. This is because QSFC is calculated conditionally and at some grid points, it takes the input value. Details

WRFDA

Update the VAR namelist readme to include change of use_clddet option since V4.2. Details
Update bufrseviri.inc so that if seviri is turned on in WRFDA after Feb 2018, many error messages in the log file will be suppressed and meteosat 11 data is ingested into WRFDA to be used by the CRTM; correct the CRTM sensor information, add the HIRS4 information file to recognize NOAA 19 and not just NOAA 18 in HIRS4 and AMSUS files; update VARBC.in to include templates of NOAA 15-17 HIRS3 and NOAA 18-19 HIRS4. Details

Compiling

Fix undefined behavior in RSL_LITE.Detail

Miscellaneous

Add check to prevent CLM land surface scheme from being used with any urban scheme. Details

If downloading the files from the section below, please choose either the v4.4.2.tar.gz file, or the v4.4.2.zip file. DO NOT choose those named "Source Code." They do not include the mandatory NoahMP submodule - needed for compiling WRF.

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by kkeene44 over 2 years ago

WRF-ARW - WRF Version 4.4.1 (Bug-fix Release)

The WRF model was updated to Version 4.4.1 on August 25, 2021.

Acknowledgements: We would like to thank Jatin Kala (Murdoch Univeristy, Australia), Louis Marelle (LATMOS, CNRS; Paris, France), J. Shpund (PNNL), Jordan Schnell (NOAA), Robert Gilliam (US EPA), Tim Juliano and Maria Frediani (NCAR/RAL) for their contributions to this release.

This is a bug fix release. Following are fixes associated with this release.

Compiling

A modification addresses an issue that prevented WRF from compiling when choosing either the specified move (option 2) or vortex-following (option 3) nesting options, while using an Intel compiler. This issue began with WRFV4.4. Details
Remove duplicate reference to use module_dm in module_stoch.F (added in v4.4), which causes a compilation error when using an older version of the Intel compiler (12.1.3). Details
Add a check for status code -4001 from Cray CCE Fortran compiler related to reading CAMtr file. Details

Physics

Fix a minor bug for logical checking of time-varying aerosol option when climatological aerosol option selected. Details.
Recommend not to use the full spectral bin microphysics until it is updated. Use the fast version instead. Details
Update VEGPARM.TBL for NLCD40 landuse case to work with recent updates of the LSM initialization that reads the file. Details
Fix a bug associated with ozone input option o3input = 0. When o3input is set to 0, ozone is effectively not used in the RRTMG radiation option before this fix. Details.

Chemistry

Add gocart_dust dependency to mosaic_addemiss, which is necessary for WRF-CHEM to compile succeessfully. Details

If downloading the files from the section below, please choose either the v4.4.1.tar.gz file, or the v4.4.1.zip file. DO NOT choose those named "Source Code." They do not include the mandatory NoahMP submodule - needed for compiling WRF.

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by kkeene44 over 2 years ago

WRF-ARW - WRF Version 4.4

The WRF Model has been updated to Version 4.4 on April 26, 2022.

The WRF Pre-processing System (WPS) has been updated to Version 4.4.

Acknowledgements: We would like to thank Do Ngoc Khanh (Tokyo Institute of Technology), Ted Mansell (NOAA/NSSL), Joseph Olson (NOAA/GSL), Mike Iacono and John Henderson (AER), Dongmei Xu (Nanjing University of Information Science and Technology), Wen-Jou Chen,Ying-Jhen Chen, and Yi-Chuan Lo (Central Weather Bureau, Taiwan), Alexander Ukhov (KAUST), Marc Honnorat (EXWEXs), David Wong, (US EPA), Jaemo Yang (NREL), Patricia Balle (HPE), Tzu-Chin Tsai (National Taiwan University), Mathieu Le Corre (University of Toulouse, France, and University of Brest, France), Jordan Schnell (NOAA), Rebecca Schwantes (CIRES/NOAA CSL), and Prasanth Valayamkunnath, Alma Hodzic, Davide Del Vento, Masih Eghdami, Soyoung Ha, Jamie Bresch, Hugh Morrison, Craig Schwartz, Cenlin He, Rajesh Kumar, Ryan Cabell, Ju Hye Kim, Greg Thompson, Maria Frediani, Patrick Hawbecker, Timothy Juliano, and Pedro Jimenez (NCAR), and David Gill.

New in Version 4.4

Physics

The NoahMP code has been moved to an external repository and is linked to WRF repository via a submodule. This submodule is automatically pulled in to build the code during compiling, as long as the compilation is done with internet access. If compiling without internet access, or you are not sure if the computer you are using to compile supports git (git is used during the compiling process to access code in submodule), or in the rare case that one checks out another branch, other than ‘master,’ it may be necessary to issue the following commands while connected to the internet to ensure the submodule is available for compiling.
git clone https://github.com/wrf-model/WRF.git
git submodule update --init --recursive
Or in one command:
git clone --recurse-submodule https://github.com/wrf-model/WRF.git

After this, it is okay to go offline to compile. In the case you would like to remove code in the submodule, use ‘clean -aa’ instead of ‘clean -a’. Details

The shallow-water roughness scheme from Jiménez and Dudhia (2018) is included for offshore roughness adjustment in water depths less than 100 m and activated by namelist option shalwater_z0 = 1. The option works with a specified depth or real bathymetry input (namelist list option shalwater_depth) (Details). The bathymetry data is available from WPS/geogrid. Please acknowledge the following presentations and publications:
[1] GEBCO Compilation Group (2021) GEBCO 2021 Grid, DOI: 10.5285/c6612cbe-50b3-0cff-e053-6c86abc09faf
[2] Jiménez, P. A., and J. Dudhia (2018). On the need to modify the sea surface roughness formulation over shallow waters. J. App. Meteor. and Clim., 57(5), 1101-1110, DOI: 10.1175/JAMC-D-17-0137.1
WRF-Solar is expanded to have a stochastic ensemble prediction system (WRF-Solar EPS) tailored for solar energy applications (Yang et al., 2021, Kim et al., 2022). The stochastic perturbations can be introduced into variables of six parameterizations controlling cloud and radiation processes (Details). A more detailed description of the model is provided on the WRF-Solar EPS website.
References:
[1] Kim, J.-H., P.A. Jimenez, M. Sengupta, J. Yang, J. Dudhia, S. Alessandrini, and Y. Xie, 2022: The WRF-Solar Ensemble Prediction System to provide solar irradiance probabilistic forecasts. IEEE J. of Photovoltaics, 12, 141-142.
[2] Yang, J., J.-H. Kim, P.A. Jimenez, M. Sengupta, J. Dudhia, Y. Xie, A. Golnas, and R. Giering, 2021. An efficient method to identify uncertainties of WRF-Solar variables in forecasting solar irradiance using a tangent linear sensitivity analysis. Solar Energy, 220, 509-522.
Greenhouse gas (GHG) concentrations from a number of RCPs and newer SSPs are now a default run-time option (previously they were a compile-time option). This serves two purposes: 1) Since data files provide compiled global climatological values for co2, n2o, ch4, cfc11, and cfc12, up to 2006 for RCPs and up to 2014 for SSPs, they are better estimates for historical and current runs, 2) If users have values of their own, they can be easily added to the data file. The user specifies “ghg_input=1” in the &physics namelist record for climatology, which is now the default, or “ghg_input=0” for constant values for backward compatibility. The default file used is CAMtr_volume_mixing_ratio.SSP245. A simple function for CO2 is now the default when choosing not to use the climo GHG files for RRTM - previously this function was only in RRTMG schemes. This option is only available for radiation schemes CAM, RRTMG, RRTMG-Fast, and RRTM. Details
The option is added to output 16 accumulated physics tendencies for potential temperature, water vapor mixing ratio, and u and v components of wind. This is turned on with the namelist option “acc_phy_tend=1” in &physics. Details

Other New Options

The capability to couple WRF and CMAQ (an air quality model developed at the US EPA) is added. As a coupled model, WRF-CMAQ can be run in a one-way fashion (i.e., meteorological information will be transferred to CMAQ directly to drive its calculation), or in a two-way fashion, (i.e., the aerosol information from CMAQ will be fed back to the RRTMG radiation physics in WRF). This is available beginning in WRFV4.4 and CMAQ v5.3.3.3. Refer to the EPA’s WRF-CMAQ Model webpage for instructions. Details
The ability to write compressed NetCDF4 files in parallel, via NetCDF V4.7.4 (and later) is added. Performance is slower than pnetcdf, but can be notably faster than regular NetCDF on parallel file systems. As expected, compression provides files significantly smaller than pnetcdf generates. Details

Improvements and Bug Fixes

Physics

Microphysics

Thompson Aerosol-aware microphysics: Code enhancements are introduced to the processing of aerosol information in real.exe; specifically, legacy code to support the old format aerosol data before 3.9 is removed and logic in the real program is improved. Details
Thompson Aerosol-aware microphysics: A black carbon aerosol category is added. To include black carbon in this microphysics as well as its radiative effect, an updated data file, QNWFA_QNIFA_QNBCA_SIGMA_MONTHLY.dat, is required and set wif_input_opt = 2 (option 1 works with water/ice friendly aerosol data only). Support is also added for ingesting real-time data that contains these aerosols (such as from GEOS-5) by setting use_rap_aero_icbc = .true. (Details). Time-varying surface aerosol emission from either climatology or real-time data is now handled using namelist option qna_update = 1 with input file name wrfqnainp_d0* (must set auxinput17_interval and io_form_auxinput17). In addition, biomass burning organic and black carbon aerosols can be added to the model by setting namelist wif_fire_emit = .true.. This option will work with data with and without black carbon, but it is particularly important during active wildfire hence when using a first-guess aerosol source that has information about biomass burning emissions (e.g., GEOS-5). When using this option, it is assumed this aerosol is evenly distributed in the PBL (controlled by namelist wif_fire_inj = 1). Details
Thompson microphysics: An update to this, plus the cloud-fraction scheme “icloud=3” to better match observations. Modifications include updates to RRTMG LW and SW, and RRTMG fast LW and SW. Details
P3 microphysics: Four options (50, 51, 52, and 53) of the P3 scheme are updated. The 3-moment lookup table file size is reduced by an order of magnitude (file names are changed for the lookup tables). Improved flexibility has been added by including the Kogan (2013, JAS) option for cloud autoconversion/accretion. A few minor bug fixes to harden the code have been added (Details).
References:
[1] Morrison, H., and J. A. Milbrandt, 2015: Parameterization of cloud microphysics based on the prediction of bulk particle properties. Part 1: Scheme description and idealized tests. J. Atmos. Sci., 72, 287-311.
[2] Milbrandt, J. A., and H. Morrison, 2016: Parameterization of cloud microphysics based on the prediction of bulk particle properties. Part 3: Introduction to multiple free categories. J. Atmos. Sci., 73, 975-995.
[3] Milbrandt, J. A., H. Morrison, D. T. Dawson II, and M. Paukert, 2021: A triple-moment representation of ice in the Predicted Particle Properties (P3) microphysics scheme. J. Atmos. Sci., 439-458.
NTU microphysics: Fixed uninitialized variable BHTMP. Details.
NSSL microphysics: Fixed a few uninitialized variables. Details

Radiation

A new cloud overlap option is available for RRTMG radiation schemes. The option is set by using the exponential-random namelist parameter “cldovrlp=5.” A decorrelation length option, “idcor,” is also added. It can be set as a fixed value at 2500 m (idcor=0), or a function of latitude (idcor=1). The decorrelation option also works for cloud overlapping option 4 (exponential). Details
KIAPS Longwave Radiation: In V4.2, an error was introduced when attempting to fix a compilation failure, but ended up introducing yet another compile-time failure. The solution is to revert the entire KIAPS longwave radiation scheme to the original pristine condition. Details

PBL

Shin-Hong PBL: A minor bug fix to avoid a rare divide by zero issue. Details

LSM

NoahMP LSM: Updated with the following new development: Details

Added three runoff schemes (opt_run = 6,7,8).
Added a new drainage scheme (opt_tdrn)
Added new soil step control capability (namelist soiltstep)
Added water and energy budget terms for output and diagnosis (namelist noahmp_output)
Updated flood irrigation parameters
Added constraints for snow thickness during compaction
Updated snow viscosity coefficient in compaction process
Added canopy heat storage in canopy temperature calculation
Updated formulation and parameter values for TOPMODEL groundwater scheme (opt_run=1)
Bug fixes for SWE cap value in initialization, NoahMP snow depth update, snow layer index update, floating issues in runoff scheme, and BATS albedo scheme.
NoahMP Glacier: An issue is resolved that occurred when there was no sunlight (cosz<0), causing the NoahMP BATS snow albedo calculation using cosz to give a “divided by zero” error in glacier code. Details
NoahMP: Corrected declaration for two variables, FV and USTAR, from ‘out’ to ‘inout’ ones. Details

A fix addresses a runtime segmentation fault when using “sf_surface_mosaic=1” with “use_wudapt_lcz=0.” Details
Noah LSM: Fixed memory error that may cause occasional model failure when WRF-Hydro is not used. Details
Noah LSM: Fixed two input variables of absorbed longwave radiation LWDN and net shortwave radiation SOLNET for non-urban fraction of the LSM call. Details
RUC LSM: The default method to diagnose 2 m temperature and mixing ratio is changed to the flux method in RUC LSM, the same as used in NCEP operational RAPv5/HRRRv4. Details

Surface Layer

Updates and bug fixes to the MYNN surface layer scheme to bring this scheme closer to the CCPP version (Details), described in Olson et al. (2021):
Olson, Joseph B., Tanya Smirnova, Jaymes S. Kenyon, David D. Turner, John M. Brown, Weizhong Zheng, Benjamin Green (June 2021): A description of the MYNN surface layer scheme. NOAA Technical Memorandum. OAR GSL-67. pp. 33., DOI:10.25923/f6a8-bc75

Urban

A minor bug fix for urban wind directions when using detailed urban morphology maps. Details
An issue that caused an underestimation of the rain rate used in urban physics is resolved. Details
Fixed an uninitialized emissivity variable in Noah driver for urban option. Details

Orographic Gravity Wave Drag

The description for the 10 orographic gravity wave drag variables (gwd_opt=1 only) are updated. Details

Others

Occasionally when running real.exe on multiple processes with MPI, one or more processes crashed with error “Unrecognized physics suite.” A fix for the MPI synchronization bug related to (not used) split communicators in the real program provides a solution to this issue. For users that had no trouble with this problem, this fix will have no impact. Details
To address an issue with compiling WRF on AWS cloud computing platforms (and perhaps other OSs), with Intel, a work-around is introduced to include the unix command “true” as the last command in a string of registry-based commands. This allows WRF to build on AWS with the Intel compiler. For other users, this will have no impact. Details
The configure file is modified to address a compiling error when compiling the vortex-following “hires” option related to the landread.c file. Details
Beginning with V4.2.2, a bug was introduced that prevented users from configuring the model with the “INTEL (ifort/icc): HSW/BDW” option. This is now resolved. Details
The NMM dynamical core and some NMM-specific source code are removed from the WRF repository. Details
A new variable that accumulates GHI over time is added to the WRF-Solar diagnostics package. Cloud base and cloud height units (AGL) are clarified in the registry. Details

WRF-Fire

Added the ability to customize vertical levels for fire runs. Vertical nesting is still not supported. Details
A new module to parameterize firebrand spotting for WRF-Fire is added. This is a passive Lagrangian transport scheme to transport and burnout firebrands generated at the fire front. The scheme is activated when ifire==2 by setting the namelist option fs_firebrand_gen_lim to an integer greater than zero (default is 0, i.e., scheme is off). It runs with dmpar- and serially-compiled code and in the inner nest (grid_id==max_dom). It was designed and tested using a mesoscale to LES domain configuration. Details

WRFDA

WRFDA is now able to ingest NCEP’s satwnd.bufr (gdas.satwnd.tHHz.YYYYMMDD.burf) to assimilate more AMVs that are not included in NCEP’s prepbufr files. More flexible thinning/superobbing options (thin_conv_opt) and obs error specification (uv_error_opt) are implemented for using satwnd.bufr. See WRFDA/var/README.namelist for details. Details
A new capability to assimilate GPM-GMI radiance data is introduced.
Reference:
Shen, et al., 2021: Assimilation of GPM Microwave Imager Radiance data with the WRF Hybrid 3DEnVar System for the Prediction of Typhoon Chan-hom (2015), Atmospheric Research. 251, 105422. Details
WRFDA is newly-expanded for “chem_cv_options=108” for surface chemical data assimilation (PM2.5, PM10, CO, SO2, NO2, O3) using the RACM/MADE-VBS chemistry option in WRF-Chem. Details
Reference:
Ha, Soyoung, 2021: Implementation of aerosol data assimilation in WRFDA (V4.0.3) for WRF-Chem (V3.9.1) using the MADE/VBS scheme. Geoscientific Model Development. gmd-2021-74.
The WRFDA code is corrected to reduce memory usage for “ob_format=1” (i.e., bufr for conventional obs) and “thin_conv=.true.” Details
WRFDA is fixed and enhanced for applications related to offline varbc (use_background_error=.false.), 3DEns-FGAT (num_fgat_time>0 and ensdim_alpha>0 and je_factor=1.0), precondition_cg=.true. and orthonorm_gradient=.true. when ensdim_alpha>0, full-ensemble mode (ensdim_alpha>0 and je_factor=1.0). Details
WRFDA gfsref code is restructured to move QC steps in da_get_innov_vector_gpsref.inc to a separate subroutine in a new file, “da_qc_gfsref.inc.” New namelist variables are added to allow more flexibility for testing gpsref QC impact. Several issues with gpsref diagnostic output have been fixed and qcstat information for gpsref is now in new qcstat_gpsref_xx. Details
Modifications to address compilation errors due to aqueous variables in da_get_innov_vector_chem_sfc.inc. Details
Fixed an issue that occurred when assimilating with pm2.5 (e.g., chemicda_opt = 3 or 5), using chem_opt=108. pm10 is treated as (pm10 - pm2.5) residuals in both observations and priors. PM10 obs errors are also adjusted for such residuals. Details
Bug fixes for direct assimilation of radar reflectivity. The namelist parameter “qthres,” which is used for the hydrometeor base state of TL/AD for no-rain background, should not be used in the reflectivity forward operator. QC for “radar_rf_opt=2” should not be confused with that of retrieval radar DA. Details
An issue is resolved related to inconsistency between TAMDAR and Radiance VarBC code. Details

WRF-Chem

Subgrid boundary layer clouds from the MYNN PBL scheme are now allowed to be accounted for in photolysis, using “phot_opt=4” (TUV). This feature is activated by setting the &chem namelist option “phot_blcld=.true.” “icloud_bl=1” must also be set. Details
Minor bug fix for aerosol optics Mie extrapolation cases in chemistry. When aerosol refractive indices are outside the lookup table ranges, extrapolation in the parameterized Mie calculation will create unrealistic values for the extrapolated aerosol optical parameters. These unphysical fluctuations cause unrealistically high AOD values and the model crashes. The solution is to force the out-of-range aerosol refractive index to be the minimum or maximum boundary values of the lookup table. Details
Modifications address an issue in which HONO, TERP and CO2 were not automatically pulled into the RACM mechanisms, even though they were in the anthropogenic emissions file. Details
An issue is corrected related to inconsistent data units for gas species (chem_cv_options). Details

WRF-Hydro

WRF-Hydro has been updated to version v5.3. Details

If downloading the files from the section below, please choose either the v4.4.tar.gz file, or the v4.4.zip file. DO NOT choose those named "Source Code." They do not include the mandatory NoahMP submodule - needed for compiling WRF.

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by kkeene44 about 3 years ago

WRF-ARW - WRF Version 4.3.3 (Bug-fix Release)

The WRF model has been updated to Version 4.3.3 on January 11, 2022.

This is a bug fix release. Following are fixes associated with this release.

Physics

The v4.3.2 release contained an error in the commit that addressed the situation of segmentation faults when running a PBL on a coarse grid and no PBL on the fine grid. If a user is running the option to get the time series output, then the model will always report that there are PBL + LES issues if the user is choosing the MYNN or EEPS PBL schemes. This bug is fixed in the v4.3.3 release.

Chemistry

For WRF Chem, ISORROPIA in chem_opts = 100,108,109 is now controlled by a logical namelist option do_isorropia
Output photolysis rates for NO2 and O1D are included in history stream by default
An indexing bug that incorrectly enhanced vertical mixing at the surface layer in the WRF Chem dry_dep_driver.F filei is fixed.

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by smileMchen over 3 years ago

WRF-ARW - WRF Version 4.3.2 (Bug-fix Release)

The WRF model has been updated to Version 4.3.2 on 14 December, 2021.

Acknowledgements: We would like to thank Patricia Balle (HPE), Jamie Bresch (NCAR), Jordan Schnell (NOAA), Yaping Shao (Universität zu Köln), Martina Klose (Karlsruhe Institute of Technology), Piotr Kasprzyk (IETU Katowice), Theodore M. Giannaros (National Observatory of Athens, Greece)

This is a bug fix release. Following are fixes associated with this release.

Compiling

Added a new stanza that is only for ARM processors (right now aarch64 and armv7l) with GNU.

Physics

Added time-varying greenhouse gases to RRTMG shortwave options: RRTMG and RRTMG-fast.
A test was introduced to stop the incorrect combination of particular PBL schemes on coarser grids with the LES PBL option selected on a finer grid. Previously, this problem caused an inconsistent number of variables on the CG and FG, which caused segmentation faults when trying to do feedback or advection of unavailable fields.
A few syntax errors in NSSL, NTU, an NoahMP were fixed. Most compilers skipped over them. Either "small impact" or "no impact" for users for whom the code already compiled.

WRF-Chem

In WRF Chem, remove photolysis rates from the default history stream (these photolysis rates are not used by most modules).

WRF-fire

Bug fix for fire module related to the initialization of 1000-h fuel moisture content due to an incorrect indexing assignment from the namelist entries for fuel classes.

Registry

Starting with release-v4.3, the CLM specific variable PCT_PFT appeared in the standard wrfout file, even when the CLM option was not selected. A small modification to the standard Registry file removed this incorrectly output 3d array. Other than removing this innocuous zero-valued 3d array from the output stream, there is no other impact for users.

Miscellaneous

A fatal error in the real program is now issued when a user requests lateral boundary topography smoothing, but lacks the first-guess soil elevation data to do the weighted averaging of topography. Previously, this discrepancy was ignored, and the d01 high-resolution elevation data was quietly left undisturbed. We have seen this problem when a user chooses to initialize the model with ERA5 data, but when insufficient time invariant files are used.

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by kkeene44 over 3 years ago

WRF-ARW -

The WRF model has been updated to Version 4.3.1 on October 28, 2021.

Acknowledgements: We would like to thank Jesus Fernandez (Instituto de Fisica de Cantabria CSIC-UC, Santander, Spain), Andrea Zonato (University of Trento), ed Mansell (NOAA/National Severe Storms Lab), Hugh Morrison (NCAR), Matthias Demuzere (Ruhr University Bochum) and Piotr Kasprzyk (IETU Katowice) for their contributions to this release.

This is a bug fix release. Following are fixes associated with this release.

Compiling

module_ra_eclipse.o is added to radiation driver's dependency list to avoid occasional compilation failures.

Physics

Bug fix in Noah LSM routine to initialize several local variables. Without their initialization, the same case run with different number of processors cannot yield identical results.
Fix in initialization for rand_perturb=1 option for the stochastic capabilities. The error is caused by an incorrect argument list.
New CAMtr_volume_mixing_ratio data files for CMIP6 SSPs are added to acount for yearly global average concentrations from input4MIPs (v1.2.0 for historical part and v1.2.1 for scenario evolution).
The photovoltaic panel (PVP) parameterization scheme, new in version 4.3, has been updated, since the previous one provided unreasonable results due to wrong calculation of the radiation balance. While before a simple energy balance was adopted to calculate the PVP temperature and related heat fluxes, now a prognostic equation is solved, accounting for the thermal coefficients of the photovoltaic panels.
For the NSSL 2 moment MP scheme (mp_physics=17, 18, and 22), the formulation of snow radar reflectivity is reverted back to version (pre-4.1.3) using size-dependent snow density. The diagnosed bright band is still turned on by default: iusewetsnow=1. (Setting iusewetsnow=3 will switch back to old formulation.)
Fix an initialization problem in Thompson scheme (mp=28). Since V4.0 the aerosol aware Thompson (mp=28) option has not been initializing the surface aerosol source array consistently between a parent and nested domain due to grid-size dependence. This has been removed and the results will change if dx is not 20 km, but now nests will have the same source.
An error in longwave radiation calculation in the BEP-BEM urban module is fixed. This error affects the calculation of TSK, which now is correctly calculated.
The P3 microphysics scheme is updated with bug fixes and code cleanup.
Units specified in Registry are corrected for several NTU microphysics variables.
Dimensional errors in subroutine spec_bdytend_perburb and spec_bdytend_perturb_chem are fixed (in module_bc.F). This routine is used to perturb fields on the model lateral boundaries for stochastic processing, and only affects results if perturb_bdy or perturb_chem_bdy are not 0 and they are 0 by default.
Some of the urban parameter values for LCZ are changed to make them more generic.

Data Assimilation

Logical variable got_var_sso is added to the history file for cycling DA application for use with the topo_wind option.

Registry

A temporary fix has been put in place for the moving nest capability when using the Intel compiler. Previously, the WRF source code did not build due to a segmentation fault in the execution of the tools/registry program.

Miscellaneous

Increase total time periods available to process for climate simulations. The restriction of 10000 time periods in a single file has been increased to 60000 (a value that allows 3-hourly data for 20 years).
Revert print and when to turn on damping to pre 4.3 code. In WRFV4.3, w_damping is turned on at CFL = 1.2, and CFL is printed when it exceeds 1.2. This results in many CFL prints. In the release of WRFV4.3.1, w_damping starts at CFL = 1.0, and CFL prints appear when CFL exceeds 2.0. When IEVA is on, w_damping starts at CFL = w_crit_cfl set in the namelist, and CFL prints only appear for those exceeding 2.0.
The description of option MAD-WRF is improved

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by smileMchen over 3 years ago

WRF-ARW - WRF Version 4.3

The WRF Model has been updated to Version 4.3 on May 10, 2021.

Acknowledgements: We would like to thank Matthias Göbel (University of Innsbruck), Jaemo Yang and Yu Xie (NREL), Will Hatheway, Douglas Lowe (University of Manchester), Ted Mansell and Louis Wicker (NOAA/NSSL), Sam Elliott (TQI), Megan Bela and Stu McKeen (CU Boulder CIRES/NOAA CSL), Xinzhong Liang (University of Maryland), Tzu-Chin Tsai and Jen-Ping Chen (National Taiwan University), Robert Gilliam and Jonathan Pleim (US EPA), Alexander Ukhov (KAUST), Zhixiao Zhang (University of Utah), Adam Varble (PNNL and University of Utah), Katelyn Barber and Brian Gaudet (PNNL), Robert Arthur, Katie Lundquist, and Jeff Mirocha (LLNL), Sam Levis (SLevis Consulting), Andrea Zonato (University of Trento, Italy), Michael Toy (NOAA/GSL), Chunxi Zhang (NCEP) and Yuqing Wang (University of Hawaii), Eric A. Hendricks (NCAR/NSAP), James Ruppert (Penn State University), Isaac Rowe (University of Kentucky), Alex Montornès (University of Barcelona), Yunyao Li and Xin-Liang Zhong (University of Maryland), Jan Mandel (University of Colorado, Denver), Stacy Walters (formerly NCAR), Luke Conibear (University of Leeds), Sonia Lasher-Trapp (University of Illinois), T. Iguchi (NASA/Goddard), Wei Sun (Chinese Academy of Sciences), Greg Thompson (JCSDA), Prasanth Valayamkunnath, Pedro Jimenez, Jamie Bresch, Craig Schwartz, Jim Bresch, Cenlin He, Hugh Morrison, Masih Eghdami, Timothy W. Juliano, Negin Sobhani, Dave Lawrence, Bill Sacks, Jared Lee, Laura Fowler, and Mary Barth (NCAR)

Note: WRF v4.3 is the last release of the WRF model that will support the NMM dynamical core and HWRF capability. Beginning with WRF v4.3.1 and with all subsequent releases, the NMM and HWRF capabilities will incrementally be removed from the model source code.

New in Version 4.3

Physics

A turbulence kinetic energy (TKE) and TKE dissipation rate (ε) based 1.5-order closure PBL parameterization (E–ε, EEPS) is added (Zhang et al. 2020, MWR). Works with surface layer options, 1, 91, and 5.
An updated version of P3 now includes a one-ice category, 3-moment ice option. Milbrandt et al. (2021) [Milbrandt, J. A., H. Morrison, D. T. Dawson II, and M. Paukert, 2021: A triple-moment representation of ice in the Predicted Particle Properties (P3) microphysics scheme, J. Atmos. Sci., 78(2), 439-458, https://doi.org/10.1175/JAS-D-20-0084.1]
The NTU (National Taiwan University) scheme (mp_physics = 56). It applies double moments for the liquid-phase and triple moments for the ice-phase hydrometeors together with the consideration for ice crystal shape and density variations. There are five major features to NTU scheme: condensation nuclei (CN) and ice nuclei (IN) are tracked separately for the processes of cloud/rain activation and ice deposition-nucleation using predicted supersaturation; applying the triple-moment (the zeroth, second, and third moments) closure method to describe the evolution of ice particle’s spectrum; solid-phase hydrometeors’ classification (pristine ice, snow aggregate, rimed ice, and hailstone) is redefined according to their key formation mechanisms; ice crystals’ shape and apparent density can evolve gradually according to the growth conditions; and fall speed of each moment for frozen particles depends on shape and density.

For the detailed parameterizations in the NTU scheme, the liquid-based formulae of the bulk conversion rates are adopted from [1], and the technique for aerosol activation to rain/cloud is introduced in [2] and [3]. The triple-moment approach, representation of ice properties, mixed- and solid-based equation/kernels are illustrated in [4] and [5].

References:
[1] Chen, J.-P. and S.-T. Liu, 2004: Physically based two-moment bulkwater parameterization for warm-cloud microphysics. Quart. J. Roy. Meteor. Soc., 130, 51–78, doi:10.1256/qj.03.41.
[2] Cheng, C.-T., W.-C. Wang, and J.-P. Chen, 2007: A modeling study of aerosol impacts on cloud microphysics and radiative properties. Quart. J. Roy. Meteor. Soc., 133, 283–297, doi:10.1002/qj.25.
[3] Cheng, C.-T., W.-C. Wang, and J.-P. Chen, 2010: Simulation of the effects of increasing cloud condensation nuclei on mixed-phase clouds and precipitation of a front system. Atmos. Res., 96, 461–476, doi:10.1016/j.atmosres.2010.02.005.
[4] Chen, J.-P. and T.-C. Tsai, 2016: Triple-moment modal parameterization for the adaptive growth habit of pristine ice crystals. J. Atmos. Sci., 73, 2105-2122, doi:10.1175/JAS-D-15-0220.1.
[5] Tsai, T.-C. and J.-P. Chen, 2020: Multi-moment ice bulk microphysics scheme with consideration for particle shape and apparent density. Part I: Methodology and idealized simulation. J. Atmos. Sci., 77, 1821–1850, doi:10.1175/JAS-D-19-0125.1.
Added a new option (ra_sw_eclipse) to model the effect of eclipses based on the Bessel's method in four radiation schemes: RRTMG, Dudhia, Goddard, and old Goddard. Additional output includes the location of the eclipse: ELAT_TRACK and ELON_TRACK. For reference, see this paper: Montornès, A., Codina, B., Zack, J. W., and Sola, Y.: Implementation of Bessel's method for solar eclipses prediction in the WRF-ARW model, Atmos. Chem. Phys., 16, 5949–5967, https://doi.org/10.5194/acp-16-5949-2016, 2016.
A new orographic gravity wave drag option is added (gwd_opt = 3). The scheme includes large-scale orographic gravity wave drag in a similar way as in gwd_opt = 1. In addition, two subgrid-scale sources of orographic drag are added in this option. The additional schemes are applicable at coarse horizontal grid resolution (~100s km) down to fine resolutions on the order of 1km, and they account for small-scale ~1km topographic variations. Setting namelist option gwd_diags = 1 (default is "0") causes diagnostic data to be output to the wrfout_d0x files. The data include: 1) 2D vertically integrated momentum fluxes, i.e., surface stress, DUSFC_xx (x-direction) and DVSFC_xx (y-direction), where xx={LS,BL,SS,FD}, which are "large-scale", "blocking", "small-scale" and "form drag" contributions, respectively; and 2) 3D drag forces DTAUX3D_xx (x-direction) and DTAUY3D_xx (y-direction), where xx={LS,BL,SS,FD}. This scheme is contributed by NOAA/GSL.
The multilayer BEP (Building Effects Parameterization) and BEP+BEM (BEP with the Building Energy Model) urban canopy models (UCMs) are added to the Yonsei University (YSU) planetary boundary layer parameterization. Reference: Hendricks, E. A., J. C. Knievel, and Y. Wang, 2020: Addition of multilayer urban canopy models to a nonlocal planetary boundary layer parameterization and evaluation using ideal and real cases, J. Appl. Met. Clim., 59, 1369-1392.
Implementation of Local Climate Zone (LCZ using WUDAPT (31-41) landuse classes, along with standard urban classes (31-33). (Ref: Stewart, I.D. and Oke, T.R. (2012) Local Climate Zones for Urban Temperature Studies. Bulletin of the American Meteorological Society, 93, 1879-1900. http://dx.doi.org/10.1175/BAMS-D-11-00019.1). More information can be found at https://ral.ucar.edu/sites/default/files/public/product-tool/urban-canopy-model/WRF_urban_update_Readme_file_WRF4.3.pdf. Data from https://www.wudapt.org/ is required. Note that the parameters in the new urban table, URBPARM_LCZ.TBL, may vary greatly from city to city. The default values are probably not appropriate for any given city. Users should adapt these values based on the city they are working with.
WRF-urban updates for green roof, solar panel, and new building drag coefficient for BEP+BEM. The paper describing the green roof update is under development. The new building drag coefficient is based on Santiago and Martilli (2010) and Gutierrez et al. (2015).
[1] Santiago, J. L. and Martilli, A. (2010). A Dynamic Urban Canopy Parameterization for Mesoscale Models Based on Computational Fluid Dynamics Reynolds-Averaged Navier-Stokes Microscale Simulations. Boundary-Layer Meteorology, 137(3):417-439.
[2] Gutí errez. E., Martilli, A., Santiago, J. L., and González, J. E. (2015). A Mechanical Drag Coefficient Formulation and Urban Canopy Parameter Assimilation Technique for Complex Urban Environments. Boundary-Layer Meteorology, 157(2):333-341.
[3] A. Zonato, A. Martilli, E. Gutierrez, F. Chen, C. He, M. Barlage, D. Zardi, and L. Giovannini (2021): Exploring the effects of rooftop mitigation strategies on urban temperatures and energy consumption, Atmospheric Chemistry and Physics (under review).
A new dynamic irrigation scheme is implemented in NoahMP to estimate irrigation water requirements and apply water in the irrigated croplands. This new scheme irrigates the croplands based on sprinkler, micro, or surface flooding methods. The various options are selected from the namelist, and are described in the run/README.namelist file.
The capability to couple Community Terrestrial Systems Model (CTSM) with WRF via LIghtweight Land Atmosphere Coupler (LILAC) is added and activated by setting namelist option sf_surface_physics to 6. This is the initial beta release of WRF-CTSM coupling capability. For instructions on how to run WRF with CTSM please check instructions on using CTSM with WRF. See the WRF-CTSM User’s Guide (https://escomp.github.io/ctsm-docs/versions/master/html/lilac/specific-atm-models/wrf.html). Questions regarding this capability can be addressed to the CTSM Forum (https://bb.cgd.ucar.edu/cesm/forums/ctsm-clm-mosart-rtm.134/).

Other New Options

The Implicit-Explicit Vertical Advection (IEVA) scheme has been implemented that permits a larger time step by partitioning the vertical transport into an explicit piece, which uses the normal vertical schemes present in WRF, and an implicit piece which uses implicit transport (which is unconditionally stable). The combined scheme permits a larger time step than has previously been used, and reduces w-filtering. The scheme will be useful for CONUS-scale CAM (convection allowing model) simulations (dx ~ 2-3 km) when the number of vertical levels > 50. In these cases, time steps can increase to as large as 25 s, depending on the problem. Overall integration efficiency increases ~ 15%, and the IEVA solutions are closer to a benchmark run using a smaller time step. The option is activated by setting zadvect_implicit = 1 in &dynamics.
Note: Must set reasonable_time_step_ratio to enable the use of larger time steps. ( Wicker, L. J., and W. C. Skamarock, 2020: An Implicit–Explicit Vertical Transport Scheme for Convection-Allowing Models. Mon. Wea. Rev., 148, 3893–3910)
A new namelist, w_crit_cfl, is added in &dynamics. This is the threshold CFL when W damping starts if w_damping option is on. When IEVA is activated, w_crit_cfl may be set to 2.0.
MAD-WRF: A solar irradiance nowcasting system that combines the benefits of data assimilation and the philosophy of advecting cloud properties of MADCast, and the cloud physics of WRF-Solar (https://ral.ucar.edu/projects/mad-wrf)
A new README file (README.physics_files) is now available in the run directory (and is linked to the test/em_real directory during the compile). This file lists all physics options that require additional files for running, and the files necessary for each scheme. These additional files are already in the run/ directory and are linked automatically to the test/em_real directory during compile.
In combination with updates to the cloud fraction scheme (icloud=3), a new feature permits the usage of the cloud fraction scheme as part of real.exe, specifically designed to use with "cold start" simulations to reduce the spin-up problem of clouds and associated radiation.
The git commit information is now part of the WRF modeling system standard output, one of the strings within the executables, and also in the output of the ./compile command.
Air pressure was added as a new height-level diagnostic, which is useful for a model validation against barometers mounted above ground (e.g., on a wind turbine nacelle). Also, the vertical interpolation of pressure to height surfaces was switched to linear in ln(p). The resultant pressure and relative humidity changes are negligible in the lowest 1 km AGL.
A new namelist option is included that allows a user to avoid double stagger averaging of omega (coordinate velocity) in the vertical advection of geopotential. The new namelist option phi_adv_z allows switching between the two formulations: phi_adv_z = 1 (the old formulation, which remains the default), and phi_adv_z = 2 (the new formulation).

Improvements and Bug Fixes

Physics

Microphysics

NSSL microphysics updates:
- Enabled regeneration of CCN by droplet evaporation and by nudged restore in cloud-free air (mp_physics=18; default 1-hr time constant for nudging)
- Enabled radar reflectivity from cloud ice
- Added internal option for ice crystal nucleation by DeMott et al. (2010, PNAS)
- Allow greater fraction of hail to melt in one time step
- Reduced minimum number concentration (based on CAPS input)
- Increased resolution of lookup table for incomplete gamma functions
- Bug fix to limit fall speed air density factor for high model lids
- Fixed issue of spurious creation of large concentrations of very small droplets and transient large condensation
- Fixed issue of negligible "seed" values of graupel
- Fixed effective radius calculation of snow
Thompson MP scheme: Fixed a problem with a single constant that caused approximately 4x too much fake surface aerosol emissions. Relatively small tunings for one-moment graupel Y-intercept parameter diagnosis, creating a better fit to Paul Field et al (2019) analysis of T-28 aircraft hail data.

Cumulus

Deng shallow cumulus scheme: Check is added for water vapor mixing ratio (Qv), if negative, setting it to a small positive number. Previously, the scheme died when taking the log of negative Qv.
Deng shallow cumulus scheme: When the updraft top reached the domain top, the model would crash due to accessing out of bounds memory. The modification is to limit the vertical index of the updraft top.
Fixed staggering of shallow cumulus tendencies of u and v at periodic boundaries with a PERIOD communication.

Radiation

A fix is introduced to the RRTMG SW parameterization. When the model top is relatively low, and the threshold for tropopause pressure is not reached, laytrop = nlayers. In this case, some indices in sfluxzen can be filled with junk, ultimately leading to NaN values for SWDOWN and failure in the radiation driver. The modification is to initialize sfluxzen to zero before being filled.
Three-dimensional clear sky radiative heat tendencies are now available as output variables when using the following radiation schemes: RRTM longwave (ra_lw_phys=1), CAM shortwave and longwave (ra_sw_phys, ra_lw_phys=3), RRTMG shortwave and longwave (ra_sw_phys, ra_lw_phys=4), RRTMG-fast version shortwave and longwave (ra_sw_phys, ra_lw_phys=24), and RRTMG-K shortwave and longwave (ra_sw_phys, ra_lw_phys=14). Similar to the all-sky longwave and shortwave radiative heat tendencies (RTHRATLW and RTHRATSW), the clear-sky tendencies are not included in the default history output stream, but can now be added to it. They can be found in Registry.EM_COMMON as RTHRATLWC for longwave clear-sky heating rate and RTHRATSWC for shortwave clear-sky heating rate.
Fixed a bug in surface downward diagnostic output of long- and short-wave fluxes and two other bugs involving wrong numerical values used in the RRTMG-K code.

PBL

UW PBL: input radiative tendency corrected to be T tendency instead of theta (small effect of temperature vs theta, up through the top of the PBL).
The namelist option ysu_topdown_pblmix for the YSU PBL scheme has been turned on by default in the Registry files. To turn it off, set ysu_topdown_pblmix=0 in the &physics namelist record (single entry).
A minor fix to BouLac PBL tendency cleans up the assignment, no longer assuming that the tendency is cumulative (no effect on model results).
Modifications to the ACM2 PBL scheme for rare crashes. The overall effects of these changes are generally small although the indexing fix will prevent model crashes in extreme terrain (e.g. Tibet).
Fixed a few bugs for the 3DTKE scheme: 1) Fixed a bug when converting surface heat flux from dry to moist theta when use_theta_m is enabled, which is the default; 2) Fixed divide by zero; 3) Fixed computed vertical index that is out of range.

LSM

NoahMP:
- Corrected code that forced users to turn on the crop scheme when using the irrigation scheme (iopt_irr=2). Now the irrigation scheme can be used without turning on the crop scheme.
- Change local variable QRAINXY (rain rate on the ground after canopy interception) to output variable.
- Several snow-relevant enhancements and new features to improve snow simulations based on NCAR Noah-MP LSM team tests. These are: (1) new snow retention process improves streamflow modeling (tested with Noah-MP in National Water Model); (2) optimized parameter values for snow cover parameters (SCFFAC & MFSNO) improve snow water equivalent (SWE), snow depth, and surface albedo simulation; (3) new snow emissivity value improves WRF surface temperature simulation; (4) new wet-bulb temperature snow-rain partitioning scheme (opt_snf=5) improves SWE simulation see Wang et al. 2019 GRL; (5) enhanced flexibility of tuning model parameters by bringing the above-mentioned parameters from hard-coded values to MPTABLE; (6) Fix the incorrect values for wind-canopy absorption coefficient (CWPVT) parameter and the unrealistic value range of bulk leaf boundary layer resistance (RB). (He et al., 2021, in preparation)
Modifications to the Pleim-Xiu LSM. The overall effects of these changes are generally small, where the biggest effects on surface statistics was the added evaporation from the ground in vegetative fraction of the grid cell. This, by nature, increases near surface water vapor mixing ratio. In some of our tests this increased the bias in some areas.

Surface Layer

Several surface layer schemes are made thread-safe if the fractional seaice option is used. These are revised MM5, original MM5, MYJ, QNSE and PX surface layer scheme (sf_sfclay_physics = 1, 91, 2, 4, and 7).
Fixed a problem that caused the input field “wspd” to not be defined when the combination of MYJ surface layer, MYNN PBL, and the fractional seaice option are used.

Dynamics

Corrected the tke seed value for km_opt 2 and 5 to be non-zero if and only if the surface heat and momentum fluxes are zero depending on isfflx, diff_opt, and bl_pbl_physics.
Fixed a bug in hybrid vertical coordinate option 3. Those using the default with the hybrid option activated (hyb_opt=2), or those reverting back to a terrain following coordinate (hyb_opt=0) are not impacted.
Vertical nesting now works with the hybrid vertical coordinate.
Vertical nesting now works with the moist theta option use_theta_m=1. The redundant "rebalance" option 2 was removed. The rebalance option must now be set to 1 when vertical refinement is turned on.
Bug fix to use correct rho value in cal_titau subroutines. In the cal_titau subroutines in module_diffusion_em, rho was not always interpolated to the correct location on the staggered grid. This issue affected the subgrid stress terms when sfs_opt .gt. 1. Additionally, when sfs_opt .eq. 0 and m_opt .eq. 1, the actual subgrid stresses applied in the code were correct, but the output was wrong. A rhoavg variable was added to cal_titau_12_21, cal_titau_13_31, and cal_titau_23_32 and used instead of the cell-centered rho variable.
Missing corner points are set when setting boundary conditions of 3D variables and all staggered variables are treated equally in symmetric BC. Thereby the periodic BC becomes cleaner and more consistent with set_physical_bc2d and the symmetric BC are fixed as horizontal deformation is treated now as a staggered variable which leads to correct BC for SGS momentum fluxes at the northern boundaries.

Others

Corrected the incorrect treatment of cloud transmittances of diffuse radiation by clouds in the FARMS system. The impact of these modifications is small.
A fix is introduced to the FARMS parameterization. When the cosine of solar zenith angle is small, the ice cloud optical thickness can be negative, and the model can crash due to log10 of a negative number. The modification is to limit the value of cloud optical thickness.
Check is put in to avoid using Deng shallow convection scheme with icloud_bl (MYNN PBL).
Fixed a restart problem when using 3 or more total domains with feedback turned on, and using RRTMG with o3input=2 (default). Restart now works at all times.
Fixed a restart problem for nested runs when feedback and smooth_option are both on.
The feedback and smoothing are removed for all GWD input fields.
Fix for instances of dt = 0 when using adaptive time stepping together and when step_to_output_time is set to true.
A few updates to the diagnostic schemes were added. The most unstable (MU) layer and precipitable water (PWAT) were inconsistent with the conventional definitions used by the meteorological communities. LFC, CAPE and CIN calculations previously did not work well while dealing with multiple inversion layers.
It is no longer necessary to execute a 'run_me_first.csh' script for any of the idealized cases. All necessary physics input files and look-up tables are linked to the corresponding directory during the compilation step
A bug fix was added to handle the incorrect handling file I/O by the ARW real program when both the following are activated: all_ic_times and multi_bdy_files. The real program effectively skipped outputting some of the processed metgrid time periods by putting the data into an incorrect time. If a user has not seen this failure mode, then this PR does not impact them.
Added missing OMP directives in ARW solve for the option do_avgflx_em=1. Many OpenMP implementations assume the "END PARALLEL DO", but it is more conventional within WRF to have it.
Bug fix for idealized fire case that caused a fatal error when simulations were run with a nested domain.
The netcdf diffwrf program no longer by-passes integer data, but considers these types of gridded fields as eligible for comparison.
The files HLC.TBL, and wind-turbine-1.tbl are now removed from the test/ directories when doing 'clean -a'.
Namelist variable spec_bdy_final_mu is removed. It is no longer used.
The namelist.input file in test/em_real directory has been modified to update several variables to bring case location/date/time to a more recent event (2019 Hurricane Dorian), to include settings for current recommendations, and to be consistent with the updated namelist.wps file.

WRFDA

Added new capability for assimilating surface pm2.5, pm10, O3, CO, NO2, SO2 using 3DVar. (Sun, W., Liu, Z., Chen, D., Zhao, P., and Chen, M., 2020: Development and application of the WRFDA-Chem three-dimensional variational (3DVAR) system: aiming to improve air quality forecasting and diagnose model deficiencies, Atmos. Chem. Phys., 20, 9311-9329.)
Added new capability for Multi-Resolution Incremental 4DVar (MRI-4DVar).
(Liu, Z., J. Ban, J.-S, Hong, and Y.-H. Kuo, 2020: Multi-resolution incremental 4D-Var for WRF: Implementation and application at convective scale, Q. J. R. Meteorol. Soc., 146, 3661-3674.)
WRFDA gen_be_v3 is updated for cv_options=5 and ensemble applications. Please see var/gen_be_v3/README.gen_be_v3.

WRF-Chem

KPP is now compiled unoptimized. This will help with compiling WRF-Chem on ubuntu (and other smaller memory) systems. Timing tests indicate that this does not cause a significant timing impact on the WRF build time (less than 5%).
A fix to ensure NH3, CH4, CO2 fire emissions are correctly read-in to WRF-Chem when using RACM chemical mechanisms. Prior to this fix, values for those parameters were low. Users should now expect significantly higher values when using RACM chemical mechanisms.
The CLM surface scheme and associated subroutines in the physics and chemistry packages have been modified to be made consistent with the MEGANv2.1 biogenic emission model, as opposed to earlier code, which used MEGANv2.0. This upgrade to WRF-CLM parallels the use of MEGANv2.1 in CLM4.0, as described in Guenther et al. (2012, GMD). Within module_sf_clm.F, the emissions for 150 species are computed, and then mapped to the chemistry packages (e.g., SARPC99) in WRF-Chem. Emissions are a function of plant functional type (PFT) and emission factors that depend on temperature, light, and LAI.
Improvements related to convective transport of trace gases. Addition of the trace gases and aerosol subgrid-scale convective transport (Grell and Freitas, 2014; Li et al., 2018), subgrid-scale wet scavenging (including the improvement of ice retention factors and the conversion ratio of cloud water to rainwater, Li et al., 2019), and aqueous chemistry to the GF cumulus parameterization.
pH diagnostics added for mozart_mosaic_4bin_aq chemistry.

WRF-Hydro

WRF-Hydro has been updated to version v5.2.0.

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by kkeene44 almost 4 years ago

WRF-ARW - WRF Version 4.2.2

The WRF model has been updated to Version 4.2.2 on January 15, 2021.

Acknowledgements: We would like to thank Calvin Howes (UCLA), Richard Easter (PNNL), Robert Arthur and Jeff Mirocha (LLNL), Matthias Göbel (University of Innsbruck), Stacy Walters (formerly NCAR), Ryan Cabell (NCAR), Soonyoung Roh and Hwan-Jin Song (National Institute of Meteorological Science, Korea), Arianna Valmassoi (Uni-Bonn), Ted Mansell (NOAA/NuSL), Marc Honnorat (EXWEXs), Lukas Pilz (Heidelberg University), Timothy W. Juliano, Ju-Hye Kim, Pedro A. Jimenez, Jared Lee, Thomas Brummet (NCAR/RAL), Tim Raupach (UNSW CCRC) for their contributions to this release.

This is a bug fix release. Following are fixes associated with this release.

Compiling

Changes are made so that WRF-Hydro can be built successfully with the INTEL (ifort/icc)
New features are added to the WRF code so that the model can be built on a Raspberry Pi with the GNU/8 compiler.
Due to subroutine and function argument inconsistencies, the WRF code will not compile with gfortran/10. The fix introduces a new compiler flag for all gfortran stanzas (-fallow-argument-mismatch -fallow-invalid-boz).

Dynamics

To close the dry theta budget at the lowest grid point when use_theta_m is enabled, the surface heat fluxes, and in the case of km_opt=5 also the non-local heat flux, are multiplied with 1+Rv/Rd*qv. This should be done only for qv at the surface.
Fix incorrect index limits asscoiated with advection. When using advection order < 5, and either open or specified boundary conditions, the upper limit of the indexing in the horizontal advection loops of geopotential was previously incorrect. For horizontal advection order >= 5, the index limits also were incorrect.
Missing corner points are set when setting boundary conditions of 3D variables and all staggered variables are treated equally in symmetric BC. Thereby the periodic BC become cleaner and more consistent with set_physical_bc2d, which leads to correct BC for SGS momentum fluxes at the northern boundaries.
Fix divide-by-zero errors in the km_opt=5 option.
Bug fix in the ideal case initialization routines, in which incorrect indices were used in the calculation of ph_1 after a perturbation bubble was added. For the idealized cases that do not use the hybrid vertical coordinate by default, there is no impact at all.

Physics

Bug fix to use correct rho value in cal_titau subroutines. In the cal_titau subroutines in module_diffusion_em, rho was not always interpolated to the correct location on the staggered grid. A rhoavg variable was added to cal_titau_12_21, cal_titau_13_31, and cal_titau_23_32 and used instead of the cell-centered rho variable.
Bug fix to correctly calcuate surface downward diagnostic output of long- and short-wave fluxes and correct two wrong numerical values used in the code. This bug fix only affects RRTMG-K radiation option.
Fix incorrect parentheses in weighted average for QSFC when sea-ice fraction is used in myjsfc_seaice_wrapper and qnsesfc_seaice_wrapper.
Correct the tke seed value to be non-zero if and only if the surface heat and momentum fluxes are zero.
Bug fix for temporal selection that spans across different years and different local days. This fix affects the irrigation scheme. Previously, the irrigation schemes weren't working when the date selection was across two different years.
A minor fix to avoid occasional instabilities in the surface layer option 1.
Bug fix to account for the effective radius of the unresolved hydrometeors by FARMS. The fix is valid for WSM6, Thompson and all other microphysics options that do not provide effective radius.
Bug fixes for the channel method to ensure that the water used for irrigation in the channel method can be passed to the land surface model.
Several fixes are introduced to the Deng's shallow cumulus parameterization: add no resolved hydrometeors to the unresolved ones to avoid double counting; remove clouds if the updraft is not active; correct conversion factor for precipitation; correct initialization of radiusc variable; limit variable NT to avoid array out of bounds.
Incorporate specific-crop parameters for photosynthesis and stomata subroutine in Noah-MP-Crop. The current release WRF-Crop option 1 uses generic-crop parameters in stomata subroutine to calculate photosynthesis and stomata resistance, regardless of crop type. Also some crop parameters in MPTABLE.TBL for corn and soybean are not consistent with the values used in the publication. This update fixes these parameter values in the crop parameter section in MPTABLE.TBL, as well as adds specific-crop photosynthesis-stomata parameters to this section.

Data Assimilation

Bug fix for WRFDA incorrect application of errfac.dat at each outer loop when use_obs_errfac=true and max_ext_its > 1.
Big fix for WRFDA incorrect ensemble perturbation values for qice, qsnow and qgraup when alpha_hydrometeors=true and ep_format=2.

Chemistry

Fix improper indexing in a conversion factor, which caused addressing (seg fault) or floating point errors when using aircraft emissions for WRF Chem. This has been corrected by matching the indexing on the LHS with the RHS.
Fix a namelist option incompatibility between chem_opt (16) and gas_drydep_opt (1) , which leads to trying to access a variable that is not available in memory. A fatal error was added when the Wesely scheme initialization is called with numgas = 0.

Diagnostics

Fix lh_urb2d diagnostic to be urban area average instead of grid cell average consistent with other urb2d diagnostics associated with NoahMP. No effect on run results.
Fix mixactivate for diagnostic CCN in sectional schemes. This fix impacts Morrison and Lin et al. microphysics schemes.
Correct solar zenith angle calculation and check presence of hydrometeor type before calculating cloud parameters for the solar diagnostic package. In addition, shortwave irradiance variables were added to the tslist when activating the RRTMG or RRTMG FAST shortwave radiation scheme or FARMS.

Registry

Fix the problem of non-reproducibility due to the advection order. Previously, the halo exchange (used during distributed memory communication) had inconsistent use of h_mom_adv_order. Consequently, the model results were not bit-identical with different number of processors (when the advection order differed between the momentum and scalar variables).

Miscellaneous

Fix an inconsistency in the declaration of integer kinds in the routine that processes binary formatted data.
Examples.namelist is modified to provide sufficient information for trajectory run.
Bug fix to make the stand-alone vertical interpolation program to work properly. This utility is only used for internal testing and thus has no impact on users.
An MPI barrier is added at the end of wrf_dm_initialize to force all of the processes to be synchronized before checking the namelist consistency. This barrier avoids occasional crash in setup_physics_suite (in share/module_check_a_mundo.F) when running real.exe.
Remove or fix some obsolescent and non-standard features in several physics modules which were preventing WRF to be built with IBM XL v16.1.1 on Power9.

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by smileMchen over 4 years ago

WRF-ARW - WRF Version 4.2.1

The WRF model has been updated to Version 4.2.1 on July 22, 2020.

Acknowledgements: We would like to thank John Collins (Edge Hill University), James Ruppert (Penn State University), Nicolas Baldeck (OpenMeteoData),Cristina L. Archer (University of Delaware), Sicheng Wu and Yulong Ma (University of Delaware, CReW), and Pedro A. Jimenez (NCAR/RAL), Franciano Puhales (Federal University of Santa Maria, Brazil), Ted Mansell (NSSL), Manuel Luis Aznar (University of La Laguna), Isaac Rowe (University of Kentucky), Matthias Göbel (University of Innsbruck), Matthias Göbel (University of Innsbruck), Yulong Ma (University of Delaware), Christopher Thomas (Climate Change Research Centre, UNSW Australia), Michael Newman (Woodard & Curran), Bart Brashers (Ramboll), Peng Zimu (Peking University) for their contributions to this release.

This is a bug fix release. Following are fixes associated with this release.

Compiling

Changes are made to implement static analysis regression testing of WRF (#1252, #1253))
A compiler error raised by ifort 13 for a fire module routine due to omission of the keyword "SUBROUTINE" in the "END SUBROUTINE" statement has been addressed. (#1223)
A minor modification has been introduced into the SBM Polar Radar microphysics code to allow the WRF code to build with PGI compilers. The nonstandard derf() intrinsic has been replaced by the 2008 standard compliant erf() intrinsic.(#1183, #1246)
The non-standard spacing in macro is removed, which makes the WRF system able to be built with LLVM clang. (#1221)

Dynamics

Fix the 1D advection testing program broken with the hybrid vertical coordinate. (#1247)
For periodic boundary conditions, the slope calculation for the lateral boundaries is made equivalent to the calculation for the interior of the domain (#1217)
Fixed units of surface momentum flux in dyn_em/module_diffusion_em.F when using m_opt. This bug fix divides the surface fluxes by the density to obtain the correct units. This change only affects the output of the SGS fluxes (nba_mij). The tendencies and thus the results are not changed.

Physics

Fix a problem in the advection of TKE generated by the wind turbines if the user instructs MYNN to advect the TKE (setting mynn_tkeadvect = .true. for each domain). There is also an empirical reduction of the wind turbines' TKE coefficient based on LES results (Archer et al. 2020). (#1235). Note that TKE generated by the wind turbines was advected in WRF V3.3 - V3.4. In WRF V3.5, the TKE was advected if mynn_tkeadvect option is on. The error was introduced in WRF V3.6 when the wind turbine scheme was significantly updated.
Fix a divided by zero problem in Shin-Hong PBL (#1219)
Fix an error in aer_opt =2 and Goddard SW radiation to fill in the 2D arrays with the aerosol optical properties (SSA and g) from namelist.input. (#1190)
Fix a floating point log(0) for the NSSL microphysics. Also, minor changes are made to lookup table. (#1216)
Fix a bug in RRTMK radiation shceme. This bug leads to the output shortwave and longwave tendencies being in swapped positions. However, note that this bug does not affect the sum of shortwave and longwave heat tendency, so it does not affect model results (#1257).

Data Assimilation

Fix a problem related to VarBC when num_fgat_time > 1 (#1222)

Chemistry

Enable TUV diagnostics only for debugging (debug_level >=100). This is because when using WRF-Chem with the new Photolysis option (phot_opt = 4) activated, the model spent too much time on looping and writing the TUV.diags which is only used for debugging photolysis rates.

Diagnostics

Fix the segmentation fault problem when time series is output with the solar diagnostic option activated (#1228)
Fix occasionally incorrect calculation of the times of minimum and maximum of fields within the clWRF diagnostics module (#1208)

Registry

The urban module caused segmentation fault when building WRF with real*8, which was traced to a default assignment to the wrong horizontal interpolation and feedback routines. This problem is fixed. (#1213)

Miscellaneous

README.grid_fdda and examples.namelist are modified for correct setting of rinblw when observational nudging is applied to nested domains. (#1248)

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by smileMchen almost 5 years ago

WRF-ARW - WRF Version 4.2

The WRF Model has been updated to Version 4.2 on April 23, 2020.

The WRF Pre-processing System (WPS) has been updated to Version 4.2.

Acknowledgements: We would like to thank Tieh-Yong Koh, Ricardo Fonseca, and Tengfei Zhang (Earth Observatory of Singapore, Nanyang Technological University, Singapore), Chee-Kiat Teo (Temasek Laboratories, Nanyang Technological University, Singapore), Kiran Alapaty and Rob Gilliam (EPA), Robert Arthur, Katie Lundquist (LLNL), Tina Chow and Jason Simon (U.C. Berkeley), Jacob Shpund, Alexander Khain, and Barry Lynn (The Hebrew University, Jerusalem, Israel), Arianna Valmassoi (University College Dublin/NCAR), Jan Mandel and Angel Farguell (University of Colorado Denver), Adam Kochanski (University of Utah), Anders Jensen, Changhai Liu, Greg Thompson, Ju Hye Kim, Pedro A. Jimenez, and Timothy Juliano (NCAR/RAL), Milan Curcic (University of Miami), Negin Sobhani and William Sacks (NCAR/CISL), Mike Iacono (AER), Yago Riveiro (Air Quality and Odor Management), Brian Reen (Army Research Laboratory), Piotr Kasprzyk (CIRI), Robert Rozumalski (COMET), Theodoros Christoudias (Cyprus Institute), Emily Collier (FAU Erlangen-Nürnberg), Han Lang (Fujitsu America, Inc.), Kezhen Chong (Georgia Institute of Technology), Huang Wei (HPE), Thomas Auligné (JCSDA), Alexander Ukhov (KAUST), Sooya Bae (KIAPS), Kyo-Sun Lim (Kyungpook National University), Hugh Morrison, Jamie Bresch, Jim Bresch, Judith Berner, Kevin Manning, Michael Duda, and Yali Wu (NCAR/MMM), Michael Kavulich (NCAR/DTC), Stacy Walters (NCAR/ACOM), James Kenyon (NOAA/CIRES), Joseph Olson, Ravan Ahmadov, and Tanya Smirnova (NOAA), Yu Xie (NREL), Ted Mansell (NSSL), Dongmei Xu (NUIST/NCAR), Shizhang Wang, Xin Zhang, and Feng Gao (NUIST), Daniel Wesloh (Penn State University), Xu Zhang (Shanghai Typhoon Institute/CMA, China), Rob Fovell (SUNY Albany), Ivan Toman (University of Zadar, Croatia), Maik Reichert, Richard Carpenter, and Tao Sun (NUIST).

New in Version 4.2

Physics

Planetary Boundary Layer

A three dimensional scale-adaptive TKE subgrid mixing parameterization scheme (SMS-3DTKE, km_opt=5) developed by Shanghai Meteorological Service (SMS), China. (Zhang, X., Bao, J., Chen, B., and Grell, E., 2018: A Three-Dimensional Scale-Adaptive Turbulent Kinetic Energy Scheme in the WRF-ARW Model. Mon. Wea. Rev., 146(7), 2023-2045, https://doi.org/10.1175/MWR-D-17-0356.1) . The subgrid mixing parameterization extends the original TKE model (km_opt=2, Deardorff 1980) that is usually used as an LES subgrid model to the mesoscale limit in the framework of the WRF-ARW. The scheme can be used in LES, mesoscale and the gray zone resolutions in between. The option must be used with diff_opt=2. The scheme can replace LES subgrid model (km_opt=2) and conventional PBL schemes. When the scheme is turned on (diff_opt=2 and km_opt=5), PBL schemes must be turned off (bl_pbl_physics = 0). This scheme can be used with sf_sfclay_physics = 1, 5, 91. (Thanks to Xu Zhang of SMS, CMA, China)

Surface Irrigation

Three surface irrigation options (Valmassoi et al., 2020, Geosci. Model Dev., to be published) are added to the Noah LSM to parameterize irrigation as a function of three evaporative processes that the water applied undergoes. The irrigation water can evaporate (i) from the surface (Channel method), (ii) from the canopy and then drip on the surface (Drip method), or (iii) from the atmospheric lowest mass-model layer (Sprinkler method) and then follow the previous process. The schemes have an explicit treatment regarding the water amount, location, and timing. (Thanks to Arianna Valmassoi (University College Dublin/NCAR.)

Radiation

A new parameterization, FARMS (Fast All-sky Radiation Model for Solar applications, Xie et al., 2016, Solar Energy), provides fast calculation of the surface shortwave irradiance allowing the WRF model to call the scheme every model time step, resulting in an irradiance that responds to changing clouds. The scheme is activated by setting swint_opt = 2 in the WRF namelist.input file and is added to complete the WRF-Solar suite (Jimenez et al., 2016, Bull. Amer. Meteor. Soc.). (Thanks to Yu Xie of NREL, and P. Jimenez of NCAR/RAL)

Other New Options

A new diagnostic package for users with interests in parameters that are potentially relevant to solar applications is added. In the new package, The following variables are calculated: solar zenith angle, clearness index, 2-D max cloud fraction, water vapor path, liquid water path, ice water path, snow water path, total water path (liquid + ice + snow), liquid cloud effective radius, ice effective radius, snow effective radius, liquid cloud optical thickness, ice optical thickness, snow optical thickness, cloud base height, and cloud top height. Also, for liquid and ice variables, the 'total' water path, effective radius, and optical thickness are computed, where the 'total' variables also account for subgrid (unresolved) hydrometeors. The diagnostics is activated by namelist variable solar_diagnostics (=1) in namelist record @diags. (Thanks to Timothy Juliano, Pedro Jimenez, and Ju-Hye Kim of NCAR/RAL)
Introducing a Windows-based WRF build and run via the free Unix-like CYGWIN environment. (Thanks to Daniel Wesloh of The Pennsylvania State University Department of Meteorology and Atmospheric Science)

Improvements and Bug Fixes

Physics

Microphysics

Updates to the Thompson MP scheme:
A significant alteration is introduced for diagnosing the Y-intercept parameter for the one-moment graupel. These changes affect results for any graupel area. The changes are based on observations of graupel/hail size spectra from aircraft observations. The older diagnosis used an inverse relationship for N0 as a function of graupel mass mixing ratio (Qg) plus a functional dependence on supercooled liquid water (to approximate favorable wet growth). A recent analysis of T-28 aircraft-observed data by Field et al. (2019) shows observations in direct contradiction; particularly that N0 is more directly (not inversely) proportional to Qg. After 5 slight modifications to the constants in the formula, this change produces the best match to observed radar reflectivities, particularly with the upper tail of the distribution.
Updates to the Fast (33 bins) Spectral-bin Microphysics (FSBM):
- The update includes: A switch to use either graupel or hail; Cloud base nucleation; 3 log-normal user-defined aerosol distribution; Updated collision-coalescence; spontaneous rain and snow breakup; and code structure changes.
- A forward polarimetric operator is coupled to the FSBM scheme. The user can see the total reflectivity field, as well as the per hydrometeor total reflectivity (rain, snow, graupel/hail). This can be activated by the namelist option sbm_diagnostics = 1 in &physics.
- In order to run the new FSBM scheme, users need to download an external directory named "SBM_input_33" consisting of mandatory input tables and place it in the 'run' directory. In case the coupled polarimetric forward operator is to be used (e.g., 'sbm_diagnostics = 1'), a second directory consisting of scattering amplitudes named "SBM_scatter_amplit.tgz" needs to be placed in the 'run' directory. Both directories are compressed and can be downloaded at the following link: https://www2.mmm.ucar.edu/wrf/src/wrf_files/
- When using SBM Fast microphysics scheme and simulating deep continental convective systems (as opposed to maritime convection), it is advised to use a 43-bins version of the updated FSBM scheme. In order to run the code using 43 bins, one needs to modify: 1) The registry (to include more bins for the advection), 2) the corresponding loop indices inside the scheme, and the bin-space parameter (NKR) placed at the top of the main module, and 3) switch between the 33-bins input tables with 43-bins.
- The FAST SBM code is not currently configured to work with default 8-byte reals.
- More info about the updated scheme and corresponding results can be found in Shpund, J., Khain, A., Lynn, B., Fan, J., Han, B., Ryzhkov, A., Snyder, J., Dudhia, J. and Gill, D., 2019. Simulating a Mesoscale Convective System Using WRF with a New Spectral Bin Microphysics: 1: Hail vs Graupel. Journal of Geophysical Research: Atmospheres. https://doi.org/10.1029/2019JD030576.
  The 43 lookup tables can be downloaded here:
  https://www2.mmm.ucar.edu/wrf/src/wrf_files/

Cumulus

Updates to Multi-scale Kain-Fritcsh cumulus parameterization scheme:
*Added Zhang and McFarlane (JGR, 100, 1995) convective momentum transport (CMT) to the MSKF scheme (improves speed of propagation of storms and squall lines and associated changes in horizontal winds components). The option is controlled by cmt_opt_flag inside module_cu_mskf.F, and turned on.
*The MSKF cumulus parameterization scheme is now working with other PBL schemes, in addition to YSU. But the newly added 3DTKE scheme cannot yet work with MSKF.
A bug fix for the BMJ cumulus scheme corrects an error with conservation of enthalpy. The fixed code produces less surface rain. More detail can be found in the following paper: https://www.geosci-model-dev.net/8/2915/2015/
A subgrid-scale cloud-radiation feedback scheme for Betts-Miller-Janjic (BMJ) cumulus is added and activated by namelist option bmj_rad_feedback option in &physics. The scheme is based on precipitation rate. See the following paper for science details: Koh, T.-Y. and R. Fonseca (2016), "Subgrid-scale cloud-radiation feedback for the Betts-Miller-Janjic convection scheme", Quarterly Journal of the Royal Meteorological Society, 142(695), 989-1006. doi: 10.1002/qj.2702. For tropical applications, the authors recommend this scheme to be used with the WDM5, WDM6 or WDM7 cloud microphysics schemes, and with the RRTMG / RRTMG fast radiation scheme. The BMJ scheme should also be modified according to Fonseca et al. (2015) by reducing the "Fs" parameter from its default value of 0.85 to 0.6 so that it produces a good estimate of the observed tropical precipitation on pentad time-scales as given by the TRMM satellite data. (Thanks to R. Fonseca of Earth Observatory of Singapore, Nanyang Technological University, Singapore)
One of the NCEP SAS cumulus schemes, previously defined as cu_physics=4, has been changed to cu_physics=95 many versions ago. Now if an ARW user requests cu_physics=4 (which doesn’t work for ARW), the code gracefully stops and suggests using cu_physics=95 instead.
A check is added if Deng shallow scheme (shcu_physics=5) and Thompson cloud fraction icloud=3 are used together. Since Deng shallow scheme provides subgrid cloud, it should not be used with Thompson cloud fraction. The PBL check for use with Deng shallow scheme is modified: Deng scheme may be used with MYNN, but EDMF and icloud_bl have to be off. But making these options off will change the behavior of the MYNN PBL option. Hence this is not recommended.

Radiation

Updates to RRTMG radiation:
Code fixes for the single-scattering albedo and asymmetry parameter inputs to the “denom” equation. This is a modified fix to the one provided in 4.1.3.
RRTMG (and fast version) longwave and shortwave CO2 value is now a function of year, which replaces the outdated constant value of 379 ppm (applies to 2005).
The vertical refinement option now works with the RRTMG radiation scheme.

PBL

Updates to the MYNN-EDMF (Thanks to Joe Olson and Jaymes Kenyon of NOAA/GSL and Wayne Angevine NOAA/PSL):
- Modifications to better balance the new heating terms (subsidence and TKE dissipation) with strength of detrainment. Helps reduce daytime PBL cold bias in RAP/HRRR retrospective tests.
- Changes which impact the resolved and subgrid-scale (SGS) clouds: addition of QI_BL instead of both qc and qi lumped into a single array. This helps reduce the excessive SGS clouds in the polar region.
- Bug fix to for diagnostic-decay of SGS clouds
- Top-down diffusion linked to cloud-top radiative cooling is now active by default.
- Changed the scale-aware tapering functions to allow for more mixing/transport at dx < 1 km, since the MYNN-EDMF seems to naturally taper itself off below dx = 1km.
- Added namelist control, using bl_mynn_output, for allocation/output of 10 3D diagnostic output, and removed TKE_PBL in output (only QKE).
- This code is similar to what's in the upcoming HRRRv4, which will become operational in June 2020.

LSM

The default VEGPARM.TBL file has been modified for the RUC LSM scheme. Several roughness length categories were increased for use with both the USGS and the MODIS data sets. Similarly, the STAS-RUC section of the SOILPARM.TBL file was modified to be consistent with the RUC developer's requests. The impact is only to RUC LSM.
The U10/V10 diagnostic variables output from the MYJ and QNSE surface layer schemes have been replaced. The standard fields from MO theory are now output as U10/V10, and the original MYJ SFC scheme diagnostic values of U10/V10 are output with the names U10E/V10E. For those using the MYJ surface layer scheme, this change will modify the resultant U10/V10 values that are used inside of the diagnostics and the time series options.

Urban

Fixed a restart problem with all three urban options by adding a few initial urban fields (building height, building fraction, etc.) to restart file. This is a problem if those urban initial fields are non zero. A user could get these fields from either the NUDAPT44 data directly or have URB_PARAM in the GEOGRID.TBL which uses a pre-processed dataset. Since NUDAPT data is US only, this fix only affects US domain applications.

Surface Layer

Updates to MYNN surface layer: Uses z/L calculation from Li et al. (2010 BLM) for the first time step. After that uses the method similar to that in module_sf_sfclayrev.F. Uses the blended Monin-Obukhov flux-profile relationships COARE (Fairall et al 2003) for the unstable regime. Tweaks to Zilitinkevich (1995) thermal roughness length; Czil now set to 0.085. Removed the Zilitinkevich thermal roughness lengths for water. Removed several 3D variables (qc_bl,cldfra_bl,tsq,qsq,cov,Sh3D,el_pbl). Updated MYNN sea-ice code. This MYNN surface layer scheme is very similar to what's in the upcoming HRRRv4. This new version of the MYNN surface layer scheme will become operational in June 2020. (Thanks to Joe Olson of NOAA/GSD)

Dynamics

Fix based on theory for horizontal diffusion of w in diff_opt=2 (Thanks to Lundquist of LLNL and Chow of U.C. Berkeley).
Defaults for damp_opt have been changed from 0 to 3, or w-Rayleigh damping, and its coefficient dampcoef to 0.2. (Ref: Klemp, J. B., J. Dudhia, and A. D. Hassiotis, 2008, MWR)

Others

All of the odd nest ratios (excluding 3:1) had an incorrect feedback. The indexing was off by (ratio/2 -1). For example, a 5:1 ratio would be off by a single fine grid point. Even ratio feedback was not impacted.
The gwd_opt namelist has been changed to domain-dependent namelist. Hence one may turn off the option in a high-resolution domain.
The ARW real and WRF code now ignore the provided grid distance for any domain other than the most-coarse grid (using the parent grid distance and the grid distance ratio to compute the child grid distance resolution). This update to the namelist processing for grid distance provides user-level behavior in the real/WRF namelist.input file similar to the WPS geogrid and metgrid programs’ usage of the namelist.wps file.
A run-time capability for multiple LBC files is now available. This supersedes and replaces the compile-time option. This is a run-time option accessed through &bdy_control namelist logical variable multi_bdy_files. This option requires that the lateral boundary file names in &time_control include a date: bdy_inname = "wrfbdy_d<domain>_<date>", and each of the wrfbdy file must contain boundary conditions for one time period. All other functionality of the capability is identical to the previously (and no longer) existing compile-time option.
Added vertical motion (suffix = WW) and pressure (suffix = PR) to time series profile output. The simulation start date is now in the title or all time series output. The geopotential height in the time series output is now at half levels (similar to as is done for the vertical velocity WW).
The default setting for the WRF model is to NOT compile either of the two derivative RRTMG schemes: FAST or KIAPS. These codes break the GNU compiler (v6.3.0 through v8.*). To re-enable the schemes, modify the configure.wrf file after running the ./configure command. Replace the "zero" values with "ones" in the lines -DBUILD_RRTMG_FAST=0 or -DBUILD_RRTMK=0. This is a compile-time modification.
nwp_diagnostics are now compatible with adaptive time steps.
By changing the names of the older ESMF modules in the WRF source code (for example, from ESMF_AlarmMod to WRF_ESMF_AlarmMod), the WRF model may now be built with released ESMF library versions.
The WRF model is now able to use > 9999 MPI processes. Users will notice that the file names are modified for jobs requesting > 9999 MPI processes. For example, the command mpirun -np 10001 wrf.exe will produce file names rsl.out.00000000 through rsl.out.00010000. For convenience, to only have output go to the master rsl files, a run-time commented-out flag is provided in the configure.wrf file, as part of the CFLAGS_LOCAL makefile macro.
Most of the WRF modeling system uses the netcdf.inc file instead of use netcdf, and most of the WRF modeling system uses mpif.h instead of use mpi. The Fortran module constructs are removed for consistency with the rest of the WRF system and also to lower the bar for porting the WRF code to small Linux boxes. With this modification, the WRF code now builds and runs with newer GNU compilers (starting to become important with Fortran 2003 and Fortran 2008 constructs used inside the model) but still using the default installed libraries on the small desktop Linux machines (where the libs are built with GNU 4.8-ish). In addition to assisting the user base that has access to user-managed Linux systems, these mods simplify building containers for WRF which can be used for automated testing or for cloud based purposes. Lastly, this modification allows the Windows-based CYGWIN environment to use a fairly recent version of the GNU compiler that was not necessarily used to build the automatically installed netcdf or mpi libraries.
A fix was implemented for observation nudging when nudging temperature observations on the boundary edge of a tile that used an incorrect nudging term.

WRFDA

New Features

New capability for directly assimilating radar reflectivity using a new observation operator and its TL/AD operator considering snow and graupel.
Wang, S. and Liu, Z.: A radar reflectivity operator with ice-phase hydrometeors for variational data assimilation (version 1.0) and its evaluation with real radar data, Geosci. Model Dev., 2019, 12, 4031–4051.
New capability for variational bias correction of TAMDAR T observations.
Gao Feng, Liu Zhiquan, Ma Juhui, A. Jacobs Neil, P. Childs Peter, & Wang, Hongli. (2019). Variational Bias Correction of TAMDAR Temperature Observations in the WRF Data Assimilation System. Monthly Weather Review. 147. 10.1175/MWR-D-18-0025.1.
New DA ep_format option to read in ensemble perturbation (ep) generated by existing gen_be_ep2 and new gen_be_v3 utilities.
New stand-alone gen_be_v3 for generating ep (ensemble perturbation) and be (background error) for WRFDA alphacv_method=2, alpha_hydrometeors=true, cv_options=7 and cloud_cv_options=2 applications. The major advantage of gen_be_v3 is that it is much more efficient than the existing gen_be package.
A new channel-based cloud detection scheme for Infrared sensors is added based on the Particle Filter method.
Xu D., T. Auligné, G. Descombes, and C. Snyder, 2016: A method for retrieving clouds with satellite infrared radiances using the particle filter. Geosci. Model Dev., 9, 3919–3932.

Updated Features

WRFDA analysis_type="RANDOMCV" is improved. To get multiple perturbed output in one WRFDA run, set a new namelist n_randomcv to the desired number.
For EnVar DA applications, logical alpha_vertloc namelist switch is replaced by new integer option alpha_vertloc_opt. See var/README.namelist.
For additional details regarding WRFDA release information, see https://www2.mmm.ucar.edu/wrf/users/wrfda/updates-4.2.html

##WRF-Chem

"ozone" is added to the list of profile fields for time series, but only available when building with WRF Chem.
Prevent 0 index in GOCART settling array.
Removing duplicate lightning option entries from WRF-Chem. They are set under &physics.
Fixing extra biomass burning mass aloft.

WRF-Hydro

WRF-Hydro v5.1.2 includes the following new features and fixes: Water class value for output masking is now set to the appropriate water class code instead of being hard-coded to 16. WRF-Hydro can now be linked to versions of netCDF newer than v4.6 without crashing. Overland routing has been refactored into a new object-oriented Fortran module. The soil moisture smcmax check threshold was reduced from 10E-6 to 10E-5 due to some conflicts still triggering at the lower threshold. Surface runoff that was mapped through groundwater (GW) buckets in non-UDMP configurations now behaves as pass-through only. Bug fix to snow liquid water movement. NetCDF 3 outputs were converted to netCDF 4 and the compile-time option for large file support was removed. Bugfix to reservoir physics: the uniform lake surface area was set to be consistent with the vertical wall assumption. Makefiles no longer generate intermediate preprocessed source files. Please see https://github.com/NCAR/wrf_hydro_nwm_public/releases for more detailed release notes.

WRF-Fire

Fuel moisture content (FMC) model is coupled with the coupled fire-atmosphere model. FMC has a strong effect on fire spread, including the diurnal cycle, and previously was a constant in the fuel description. Illustration of the effect of the coupled FMC model on diurnal fire spread variability in the 2012 Barker Canyon Fire, from the following: https://user-images.githubusercontent.com/6204560/53294946-ed342c00-37ad-11e9-8d0b-d50d268dc027.png Reference: J. Mandel, S. Amram, J.D. Beezley, G. Kelman, A.K. Kochanski, V.Y. Kondratenko, B.H. Lynn, B. Regev, M. Vejmelka, Recent advances and applications of WRF-SFIRE. Natural Hazards and Earth System Science, 14, 829-2845, 2014, https://doi.org/10.5194/nhess-14-2829-2014 (Thanks to Jan Mandel and Angel Farguell (University of Colorado Denver), Adam Kochanski (University of Utah))

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by kkeene44 about 5 years ago

WRF-ARW - WRF Version 4.1.5

The WRF model has been updated to Version 4.1.5 on March 10, 2020.

This is a bug fix release. Following are fixes associated with this release.

Physics

Problem:
Newer versions of the GNU compiler more tightly conform to the Fortran standard and do not immediately short circuit IF test processing once a statement is known to be FALSE. This newer and legal behavior recently caused troubles with NoahMP regarding the frequency to process ground water.

When running the NoahMP scheme with the ground water option (opt_run=5), the STEPWTD variable is initialized to a non-zero value. However when opt_run!=5, then the STEPWTD variable remains an uninitialized zero value. In the WRF source code, there are a couple of places where mod(x,STEPWTD) occurs. This causes an immediate "divide by zero" error in the first call to the surface driver.

Solution:
The IF test logic is re-arranged so that only when opt_run=5 (which means that STEPWTD is non-zero) is there, a second IF test (using mod(x,STEPWTD)) to compute whether additional processing takes place.

The purpose of this modification is to allow users access to the latest v4.1.* branch that works with NoahMP. This release fixes that trouble with modified logic, but the still existing problem could extend (depending on the compiler) back to version 3.6.

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by smileMchen about 5 years ago

WRF-ARW - WRF Version 4.1.4

The WRF model has been updated to Version 4.1.4 on February 12, 2020.

Acknowledgements: We would like to thank Yali Wu (NCAR MMM), Kyo-Sun Lim (Kyungpook National University), Sooya Bae (KIAPS) and Yago Riveiro (Air Quality and Odor Management - AQOM) for their contributions to this release.

This is a bug fix release. Following are fixes associated with this release.

Compiling

The replicated symbol "SMOOTHER" in subroutine feedback_domain_nmm_part2 is renamed to be "SMOOTHR" to get HWRF to build and run the regresion test cases
An ifdef is put around the routine (module_dm.F) that handles pops and pushes of the the communicator stack to make sure the communicators are only required for distributed memory MPI jobs.
The incorrect use of whereis in configure to locate nc-config is fixed.

Physics

Changes in the snow/graupel/rain mass and rain number concentration caused by melting of snow/graupel are made to occur at the same time for WDM5, WDM6 and WDM7 microphysics. Also, the cloud water number concentration is added to the CCN number concentration. These changes are more physically reasonable and would leads to decreases in the generation of rain number concentrations and increases in CCN number concentration .
Fix a problem in Deng shallow scheme when it is combined with other sub-grid cloud schemes (e.g. non-microphysics options that have radiation feedback).
Fix an occasional divided-by-zero error in SFCLAYREV scheme

Data Assimilation

Fix a problem in the calculation of effective radius of rain, snow and graupel.

Computing

The maximum number of eta levels is increased from 1001 to 10001 to avoid REAL porgram failure when the number of vertical levels is more than a thousand.

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by smileMchen about 5 years ago

WRF-ARW - WRF Version 4.1.3

The WRF model has been updated to Version 4.1.3 on November 25, 2019.

Acknowledgements: We would like to thank Alexander Ukhov (KAUST) , Brian Reen (Army Research Lab), Pedro Jimenez and Greg Thompson of NCAR/RAL, Ted Mansell of NSSL, Alexander Ukhov, KAUST, Gerardo Cisneros-Stoianowski (Mellanox Technologies), Kevin Manning and Jamie Bresch of NCAR/MMM for their contributions to this release.

This is a bug fix release. Following are fixes associated with this release.

Compiling

A modification was added to the compile script to print out the correct version number of the WRF model in the compile log.

Physics

Deng shallow scheme fixed for occasional blow-ups due to uninitialized RADIUSC or a divide by zero.
Correction to NSSL 2-moment scheme. Snow reflectivities were too low and were reverted to a previous version. Effective radius calculations for droplets/ice (used only by RRTMG) were only calculated at history output times and are now done every time step.
In 3.5.1, a bug fix was introduced to fix the possibility of dividing by zero in subroutine reftra_sw, in module_ra_rrtmg_sw.F. The fix was to limit the denominator to a small positive number. However, this denominator can be both positive and negative. That V3.5.1 fix was not correct, and introduced cold bias in areas of cloud. This PR fixes the nearly identical problem in all three module_ra_rrtmg_sw*.F files.
Corrected a problem that caused model to stop when using sf_ocean_physics = 1, with WRF code compiled for bounds checking.
Using the YSU PBL scheme in combination with either the WSM3 or Kessler MP schemes (no ice) is now fixed to not use ice array.
When using the MYNN PBL scheme, with icloud_bl = 1 (which is default), restarts did not give bit-for-bit results when compared to a non-restart run. This has been a problem since the option was introduced in V3.8, but is now corrected.
For Thompson MP, relatively small bug fixes and safety checks have been implemented. Specifically: avoid exceeding max amounts of snow, graupel, or cloud water; adjust snow terminal fall speed; initialize snow moment variables; permit larger max values of ice crystal concentrations.
Corrected vertical index error in fire module to avoid accessing values of a half-level variable with full-level index.
Limited max_ts_level to be less than the number of half levels.
namelist checks are now conducted only on used or allowed domains.

Nudging

Fixed a problem with observation nudging diagnostic prints. Sometimes the information from multiple observations were mixed together, and sometimes missing observations were printed out.
Fixed conflict in file units used in obs nudging and Thompson microphysics.

Idealized Cases

If a user chooses to explicitly list the eta levels for an idealized case, a check has been added to verify that valid eta levels are supplied for idealized cases in the namelist.input file.

Data Assimilation

A bug fix in the trajectory calculation of specific humidity for CRTM may slightly improve radiance DA results that use water-vapor-sensitive channels. There is not expected to be a large impact.
Fix WRFDA implementation of alpha_std_dev for tuning flow-dependent B. Current namelist alpha_std_dev cannot be set to values other than 1.0 (default). After the fix, setting alpha_std_dev to a value larger than 1.0 inflates the spread, and setting alpha_std_dev to a value smaller than 1.0 deflates the spread.

Chemistry

Corrected a typo in diagnostic dust gravitational settling flux.

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by kkeene44 over 5 years ago

WRF-ARW - WRF Version 4.1.2

WRF Model Version 4.1.2

The WRF model has been updated to Version 4.1.2 on July 12, 2019.

Acknowledgements: We would like to thank Changhai Liu (NCAR), Nicolas Baldeck (OpenMeteoData), Emily Collier (FAU Erlangen-Nürnberg), Han Lung (Fujitsu America, Inc.), Jamie Bresch (NCAR) and Anders A. Jensen (NCAR-RAL) for their contributions to this release

This is a bug fix release. Following are fixes associated with this release.

Compiling

A modification was added to the compile script to print out the correct version number of the WRF model in the compile log.

Diagnostics

Bug fix in temperature and winds computation when the “traditional fields” diagnostics option is turned on.
- For temperature, the moisture is incorrectly removed from the already dry potential temperature perturbation, which leads to a cooler temperature that can be 8-10 K lower in locations near the surface with high moisture content.
- For winds, when computing the earth-relative winds from the model’s projection-relative winds, both zonal and meridional winds were too large by a factor of two.

Physics

Logic was added to include the buckets for accumulated radiation fields when the user selects bucket_J > 0. A bug was introduced since WRFV3.9, when the user selected the bucket_J >0, no additional fields were computed or output.
Bug fix to ensure restart runs are reproducible in nested domains when topographic shading is turned on.
Bug fixes in Jensen ISHMAEL microphysics
- Add allocation for single-precision lookup tables which are defined from the double-precision ones. Without this, certain versions of PGI throw an error.
- Move the calculation of the inherent growth ratio to prevent it from being undefined in certain instances.
- Store shape parameter, use it to diagnose shape during nucleation and aggregation (transfer between ice categories), and ensure that this value stays in bounds. This prevents the updated shape from going out of bounds in certain instances.
- Move the mass check for aggregation to occur to around the aggregation subroutine call to prevent the aggregate diagnostics from being undefined.
- Add a fix to ensure that nucleation won't turn planar particles into columnar ones by limiting this shape to spherical.
- Set a low limit on nucleation size to 2 microns.

Data Assimilation:

Fix to allow proper QVAPOR analysis when mp_physics is not set in namelist.input.
Fix WRFDA V4.1.1 serial compilation failure.
Fix WRFDA 4DVAR V4.1.1 compilation failure.
Fix Radar DA memory leak. Some allocated arrays for radar were not deallocated.

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by smileMchen almost 6 years ago

WRF-ARW - WRF Version 4.1.1

WRF Model Version 4.1.1

The WRF model has been updated to Version 4.1.1 on June 4 , 2019.

Acknowledgements: We would like to thank Katie Lundquist (LLNL), Maik Reichert,
RCarpenter and Robert Rozumalski (COMET) for their contributions to this release

This is a bug fix release. Following are fixes associated with this release.

Compiling

Following the correct syntax of using CPP #includes, angle brackets are replaced by quotes (frame/module_alloc_space.h)
The “CHUNK” variable in the Goddard MP and radiation schemes are changed from “CHUNK=16” to “chunk=16” to avoid confusion caused by the arch/configure.default -DCHUNK compiler flag.
Enforce conditional use of IEEE Fortran 2003 capabilities in v4.1 Goddard schemes, which avoids the failure of WRF compilation using GNU on Linux.

Vertical refinement

Bug fix for the vertical refinement options. From version 4.0 through 4.1, both vertical refinement options were broken due to the introduction of hybrid vertical coordinate option. For online vertical refinement (within the WRF model using namelist vertical_refine_method) with only a single input domain, the incorrect WRF code had uninitialized arrays that defined the vertical coordinate in the child domain. However, any users with vertical refinement activated AND with input_from_file = t,t were never impacted. For offline vertical refinement (using ndown with vert_refine_fact), the initial conditions also only had the lower portion of the 1d arrays initialized (and those were initialized incorrectly). Modifications were introduced to fix both refinement options, and the modifications now support the hybrid and the terrain-following vertical coordinate.

Dynamics

Bug fix for the option diff_6th_opt. A logical variable is not specified before it is used, which can lead to model failure when running with this option on some machines and using some compilers.

Physics

Force positive definite with small mass and number from RRTMK for cloud, snow and graupel in the computation of effective radius to avoid floating exception

Miscellaneous

Fix misleading error message when an auxiliary file was missing
Correct the typo in the Ishmael DFI package
Correct some instructions in README.namelist
Add missing commas to LANDUSE.TBL
Fix compiler-sensitive looping for random seeds initialization in stochastic processing
Add controlled stop when timing problems with the boundary data are detected

Data Assimilation

Bug fix for missing values in bufr files. This fix prevent the errors from occurring in the debug build
Bug fix for dependencies in DA dm_bcast_interfaces. All known missing use statements related to wrf_dm_bcast_* calls are corrected.
Avoid DA EnVAR unnecessary ep allocation to reduce memory requirement
Fix DA dual-res hybrid significantly different results with different number of processors.

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by smileMchen almost 6 years ago

WRF-ARW - WRF Version 4.1

The WRF Model has been updated to Version 4.1 on April 12, 2019.

The WRF Pre-Processing System (WPS) has been updated to Version 4.1.

Acknowledgements: We would like to thank A.-J. Deng (Penn State U.), John Henderson and Mike Iacono (AER Inc.), Patrick Hawbecker (National Renewable Energy Laboratory), HWRF team at (NCEP-EMC), Xiaoyu Xu (Inst. Of Urban Meteorology, China), Patrick Campbell, Jonathan Pleim, Limei Ran, and Robert Gilliam (US Environmental Protection Agency), Toshi Matsui and Carlos Cruz (NASA-GSFC), Jeremiah Johnson and Bart Brashers (Ramboll), Yuanbing Wang (Nanjing U. of Information Science & Technology, China), Joe Olson and Tanya Smirnova (NOAA-ESRL-GSD), Kevin Thomas (U. of Oklahoma-CAPS), Siou-Ying Jiang (Central Weather Bureau, Taiwan), Hamada S. Badr (Johns Hopkins U.), Brett Wilt of IBM, Sooya Bae and Songyou Hong (KIAPS, S. Korea.), and Anders Jensen, Michael Barlage, Greg Thompson, Mei Xu and Jamie Bresch (NCAR) for their contributions to this release.

New in Version 4.1:

Physics

Microphysics

Jensen Ice-Spheroids Habit Model with Aspect-ratio Evolution (ISHMAEL) microphysics (mp_physics = 55). (Jensen, A.A., J.Y. Harrington, H. Morrison, and J.A. Milbrandt, 2017: Predicting Ice Shape Evolution in a Bulk Microphysics Model. J. Atmos. Sci., 74, 2081–2104,https://doi.org/10.1175/JAS-D-16-0350.1)
The parameterization predicts ice particle shape evolution for planar-nucleated (QICE) and columnar-nucleated (QICE2) ice. A third ice species, aggregates (QICE3), is also prognosed. Ice shape can change rapidly when ice grows by vapor deposition or riming, and thus, short microphysical timesteps are preferred. It is recommended that the option is used with 3 km grid spacing or less, and 6 sec timestep or less, but the scheme can be run at larger grid sizes and longer time steps. It’s also recommended to use positive-definite advection scheme (i.e. moist_adv_opt=2 and scalar_adv_opt = 2) with this scheme. The scheme is coupled to RRTMG radiation (ra_sw_physics=4, ra_lw_physics=4). The scheme is not officially coupled to cumulus parameterizations, but QICE is used for ice initialized as planar. The other two ice species (QICE2, QICE3) currently do not interact with cumulus parameterizations.
WRF Double Moment 7-Class Scheme (WDM7, mp_physics=26): This is a new double-moment microphysics with hail as a new prognostic variable. The added hail category helps to reduce the widespread light precipitation (by reducing the formation of graupel) and enhance the heavy precipitation. Both aspects are viewed as improvement over the WDM6 scheme.
WRF Single Moment 7-Class Scheme (WSM7, mp_physics=24): This is a new microphysics with hail as a new prognostic variable. The added hail category helps to reduce the light precipitation (by reducing the amount of graupel) in the stratiform region, and enhance the heavy precipitation in the leading edge of a simulated 3D idealized squall line as well as a real-data case. Both aspects are viewed as improvement over the WSM6 scheme.
Goddard 4-Ice microphysics scheme (mp_physics = 7) (Tao et al. 2016): This scheme replaces the old Goddard microphysics scheme. It has both graupel and hail. It outputs effective radii for microphysical species to interact with both the updated Goddard long- and short-wave radiation and RRTMG radiation. The 4ICE scheme was built upon improved versions of the Goddard 3ICE scheme developed for the Goddard Cumulus Ensemble model and includes dozens of new/modified functions regarding ice microphysics and unique particle-size mapping and observation-driven snow density.

Cumulus

Deng shallow convection option (shcu_physics=5) (Deng et al. 2003): It’s a mass-flux scheme. It considers both buoyant updraft and cloud with nearly neutral buoyancy, and predicts cloud fraction and cloud liquid content which then interact with radiation physics. Even though it is named a ‘shallow convection scheme’, it does have a deep convection component. In conditionally unstable environment, the scheme transitions to a version of Kain-Fritsch scheme. However, as the code is implemented in WRF now, this deep component is not very active, but it has been tested with a deep scheme on too with reasonable results.

Radiation

Goddard Radiation (ra_lw/sw_physics=5) (Matsui et al. 2018): The newest version of the Goddard radiation scheme involves three major improvements. First, a size-, shape-, and radiation-spectrum-consistent single scattering database (Yang et al., 2013) has been incorporated to represent hydrometeor-consistent optical properties (including rain, snow, graupel, and hail). Second, the molecular absorption database has been updated from HITRAN1996 to HITRAN2012, which reduces the biases for clear-sky radiation flux. Third, the radiation code has been vectorized for improved computational efficiency. This scheme is best used with Goddard 4ICE microphysics scheme.

Improvements and Bug Fixes

Physics

Land Surface Model

Noah-MP:
- Updated with code modifications that are present in WRF-Hydro/Noah-MP. Non-answer changing: parameters data structure expanded, read/write statement formatting; Answer changing: remove normalized LAI, snowpack heat conservation, increase maximum SWE limit.
- When running Noah-MP with no urban scheme (bulk method), the roughness length in urban areas used a bare soil value. This resulted in a high temperature, high wind speed, and low sensible heat flux over cities. The correction uses z0 from the Noah-MP look-up table.
PX LSM: The new option (pxlsm_modis_veg = 1) was added in the physics section of the namelist file. When activated, this option uses the time-varying VEGFRA and LAI from the wrflowinp_d01 file instead of the look-up values in the PX data table. Users are encouraged to use this option if they are using WPS v4.0 and later that have the higher resolution MODIS greenfrac dataset. Also, the soil calculation in the PX LSM were modified to use analytical functions from Noilhan and Mahfouf (1996) for field capacity, saturation and wilting point based on fractional soil data. Soil categories were updated in PX to 16 class consistent with the WRF system and other LSMs.
RUC LSM: flag_sm_adj is added to turn on/off soil moisture adjustment for RUC LSM. This flag should be turned on when RUC soil moisture is initialized from the Noah soil moisture (such as data from GFS, NAM), and turned off if it is initialized from other than Noah LSM. The default value is 0 => do not do soil moisture adjustment.

Cumulus

MSKF scheme: Set a hard limit for the convective time-scale to not go beyond 24 hours in the scheme. This helps with some occasional model blowups.

PBL

MYNN Updates: The biggest improvement is the reduction in the downward shortwave radiation bias through better cloud fraction and subgrid scale mixing ratios. Improved ensemble spread from changes to SPP in MYNN. Added a new option to mix 2nd moments and aerosols in MYNN as opposed to the scalar_pblmix option. This is activated with bl_mynn_mixscalars = 1; automatically sets scalar_pblmix = 0. Added the non-Gaussian buoyancy flux functions of Bechtold and Siebesma (1998, JAS) to improve the turbulent mixing in cloudy environments (only activated when bl_mynn_cloudpdf = 2). Added heating due to dissipation of TKE (activated by default - small impact). Important bug fix for wrf chem when transporting chemical species in the MYNN mass-flux scheme. Removed a lot of unused code, i.e., 2nd mass-flux scheme (now only bl_mynn_edmf = 1, no longer an option 2).
YSU, and Shin-Hong PBL schemes: fix the constant h1 in its 8th decimal place from 0.33333335 to 0.33333333. YSU now outputs exchange coefficient for momentum, exch_m.

Urban

BEM: Currently, the option applies air conditioning to all floors of all buildings on a model grid. Two additional urban parameter table options are added to prescribe the fraction of building (BLDAC_FRC) and fraction of floors in a building (COOLED_FRC) that have air conditioning. This reduced the air conditioning energy load and compares better with observations.
BEP and BEM: These multi-layer urban models pack 4D, 5D, and 6D arrays into 3D. The original code used one dimension (num_urban_layers) which depends on the largest 6D array. With no answer changes, the new code creates separate mapping for the different urban arrays so that the arrays are maximally filled for each urban array. Compared to the original code, the memory cost (relative to a non-urban run) of using UCM is reduced from 31% to 5%, BEP from 31% to 15% and BEM from 707% to 274%, equivalent to a 20%, 45% and 64% reduction in memory used, respectively. When users run BEP or BEM, new namelist values for the urban dimensions must be set (num_urban_ndm, num_urban_ng, num_urban_nwr, num_urban_ngb, num_urban_nf, num_urban_nz, num_urban_nbui) and equal the values in the BEP and BEM modules (see run/README.namelist).

Radiation

An error in CAM ozone time interpolation exists in WRF prior to version V4.1. The ozone interpolation is off by a month; zero value (which is wrong) and value valid at Jan 16 are used for interpolating dates from Dec 15 - Jan 16. This fix affects RRTMG radiation with o3input=2.

Microphysics

Thompson MP schemes (mp_physics= 8, 28): A few adjustments made for rain number concentration tendencies, two of which are bug fixes. Also a minor adjustment for consistency to minimum cloud ice size in event a balance of ice number is needed. Lastly a bug fix for partly-melted snow. The most notable changes are: (1) a larger mean rain size from cloud water to rain autoconversion that helps to produce a stronger leading edge of convective squall lines; (2) partly-melted snow was falling too rapidly while this fix produces very realistic progression of melting snow falling gradually toward rain velocity as it melts.
Thompson MP schemes (mp_physics = 8, 28): Most GRIB files contain mass mixing ratios for cloud/precip variables, but not all of them contain number concentrations. A good example is pressure level RAP/HRRR files whereas the native level model files contain both mass and number. Artificially initializing a model without the associated number concentration fields will cause greatly inferior assumptions at run time versus fixing the number concentration variables with something sensible during real.exe to get both initial and boundary condition 2-moment variables reasonably in sync. The solution here is creating number of cloud droplets, cloud ice, and rain for any of those that have mass but not number.
P3 microphysics has had several updates and a few bug-fixes since v4.0 including an important one for providing snow/ice as part of RAINNC. This brings WRF up to the P3 v3.1.11.
WSM6 and WDM6 schemes have been fixed to allow more freezing rain.

Cloud Fraction

Cloud fraction scheme (icloud=3): Compared to prior version generally produces lower cloud fraction amount at high altitude (cold tropopause level temperatures for example). Full year (every 5th day) test gives results with better shortwave radiation at ground comparison to observations.

Others

The time-series capability has been upgraded.
The output option has been modified to allow output at i,j coordinate locations. One must specify i/j coordinates in the tslist and change the column headers to i/j (instead of lat/lon).
Both real and ideal cases are capable of outputting all variables (including U/V) at the cell centers by setting the namelist option tslist_unstagger_winds = .true.
Time series output now includes vertical velocity and pressure.
When nwp_diagnostics = 1, the max radar reflectivity output is not the maximum value between history output, because it is only computed once at the history output time. This change fixes it and will also turn on radar reflectivity calculation (do_radar_ref) if it isn't turned on. Note that this change has increased computational time significantly due to reflectivity being computed at every model time step.
In V3.8, we changed default surface_input_source option value from 1 to 3, which uses the dominant categories computed in geogrid. But we neglected that when dominant categories for land and soil are recomputed in real, the real program does checks for mismatches, making sure the lower resolution isltyp matches with landmask data. This has two effects: one is the actual isltyp, ivgtyp, and xland can be different from that from surface_input_source = 1 (a check later in real program effectively matched ivgtyp based on isltyp); and the other effect is that for certain cases, the WRF model would stop due to un-matched isltyp and ivgtyp. This change fixes this problem, and isltyp, ivgtyp, and xland are now identical to those coming from surface_input_source = 1.
A new soil parameter table is added (SOILPARM.TBL_Kishne_2017). The new table is based on the Kishne et al. paper. The new values reduce dry soil moisture threshold, porosity, and field capacity for some categories, and increases wilting point for all categories. In general it favors more evaporation from the soil. To test the model using this table, copy this table to the default name SOILPARM.TBL. (Kishné, Andrea Sz, Yimam, Y. T., Morgan, C. L. S., and Dornblaser, B. C., 2017: Evaluation and improvement of the default soil hydraulic parameters for the noah land surface model. Geoderma, 285, 247-259.)
For the initial releases of major version 4 (WRF v4.0, v4.0.1, v4.0.2, and v4.0.3) , when the option use_theta_m=1 (the default) was activated, a number of diagnostic utilities in the WRF model incorrectly used a moist potential temperature (when a dry potential temperature was required). These include: DFI, nest adjustment for T and Qv, time series (the *.TH file), the p- and z-level vertically interpolated diagnostics (impact T and Qv), and the RASM mean_output_calc and diurnalcycle_output_calc utilities. In areas of tropical moisture, the temperature errors could be in excess of 8 K (while winter-time Canadian shield errors are less than 0.5 K).
Since the height level interpolation routine was introduced in v3.8, it has used an incorrect formulation for the AGL computations.The AGL computation of the height-based diagnostics instead computed a height above mean sea level. In areas of high terrain, even with a 6.5 K/km lapse rate, large errors near the intended surface would have occurred. This has now been corrected.
A few commonly used variables (i.e., precipitable water, relative humidity, total rain, total liquid rain, potential temperature) are included in the output of traditional fields. By setting diag_nwp2=1 (&diags), one can get these variables.
Fixed a bug that causes AFWA diagnostics option to fail or hang at the time of writing output. Note also that the full set of AFWA diagnostics only works with a small number of microphysics options at the moment.
After hrrr_cycling was introduced to WRF, the stochastic solutions were not bit reproducible when compared to the solutions that used a restart. Also, with the newer GNU compilers, the size of the seed array exceeded the vertical dimensions in WRF (the vertical dimension was originally used as an easy and safe value). Both problems have been addressed.
For the real program, the order of the vertical interpolation option is now always set to a linear interpolation for meteorological fields that are supposed to be positive definite (for example, moist process constituents and number concentration). However, the RH field in the real program continues to use a user-defined vertical interpolation option.
OBSGRID input data may now be used with V4.1 real.exe, regardless of the version of OBSGRID.

WRFDA

New Features:

Himawari AHI radiance data assimilation
Wang, Y., Z. Liu, S. Yang, J. Min, L. Chen, Y. Chen, and T. Zhang, 2018:
Added value of assimilating Himawari-8 AHI water vapor radiances on analyses
and forecasts for "7.19" severe storm over north China. J. Geophys. Res. Atmos.,
123, https://doi.org/10.1002/2017JD027697.
3DVAR and hybrid-3D/4DEnVar now work with moist potential temperature.

Updated Features:

Update WRFDA to handle RTTOV v12.1. New emissivity formulations are enabled as well as HDF-only emissivity atlas files. Compiling WRFDA for RTTOV now requires an HDF5 library.
The code-base of CRTM carried with WRFDA is updated from v2.2.3 to v2.3.0.
Starting with V4.0, CRTM coefficient files are NOT included in any of the WRF or WRFDA tar files. There are two ways to get this data.
1. A subset of coefficient files can be download from
http://www2.mmm.ucar.edu/wrf/users/wrfda/download/crtm_coeffs.html
2. The full set of CRTM coefficients may be downloaded from:
ftp://ftp.emc.ncep.noaa.gov/jcsda/CRTM/REL-2.3.0/crtm_v2.3.0.tar.gz
For external users, replace the link ‘var/run/crtm_coeffs’ with the directory
contained in the appropriate tar file. A copy is also available under ~wrfhelp
on the NCAR cheyenne supercomputer (for internal use).
Enhance the gts_omb_oma diagnostics by appending time slot information after the level information.
Allow pseudo ob at not-first time slot for 4DVAR.
New WRFDA namelist options for specifying errors for radar rhv(rrn/rsn/rgr)
New WRFDA namelist options for turning on/off some diagnostic output
Fix to allow assimilation of SSMI TPW contained (prior to Nov 2009) in prepbufr obs files when ob_format=1 and use_ssmiretrievalobs=.true.

Note: See http://www2.mmm.ucar.edu/wrf/users/wrfda/updates-4.1.html for a full list of updates. For more information about WRFDA, visit the WRFDA Users home page
http://www2.mmm.ucar.edu/wrf/users/wrfda/index.html

WRF-Chem

Correction to the mapping of MEGAN emissions of A-PINENE and B-PINENE for the CRIv2R5 gas-phase chemical scheme (those constituents were swapped).
It was found that the dust (and sea salt) mass balance is violated in routine settling (module module_gocart_settling.F). In particular, the total mass of the dust (and sea salt) is increasing, when the dust (and sea salt) settles from the upper to the lower layers. In the code the dust (and sea salt) mixing ratios are transferred between the layers, which is not correct, since the air density is changing. The discretization scheme was replaced to the one, which accounts for air density change and allows to conserve the dust (and sea salt) total mass.
Added CH4 and CO2 initialization to RACM_SOA_VBS_KPP (chem_opt=108) and RACM_SOA_VBS_KPP_AQCHEM (chem_opt=109). Added NH3 and CH4 to plumerise transport for biomass_burn_opt=1. Updated CO2 to 400 ppm.
The sub-grid cloud effect to the optical depth in radiation (cu_rad_feedback = .true.) requires setting additional time-averaged variables for cumulus physics (cu_diag = 1). This option now works with the Grell-Devenyi (GD) cumulus scheme, joining Grell-Freitas (GF), Grell 3D (G3), and Kain–Fritsch Cumulus Potential (KF-CuP). These changes only impact WRF-Chem runs.

HWRF

Operational 2018 HWRF updates to the WRF model, including:
(i) RRTMG cloud overlap options
(ii) Ensemble perturbations to PBL scheme (GFSEDMF)

Note: The active HWRF developers are using a fork in the NCAR organization (gitHub.com/NCAR/HWRFdev). Anyone interested in HWRF development should contact hwrf-help at ucar dot edu.

WRF-Hydro

Updated WRF-Hydro code to v5.0.3 - associated WRF-Hydro release notes can be found here: https://github.com/NCAR/wrf_hydro_nwm_public/blob/bugfix-5.0.x/NEWS.md

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by kkeene44 about 6 years ago

WRF-ARW - WRF Version 4.0.3

The WRF model has been updated to Version 4.0.3 on December 18, 2018.

Acknowledgements: We would like to thank Xu Zhang (Shanghai Typhoon Institute), J. Fidel González-Rouco (UCM & IGEO-CSIC, Spain), Jorge Navarro and Elena García-Bustamante (CIEMAT, Spain), and Carlos Ross (DFM Consultores) for their contributions to this release.

This is a bug fix release. Following are fixes associated with this release.

Dynamics:

LES options in V4.0 with use_theta_m=1 (default) contained an error with surface fluxes leading to a cooler PBL. QFX term has been added to surface theta_m flux to correct this problem.
Fixed a non-bit-for-bit restart problem with the diff_6th_opt option. The error was due to the start and end indices too close to the domain boundaries when doing 6th order horizontal diffusions. Because the loop indices are changes, the results from any run that uses this option will also change. This was a bug since this option was added to the model in V3.1. Thanks to Jorge Navarro for the initial report that led to the discovery of the bug.

Physics:

Fix for a problem that occurred if a user were to use NLCD and MODIS landuse static data together for an urban application. All urban land use categories were assigned type 13, which is correct for MODIS, but NLCD urban high-resolution categories should be assigned 24, 25, and 26 (and should be preferentially used, when available).

Ideal cases:

Idealized cases will now utilize the “set_physics_rconfigs” and “check_nml_consistency” subroutines from share/module_check_a_mundo.F. This will ensure that namelist requirements are not overlooked, and that the correct number of soil layers is assigned to cases, depending on what (if any) LSM scheme is used.

Data Assimilation:

A syntax error was introduced with the v4.0 release that breaks the build of the WRFPLUS code (and only the WRFPLUS code) v19.0.1 of the Intel compiler. An OpenMP parallel loop was incorrectly constructed. That syntax error has been removed, though it is important to note that WRFPLUS never was intended to work with OpenMP, and the WRFPLUS code continues to only work with serial or DM build options.

Ndown:

A fix to the Registry.EM_COMMON file to correct a problem causing a few variables to be identically zero in the fine-grid ndown run or nesting run when input_from_file = .false.

Configuration:

WRF now supports the correct compiler flag for OpenMP builds with newer Intel compilers, v18.0 and later. The Intel-recommended flag (-qopenmp) was tested successfully with Intel 16.0.1, 17.0.1, and 18.0.1.

Obs-Nudging:

After removing a hard-coded limit for an assumed maximum number of MPI tasks, the observation nudging code for WRF now supports more than 1024 MPI tasks. If users previously ran the obs-nudging code with 1024 or fewer MPI tasks, the original code is OK; however, if users attempted to run obs-nudging with >1024 MPI tasks, likely the code died with a segmentation fault while trying to access an address for an array index that was not available.

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by kkeene44 over 6 years ago

WRF-ARW - WRF Version V4.0.2

WRF Model Version 4.0.2: UPDATES

The WRF model has been updated to Version 4.0.2 on November 9, 2018.

The WRF Preprocessing System (WPS) has also been updated to Version 4.0.2.

Acknowledgements: We would like to thank Tanya Smirnova (NOAA), Alexander Davies (US Naval Academy), Jeremy Silver (U. Melbourne), @sodoesburningbus, and Alexander Ukhov (KAUST) for their contributions to this release.

This is a bug fix release. Following are fixes associated with this release.

Dynamics:
- Fix zero-valued d02 dry potential temperature perturbation from horizontal interpolation.This fix does not change the WRF model simulation results. However, if users grab the initial condition dry potential temperature perturbation field, "T" in the netcdf file, that value will be identically zero (the field is initialized within the solver routine for all subsequent time steps). Most WRF post-processors will assume that the field is a horizontally and vertically uniform dry potential temperature = 300 K, which would provide almost reasonable-looking values for temperature near the surface. After the first time step, the model results will be correct. This bug is in v4.0 and v4.0.1.
Microphysics:
- A loop counter and an array index are changed from real to integers in the Morrison 2-moment scheme with consideration of CESM-NCSU RCP4.5 climatological aerosol. This change makes the code able to compile on Cray.
- Explicitly name a small value as a PARAMETER so that the small value can be represented with single precision on Cray. This is also for the Morrison 2-moment scheme with consideration of CESM-NCSU RCP4.5 climatological aerosol.
Cumulus parameterization:
- A minor modification was made to the Grell 3 scheme to prevent random seed array bounds errors in the cup_forcing_ens_3d subroutine
- Correct a few syntaxes in multi-scale Kain-Fritsch to make the code able to compile on Cray
Land surface scheme
- Bug fix in RUC LSM for the computation of LAI for particular day when LAI parameter is read in from the VEGPARM.TBL. The bug makes the value of LAI smaller than it should be.
- Bug fix for Noah LSM when rdalai2d is set to true. In areas where LAI is zero (and vegetation fraction isn’t), unrealistic sensible and latent heat fluxes are produced due to this bug, leading to model crash.
PBL scheme
- The MYJ and QNSE PBLs use surface water vapor mixing ratio as low boundary condition for moisture and compute it from the latent heat flux, while other PBL schemes use directly latent heat flux as low boundary condition. Therefore, there is a special computation of latent heat flux inside RUC LSM to get correct surface Qv. In original code the check in the surface driver is made on the option of surface layer scheme, and not PBL scheme. Therefore, if MYJ PBL is not used together with MYJ surface layer scheme, there will be inaccuracy in low boundary condition for moisture.
- A check was implemented to stop the real or wrf program if bl_mynn_edmf (which is set to 1 by default in the MYNN scheme) is being used in combination with either the ishallow or shcu_physics option turned on. A message informs the user that bl_mynn_edmf cannot be used with the other 2 options and the namelist must be modified.
Ideal case:
- The default value for vegfra in idealized cases was modified from a fraction to a percentage, consistent with the Noah LSM driver
Data Assimilation:
- Bug fix for WRFDA RANDOMCV when put_rand_seed=.false. and the compiler has a minimum size of random seed array larger than 10.
- Fix WRFDA compilation when RTTOV-11 is built with HDF5-enabled.
Chemistry:
- Correct the chemical reaction rate constant for reaction: SO2 + OH -> SO4. This change applies only to the following chem_opt options: MOZART_KPP,MOZCART_KPP,T1_MOZCART_KPP,MOZART_MOSAIC_4BIN_KPP,MOZART_MOSAIC_4BIN_AQ_KPP.
- Bug fix in optical_prep_gocart to better represent dust particles that have radii less than 0.46 microns. With this fix, the underestimation of AOD is overcome.
- Bug fix for the contribution calculation of dust and sea salt bins into the MOZAIC bins. This fix leads to increased contribution of small dust and sea salt particles to AOD.

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by smileMchen over 6 years ago

WRF-ARW - WRF Version 4.0.1

WRF Model Version 4.0.1: UPDATES

The WRF model has been updated to Version 4.0.1 on October 2, 2018.
The WRF Preprocessing System (WPS) has also been updated to Version 4.0.1.

Acknowledgements: We would like to thank Mike Barlage, Jim Bresch, Greg Thompson, and May Wong for their contributions to this release.

This is a bug fix release. Following are fixes associated with this release.

Dynamics:

Hybrid vertical coordinate - In the original implementation of the hybrid vertical coordinate (v3.9), the cpp directives utilized capitalization to differentiate between the use of full-level (c1f, for example) and half-level (c1h, for example) 1d coefficient arrays. For several locations in module_small_step_em.F, the use of the "mu" array with the "t_2" field incorrectly was assigned to full levels. These errors are minor.
Moist potential temperature - In the porting of the theta_m mods from being ONLY inside of the WRF model to being part of the initialization step in the real and ideal programs, the start_em file was modified. Part of the computation of the pressure involves dry_theta * ( 1 + Rv/Rd Qv ). When the incoming data from the real or ideal program has the namelist option use_theta_m==1, then the factor ( 1 + Rv/Rd Qv ) was already included in the theta field. This error is dramatic for the initial conditions (more than 30 hPa in some locations), but the error dies out over time.

Physics:

Microphysics, Thompson:
- Bug has been in the code for approximately 10 years, skipping a reduced number of rain.
- Increased maximum ice concentration permitting 10000/Liter whereas previously 500/Liter upper limit.
- A series of new checks using "EQV" construct to provide some loops to be skipped.
- Minor changes to lower limits of precip to eliminate really tiny values (~1E-9).
- Removed extraneous ELSE blocks in reading lookup tables.
- Eliminate air density - fix for extreme model top heights by Tanya Smirnovich for FIM model (and MPAS).
- Set a higher value for min dust concentration when nucleating ice.
Land Surface Model, NoahMP:
- Noah-MP frozen soil initialization- An incorrect sign change was introduced in v4.0, impacting soil moisture and soil temperature initialization.
- Array out of bounds Noah-MP - Fix possible/likely array out of bounds by assuming homogeneous soil with depth.Only applies to opt_run=2.
- Noah-MP snow liquid water movement - prevent excessive gravitational water movement. Fixes unrealistic snow density values during melt season.
- Noah-MP divide by zero - Bug fix in v4.0 introduced a possible divide by zero when LAI is zero.
- Noah-MP leaf aerodynamic resistance - limit leaf aerodynamic resistance to prevent very large canopy exchange coefficients with high wind speed.

Fire:

Wrong namelist entries given in test/em_real/namelist.input.fire

Data Assimilation:

Cheyenne path change for crtm_coeffs_2.2.3 corrected for new file system directory structure.
An array out-of-bounds error is fixed when cloud_cv_options is set (0, 1, or 2) and mp_physics is set (any value besides 0 or 98) such that any of xa%qrn, xa%qcw, xa%qci, xa%qsn, or xa%qgr are disabled and such that moist scalars qc, qr, qi, qs, and qg are enabled. See registry.var for more info. Error mode without fix: WRFDA will have a fatal error during the analysis update when bounds checking is turned on (debug mode). It is also possible that those moist scalars and correlated variables will be incorrectly updated when bounds-checking is turned off.
Bug fix for max_ext_its>=2 and cloud/w BE scaling factors (i.e. var_scaling6..11 and len_scaling6..11) that vary with outerloops.
Bug fix for ob.radar that contains rv=-999.99, rv_err=-888888.0, rf_err=-888888.0 to avoid excess check_max_iv prints.

Chemistry:

convert_emiss compilation - Minor bug fix for convert_emiss compilation issue. The build target for convert_emiss was accidentally removed when error detection (not recognizing an argument for the ./compile script) was introduced with v4.0. Users will not be able to build the WRF Chem executable convert_emiss without this fix.

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by kkeene44 over 6 years ago

WRF-ARW - WRF Version 4.0

WRF Model Version 4.0: UPDATES

WRF Model Update
The WRF model tar file has been updated to Version 4.0 on June 8, 2018. For known problems in V4.0, please click here.

The WRF Pre-Processing System (WPS) has been updated to Version 4.0.

Note: V4.0 WRF can use wrfinput and wrfbdy files from earlier versions of the model. But one must use appropriate options to do so (see below). If one would like to try the latest options, one must re-run program real.exe from this version of the code.

As always, be cautious when using new options.

If you are interested in seeing how V4.0 has been tested, click here.

New in Version 4.0:

Dynamics:

The hybrid sigma-pressure vertical coordinate introduced in V3.9 is now the default. It can be turned of by using namelist option hybrid_opt=0. See doc/README.hybrid_vert_coord for details.
Temperature variable is now moist theta by default. This has been an option since V3.7, controlled by namelist use_theta_m.

Physics

Microphysics:

P3 (Predicted Particle Property scheme) microphysics: Added an option for including the second free array for ice. Added a new data file p3_lookup_table_2.dat. Activated by using mp_physics = 52. (Thanks to Morrison of NCAR)
Thompson-Eidhammer microphysics (mp_physics = 28): Added a surface dust emissions scheme (near copy of GOCART) to link into the ice-friendly aerosol variable. This option is activated by namelist option dust_emis = 1. (Thanks to Thompson, M. Xu and T. Eidhammer of NCAR)

Cumulus/Microphysics:

An option is added use CESM-NCSU (North Carolina State University) prescribed RCP4.5 climatological aerosol data (spatially and temporally varying) to interact with cloud simulated by MSKF (cu_physics=11) and a special aerosol Morrison scheme (mp_physics=40). A subgrid-scale cloud microphysics scheme of Song and Zhang (2011) is added in MSKF, and the aerosol activation and ice nucleation parameterizations of Song and Zhang (2011) are added into a version of Morrison scheme (mp_physics = 40). The cloud produced by the MSKF is passed onto radiation. Radar reflectivity from the microphysics in MSKF can be output as mskf_refl_10cm. The options are activated by namelist aercu_opt (= 0, no aerosol interaction; = 1, interaction with MSKF only; = 2, interaction with both MSKF and Morrison). The aerosol concentration can be changed by a namelist aercu_fct (default value is 1.). The aerosol comes in on 30 pressure levels, and has 10 types. The input data can be found on Cheyenne in /glade/p/work/wrfhelp/WRF_files, and on the web at http://www2.mmm.ucar.edu/wrf/src/wrf_files/CESM_RCP4.5_Aerosol_Data.tar.gz. A user needs to link one of these files to CESM_RCP4.5_Aerosol_Data.dat before running. (Thanks to Timothy Glotfelty, Patrick Hawbecker, and Kiran Alapaty of EPA)

Cumulus:

KSAS (cu_physics=14): A scale-aware verson of NSAS. Should be run with shcu_physics = 4, the NSAS shallow scheme, which is part of the NSAS scheme. The old NSAS scheme has been moved to cu_physics = 96. (Thanks to Korea Instituted of Atmospheric Prediction Systems or KIAPS and Jihyeon Jang of NCAR).

Land-Surface Model:

Noah-MP:

Gecros crop model (Genotype-by-Environment interaction on CROp growth Simulator, Germany) is added to NoahMP and can be activated by option opt_crop = 2. (Thanks to J. Ingwersen of U. Hohenheim and M. Barlage of NCAR)
Added capability for Noah-MP to read soil composition data and use pedotransfer (PT) function to determine soil properties (opt_soil in namelist). PT function coefficients are controlled from MPTABLE. Users can add additional PT functions controlled by opt_pedo in namelist. Special version of GEOGRID.TBL and datasets added in WPS. (Thanks to M. Barlage of NCAR)
Variables for NoahMP is now in a separate registry file, registry.noahmp (Thanks to M. Barlage of NCAR)

Radiation:

RRTMK (ra_sw_physics=14, ra_lw_physics=14): A version of RRTMG improved by KIAPS: G-packed McICA was sped up without losing accuracy, and two-stream approximation (TSA) for scattering process is revised. Overall speed of the code isn't changed. (Thanks to Korea Instituted of Atmospheric Prediction Systems. Reference: Baek, S. (2017), A revised radiation package of G-packed McICA and two-stream approximation: Performance evaluation in a global weather forecasting model, J. Adv. Model. Earth Syst., 9, doi:10.1002/2017MS000994)

Improvements and Bug fixes:

Physics:

RAP/HRRR Physics: (Thanks to Joe Olson, Tanya Smirnova, and others at GSD/NOAA)

MYNN PBL (bl_pbl_physics=5) and MYNN surface layer (sf_sfclay_physics=5):

(a) Enhancement of mass-flux scheme, tweak of mixing length and stratus component of subgrid clouds. Changed default option to EDMF and mixing length option 2.
(b) SPP fixes: SPP turned off for subgrid cloud decay; and missing height information for length scale is added.
(c) MYNN surface layer: Changed 2 m exchange coefficient calculation which affects 2-m diagnostics, but not model solution;
(d) Switched to use THv for TH* calculation to be consistent with RIB calculation and further relaxed some limits on PSIT, PSIX, and UST following Jimenez et al. (2012). Small impact. Using Li et al. (2010, BLM) for RIB -> z/L mapping. Changed default z0,zt,zq over water to be COARE 3.0 (was 3.5 in v3.9).

RUC LSM (sf_surface_physics = 3):

(a) Changed empirical formulation of frozen precipitation density to match better observed snow accumulation on the ground;
(b) Applied snow "mosaic" approach only when temperature at the first atmospheric level is above 271K. This change prevents from two warm skin temperature over shallow snow;
(c) Modified 2-m diagnostics of temp and mixing ratio to not use 'flux' method.

Microphysics

NSSL scheme: (Thanks to Ted Mansell of NSSL)

Added a new default ice crystal fall speed formula (icefallopt = 3), which has faster speeds for small ice particles. The main effect is on anvil clouds to help them settle a bit faster. The old behavior can be recovered with icefallopt=1;
Added an option to change ice crystal and snow fall speeds by an arbitrary factor, but default factors are 1.0 (no impact with default settings);
Added a calculation of reflectivity-weighted snow fall speed instead of using mass-weighted (for size-sorting limiter);
Added a check for large supersaturation w.r.t. liquid at ice-only temperatures and nucleate/condense if needed;
Added option to limit accretion of droplets smaller than "exwmindiam" by setting a "reserve value" that is excluded from accretion so that a minimum value of mass and number remains. (Default is no limit, i.e., exwmindiam = 0.0 );
Added two options for rime density on graupel/hail (irimdenopt). Default is original Heymsfield and Pflaum formula. Extra options are for Cober and List (1993) and Macklin;
Added a rime density calculation for rain accreted by graupel instead of always assuming high density;
Added a calculation of a separate rate for number of rain drops produced by snow melting rather than using the snow number loss rate;
Separated rain-ice freezing from flag ifrzg, and is now controlled by ifiacrg (no impact -- for sensitivity study only)
Snow aggregation coefficient set back to 1.0 (increases aggregation rate);
Made consistent use of dtpinv to replace divisions by dtp (no impact).

WDM6 (mp_physics = 16): Added dependency of qcr (cloud and rain water conversion rate) on land and ocean. The lower limit of nccn is set to clean air value over ocean, and this improves effective radius calculation and radiation feedback (From KIAPS).
Morrison schemes (mp_physics = 10, 40): A new namelist, morr_rimed_ice, is added to control graupel or hail option for these schemes. This was controlled by hail_opt. The default value is set to use hail (1). (From H. Morrison)

Planetary boundary layer and surface layer physics:

Q2 diagnostics: limited to 5% over lowest model level QV to avoid unrealistic spikes at sunrise and sunset times.
ACM2 (bl_pbl_physics=7): used height instead eta values to compute layer structure. This is more consistent with MPAS, but also gives some improvement for 2 m T in WRF (Thanks to R. Gilliam and J. Pleim of EPA)
MM5 surface layer (sf_sfclay_physics = 1, 91): added a option to use Czil = 0.1 by setting iz0tlnd = 2. This provides a way to change the sensitivity of exchange coefficient.
MYJ suface layer (sf_sfclay_physics = 2): Allowed to work with LES by providing frictional velocity for momentum from the scheme. ARW only.

Land Surface physics:

NoahMP: a patch communication is needed for the groundwater model with generic dataset. The bug was introduced in 3.9. (Thanks to M. Barlage of RAL/NCAR).
Noah: output variable CANWAT is incorrectly scaled on output by 1000. No effect to model solutions (Thanks to P. Campbell of EPA).
Several variables, RA (aerodynamic resistance), RS (stomatal resistance), WSPD, LAI_mosaic, RS_mosaic, are diagnosed from Noah, Noah-mosaic and NoahMP and can be added to the output for CMAQ applications (EPA). (Thanks to P. Campbell of EPA)
Urban physics: Added variable TSK_RURAL for BEP/BEM with Noah to prevent possible model failure when aggregated TSK is disaggregated (Thanks to Alberto Martilli of CIEMAT, Spain)
PX LSM: removed arbitrarily high values of vegetation fraction and leaf area index computed internally for rangeland, scrubland etc. in western US. This improves the moist and cool biases and high precipitation biases. Also added data for 21 category MODIS. (Thanks to R. Gilliam of EPA)
SSiB LSM: fixed a bug in soil layer initialization.

Ocean:

sf_ocean_physics = 1: Allowed simple mixed layer ocean model to be turned off in DFI, and turned back on in forward integration.

Radiation:

Two fixes for using CAM climatological ozone (o3input=2) with RRTMG: 1) if one uses physics suite specification in V3.9, then o3input is mistakenly reset to 0 for all domains (This means the constant ozone profile is used); 2) if one uses a nest, the nest didn't use any ozone data. This bug was introduced in V3.9 as well.
Clear-sky surface diagnostics for direct and direct normal surface solar radiation are added (Thanks to Pedro Jimenez of NCAR).
Use of GHG: The model will now stop if input file CAMtr_volume_mixing_ratio is not found, rather than reverting to default values. (Thanks to Claire Carouge of ARC Centre of Excellence for Climate Systems Science, Australia)

Others:

Gravity wave drag (gwd_opt=1): 1) Gravity wave drag generalized for all projections. Winds are rotated to lat/long at beginning of GWD code so that orographic data from geogrid can be used accurately. Tendencies are rotated back afterwards. The stop when the grid is rotated too much relative to lat/lon, which was introduced in V3.9, is removed. 2) Fixed a problem with excessive orographic stress in unstable conditions (Thanks to M. Toy of NOAA and Hyun-Joo Choi of KIAPS).
WRF-Fire updates. (Thanks to Domingo Munoz-Esparza of RAL/NCAR)

Proper parallelization of the entire level-set algorithm using MPI
3rd-order Runge-Kutta time integration scheme for the level-set equations
3rd- and 5th-order weighted essentially non-oscillatory (WENO) schemes for the discretization of the advection term in the level-set equations [also 2nd- and 4th-order schemes are available]
A hybrid-order level-set method with locally reduced artificial viscosity (high-order over a local band
Reinitialization PDE for the level set (also using hybrid-order discretization)
Transport and dispersion of fire-generated smoke as a passive tracer
Diagnosed flame length
Option to extrapolate mid-flame height winds to reduce over-coupling of atmospheric fields impacted by the fire
Initialization from observed fire area
Scott & Burgan (40) fuel model crosswalk to Anderson (13)
Several bug fixes
Modified test/em_fire/namelist.input* files.

Hailcast option (hailcast_opt = 1): Added haildt to avoid computing at every model time step; Added the use of the adiabatic liquid cloud water profile to account for typical convective-allowing resolutions not fully resolving the precipitation-free bounded weak echo region; the output variable is changed to maximum hail size rather than the mean and stand deviation; The model will stop if this option is used in a domain where a cumulus scheme is activated (Thanks to Becky Adams-Selin of AER and Conrad Ziegler of NSSL/NOAA).
Spectral nudging of water vapor: added and controlled by namelist. A capability to place an upper-layer limit on spectral nudging of QV and T via a new namelist variable "ktrop" is added too. (Thanks to Tanya Spero of National Exposure Research Lab).
Obs-nudging: Fixed a bug for variable reserf when model levels are greater than 100 (Thanks to B. Reen of Army Research Laboratory)
output_diagnostics = 1: Fixed a potential issue in climate extreme computation where roundoff errors may result in negative value for square root computation. (Thanks to C. Carouge of ARC Centre of Excellence for Climate Systems Science, Australia)
Physics sharing with MPAS: MPAS directives are added to a few physics options (new Tiedtke CPS, YSU PBL, Noah LSM, WSM6) so that potentially it will allow MPAS to get the latest versions of WRF physics.

Initialization/Input/Output:

Added a new way to computed vertical levels, which has a smooth variation of dz (see User's Guide, Chapter 5 for details). This method is the default and can be changed by the new namelist auto_levels_opt (= 1 gives the old set of levels). Together with namelists dzbot, dzstretch_s, dzstretch_u (new) and max_dz, one can adjust to give you different number of vertical levels.
force_use_old_data: Set to true, if using wrfinput data from older versions of the model. This is a new namelist in &time_control. When using old input data, one should also set hybrid_opt =0 and use_theta_m = 0 in &dynamics, since the defaults for V4 are hybrid_opt =2 and use_theta_m = 1. If this option isn't active, checks in the model will make suggestions to you.
A general module_initialize_ideal.F now handles ideal cases quarter_ss, squall2d_x, squall2d_y, b_wave, les, seabreeze2d_x, convrad, hill2d_x and grav2d_x. The case is determined by namelist variable ideal_case (added to namelist template in each ideal case directories). Compile-time selection is not changed, since it still needs to know which test directory to link executables to.
force_ideal.nc in test/em_scm_xy/: String 'V4' is added to the TITLE of global attribute so that it can be used by V4 model.
real and wrf: in V3.9, checks were added for mixed layer ocean and input, gwd and input, urban and input to catch problems with using those options but no input data. But checks also stop the model if there is input for a particular option, but the option was turned off. This is fixed in V4.0.
Urban fraction initialization: To use new NLCD FRC_URB2D input file (see WPS updates) from WPS/geogrid program properly in the model, a match between this data and other datasets (such as MODIS or USGS landuse, which may cover the non-urban areas) is put in.
The default input SNOALB is changed from old 0.15 deg NCEP data to new 0.05 deg MODIS snow albedo data. If one uses usemonalb = .true. (default value is false), the default dataset from geogrid is changed from 0.15 deg NCEP data to new 0.05 deg MODIS data.
RSL files: Removed non-fatal messages from these files.
A mechanism to output a small selected fields is added. This writes to auxiliary stream 1. It writes total pressure, total geopotential height, temperature, pressure, rotated u and v (earth relative), wind speed, mslp. This can be activated by using namelist option diag_nwp2 = 1. This option will be developed further in the coming months.

Dynamics:

diff_6th_opt was fixed for conservation. The option can be turned off if the terrain slope exceeds a namelist specified threshold value (diff_6th_slopeopt=1, and diff_6th_thresh=0.1).
Options are added to turn on/off 2nd and 6th order diffusion for all 4D arrays.
Diffusion for lat/lon grid: corrected map-scale factors for two terms. (Thanks to GSD/NOAA).

Namelist Changes:

Added a new namelist template, namelist.input.pbl-les, in test/em_real directory for running model from km to LES scales.
Simplied to reduce user errors: removed many variables those values cannot be changed. Examples are num_soil_layers, spec_zone, relax_zone (can be derived from spec_zone and spec_bdy_width), {start, end}_{minute,second}, debug_level,
e_vert in all real-data namelist templates are increased to 33 based on the new default method to compute eta levels. If you set eta_levels, there is no restrictions. But caution must be used to ensure there are adequate number of vertical levels.
damp_opt for real-data templates are changed from 0 to 3.
num_metgrid_levels in all real-data namelist templates are changed from 27 (number of old GFS data levels) to 32 (new GFS data levels).

Removed:

SLOPECAT field: removed from WPS and model.
directory dyn_exp/: removed.
run_restart.tar in test/em_b_wave/: removed.

Software:

The top model directory is now WRF/.
The default build is now netCDF4 with compression if netCDF4 is built with HDF5 library. This allows WRF netCDF input and output files written with compression. If the desired library doesn't exist, it reverts to classic netCDF. It is noted that the compression may be slower than non-compression writes, but the advantage is it writes files with smaller sizes.
Determine the resultant x- and y-direction computational patch sizes based on the full domain size and the number of MPI ranks used in each direction. If the patch size < 10 grid cells in either direction, the model will stop.
Support for OpenMPI version 1 is removed..
Cleaned up configure and compile scripts.

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by davegill almost 7 years ago

WRF-ARW - Version 3.9

This is the version 3.9 release of the WRF code. The release includes all components of the system:

WRF model (including ARW and NMM/HWRF)
WRF-Chem
WRFDA

WPS and WRFPLUS have also been updated, but these codes are kept separately from this repository. Click those respective links for further details.

The highlights of the new release are:

WRF ARW model:

Hybrid sigma-pressure vertical coordinate as an alternative vertical coordinate
P3 microphysics
Stochastically perturbed parameter scheme and stochastically perturbed physics tendencies
Improvement to RAP/HRRR physics and various other physics
Runtime performance improvement through optimization of advection routines

See http://www2.mmm.ucar.edu/wrf/users/wrfv3.9/updates-3.9.html for full details

WRFDA:

New assimilation capability: 4DEnVar
Cloudy radiance assimilation capability for AMSR2
Radar "null-echo" assimilation

See http://www2.mmm.ucar.edu/wrf/users/wrfda/updates-3.9.html for more information

WRF-Chem:

TUV photolysis option
Coupled WSM6 microphysics with MOZART gas wet scavenging
Modified GOCART dust scheme for the SORGAM and MOSAIC chemistry options
A new trajectory option that calculates meteorological and chemical properties
along air trajectories

WRF-NMM-HWRF:

The 2016 operational update, which includes surface layer and ocean coupling updates; the GFS EDMF PBL scheme, and a scale aware SAS cumulus scheme.

For detailed code changes, updates and to download the latest version of the code, please visit http://www2.mmm.ucar.edu/wrf/users/downloads.html and http://www.dtcenter.org/HurrWRF/users

We would like to thank many people who contributed to the v3.9 release:

Hugh Morrison, Jason Milbrandt, Judith Berner, Joe Olson, Tanya Smirnova,
Georg Grell, Francisco Salamanca, Yizhou Zhang, Xing Liu, Gonzalo Miguez-Macho,
Mike Barlage, Chunxi Zhang, Sebastien Masson, Jihyeon Jang, Kun Liu, Becky Adams-Selin,
Ted Mansell, Greg Thompson, Raquel Lorente, Pedro Jimenez, Jiachun Yang, Yu Gu,
Mei Xu, Tim Glotfelty, Kiran Alapaty, Josh Laughner, Brian Reen, Milan Curcic,
Negin Sobhani, Davide Del Vento, Samuel Elliot, Patricia Balle, Pete Johnsen,
Laura Fowler, Robert Gilliam, Katie Lundquist, William Gustafson, Jim Bresch,
Kevin Manning, Mark Seefeldt, Jose Renteria, Wei Yu, Saulo Freitas, John Cazes,
Balwinder Singh, Megan Bela, Chun Zhao, D. Lowe, K. Wang, P. Tuccella, Branko Kosovic,
S. McKeen, R. McVay, Mary Barth, P. Saide, M. Hirtl, Y. Zhang, G. Pfister, Y-H Ryu,
David Avisar, Hongli Wang, Yu-Shin Kim, Ki-Hong Min, Yuanbing Wang, Xiaowen Tang,
Thomas Schwitalla, Yuanbing Wang, Jeffrey Willison, Chun Yang, Heimo Truhetz,
Nils Gustafsson, Feng Gao, National Water team, NCEP HWRF team, Fujitsu Japan and Taiwan

Atmosphere - Meteorological Observation and Forecast - Fortran
Published by mkavulich almost 8 years ago