Output options

Whilst running, the program keeps the user informed of how the run is going by outputting an information line at the end of each iteration. An example would be

Iteration number 4 is complete. 16 shells converged.( 107 levels didnot)

The iteration number is shown together with the number of shells that have converged. For a shell to converge the change in the source function from one iteration to the next for each level must satisfy equation 4.4. For overall convergence all shells are required to converge in two consecutive iterations. The number shown in brackets is the total number of levels in all shells that did not converge (eg. if a model has 25 shells and is modelling 10 levels then this number will start off with 250 levels not converging and reduce to zero when overall convergence occurs).

After convergence the GKS window is opened and the program calculates the output for each transition at all required positions. As soon as each output has been calculated it is plotted in the GKS window. Note that since these outputs occur in the order that the program calculates them the order on the screen is not as might be expected so attention needs to be paid to the titles on each individual graph.

Once the program has completed its run an option to save is presented (this option is only presented if model run has been changed in some way from a previous model - there is no point saving 2 identical copies of one model). If chosen saving takes place without any need for further input from the user.

Finally an option to have further outputs is presented. If no further outputs are needed the program is terminated otherwise the program calls the output subroutine and is thereafter the same as if option 1 had been chosen from the initial menu (ie. View some old runs).

The output subroutine allows output to Xwindows or Postscript. Xwindows displays to GKS windows and immediately proceeds to draw the output. If Postscript is selected, however, a further set of options is displayed:

Which output size do you want?
1. Full Size - landscape.
2. Full Size - portrait.
3. User defined size.

Options 1 and 2 are self-explanatory. However, option 3 gives output as encapsulated postscript and asks for input as to the desired size of the plot.

Both the Xwindows and Postscript options produce 2 pages of output with the first showing the modelled and observed (if the data file was present in the  /Sten/Mcs directory) antenna temperatures at the appropriate velocities for each transition and each position. The second page shows the optical depth at the same velocities as used on the first page, again for each transition and position. When using GKS windows the user is prompted to press return to get the second page.

Finally, the option is presented to produce a standard printout. If anything other than 'y' is selected the subroutine returns to wherever it was called from otherwise a A4 page of postscript is produced containing both of the above mentioned 2 pages of output and also listing all the model parameters used. This is stored as a gks.ps file in the 'Sten' directory.