The subroutines save and restor can be used to save the state of common blocks immediately after a return from subroutine mgghat, and restore a saved state, respectively.
In Version 1.1, save is useful for programs where more than one elliptic PDE is solved, and the solution of one is used in defining the other (see the `system of equations' and `time dependent' examples in section 11). The solution contained in the saved state can be evaluated with the function ssolut (see section 7.2).
In Version 1.1 of MGGHAT , restor is not very useful since MGGHAT reinitializes the data structures on every call. It is intended to be used in a later version that will have additional capability.
subroutine save(rwrk,iwrk,lwrk) copies all data structures in the common blocks into the three arrays real rwrk, integer iwrk, and logical lwrk. These arrays should be dimensioned with the parameters ndrwrk, ndiwrk, and ndlwrk, respectively. These parameters are in `commons'.