The programming language python in earth system simulations

-scale vary from a planetary scale and million years for convection problems to 100km and 10 years for fault systems simulations. Various techniques are in use to deal with the time dependency (e.g. Crank-Nicholson), with the non-linearity (e.g. Newton-Raphson) and weakly coupled equations (e.g. non-linear Gauss-Seidel). Besides these high-level solution algorithms discretization methods (e.g. finite element method (FEM), boundary element method (BEM)) are used to deal with spatial derivatives. Typically, large-scale, three dimensional meshes are required to resolve geometrical complexity (e.g. in the case of fault systems) or features in the solution (e.g. in mantel convection simulations). The modelling environment escript allows the rapid implementation of new physics as required for the development of simulation codes in earth sciences. Its main object is to provide a programming language, where the user can define new models and rapidly develop high-level solution algorithms. The current implementation is linked with the finite element package finley as a PDE solver. However, the design is open and other discretization technologies such as finite differences and boundary element methods could be included. escript is implemented as an extension of the interactive programming environment python (see www.python.org). Key concepts introduced are Data objects, which are holding values on nodes or elements of the finite element mesh, and linearPDE objects, which are defining linear partial differential equations to be solved by the underlying discretization technology. In this paper we will show the basic concepts of escript and will show how escript is used to implement a simulation code for interacting fault systems. We will show some results of large-scale, parallel simulations on an SGI Altix system. Acknowledgements: Project work is supported by Australian Commonwealth Government through the Australian Computational Earth Systems Simulator Major National Research Facility, Queensland State Government Smart State Research Facility Fund, The University of Queensland and SGI.

Rough Notes on Programming AVR Microcontrollers in C

These notes follow on from the material that you studied in CSSE1000 Introduction to Computer Systems. There you studied details of logic gates, binary numbers and instruction set architectures using the Atmel AVR microcontroller family as an example. In your present course (METR2800 Team Project I), you need to get on to designing and building an application which will include such a microcontroller. These notes focus on programming an AVR microcontroller in C and provide a number of example programs to illustrate the use of some of the AVR peripheral devices.

A scheme for efficient quantum computation wtih linear optics

Quantum computers promise to increase greatly the efficiency of solving problems such as factoring large integers, combinatorial optimization and quantum physics simulation. One of the greatest challenges now is to implement the basic quantum-computational elements in a physical system and to demonstrate that they can be reliably and scalably controlled. One of the earliest proposals for quantum computation is based on implementing a quantum bit with two optical modes containing one photon. The proposal is appealing because of the ease with which photon interference can be observed. Until now, it suffered from the requirement for non-linear couplings between optical modes containing few photons. Here we show that efficient quantum computation is possible using only beam splitters, phase shifters, single photon sources and photo-detectors. Our methods exploit feedback from photo-detectors and are robust against errors from photon loss and detector inefficiency. The basic elements are accessible to experimental investigation with current technology.

Turbulence in Small Sub-Tropical Estuary with Semi-Diurnal Tides

In natural estuaries, contaminant transport is driven by the turbulent momentum mixing. The predictions of scalar dispersion can rarely be predicted accurately because of a lack of fundamental understanding of the turbulence structure in estuaries. Herein detailed turbulence field measurements were conducted at high frequency and continuously for up to 50 hours per investigation in a small subtropical estuary with semi-diurnal tides. Acoustic Doppler velocimetry was deemed the most appropriate measurement technique for such small estuarine systems with shallow water depths (less than 0.5 m at low tides), and a thorough post-processing technique was applied. The estuarine flow is always a fluctuating process. The bulk flow parameters fluctuated with periods comparable to tidal cycles and other large-scale processes. But turbulence properties depended upon the instantaneous local flow properties. They were little affected by the flow history, but their structure and temporal variability were influenced by a variety of mechanisms. This resulted in behaviour which deviated from that for equilibrium turbulent boundary layer induced by velocity shear only. A striking feature of the data sets is the large fluctuations in all turbulence characteristics during the tidal cycle. This feature was rarely documented, but an important difference between the data sets used in this study from earlier reported measurements is that the present data were collected continuously at high frequency during relatively long periods. The findings bring new lights in the fluctuating nature of momentum exchange coefficients and integral time and length scales. These turbulent properties should not be assumed constant.

Molecular cocrystals of carboxylic acids. XXVII - An investigation into the use of triphenylphosphine oxide cocrystals for non-linear optics: the crystal structure of the adduct of triphenylphosphine oxide with N-methylpyrrole-2-carboxylic acid

Four adducts of triphenylphosphine oxide with aromatic carboxylic acids have been synthesized and tested for second-order non-linear optical properties. These were with N-methylpyrrole-2-carboxylic acid (I), indole-2-carboxylic acid (2), 3-dimethylaminobenzoic acid (3), and thiophen-2-carboxylic acid (4). Compound (1) produced clear, colourless crystals (space group P2(1)2(1)2(1) With a 9.892(1), b 14.033(1), c 15.305(1) Angstrom, Z 4) which allowed the structure to be determined by X-ray diffraction.

Ester prodrugs of a potent analgesic, morphine-6-sulfate: syntheses, spectroscopic and physicochemical properties

The aim of this work is to develop 3-acyl prodrugs of the potent analgesic morphine-6-sulfate (M6S). These are expected to have higher potency and/or exhibit longer duration of analgesic action than the parent compound. M6S and the prodrugs were synthesized, then purified either by recrystallization or by semi-preparative HPLC and the structures confirmed by mass spectrometry, IR spectrophotometry and by detailed 1- and 2-D NMR studies. The lipophilicities of the compounds were assessed by a combination of shake-flask, group contribution and HPLC retention methods. The octanol-buffer partition coefficient could only be obtained directly for 3-heptanoylmorphine-6-sulfate, using the shake-flask method. The partition coefficients (P) for the remaining prodrugs were estimated from known methylene group contributions. A good linear relationship between log P and the HPLC log capacity factors was demonstrated. Hydrolysis of the 3-acetyl prodrug, as a representative of the group, was found to occur relatively slowly in buffers (pH range 6.15-8.01), with a small buffer catalysis contribution. The rates of enzymatic hydrolysis of the 3-acyl group in 10% rat blood and in 10% rat brain homogenate were investigated. The prodrugs followed apparent first order hydrolysis kinetics, with a significantly faster hydrolysis rate found in 10% rat brain homogenate than in 10% rat blood for all compounds. (C) 1998 Elsevier Science B.V. All rights reserved.

Response (re-)programming in aging: A kinematic analysis

Background. Age-related motor slowing may reflect either motor programming deficits, poorer movement execution, or mere strategic preferences for online guidance of movement. We controlled such preferences, limiting the extent to which movements could be programmed. Methods. Twenty-four young and 24 older adults performed a line drawing task that allowed movements to he prepared in advance in one case (i.e., cue initially available indicating target location) and not in another (i.e., no cue initially available as to target location). Participants connected large or small targets illuminated by light-emitting diodes upon a graphics tablet that sampled pen tip position at 200 Hz. Results. Older adults had a disproportionate difficulty initiating movement when prevented from programming in advance. Older adults produced slower, less efficient movements, particularly when prevented from programming under greater precision requirements. Conclusions. The slower movements of older adults do not simply reflect a preference for online control, as older adults have less efficient movements when forced to reprogram their movements. Age-related motor slowing kinematically resembles that seen in patients with cerebellar dysfunction.

Temperature and sowing date affect the linear increase of sunflower harvest index

The linearity of daily linear harvest index (HI) increase can provide a simple means to predict grain growth and yield in field crops. However, the stability of the rate of increase across genotypes and environments is uncertain. Data from three field experiments were collated to investigate the phase of linear HI increase of sunflower (Helianthus annuus L,) across environments by changing genotypes, sowing time, N level, and solar irradiation level. Linear increase in HI was similar among different genotypes, N levels, and radiation treatments (mean 0.0125 d(-1)). but significant differences occurred between sowings, The linear increase in HI was not stable at very low temperatures (down to 9 degrees C) during grain filling, due to possible limitations to biomass accumulation and translocation (mean 0.0091 d(-1)). Using the linear increase in HI to predict grain yield requires predictions of the duration from anthesis to the onset of linear HI increase (lag phase) and the cessation of linear RT increase. These studies showed that the lag phase differed, and the linear HI increase ceased when 91% of the anthesis to physiological maturity period had been completed.

The Heisenberg antiferromagnet on an anisotropic triangular lattice: linear spin-wave theory

We consider the effect of quantum spin fluctuations on the ground-state properties of the Heisenberg antiferromagnet on an anisotropic triangular lattice using linear spin-wave (LSW) theory. This model should describe the magnetic properties of the insulating phase of the kappa-(BEDT-TTF)(2)X family of superconducting molecular crystals. The ground-state energy, the staggered magnetization, magnon excitation spectra, and spin-wave velocities are computed as functions of the ratio of the antiferromagnetic exchange between the second and first neighbours, J(2)/J(1). We find that near J(2)/J(1) = 0.5, i.e., in the region where the classical spin configuration changes from a Neel-ordered phase to a spiral phase, the staggered magnetization vanishes, suggesting the possibility of a quantum disordered state. in this region, the quantum correction to the magnetization is large but finite. This is in contrast to the case for the frustrated Heisenberg model on a square lattice, for which the quantum correction diverges logarithmically at the transition from the Neel to the collinear phase. For large J(2)/J(1), the model becomes a set of chains with frustrated interchain coupling. For J(2) > 4J(1), the quantum correction to the magnetization, within LSW theory, becomes comparable to the classical magnetization, suggesting the possibility of a quantum disordered state. We show that, in this regime, the quantum fluctuations are much larger than for a set of weakly coupled chains with non-frustrated interchain coupling.

Phase diagram for a class of spin-1/2 Heisenberg models interpolating between the square-lattice, the triangular-lattice, and the linear-chain limits

We study the spin-1/2 Heisenberg models on an anisotropic two-dimensional lattice which interpolates between the square lattice at one end, a set of decoupled spin chains on the other end, and the triangular-lattice Heisenberg model in between. By series expansions around two different dimer ground states and around various commensurate and incommensurate magnetically ordered states, we establish the phase diagram for this model of a frustrated antiferromagnet. We find a particularly rich phase diagram due to the interplay of magnetic frustration, quantum fluctuations, and varying dimensionality. There is a large region of the usual two-sublattice Neel phase, a three-sublattice phase for the triangular-lattice model, a region of incommensurate magnetic order around the triangular-lattice model, and regions in parameter space where there is no magnetic order. We find that the incommensurate ordering wave vector is in general altered from its classical value by quantum fluctuations. The regime of weakly coupled chains is particularly interesting and appears to be nearly critical. [S0163-1829(99)10421-1].

Programming internet based DAI applications in Qu-Prolog

This paper presents the unique collection of additional features of Qu-Prolog, a variant of the Al programming language Prolog, and illustrates how they can be used for implementing DAI applications. By this we mean applications comprising communicating information servers, expert systems, or agents, with sophisticated reasoning capabilities and internal concurrency. Such an application exploits the key features of Qu-Prolog: support for the programming of sound non-clausal inference systems, multi-threading, and high level inter-thread message communication between Qu-Prolog query threads anywhere on the internet. The inter-thread communication uses email style symbolic names for threads, allowing easy construction of distributed applications using public names for threads. How threads react to received messages is specified by a disjunction of reaction rules which the thread periodically executes. A communications API allows smooth integration of components written in C, which to Qu-Prolog, look like remote query threads.

Optimal release strategies for biological control agents: an application of stochastic dynamic programming to population management

1. Establishing biological control agents in the field is a major step in any classical biocontrol programme, yet there are few general guidelines to help the practitioner decide what factors might enhance the establishment of such agents. 2. A stochastic dynamic programming (SDP) approach, linked to a metapopulation model, was used to find optimal release strategies (number and size of releases), given constraints on time and the number of biocontrol agents available. By modelling within a decision-making framework we derived rules of thumb that will enable biocontrol workers to choose between management options, depending on the current state of the system. 3. When there are few well-established sites, making a few large releases is the optimal strategy. For other states of the system, the optimal strategy ranges from a few large releases, through a mixed strategy (a variety of release sizes), to many small releases, as the probability of establishment of smaller inocula increases. 4. Given that the probability of establishment is rarely a known entity, we also strongly recommend a mixed strategy in the early stages of a release programme, to accelerate learning and improve the chances of finding the optimal approach.