This paper describes ZEUS-MP, a multi-physics, massively parallel, message- passing implementation of the ZEUS code. ZEUS-MP differs significantly from the ZEUS-2D code, the ZEUS-3D code, and an early "version 1" of ZEUS-MP distributed publicly in 1999. ZEUS-MP offers an MHD algorithm better suited for multidimensional flows than the ZEUS-2D module by virtue of modifications to the Method of Characteristics scheme first suggested by Hawley and Stone (1995), and is shown to compare quite favorably to the TVD scheme described by Ryu et. al (1998). ZEUS-MP is the first publicly-available ZEUS code to allow the advection of multiple chemical (or nuclear) species. Radiation hydrodynamic simulations are enabled via an implicit flux-limited radiation diffusion (FLD) module. The hydrodynamic, MHD, and FLD modules may be used in one, two, or three space dimensions. Self gravity may be included either through the assumption of a GM/r potential or a solution of Poisson's equation using one of three linear solver packages (conjugate-gradient, multigrid, and FFT) provided for that purpose. Point-mass potentials are also supported. Because ZEUS-MP is designed for simulations on parallel computing platforms, considerable attention is paid to the parallel performance characteristics of each module. Strong-scaling tests involving pure hydrodynamics (with and without self-gravity), MHD, and RHD are performed in which large problems (256^3 zones) are distributed among as many as 1024 processors of an IBM SP3. Parallel efficiency is a strong function of the amount of communication required between processors in a given algorithm, but all modules are shown to scale well on up to 1024 processors for the chosen fixed problem size.

We present light curves for 149 sources monitored with the Green Bank Interferometer. The light curves are at two radio frequencies (approximately 2.5 and 8.2 GHz) and range from 3 to 15 yrs in length, covering the interval 1979--1996, and have a typical sampling of one flux density measurement every 2 days. We have used these light curves to conduct various variability analysis (rms flux density variations and autoregressive, integrated, moving average modeling) of these sources. We find suggestive, though not unambiguous evidence, that these sources have a common, broadband mechanism for intrinsic variations, in agreement with previous studies of a subset of these source. We also find that the sources generally display a short-term variability (~ 10 d) that arises from radio-wave scattering in an extended medium. These conclusions extend those of Fiedler et al. (1987) who used a sub-sample of these data. The primary motivation for this monitoring program was the identification of extreme scattering events. In an effort to identify ESEs in a systematic manner, we have taken the wavelet transform of the light curves. We find 15 events in the light curves of 12 sources that we classify as probable ESEs. However, we also find that five ESEs previously identified from these data do not survive our wavelet selection criteria. Future identification of ESEs will probably continue to rely on both visual and systematic methods. Instructions for obtaining the data are also presented.

We report multifrequency observations of the $\gamma$-ray blazar 0528+134 with the Effelsberg 100-m telescope, the IRAM 30-m telescope, and the NRL Green Bank Interferometer. 0528+134 underwent a major radio and mm outburst in 1993 a few month after a very strong outburst in high energy $\gamma$-rays. A similar behaviour is indicated for the weaker $\gamma$-ray outburst in May 1991 and the mm activity in the beginning of 1992. This result fits nicely to the general finding that blazars are bright in $\gamma$-rays preferrentially at the beginning of a radio outburst. 0528+134 was observed with a global VLBI array at 22$\,$GHz in November 1992. Within the IRIS-S and EUROPE geodetic VLBI observing campaigns the source was also observed regularly at 8.4$\,$GHz. A new VLBI component which is not present in data at 22 $\,$GHz from 1991 also exhibits superluminal motion. Backextra\-polating in time we find that the new component was expelled from the core in the first half of 1991, near the first $\gamma$-ray outburst in Mai 1991. Newest VLBI data indicate the appearance of another new component in 1994. Although this new component has yet been found only at two epochs, it appears to move superluminally and may have been expelled from the core between fall 1992 and summer 1993, i.e. near the time of the second gamma ray outburst and at the beginning of a strong radio outburst.

We report on multifrequency radio observations made with the Effelsberg 100-m telescope, the IRAM 30-m telescope and the Green Bank Interferometer between 1992 and 1994 of the $\gamma$-ray quasar 0528+134. We present a new VLBI based map of 0518+134 at 22 GHz with sub-mas angular resolution observed in November 1992, which shows a one-sided core jet structure of $\sim$ 5 mas length. One component, seen also previously at 8.4 GHz & 22 GHz, shows apparent superluminal motion with $\beta_{app}$=4.4$\pm$1.7. In summer 1993 we observed a major outburst at all frequencies higher than a few GHz, which peaked a few months after a strong outburst in high energy $\gamma$-rays and showed the canonical time evolution of a cooling and/or expanding electron distribution. Our data indicate that the outburst in the EGRET range originated very close to the central object of the AGN and that a remnant of this outburst moved further outward in the jet until it became optically thin at radio frequencies after a few months. During the flare in July 1993 we observed with the Effelsberg 100-m telescope an unusually strong decrease of the flux density by about 50% at 4.75 GHz and 10.55 GHz and slightly less at 2.695 GHz. This behaviour is also seen in the monitoring data at 2.25 GHz and 8.3 GHz taken with the Green Bank Interferometer (NRL-GBI). An extreme scattering event by a small dense plasma cloud in the line of sight is able to match the observed time lag in the lightcurves if we take into account the mas-structure of the source and different spectra of the components on the basis of their brightness in the VLBI maps.