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.

Migrating eprints.org data to a Fez repository

This document records the process of migrating eprints.org data to a Fez repository. Fez is a Web-based digital repository and workflow management system based on Fedora (http://www.fedora.info/). At the time of migration, the University of Queensland Library was using EPrints 2.2.1 [pepper] for its ePrintsUQ repository. Once we began to develop Fez, we did not upgrade to later versions of eprints.org software since we knew we would be migrating data from ePrintsUQ to the Fez-based UQ eSpace. Since this document records our experiences of migration from an earlier version of eprints.org, anyone seeking to migrate eprints.org data into a Fez repository might encounter some small differences. Moving UQ publication data from an eprints.org repository into a Fez repository (hereafter called UQ eSpace (http://espace.uq.edu.au/) was part of a plan to integrate metadata (and, in some cases, full texts) about all UQ research outputs, including theses, images, multimedia and datasets, in a single repository. This tied in with the plan to identify and capture the research output of a single institution, the main task of the eScholarshipUQ testbed for the Australian Partnership for Sustainable Repositories project (http://www.apsr.edu.au/). The migration could not occur at UQ until the functionality in Fez was at least equal to that of the existing ePrintsUQ repository. Accordingly, as Fez development occurred throughout 2006, a list of eprints.org functionality not currently supported in Fez was created so that programming of such development could be planned for and implemented.

Protein sequence threading, the alignment problem, and a two-step strategy

Conventionally, protein structure prediction via threading relies on some nonoptimal method to align a protein sequence to each member of a library of known structures. We show how a score function (force field) can be modified so as to allow the direct application of a dynamic programming algorithm to the problem. This involves an approximation whose damage can be minimized by an optimization process during score function parameter determination. The method is compared to sequence to structure alignments using a more conventional pair-wise score function and the frozen approximation. The new method produces results comparable to the frozen approximation, but is faster and has fewer adjustable parameters. It is also free of memory of the template's original amino acid sequence, and does not suffer from a problem of nonconvergence, which can be shown to occur with the frozen approximation. Alignments generated by the simplified score function can then be ranked using a second score function with the approximations removed. (C) 1999 John Wiley & Sons, Inc.

Effects of a program of intervention on parental distress following infant death

A longitudinal study of 144 patents (65 fathers, 79 mothers) was conducted to evaluate the effectiveness of a program of intervention in relieving the psychological distress of parents affected by infant death. Participants were assessed in terms of their psychiatric disturbance, depression, anxiety, physical symptoms, dyadic adjustment, and coping strategies. The experimental group (n = 84) was offered an intervention program comprising the use of specially designed resources and contact with a trained grief worker. A control group (n = 60) was given routine community care. Parental reactions were assessed at four to six weeks postloss (prior to the implementation of the intervention program), at six months postloss, and at 15 months postloss. A series of multivariate analyses of valiance revealed that the intervention was effective in reducing the distress of parents, particularly those assessed prior to the intervention as being at high-risk of developing mourning difficulties. Effects of the intervention were noted in terms of parents' overall psychiatric disturbance, marital quality, and paternal coping strategies.

Use of advance information in patients with schizophrenia

To document possible motor disturbance in schizophrenia, we examined the ability to use advance information (or cues) to plan movements in a sequential button pressing task in 12 Clozapine medicated patients. Programming of movements under various cues revealed that patients with schizophrenia, relative to controls, initiated movements slower to the right than left, providing possible evidence for right hemineglect (left hemisphere dysfunction). Additionally, patients with schizophrenia had difficulty in the initiation of movements in the absence of a cue, suggesting internal cue generation difficulty for movement related to some form of fronto-striatal disturbance. Motor abnormalities were predominantly observed at the level of movement initiation, but not execution, contrary to basal ganglia disorders such as Parkinson's and Huntington's disease.

An L-system-based plant modeling language

Cpfg is a program for simulating and visualizing plant development, based on the theory of L-systems. A special-purpose programming language, used to specify plant models, is an essential feature of cpfg. We review postulates of L-system theory that have influenced the design of this language. We then present the main constructs of this language, and evaluate it from a user's perspective.

Sperm environment affects offspring quality in broadcast spawning marine invertebrates

The provisioning of offspring can have far-reaching consequences for later life in a wide range of organisms and generally this provisioning is thought to be under maternal influence or control. In experiments with a broadcast-spawning ascidian, we found that the size of offspring was determined by egg size and the abundance of sperm present during fertilization. Larger eggs were fertilized at low sperm concentrations, whilst smaller eggs were successfully fertilized at high sperm concentrations. These differences in fertilized egg size resulted in differences in the development rate, hatching success and mean size of the subsequent larvae. Our results suggest that, in contrast to females that reproduce by other mating systems, free-spawning mothers lack some control over the provisioning of offspring. Furthermore, because males can alter the sperm environment, they can exert paternal (non-genetic) control over key offspring characteristics.

A preliminary phylogeny of the scale insects (Hemiptera : Sternorrhyncha : Coccoidea) based on nuclear small-subunit ribosomal DNA

Scale insects (Hemiptera: Sternorrhyncha: Coccoidea) are a speciose and morphologically specialized group of plant-feeding bugs in which evolutionary relationships and thus higher classification are controversial. Sequences derived from nuclear small-subunit ribosomal DNA were used to generate a preliminary molecular phylogeny for the Coccoidea based on 39 species representing 14 putative families. Monophyly of the archaeococcoids (comprising Ortheziidae, Margarodidae sensu lato, and Phenacoleachia) was equivocal, whereas monophyly of the neococcoids was supported. Putoidae, represented by Puto yuccae, was found to be outside the remainder of the neococcoid clade. These data are consistent with a single origin (in the ancestor of the neococcoid clade) of a chromosome system involving paternal genome elimination in males. Pseudococcidae (mealybugs) appear to be sister to the rest of the neococcoids and there are indications that Coccidae (soft scales) and Kerriidae (lac scales) are sister taxa. The Eriococcidae (felt scales) was not recovered as a monophyletic group and the eriococcid genus Eriococcus sensu lato was polyphyletic. (C) 2002 Elsevier Science (USA). All rights reserved.

Darwinian grandparenting: Preferential investment in more certain kin

Studies on grandparental investment have revealed that mothers fathers are emotionally closer to their grandchildren than are fathers' mothers. In the current study, it was hypothesized that this difference is caused by the fact that fathers' mothers often have the potential to invest in genetically more certain kin (children through their daughters). To test this hypothesis, 787 participants rated their emotional closeness and exposure to their grandparents and indicated whether they had cousins through paternal and maternal aunts and uncles. Results indicated that participants felt closer to mothers' fathers than fathers' mothers only when alternate investment outlets for fathers' mothers were available. Closeness ratings to fathers fathers also were reduced when they had grandchildren through their daughters. Exposure to grandparents revealed a similar pattern of findings but did not show the same sensitivity to the presence of more certain kin and did not appear to account for the closeness ratings.

Pectinase production by solid fermentation from Aspergillus niger by a new prescription experiment

The Ordos Plateau in China is covered with up to 300,000 ha of peashrub (Caragana) which is the dominant natural vegetation and ideal for fodder production. To exploit peashrub fodder, it is crucially important to optimize the culture conditions, especially culture substrate to produce pectinase complex. In this study, a new prescription process was developed. The process, based on a uniform experimental design, first optimizes the solid substrate and second, after incubation, applies two different temperature treatments (30 degrees C for the first 30 h and 23 degrees C for the second 42 h) in the fermentation process. A multivariate regression analysis is applied to a number of independent variables (water, wheat bran, rice dextrose, ammonium sulfate, and Tween 80) to develop a predictive model of pectinase activity. A second-degree polynomial model is developed which accounts for an excellent proportion of the explained variation (R-2 = 97.7%). Using unconstrained mathematical programming, an optimized substrate prescription for pectinase production is subsequently developed. The mathematical analysis revealed that the optimal formula for pectinase production from Aspergillus niger by solid fermentation under the conditions of natural aeration, natural substrate pH (about 6.5), and environmental humidity of 60% is rice dextrose 8%, wheat bran 24%, ammonium sulfate ((NH4)(2)SO4) 6%, and water 61%. Tween 80 was found to have a negative effect on the production of pectinase in solid substrate. With this substrate prescription, pectinase produced by solid fermentation of A. niger reached 36.3IU/(gDM). Goats fed on the pectinase complex obtain an incremental increase of 0.47 kg day(-1) during the initial 25 days of feeding, which is a very promising new feeding prospect for the local peashrub. It is concluded that the new formula may be very useful for the sustainable development of and and semiarid pastures such as those of the Ordos Plateau. (c) 2005 Elsevier Inc. All rights reserved.

Risk attitude, planting conditions and the value of seasonal forecasts to a dryland wheat grower

The value of a seasonal forecasting system based on phases of the Southern Oscillation was estimated for a representative dryland wheat grower in the vicinity of Goondiwindi. In particular the effects on this estimate of risk attitude and planting conditions were examined. A recursive stochastic programming approach was used to identify the grower's utility-maximising action set in the event of each of the climate patterns over the period 1894-1991 recurring In the imminent season. The approach was repeated with and without use of the forecasts. The choices examined were, at planting, nitrogen application rate and cultivar and, later in the season, choices of proceeding with or abandoning each wheat activity, The value of the forecasting system was estimated as the maximum amount the grower could afford to pay for its use without expected utility being lowered relative to its non use.

Refining specifications to logic programs

The refinement calculus provides a framework for the stepwise development of imperative programs from specifications. In this paper we study a refinement calculus for deriving logic programs. Dealing with logic programs rather than imperative programs has the dual advantages that, due to the expressive power of logic programs, the final program is closer to the original specification, and each refinement step can achieve more. Together these reduce the overall number of derivation steps. We present a logic programming language extended with specification constructs (including general predicates, assertions, and types and invariants) to form a wide-spectrum language. General predicates allow non-executable properties to be included in specifications. Assertions, types and invariants make assumptions about the intended inputs of a procedure explicit, and can be used during refinement to optimize the constructed logic program. We provide a semantics for the extended logic programming language and derive a set of refinement laws. Finally we apply these to an example derivation.