#######################
# PUBLIC RELEASE 1998 #
#######################
   JET_35171_com.dat  ITER Profile Database:  datasets originally prepared 
   by Pam Stubberfield, Bernard Balet and Knud Thomsen JET, August 1995
   and re-issued in July 1996.

   JET shot 35171 is an ELMy H-mode shot obtained in the new JET pump divertor
   configuration. This pulse ( 1T, 1MA ) is part of a series of rho* scaling 
   experiments relevant for ITER. Plasma is deuterium.  The plasma is in a single-
   null X-point configuration with the ion gradB drift pointing towards the X-point.
   Re-processed NBI data giving approximately 10% more power in the 140KV/D beams
   was used in the TRANSP simulation.  Some 2.9MW of 80KV/D and 1.9MW of 140kV/D
   neutral beam heating is turned on at 62 secs for 4 secs.

   The dataset files in this subdirectory were prepared from TRANSP run I509 
   taking into account the up-down asymmetry of the plasma, run by B.Balet on 
   23 January 1996 on the IBM JAC-2 RS6000 workstation based on JET shot 35171 
   taken on 18 May 1995. All these files have been prepared according to the
   Standard list of Variables and File Format manual by D.Boucher (7 December 1996).
   Tools to produce these files from a TRANSP run were recently updated by Stan Kaye
   of PPPL.  These have only been modified slightly to run on the RS6000 UNIX
   workstation, to add 40 secs to all time axes and to include elongation and
   triangularity.  There are 23 1-D Ufiles and 31 2-D Ufiles. Variables GRHO1 and
   GRHO2 are now stored on flux boundaries and 2-D profiles for CURTOT, NM1(D), 
   NIMP, VROT, BPOL and QOHM are provided. All beam related quantities (energy,
   particles sources, currents etc) have been smoothed to suppress Monte-Carlo noise.

   In the file JET_35171_0d.dat data is given at 65.87 secs.  At that time, there is
   4.8MW of beam heating and the plasma is in a steady ELMy H-mode phase. 

   Values for kappa and delta in the 0d, 1d and 2d files are from the TRANSP
   output derived from the boundary input definition provided by EFIT.
        kappa = (Zmax-Zmin)/(Rmax-Rmin):   delta = (RGEO - R(Zmax))/a
   XPLIM is from XLOC ; SEPLIM is given by EFIT.
   PLTH is defined by WTH/TAUTH coming directly from TRANSP output and not 
   as defined in the standard list of variables. The 1d files for IP,VSURF 
   and IBOOT have negative values as normal for JET. 
  

   Some comments on this particular TRANSP run I509:
   50ms analysis time resolution: VMEC equilibrium LEVGEO = 6  
   New up-down asymmetric boundary definition obtained from the improved EFIT program
   Sawteeth events are not taken into account in the simulation for 
   redistribution of the fast ions and of the current.
   The Ti profile is assumed to be the same as for Te which is from the LIDAR 
   measurements mapped on the outer half of the plasma midplane. 
   TRANSP simulation using this gave a poor fit to the measured neutrons.
   Ne profile vs. radius is from the LIDAR measurements.
   Radiation is taken into account in the analysis. 
   Rotation velocity profile from CX measurements is provided.
   As there was no VB data available, Zeff was assumed constant in time
   and space at 1.1 which gave a good fit to the measured neutron data. 
   Impurity species assumed to be carbon
   Edge ptcl confinement time is set to 1s.
   Bootstrap current on ; impose q(0) = 1.1 to start
   Neoclassical beam driven current
 
   The TRANSP simulation results in a predicted W_dia lower than the measured 
   one by ~0.1MJ.


   Publication:
   Balet B., Campbell D., Christiansen J.P. et al. 22nd EPS Bournemouth UK 1995