2008-07-18
Version 2.0.0
Edited to allow number of processors to exceed number of subdomains.
INCLUDE replaced by USE HSL_MP01 and source form changed.
Pointer array components replaced by allocatables.
Private communicator added.

2005-01-18
Version 1.1.0
Bug fixes:
  * VALNAM/RHSNAM allocatable (not pointer)

2004-07-12
Version 1.0.0.
Version numbering added.

2004-01-07
In the loop DO 910 changed
IF (I.EQ.3) INQUIRE (IOLENGTH = LENGTH) RCON
to
IF (I.EQ.3) INQUIRE (IOLENGTH = LENGTH) ICON

2003-03-12
edited diagnostic printing (ICNTL(4)).

2003-01-24
edited code to use ADJUSTL and LEN_TRIM (allows blanks within filenames)
on final call, deallocate data%NELTSB on all processes except the host
removed restriction that NDOM must be at least NPROC

2002-05-29
JOB = 5 call ... a couple of arrays were not deallocated and this is now done.
Also moved some deallocate/allocate statements to outside do loops (avoids
  deallocate/allocate within a loop if process has to deal with more than one
  subdomain)

2002-01-25
Reordering is now carried out in parallel to improve timings for JOB = 2 calls.

2002-01-03
Added extra options for entering matrix data (see data%ICNTL(7)).

1999-09
Original Version
This is parallel frontal code for unsymmetric element problems