A benchmark-driven modelling approach for evaluating deployment choices on a multi-core architecture

The complexity of current and emerging architectures provides users with options about how best to use the available resources, but makes predicting performance challenging. In this work a benchmark-driven model is developed for a simple shallow water code on a Cray XE6 system, to explore how deployment choices such as domain decomposition and core affinity affect performance. The resource sharing present in modern multi-core architectures adds various levels of heterogeneity to the system. Shared resources often includes cache, memory, network controllers and in some cases floating point units (as in the AMD Bulldozer), which mean that the access time depends on the mapping of application tasks, and the core's location within the system. Heterogeneity further increases with the use of hardware-accelerators such as GPUs and the Intel Xeon Phi, where many specialist cores are attached to general-purpose cores. This trend for shared resources and non-uniform cores is expected to continue into the exascale era. The complexity of these systems means that various runtime scenarios are possible, and it has been found that under-populating nodes, altering the domain decomposition and non-standard task to core mappings can dramatically alter performance. To find this out, however, is often a process of trial and error. To better inform this process, a performance model was developed for a simple regular grid-based kernel code, shallow. The code comprises two distinct types of work, loop-based array updates and nearest-neighbour halo-exchanges. Separate performance models were developed for each part, both based on a similar methodology. Application specific benchmarks were run to measure performance for different problem sizes under different execution scenarios. These results were then fed into a performance model that derives resource usage for a given deployment scenario, with interpolation between results as necessary.

4,4 '-Dipyridyl-N,N '-dioxide complexes of metal-thiocyanate/selenocyanate: pi-stacked molecular rods as three-dimensional support for two-dimensional polymeric sheets and intra/interchain S center dot center dot center dot S interaction dependent architecture of the R-2(2)(8) synthon driven assembly of one-dimensional polymeric chains

Three supramolecular complexes of Co(II) using SCN-/SeCN- in combination with 4,4'-dipyridyl-N,N'-dioxide (dpyo), i.e., {[Co(SCN)(2)(dpyo)(2)].(dpyo)}(n) ( 1), {[Co(SCN)(2)(dpyo)(H2O)(2)].(H2O)}(n) ( 2), {[Co(SeCN)(2)(dpyo)(H2O)(2)]center dot(H2O)}(n) ( 3), have been synthesized and characterized by single-crystal X-ray analysis. Complex 1 is a rare example of a dpyo bridged two-dimensional (2D) coordination polymer, and pi-stacked dpyo supramolecular rods are generated by the lattice dpyo, passing through the rhombic grid of stacked layers, resulting in a three-dimensional (3D) superstructure. Complexes 2 and 3 are isomorphous one-dimensional (1D) coordination polymers [-Co-dpyo-Co-] that undergo self-assembly leading to a bilayer architecture derived through an R-2(2)(8) H-bonding synthon between coordinated water and dpyo oxygen. A reinvestigation of coordination polymers [Mn(SCN)(2)(dpyo)( H2O)(MeOH)](n) ( 4) and {[Fe(SCN)(2)(dpyo)(H2O)(2)]center dot(H2O)}(n) ( 5) reported recently by our group [ Manna et al. Indian J. Chem. 2006, 45A, 1813] reveals brick wall topology rather than bilayer architecture is due to the decisive role of S center dot center dot center dot S/Se center dot center dot center dot Se interactions in determining the helical nature in 4 and 5 as compared to zigzag polymeric chains in 2 and 3, although the same R-2(2)(8) synthon is responsible for supramolecular assembly in these complexes.

A statistical model for sea surface diurnal warming driven by numercial weather predictions fluxes and winds

A statistical model is derived relating the diurnal variation of sea surface temperature (SST) to the net surface heat flux and surface wind speed from a numerical weather prediction (NWP) model. The model is derived using fluxes and winds from the European Centre for Medium-Range Weather Forecasting (ECMWF) NWP model and SSTs from the Spinning Enhanced Visible and Infrared Imager (SEVIRI). In the model, diurnal warming has a linear dependence on the net surface heat flux integrated since (approximately) dawn and an inverse quadratic dependence on the maximum of the surface wind speed in the same period. The model coefficients are found by matching, for a given integrated heat flux, the frequency distributions of the maximum wind speed and the observed warming. Diurnal cooling, where it occurs, is modelled as proportional to the integrated heat flux divided by the heat capacity of the seasonal mixed layer. The model reproduces the statistics (mean, standard deviation, and 95-percentile) of the diurnal variation of SST seen by SEVIRI and reproduces the geographical pattern of mean warming seen by the Advanced Microwave Scanning Radiometer (AMSR-E). We use the functional dependencies in the statistical model to test the behaviour of two physical model of diurnal warming that display contrasting systematic errors.

Experimentally driven atomistic model of 1,2 Polybutadiene

We present an efficient method of combining wide angle neutron scattering data with detailed atomistic models, allowing us to perform a quantitative and qualitative mapping of the organisation of the chain conformation in both glass and liquid phases. The structural refinement method presented in this work is based on the exploitation of the intrachain features of the diffraction pattern and its intimate linkage with atomistic models by the use of internal coordinates for bond lengths, valence angles and torsion rotations. Atomic connectivity is defined through these coordinates that are in turn assigned by pre-defined probability distributions, thus allowing for the models in question to be built stochastically. Incremental variation of these coordinates allows for the construction of models that minimise the differences between the observed and calculated structure factors. We present a series of neutron scattering data of 1,2 polybutadiene at the region 120-400K. Analysis of the experimental data yield bond lengths for C-C and C=C of 1.54Å and 1.35Å respectively. Valence angles of the backbone were found to be at 112° and the torsion distributions are characterised by five rotational states, a three-fold trans-skew± for the backbone and gauche± for the vinyl group. Rotational states of the vinyl group were found to be equally populated, indicating a largely atactic chan. The two backbone torsion angles exhibit different behaviour with respect to temperature of their trans population, with one of them adopting an almost all trans sequence. Consequently the resulting configuration leads to a rather persistent chain, something indicated by the value of the characteristic ratio extrapolated from the model. We compare our results with theoretical predictions, computer simulations, RIS models and previously reported experimental results.

Measuring Enterprise Architecture success: a tentative model for measuring success

Interest towards Enterprise Architecture (EA) has been increasing during the last few years. EA has been found to be a crucial aspect of business survival, and thus the importance of EA implementation success is also crucial. Current literature does not have a tool to be used to measure the success of EA implementation. In this paper, a tentative model for measuring success is presented and empirically validated in EA context. Results show that the success of EA implementation can be measured indirectly by measuring the achievement of the objectives set for the implementation. Results also imply that achieving individual's objectives do not necessarily mean that organisation's objectives are achieved. The presented Success Measurement Model can be used as basis for developing measurement metrics.

Mathematical model for NF-kappaB-driven proliferation of adult neural stem cells

Neural stem cells (NSCs) are early precursors of neuronal and glial cells. NSCs are capable of generating identical progeny through virtually unlimited numbers of cell divisions (cell proliferation), producing daughter cells committed to differentiation. Nuclear factor kappa B (NF-kappaB) is an inducible, ubiquitous transcription factor also expressed in neurones, glia and neural stem cells. Recently, several pieces of evidence have been provided for a central role of NF-kappaB in NSC proliferation control. Here, we propose a novel mathematical model for NF-kappaB-driven proliferation of NSCs. We have been able to reconstruct the molecular pathway of activation and inactivation of NF-kappaB and its influence on cell proliferation by a system of nonlinear ordinary differential equations. Then we use a combination of analytical and numerical techniques to study the model dynamics. The results obtained are illustrated by computer simulations and are, in general, in accordance with biological findings reported by several independent laboratories. The model is able to both explain and predict experimental data. Understanding of proliferation mechanisms in NSCs may provide a novel outlook in both potential use in therapeutic approaches, and basic research as well.

Spherical nonparametric estimators applied to the UGAMP model integration for AMIP

The aim of this paper is essentially twofold: first, to describe the use of spherical nonparametric estimators for determining statistical diagnostic fields from ensembles of feature tracks on a global domain, and second, to report the application of these techniques to data derived from a modern general circulation model. New spherical kernel functions are introduced that are more efficiently computed than the traditional exponential kernels. The data-driven techniques of cross-validation to determine the amount elf smoothing objectively, and adaptive smoothing to vary the smoothing locally, are also considered. Also introduced are techniques for combining seasonal statistical distributions to produce longer-term statistical distributions. Although all calculations are performed globally, only the results for the Northern Hemisphere winter (December, January, February) and Southern Hemisphere winter (June, July, August) cyclonic activity are presented, discussed, and compared with previous studies. Overall, results for the two hemispheric winters are in good agreement with previous studies, both for model-based studies and observational studies.

Glacier shrinkage and climatic change in the Russian Altai from the mid-20th century: An assessment using remote sensing and PRECIS regional climate model

This paper examines changes in the surface area of glaciers in the North and South Chuya Ridges, Altai Mountains in 1952-2004 and their links with regional climatic variations. The glacier surface areas for 2004 were derived from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) imagery. Data from the World Glacier Inventory (WGI)dating to 1952 and aerial photographs from 1952 were used to estimate the changes. 256 glaciers with a combined area of 253±5.1 km2 have been identified in the region in 2004. Estimation of changes in extent of 126 glaciers with the individual areas not less than 0.5 km2 in 1952 revealed a 19.7±5.8% reduction. The observed glacier retreat is primarily driven by an increase in summer temperatures since the 1980s when air temperatures were increasing at a rate of 0.10 - 0.13oC a-1 at the glacier tongue elevation. The regional climate projections for A2 and B2 CO2 emission scenarios developed using PRECIS regional climate model indicate that summer temperatures will increase in the Altai in 2071-2100 by 6-7oC and 3-5oC respectively in comparison with 1961-1990 while annual precipitation will increase by 15% and 5%. The length of the ablation season will extend from June-August to the late April – early October. The projected increases in precipitation will not compensate for the projected warming and glaciers will continue to retreat in the 21st century under both B2 and A2 scenarios.

An adaptive finite element water column model using the Mellor-Yamada level 2.5 turbulence closure scheme

A one-dimensional water column model using the Mellor and Yamada level 2.5 parameterization of vertical turbulent fluxes is presented. The model equations are discretized with a mixed finite element scheme. Details of the finite element discrete equations are given and adaptive mesh refinement strategies are presented. The refinement criterion is an "a posteriori" error estimator based on stratification, shear and distance to surface. The model performances are assessed by studying the stress driven penetration of a turbulent layer into a stratified fluid. This example illustrates the ability of the presented model to follow some internal structures of the flow and paves the way for truly generalized vertical coordinates. (c) 2005 Elsevier Ltd. All rights reserved.

Anthropogenic warming of the oceans: Observations and model results

Observations show the oceans have warmed over the past 40 yr. with appreciable regional variation and more warming at the surface than at depth. Comparing the observations with results from two coupled ocean-atmosphere climate models [the Parallel Climate Model version 1 (PCM) and the Hadley Centre Coupled Climate Model version 3 (HadCM3)] that include anthropogenic forcing shows remarkable agreement between the observed and model-estimated warming. In this comparison the models were sampled at the same locations as gridded yearly observed data. In the top 100 m of the water column the warming is well separated from natural variability, including both variability arising from internal instabilities of the coupled ocean-atmosphere climate system and that arising from volcanism and solar fluctuations. Between 125 and 200 m the agreement is not significant, but then increases again below this level, and remains significant down to 600 m. Analysis of PCM's heat budget indicates that the warming is driven by an increase in net surface heat flux that reaches 0.7 W m(-2) by the 1990s; the downward longwave flux increases bv 3.7 W m(-2). which is not fully compensated by an increase in the upward longwave flux of 2.2 W m(-2). Latent and net solar heat fluxes each decrease by about 0.6 W m(-2). The changes in the individual longwave components are distinguishable from the preindustrial mean by the 1920s, but due to cancellation of components. changes in the net surface heat flux do not become well separated from zero until the 1960s. Changes in advection can also play an important role in local ocean warming due to anthropogenic forcing, depending, on the location. The observed sampling of ocean temperature is highly variable in space and time. but sufficient to detect the anthropogenic warming signal in all basins, at least in the surface layers, bv the 1980s.

Influence of the seasonal cycle on the termination of El Nino events in a coupled general circulation model

In this study, the mechanisms leading to the El Nino peak and demise are explored through a coupled general circulation model ensemble approach evaluated against observations. The results here suggest that the timing of the peak and demise for intense El Nino events is highly predictable as the evolution of the coupled system is strongly driven by a southward shift of the intense equatorial Pacific westerly anomalies during boreal winter. In fact, this systematic late-year shift drives an intense eastern Pacific thermocline shallowing, constraining a rapid El Nino demise in the following months. This wind shift results from a southward displacement in winter of the central Pacific warmest SSTs in response to the seasonal evolution of solar insolation. In contrast, the intensity of this seasonal feedback mechanism and its impact on the coupled system are significantly weaker in moderate El Nino events, resulting in a less pronounced thermocline shallowing. This shallowing transfers the coupled system into an unstable state in spring but is not sufficient to systematically constrain the equatorial Pacific evolution toward a rapid El Nino termination. However, for some moderate events, the occurrence of intense easterly wind anomalies in the eastern Pacific during that period initiate a rapid surge of cold SSTs leading to La Nina conditions. In other cases, weaker trade winds combined with a slightly deeper thermocline allow the coupled system to maintain a broad warm phase evolving through the entire spring and summer and a delayed El Nino demise, an evolution that is similar to the prolonged 1986/87 El Nino event. La Nina events also show a similar tendency to peak in boreal winter, with characteristics and mechanisms mainly symmetric to those described for moderate El Nino cases.

Wave-driven wind jets in the marine atmospheric boundary layer

The interaction between ocean surface waves and the overlying wind leads to a transfer of momentum across the air–sea interface. Atmospheric and oceanic models typically allow for momentum transfer to be directed only downward, from the atmosphere to the ocean. Recent observations have suggested that momentum can also be transferred upward when long wavelength waves, characteristic of remotely generated swell, propagate faster than the wind speed. The effect of upward momentum transfer on the marine atmospheric boundary layer is investigated here using idealized models that solve the momentum budget above the ocean surface. A variant of the classical Ekman model that accounts for the wave-induced stress demonstrates that, although the momentum flux due to the waves penetrates only a small fraction of the depth of the boundary layer, the wind profile is profoundly changed through its whole depth. When the upward momentum transfer from surface waves sufficiently exceeds the downward turbulent momentum flux, then the near-surface wind accelerates, resulting in a low-level wave-driven wind jet. This increases the Coriolis force in the boundary layer, and so the wind turns in the opposite direction to the classical Ekman layer. Calculations of the wave-induced stress due to a wave spectrum representative of fast-moving swell demonstrate upward momentum transfer that is dominated by contributions from waves in the vicinity of the peak in the swell spectrum. This is in contrast to wind-driven waves whose wave-induced stress is dominated by very short wavelength waves. Hence the role of swell can be characterized by the inverse wave age based on the wave phase speed corresponding to the peak in the spectrum. For a spectrum of waves, the total momentum flux is found to reverse sign and become upward, from waves to wind, when the inverse wave age drops below the range 0.15–0.2, which agrees reasonably well with previously published oceanic observations.