CHEASE code

CHEASE is an open access fixed-boundary equilibrium code. It has been published in:

H. Luetjens, A. Bondeson, O. Sauter, The CHEASE code for toroidal MHD equilibria, Comput. Phys. Commun. 97 (1996) 219 (with updated comments)

The other part of the documentation is the 1992 CPC paper. There are many new features since 1996 publication. A new documentation will appear at some point, but you can contact Olivier Sauter in the meantime.

You can use the code at your wish, but we ask you to fill in the following agreement: CHEASE_pdf and send it to

CHEASE uses COCOS=2 internally, but has COCOS_in and COCOS_out namelist variable to match to any COCOS index (see below)

CHEASE is developed under svn and we publish below the various tagged version. For further details and questions, please ask the contact person:

You can get the last version from svn crpp server:

svn co https://spcsvn.epfl.ch/repos/CHEASE/chease/trunk chease

or

svn export https://spcsvn.epfl.ch/repos/CHEASE/chease/trunk chease

A new file explaining how to install and run the make commands has been added: Make_README

There is a matlab GUI within the above release (in matlab directory) which is useful to get used to running CHEASE in chease/matlab/CHEASEgui (see demos below)

To use the CHEASEgui you should have the executable "chease" in your path as well as the files o.chease_to_cols, o.chease_rz_plasma_boundary_to_cols.pl and o.chease_to_cols_v10.pl (which are in the main CHEASE/trunk/scripts_for_bin directory).

Typically you would have:

PATH=$PATH:$HOME/bin # in your .profile, .cshrc or .bashrc file (or equivalent)

cd $HOME/bin

ln –s .../chease/src-f90/chease .

ln –s .../chease/scripts_for_bin/o.chease_to_cols .

ln –s .../chease/scripts_for_bin/o.chease_rz_plasma_boundary_to_cols.pl .

ln –s .../chease/scripts_for_bin/o.chease_to_cols_v10.pl .

matlab -nodesktop

>> addpath .../chease/matlab/CHEASEgui

>> CHEASEgui

>>(then click on "plot and save" for the plasma boundary, the pprime profile and the TTprime or current profile, then choose suffix (default run_name) then click "prepare CHEASE input", then in menu "Run CHEASE")

You can create a file containing all the output 1D profiles with (from the chease directory):

chease > o.cheaseoutput

o.chease_to_cols o.cheaseoutput o.cheaseoutput.cols

matlab

>> addpath matlab/CHEASEgui

>> plotdatafile(' o.cheaseoutput.cols') % which will open a GUI to plot profiles

Useful Notes

· Notes on the CHEASE normalizations

· COCOS paper defining uniquely the multiple options for Coordinate COnventionS (by O. Sauter et S. Y Medvedev)

· CPO structure. CHEASE now uses the equilibrium structure originally proposed by the EU-ITM Task Force: CPO paper

Using CamStudio, the following demo are available (.avi large files)

· Demo 1 (42MB) for getting started with CHEASEgui

· Demo 2 (20MB) to easily specify a reverse q profile

New option for CHEASE with q profile in input: NSTTP = 5

You can test this option using:

make test_chease_q

Then looking at the files in /tmp/$USER/nsttp_tests

Related (old) codes

 (Note: Now the PROFFIT and BNDFIT options are included in the new CHEASE version with the new namelist variables: TENSPROF and TENSBND using the interpos routines )

We also developed two fortran codes which fit and smooth the plasma boundary (BNDFIT) and the profiles (PROFFIT). These are necessary usually when using experimental data for example.

BNDFIT version 1.1

Or: svn co http://spcsvn.epfl.ch/repos/CHEASE/BNDFIT/trunk/ BNDFIT

PROFFIT version 1.1

Or: svn co http://spcsvn.epfl.ch/repos/CHEASE/PROFFIT/trunk/ PROFFIT