Running Climate Models On The NERC Cluster Grid Using G-Rex

Compute grids are used widely in many areas of environmental science, but there has been limited uptake of grid computing by the climate modelling community, partly because the characteristics of many climate models make them difficult to use with popular grid middleware systems. In particular, climate models usually produce large volumes of output data, and running them usually involves complicated workflows implemented as shell scripts. For example, NEMO (Smith et al. 2008) is a state-of-the-art ocean model that is used currently for operational ocean forecasting in France, and will soon be used in the UK for both ocean forecasting and climate modelling. On a typical modern cluster, a particular one year global ocean simulation at 1-degree resolution takes about three hours when running on 40 processors, and produces roughly 20 GB of output as 50000 separate files. 50-year simulations are common, during which the model is resubmitted as a new job after each year. Running NEMO relies on a set of complicated shell scripts and command utilities for data pre-processing and post-processing prior to job resubmission. Grid Remote Execution (G-Rex) is a pure Java grid middleware system that allows scientific applications to be deployed as Web services on remote computer systems, and then launched and controlled as if they are running on the user's own computer. Although G-Rex is general purpose middleware it has two key features that make it particularly suitable for remote execution of climate models: (1) Output from the model is transferred back to the user while the run is in progress to prevent it from accumulating on the remote system and to allow the user to monitor the model; (2) The client component is a command-line program that can easily be incorporated into existing model work-flow scripts. G-Rex has a REST (Fielding, 2000) architectural style, which allows client programs to be very simple and lightweight and allows users to interact with model runs using only a basic HTTP client (such as a Web browser or the curl utility) if they wish. This design also allows for new client interfaces to be developed in other programming languages with relatively little effort. The G-Rex server is a standard Web application that runs inside a servlet container such as Apache Tomcat and is therefore easy to install and maintain by system administrators. G-Rex is employed as the middleware for the NERC1 Cluster Grid, a small grid of HPC2 clusters belonging to collaborating NERC research institutes. Currently the NEMO (Smith et al. 2008) and POLCOMS (Holt et al, 2008) ocean models are installed, and there are plans to install the Hadley Centre’s HadCM3 model for use in the decadal climate prediction project GCEP (Haines et al., 2008). The science projects involving NEMO on the Grid have a particular focus on data assimilation (Smith et al. 2008), a technique that involves constraining model simulations with observations. The POLCOMS model will play an important part in the GCOMS project (Holt et al, 2008), which aims to simulate the world’s coastal oceans. A typical use of G-Rex by a scientist to run a climate model on the NERC Cluster Grid proceeds as follows :(1) The scientist prepares input files on his or her local machine. (2) Using information provided by the Grid’s Ganglia3 monitoring system, the scientist selects an appropriate compute resource. (3) The scientist runs the relevant workflow script on his or her local machine. This is unmodified except that calls to run the model (e.g. with “mpirun”) are simply replaced with calls to "GRexRun" (4) The G-Rex middleware automatically handles the uploading of input files to the remote resource, and the downloading of output files back to the user, including their deletion from the remote system, during the run. (5) The scientist monitors the output files, using familiar analysis and visualization tools on his or her own local machine. G-Rex is well suited to climate modelling because it addresses many of the middleware usability issues that have led to limited uptake of grid computing by climate scientists. It is a lightweight, low-impact and easy-to-install solution that is currently designed for use in relatively small grids such as the NERC Cluster Grid. A current topic of research is the use of G-Rex as an easy-to-use front-end to larger-scale Grid resources such as the UK National Grid service.

A global climatology of wind–wave interaction

Generally, ocean waves are thought to act as a drag on the surface wind so that momentum is transferred downwards, from the atmosphere into the waves. Recent observations have suggested that when long wavelength waves, characteristic of remotely generated swell, propagate faster than the surface wind momentum can also be transferred upwards. This upward momentum transfer acts to accelerate the near-surface wind, resulting in a low-level wave-driven wind jet. Previous studies have suggested that the sign reversal of the momentum flux is well predicted by the inverse wave age, the ratio of the surface wind speed to the speed of the waves at the peak of the spectrum. ECMWF ERA-40 data has been used here to calculate the global distribution of the inverse wave age to determine whether there are regions of the ocean that are usually in the wind-driven wave regime and others that are generally in the wave-driven wind regime. The wind-driven wave regime is found to occur most often in the mid-latitude storm tracks where wind speeds are generally high. The wave-driven wind regime is found to be prevalent in the tropics where wind speeds are generally light and swell can propagate from storms at higher latitudes. The inverse wave age is also a useful indicator of the degree of coupling between the local wind and wave fields. The climatologies presented emphasise the non-equilibrium that exists between the local wind and wave fields and highlight the importance of swell in the global oceans.

Designing Hospital Bedside Systems for use by Patients with Visual Impairments

This paper describes the development of an interface to a hospital portal system for information, communication and entertainment such that it can be used easily and effectively by all patients regardless of their age, disability, computer experience or native language. Specifically, this paper reports on the work conducted to ensure that the interface design took into account the needs of visually impaired users.

Evidence for involvement of both IKCa and SKCa channels in hyperpolarizing responses of the rat middle cerebral artery

Endothelium-derived hyperpolarizing factor responses in the rat middle cerebral artery are blocked by inhibiting IKCa channels alone, contrasting with peripheral vessels where block of both IKCa and SKCa is required. As the contribution of IKCa and SKCa to endothelium-dependent hyperpolarization differs in peripheral arteries, depending on the level of arterial constriction, we investigated the possibility that SKCa might contribute to equivalent hyperpolarization in cerebral arteries under certain conditions. METHODS: Rat middle cerebral arteries (approximately 175 microm) were mounted in a wire myograph. The effect of KCa channel blockers on endothelium-dependent responses to the protease-activated receptor 2 agonist, SLIGRL (20 micromol/L), were then assessed as simultaneous changes in tension and membrane potential. These data were correlated with the distribution of arterial KCa channels revealed with immunohistochemistry. RESULTS: SLIGRL hyperpolarized and relaxed cerebral arteries undergoing variable levels of stretch-induced tone. The relaxation was unaffected by specific inhibitors of IKCa (TRAM-34, 1 micromol/L) or SKCa (apamin, 50 nmol/L) alone or in combination. In contrast, the associated smooth-muscle hyperpolarization was inhibited, but only with these blockers in combination. Blocking nitric oxide synthase (NOS) or guanylyl cyclase evoked smooth-muscle depolarization and constriction, with both hyperpolarization and relaxation to SLIGRL being abolished by TRAM-34 alone, whereas apamin had no effect. Immunolabeling showed SKCa and IKCa within the endothelium. CONCLUSIONS: In the absence of NO, IKCa underpins endothelium-dependent hyperpolarization and relaxation in cerebral arteries. However, when NOS is active SKCa contributes to hyperpolarization, whatever the extent of background contraction. These changes may have relevance in vascular disease states where NO release is compromised and when the levels of SKCa expression may be altered.

Diabatic processes modifying potential vorticity in a North Atlantic cyclone

Potential vorticity (PV) succinctly describes the evolution of large-scale atmospheric flow because of its material conservation and invertibility properties. However, diabatic processes in extratropical cyclones can modify PV and influence both mesoscale weather and the evolution of the synoptic-scale wave pattern. In this investigation, modification of PV by diabatic processes is diagnosed in a Met Office Unified Model (MetUM) simulation of a North Atlantic cyclone using a set of PV tracers. The structure of diabatic PV within the extratropical cyclone is investigated and linked to the processes responsible for it. On the mesoscale, a tripole of diabatic PV is generated across the tropopause fold extending down to the cold front. The structure results from a dipole in heating across the frontal interface due to condensation in the warm conveyor belt flanking the upper side of the fold and evaporation of precipitation in the dry intrusion and below. On isentropic surfaces intersecting the tropopause, positive diabatic PV is generated on the stratospheric side, while negative diabatic PV is generated on the tropospheric side. The stratospheric diabatic PV is generated primarily by long-wave cooling which peaks at the tropopause itself due to the sharp gradient in humidity there. The tropospheric diabatic PV originates locally from the long-wave radiation and non-locally by advection out of the top of heating associated with the large-scale cloud, convection and boundary layer schemes. In most locations there is no diabatic modification of PV at the tropopause itself but diabatic PV anomalies would influence the tropopause indirectly through the winds they induce and subsequent advection. The consequences of this diabatic PV dipole for the evolution of synoptic-scale wave patterns are discussed.

A continuum anisotropic model of sea-ice dynamics

We develop the essential ingredients of a new, continuum and anisotropic model of sea-ice dynamics designed for eventual use in climate simulation. These ingredients are a constitutive law for sea-ice stress, relating stress to the material properties of sea ice and to internal variables describing the sea-ice state, and equations describing the evolution of these variables. The sea-ice cover is treated as a densely flawed two-dimensional continuum consisting of a uniform field of thick ice that is uniformly permeated with narrow linear regions of thinner ice called leads. Lead orientation, thickness and width distributions are described by second-rank tensor internal variables: the structure, thickness and width tensors, whose dynamics are governed by corresponding evolution equations accounting for processes such as new lead generation and rotation as the ice cover deforms. These evolution equations contain contractions of higher-order tensor expressions that require closures. We develop a sea-ice stress constitutive law that relates sea-ice stress to the structure tensor, thickness tensor and strain rate. For the special case of empty leads (containing no ice), linear closures are adopted and we present calculations for simple shear, convergence and divergence.

Collaborative personal profiling for web service ranking and recommendation

Web service is one of the most fundamental technologies in implementing service oriented architecture (SOA) based applications. One essential challenge related to web service is to find suitable candidates with regard to web service consumer’s requests, which is normally called web service discovery. During a web service discovery protocol, it is expected that the consumer will find it hard to distinguish which ones are more suitable in the retrieval set, thereby making selection of web services a critical task. In this paper, inspired by the idea that the service composition pattern is significant hint for service selection, a personal profiling mechanism is proposed to improve ranking and recommendation performance. Since service selection is highly dependent on the composition process, personal knowledge is accumulated from previous service composition process and shared via collaborative filtering where a set of users with similar interest will be firstly identified. Afterwards a web service re-ranking mechanism is employed for personalised recommendation. Experimental studies are conduced and analysed to demonstrate the promising potential of this research.

Purification of coronavirus virions for Cryo-EM and proteomic analysis

Purification of intact enveloped virus particles can be useful as a first step in understanding the structure and function of both viral and host proteins that are incorporated into the virion. Purified preparations of virions can be used to address these questions using techniques such as mass spectrometry proteomics. Recent studies on the proteome of coronavirus virions have shown that in addition to the structural proteins, accessory and non-structural virus proteins and a wide variety of host cell proteins associate with virus particles. To further study the presence of virion proteins, high quality sample preparation is crucial to ensure reproducible analysis by the wide variety of methods available for proteomic analysis.

Resolution and discrimination - two sides of the same coin

The notions of resolution and discrimination of probability forecasts are revisited. It is argued that the common concept underlying both resolution and discrimination is the dependence (in the sense of probability theory) of forecasts and observations. More specifically, a forecast has no resolution if and only if it has no discrimination if and only if forecast and observation are stochastically independent. A statistical tests for independence is thus also a test for no resolution and, at the same time, for no discrimination. The resolution term in the decomposition of the logarithmic scoring rule, and the area under the Receiver Operating Characteristic will be investigated in this light.