scholarly journals libcloudph++ 0.2: single-moment bulk, double-moment bulk, and particle-based warm-rain microphysics library in C++

2014 ◽  
Vol 7 (6) ◽  
pp. 8275-8360
Author(s):  
S. Arabas ◽  
A. Jaruga ◽  
H. Pawlowska ◽  
W. W. Grabowski

Abstract. This paper introduces a library of algorithms for representing cloud microphysics in numerical models. The library is written in C++, hence the name libcloudph++. In the current release, the library covers three warm-rain schemes: the single- and double-moment bulk schemes, and the particle-based scheme with Monte-Carlo coalescence. The three schemes are intended for modelling frameworks of different dimensionality and complexity ranging from parcel models to multi-dimensional cloud-resolving (e.g. large-eddy) simulations. A two-dimensional prescribed-flow framework is used in example simulations presented in the paper with the aim of highlighting the library features. The libcloudph++ and all its mandatory dependencies are free and open-source software. The Boost.units library is used for zero-overhead dimensional analysis of the code at compile time. The particle-based scheme is implemented using the Thrust library that allows to leverage the power of graphics processing units (GPU), retaining the possibility to compile the unchanged code for execution on single or multiple standard processors (CPUs). The paper includes complete description of the programming interface (API) of the library and a performance analysis including comparison of GPU and CPU setups.

2015 ◽  
Vol 8 (6) ◽  
pp. 1677-1707 ◽  
Author(s):  
S. Arabas ◽  
A. Jaruga ◽  
H. Pawlowska ◽  
W. W. Grabowski

Abstract. This paper introduces a library of algorithms for representing cloud microphysics in numerical models. The library is written in C++, hence the name libcloudph++. In the current release, the library covers three warm-rain schemes: the single- and double-moment bulk schemes, and the particle-based scheme with Monte Carlo coalescence. The three schemes are intended for modelling frameworks of different dimensionalities and complexities ranging from parcel models to multi-dimensional cloud-resolving (e.g. large-eddy) simulations. A two-dimensional (2-D) prescribed-flow framework is used in the paper to illustrate the library features. The libcloudph++ and all its mandatory dependencies are free and open-source software. The Boost.units library is used for zero-overhead dimensional analysis of the code at compile time. The particle-based scheme is implemented using the Thrust library that allows one to leverage the power of graphics processing units (GPU), retaining the possibility of compiling the unchanged code for execution on single or multiple standard processors (CPUs). The paper includes a complete description of the programming interface (API) of the library and a performance analysis including comparison of GPU and CPU set-ups.


2017 ◽  
Vol 17 (17) ◽  
pp. 10195-10221 ◽  
Author(s):  
Constantino Listowski ◽  
Tom Lachlan-Cope

Abstract. The first intercomparisons of cloud microphysics schemes implemented in the Weather Research and Forecasting (WRF) mesoscale atmospheric model (version 3.5.1) are performed on the Antarctic Peninsula using the polar version of WRF (Polar WRF) at 5 km resolution, along with comparisons to the British Antarctic Survey's aircraft measurements (presented in part 1 of this work; Lachlan-Cope et al., 2016). This study follows previous works suggesting the misrepresentation of the cloud thermodynamic phase in order to explain large radiative biases derived at the surface in Polar WRF continent-wide (at 15 km or coarser horizontal resolution) and in the Polar WRF-based operational forecast model Antarctic Mesoscale Prediction System (AMPS) over the Larsen C Ice Shelf at 5 km horizontal resolution. Five cloud microphysics schemes are investigated: the WRF single-moment five-class scheme (WSM5), the WRF double-moment six-class scheme (WDM6), the Morrison double-moment scheme, the Thompson scheme, and the Milbrandt–Yau double-moment seven-class scheme. WSM5 (used in AMPS) and WDM6 (an upgrade version of WSM5) lead to the largest biases in observed supercooled liquid phase and surface radiative biases. The schemes simulating clouds in closest agreement to the observations are the Morrison, Thompson, and Milbrandt schemes for their better average prediction of occurrences of clouds and cloud phase. Interestingly, those three schemes are also the ones allowing for significant reduction of the longwave surface radiative bias over the Larsen C Ice Shelf (eastern side of the peninsula). This is important for surface energy budget consideration with Polar WRF since the cloud radiative effect is more pronounced in the infrared over icy surfaces. Overall, the Morrison scheme compares better to the cloud observation and radiation measurements. The fact that WSM5 and WDM6 are single-moment parameterizations for the ice crystals is responsible for their lesser ability to model the supercooled liquid clouds compared to the other schemes. However, our investigation shows that all the schemes fail at simulating the supercooled liquid mass at some temperatures (altitudes) where observations show evidence of its persistence. An ice nuclei parameterization relying on both temperature and aerosol content like DeMott et al. (2010) (not currently used in WRF cloud schemes) is in best agreement with the observations, at temperatures and aerosol concentration characteristic of the Antarctic Peninsula where the primary ice production occurs (part 1), compared to parameterization only relying on the atmospheric temperature (used by the WRF cloud schemes). Overall, a realistic double-moment ice microphysics implementation is needed for the correct representation of the supercooled liquid phase in Antarctic clouds. Moreover, a more realistic ice-nucleating particle alone is not enough to improve the cloud modelling, and water vapour and temperature biases also need to be further investigated and reduced.


2013 ◽  
Vol 13 (1) ◽  
pp. 1489-1526 ◽  
Author(s):  
D. Jarecka ◽  
H. Pawlowska ◽  
W. W. Grabowski ◽  
A. A. Wyszogrodzki

Abstract. This paper discusses aircraft observations and large-eddy simulation (LES) of the 15 May 2008, North Sea boundary-layer clouds from the EUCAARI-IMPACT field campaign. These clouds were advected from the north-east by the prevailing lower-tropspheric winds, and featured stratocumulus-over-cumulus cloud formations. Almost-solid stratocumulus deck in the upper part of the relatively deep weakly decoupled marine boundary layer overlaid a field of small cumuli with a cloud fraction of ~10%. The two cloud formations featured distinct microphysical characteristics that were in general agreement with numerous past observations of strongly-diluted shallow cumuli on the one hand and solid marine boundary-layer stratocumulus on the other. Macrophysical and microphysical cloud properties were reproduced well by the double-moment warm-rain microphysics large-eddy simulation. A novel feature of the model is its capability to locally predict homogeneity of the subgrid-scale mixing between the cloud and its cloud-free environment. In the double-moment warm-rain microphysics scheme, the homogeneity is controlled by a single parameter α, that ranges from 0 to 1 and limiting values representing the homogeneous and the extremely inhomogeneous mixing scenarios, respectively. Parameter α depends on the characteristic time scales of the droplet evaporation and of the turbulent homogenization. In the model, these scales are derived locally based on the subgrid-scale turbulent kinetic energy, spatial scale of cloudy filaments, the mean cloud droplet radius, and the humidity of the cloud-free air entrained into the cloud. Simulated mixing is on average quite inhomogeneous, with the mean parameter α around 0.7 across the entire depth of the cloud field, but with local variations across almost the entire range, especially near the base and the top of the cloud field.


Atmosphere ◽  
2019 ◽  
Vol 10 (3) ◽  
pp. 114 ◽  
Author(s):  
Dalton Behringer ◽  
Sen Chiao

This study investigated precipitation distribution patterns in association with atmospheric rivers (ARs). The Weather Research and Forecasting (WRF) model was employed to simulate two strong atmospheric river events. The precipitation forecasts were highly sensitive to cloud microphysics parameterization schemes. Thus, radar observed and simulated Z H and Z D R were evaluated to provide information about the drop-size distribution (DSD). Four microphysics schemes (WSM-5, WSM-6, Thompson, and WDM-6) with nested simulations (3 km, 1 km, and 1/3 km) were conducted. One of the events mostly contained bright-band (BB) rainfall and lasted less than 24 h, while the other contained both BB and non-bright-band (NBB) rainfall, and lasted about 27 h. For each event, there was no clear improvement in the 1/3 km model, over the 1 km model. Overall, the WDM-6 microphysics scheme best represented the rainfall and the DSD. It appears that this scheme performed well, due to its relative simplicity in ice and mixed-phase microphysics, while providing double-moment predictions of warm rain microphysics (i.e., cloud and rain mixing ratio and number concentration). The other schemes tested either provided single-moment predictions of all classes or double-moment predictions of ice and rain (Thompson). Considering the shallow nature of precipitation in atmospheric rivers and the high-frequency of the orographic effect enhancing the warm rain process, these assumptions appear to be applicable over the southern San Francisco Bay Area.


2019 ◽  
Vol 35 (17) ◽  
pp. 3181-3183 ◽  
Author(s):  
Patryk Orzechowski ◽  
Jason H Moore

Abstract Motivation In this paper, we present an open source package with the latest release of Evolutionary-based BIClustering (EBIC), a next-generation biclustering algorithm for mining genetic data. The major contribution of this paper is adding a full support for multiple graphics processing units (GPUs) support, which makes it possible to run efficiently large genomic data mining analyses. Multiple enhancements to the first release of the algorithm include integration with R and Bioconductor, and an option to exclude missing values from the analysis. Results Evolutionary-based BIClustering was applied to datasets of different sizes, including a large DNA methylation dataset with 436 444 rows. For the largest dataset we observed over 6.6-fold speedup in computation time on a cluster of eight GPUs compared to running the method on a single GPU. This proves high scalability of the method. Availability and implementation The latest version of EBIC could be downloaded from http://github.com/EpistasisLab/ebic. Installation and usage instructions are also available online. Supplementary information Supplementary data are available at Bioinformatics online.


2015 ◽  
Vol 30 (6) ◽  
pp. 1571-1589 ◽  
Author(s):  
Ki-Hong Min ◽  
Sunhee Choo ◽  
Daehyung Lee ◽  
Gyuwon Lee

Abstract The Korea Meteorological Administration (KMA) implemented a 10-yr project to develop its own global model (GM) by 2020. To reflect the complex topography and unique weather characteristics of the Korean Peninsula, a high-resolution model with accurate physics and input data is required. The WRF single-moment 6-class microphysics scheme (WSM6) and WRF double-moment 6-class microphysics scheme (WDM6) that will be implemented in the Korea GM (KGM) are evaluated. Comparisons of the contoured frequency by altitude diagram (CFAD), time–height cross sections, and vertical profiles of hydrometeors are utilized to assess the two schemes in simulating summer monsoon and convective precipitation cases over the Korean Peninsula during 2011. The results show that WSM6 and WDM6 overestimate the height of the melting level and bright band as compared to radar observations. However, the accuracy of WDM6 is in better agreement with radar observations. This is attributed to the difference in the sedimentation process simulated by the additional second-moment total number concentrations of liquid-phase particles in WDM6. WDM6 creates larger raindrops and higher relative humidity beneath the melting layer, allowing the scheme to simulate a more realistic reflectivity profile than WSM6 for the summer monsoon case. However, for the convective case, both schemes underestimate the precipitation and there is resolution dependence in the WRF Model’s ability to simulate convective precipitation.


2020 ◽  
Vol 36 (11) ◽  
pp. 3570-3572
Author(s):  
Mindaugas Margelevičius

Abstract Summary Searching for homology in the vast amount of sequence data has a particular emphasis on its speed. We present a completely rewritten version of the sensitive homology search method COMER based on alignment of protein sequence profiles, which is capable of searching big databases even on a lightweight laptop. By harnessing the power of CUDA-enabled graphics processing units, it is up to 20 times faster than HHsearch, a state-of-the-art method using vectorized instructions on modern CPUs. Availability and implementation COMER2 is cross-platform open-source software available at https://sourceforge.net/projects/comer2 and https://github.com/minmarg/comer2. It can be easily installed from source code or using stand-alone installers. Contact [email protected] Supplementary information Supplementary data are available at Bioinformatics online.


Sign in / Sign up

Export Citation Format

Share Document