Nature of physical problem Program obtainablefrom: CPC Program Library, Queen's Program TSP simulates the phenomena produced in a solid University of Belfast, N. Ireland (see application form in dipolax sample submitted to thermal cycles while in an this issue) external electric field. Dipolar orientation and thermal randomization are the two competing effects, the relative im-Computer: MODCOMP II, or IBM 370/158;Installation: Uni-portance of which varies with the field, the actual temperaversity of Liege ture and the thermal history of the sample (i.e. the previous thermal cycles). The polarization, and the displacement cur-Operating system: MAX!! Release ULG 0.4.1 (or VS2-R01.6) rent showing changes in the polarization, are computed by TSP to support theoretical predictions and to help experi-Program language used: FORTRAN IV mental interpretation. High speed storage required: 11k bytes on Modcomp II or 64k Method of solution bytes on IBM 370/158 The differential equation of the detailed balance is integrated to compute the polarization. The current is obtained No. of bits in a byte: 8 by taking the derivative of the polarization at the time. Overlay structure: none Typical running time Without plot feature for short test data MODCOMP II: 102 s, No. of magnetic tapes required: optional or IBM 370/158: 14 s. Other peripheralsused: papertapeor card-reader, lineprinter, Restrictions on the complexity of the problem 1 or 2 sequential scratch files The number of successive thermal cycles is limited to 6. Unless sizes of 4 vectors are changed, the job must not simulate No. of cards in combined program and test deck: 497 thermal cycles extended over more than 600 seconds. Card punching code: EBCDIC Unusual features of the program None in the Fortran statements. On Modcomp II, we use a