Major changes in this version

Major code changes have been undertaken to make this J-Solver version portable across UNIX platforms (LINUX PCs, workstations and CRAYs), dynamically allocatable, user friendly and avoid the use of proprietary software. J-Solver has been converted to fortran 90 syntax and makes use of IMPLICIT NONE statements. Dimensions are computed on the fly, which simplifies the user interface. It also greatly reduces the average memory requirement. In the process, however, some previously available features have been lost or abandoned:

• NCAR graphics. All calls to NCAR graphics have been commented out but can be switched on by executing the commands (in Bash)

  for f in `ls *.f`; do 
  cat $f | perl -ne 's/ctv80\/ncarg//o; print;' > foo
  mv $f $f_old
  mv foo $f
  done

  Rather than mingling graphics with number crunching code, the preferred approach is to custom write routines that dump data into a file for postprocessing. Examples of such routines are provided in the drive.F file.

• Restart option from dump file (IRST1=1 or IRST2=1) removed. This feature has been superseded by the warm-start capability.
• Radial and poloidal grid sizes have been merged. Previous versions allowed to specify independent radial and poloidal grid sizes through the exponent MNTHE and MNPSI. The grid size is now uniquely defined by MN = MAX(MNTHE, MNPSI), this to overcome a bug which caused the code to crash when MNTHE differed from MNPSI.
• Previous versions of J-Solver produced the fort.17 file. This file is now replaced by eqdsk.cdf file, which is NetCDF portable.