Test run contributed by Irina Voitsekhocitch :-
! RF assisted current ramp up
!---- JETtransp Update Tue May 19 17:17:11 BST 2009 ----------------------
!HOST = 'jac'
!UFILE = 1
FTIME = 17.0
KMDSPLUS = 0
NSHOT = 77694
TINIT = 7.0
! End JETtransp Update Tue May 19 17:17:11 BST 2009 ----------------------
NRIP=0 ! ripple loss =2 (Redi model), =0 off
ASRD=1.0
BSRD=1.0
! NMIX_KDSAW = 1 !standard Kadomtsev
! FPORCELLI=1.0 !q=1 island extends all the way to the magnetic axis
!and the mixing region is as big as it is for the Kadomtsev model;
NLSAW=.F
NLSAWB=.F
NLSAWE=.T
NLSAWI=.T
NLSAWIC=.F
XSWID1=0.0 !0.0=default, 1.0 turns off Kadomtsev model for T_i
XSWID2=0.0 !0.0=default, 1.0 turns off Kadomtsev model for current
DTSAWD=0.002
NZONES= 50 ! No. plasma zones
NZONE_NB=50 ! No. plasma zones
NZONE_FP=50 ! No. plasma zones
!
NLBCCW=.F ! orientation of B (.TRUE. for counter-clockwise)
NLJCCW=.F ! orientation of Ip (.TRUE. for counter-clockwise)
!
NHECFT=2 ! flag & harmonic number for ECE Te(freq,t) data
LEVTRK=2 ! switch is for the straight-line tracking code
ERRTKR=1.0E-3 ! tracking error (default=1.0E-3)
LEVGEO=6 ! 5 for symetric ; 6 for assymetric geometry
!VMCFTOLA(1)=-1.0e-9 ! weaken VMEC convergence test
DTMINT=1.0E-6
NLINTK=.T
LEVGB=3
DTMAXG=0.025 ! Maximum allowed geometry time step
NFPDEPMD=2 ! Beam deposition model 2=complete DEP 1=faster DEPO
NSIGEXC=1 ! 1 for excited states deposition model
NMSIGX=2 ! Beam deposition on impurity: 1 (default)=Olson, 2=Phaneuf
NDEPMOD=0 ! 0 (default), 1 for enhanced beam dep cross sections
XDEPMOD=1.6 ! enhancement factor if NDEPMOD=1
!
! PLASMA COMPOSITION:
AIMP=12.0
XZIMP=6.0
APLASM=2.0 !* SPECIES ATOMIC WTS.
BACKZ= 1.0
XZEFFI=2.0 !* PLASMA COMPOSITION ZEFF
NLZEFM=.FALSE. !* SET .TRUE. FOR PC ZEFF = MAGDIF ZEFF (nototherwayround)
CZEFFM=1.0 !*ZEFFM = ZEFFPC/CZEFFM
NLZFIN=.TRUE. !* SET .TRUE. TO READ ZEFF(T) DATA
NLZFI2=.FALSE. !* SET .TRUE. TO READ ZEFF(R,T)
NZEFMOD=2,2,2
TZEFMOD=10.27,19.43
DTZEFMOD=0.02
NLZFIM=.FALSE.
NLVISB=.FALSE. !* SET .TRUE. TO READ VIS. BREMS. DATA
NLZVBR=.FALSE. !* SET .TRUE. TO USE VIS. BREMS. DATA FOR ZEFF
!
! SYSTEM DIMENSIONS:
RMAJOR=0.0 !MAJOR RADIUS OF PLASMA SUPPLIED BY BDYMOM
RMINOR=0.0 !MINOR RADIUS OF PLASMA SUPPLIED BY BDYMOM
!
! TRANSP Time Limits
!
! TINIT = 17.00 JETtransp Fri Apr 23 11:08:36
! FTIME = 26.80 JETtransp Fri Apr 23 11:08:36
!
! INITIAL CONDITIONS
NG=1 !*INITIALLY ALL D
FRAC =1.0
GFRAC=1.0
RFRAC=1.0
NGMAX=1 !* MAX NO. OF HYDROGENIC SPECIES = 1 **
!
! OUTPUT
SEDIT= .025 !time step for output of profiles
STEDIT=.025 !time step for output of scalars
!
! FOR OUTPUTTING FAST ION D.F....(alphas & beam ions)
!SELOUT='XXX' !list of TRANSP common symbols to output
SELAVG='FBM BMVOL BDENS2 EBA2PL EBA2PP' !list of names to output
AVGTIM=1.000 !averaging time for SELAVG
MTHDAVG=2 !2 for averaging after each M.C. timestep
!OUTTIM=04.0,06.0,08.0,10.0 ! AC Output times (at most 5)
!
! PTCL BALANCE MODEL
NMODEL=4 !USE (1) CONST. FRAC (2) CONST V OR (3) CONST D PTCL MODEL
NDIFFI=1 ! CONST D INDEPENDANT OF T OR SPACE ,INPUT AS NAMELIST
DIFFUS=0.8E+4 ! VALUE OF CONSTANT DIFFUSIVITY
DFIMIN=0.0 ! MIN ALLOWED DIFFUSIVITY
DFIMAX=1.0E+5 ! MAX ALLOWED DIFFUSIVITY
DIFAC=1.0,1.0 ! MULTIPLIER ON D FOR EACH PLASMA SPECIES
VIFAC=1.0,10.0 ! MULTIPLIER ON FLOW VELOCITY
TAUPH=3*0.6 !* PTCL CONFINEMENT TIME
TAUPO=.1 !* IMPURITY CONFINEMENT TIME
NLRCYC=.F !T get taup from wall recycling
NLTAUP=.F !T=input Tau(P) ion vs. time
!
! ELECTRON CONDUCTIVITY MODEL
!NKEMOD=0,-2 !0 FOR ANALYSIS, -2 to use input KE2 to calculate T_e
NLXEIN=.F !T to read KE2 ufile
TKEMOD=3.05 !transition time
XKEMAX=1.0E+06 !Max bound on Chi_e
XKEMIN=0.01 !Min bound on Chi_e
!
! ION CONDUCTIVITY MODEL
NLVPHI=.T !T for rotating plasma
XVPHI=2.0 !Parabolic rot profile squared
NLTI2=.T !T to use Ti profile data
NLTI2TX=.T !T to use data as impurity Ti
NLNTMJ=.T !T to use TMJ for neutron rates
NKIMOD=4 !* SELECT KAPA(I) N.C. FIT MODEL
XKFAC=2.0 !* ANOMOLOUS MULTIPLIER FOR Chi_e
TKIMOD= 8.94,15.993 !* End times for NKIMODA
NKIMODA=5,100,5 !Use Chi_i=XKIMOD*Chi_CH until time TKIMOD(1)
XKIMOD=5.0,5.0 !* ANOMOLOUS MULTIPLIER FOR Chi_e period NKIMODA>0
NLFXKF=.F !*** ADJUST XKFAC TO FIT MEASURED TI IF .T
NLXKIE=.F !.T to set Chi_i=XKFAC*Chi_e
DLTKIE=2.0 ! to reduce Chi_i fluct
ALPH0I=.6 ! 3/2 convection
ALPH0E=.6 ! 3/2 convection
NLNTX=.T ! .T to read neutrons
NLTKA=.F ! .T FOR FEEDBACK ON K-ALPHA ION TEMPERATURE
NLTNTX=.F
XKFMAX=1000.
XKFMIN=0.001
!
!----- Ti analysis mode (single model)
!NLTIPRO=.T ! turn off Ti prediction
GIEFAC=1.0 ! Ti/Te when Ti data not avail. (default=1.0)
FIEFAC=1.0 ! Ti = (1-FIEFAC)*(GIEFAC*Te) + FIEFAC*Ti (default=0.0)
DTISAVE=0.01 ! delay for predictive/analytical model switchover
TIXLIM=1.0 ! max xi (r/a) of valid Ti data (default=1.0)
XTILIM=1.0 ! limit r/a of validity of Ti profile
TIFACX=1.0 ! edge Ti/Te factor (default=0.0)
TIDXSW=0.05 ! switch to using TIFACX*Te data over the xi range
!
!----- NCLASS controls
NCMODEL=2
NLOMGVTR=.T ! .T: rotation data inputs OMG or VP2 refer
NLBOOTW=.T ! .T: use NCLASS bootstrap current
XL1NCJBS=0.1 ! r/a for bootstrap smoothing; default XL1NCJBS=0.1
NLETAW=.T ! .T: use NCLASS resistivity
XL1NCETA=0.0 ! r/a for eta smoothing; default XL1NCJBS=0.0
NLVWNC=.T ! .T: compute NCLASS radial electrostatic potential profile
NLSQUEEZ=.F ! .T for "squeezing factor" correction
NLTINC=.T ! .T to use Tx/Tmj, i.e. impurity temperature
NLNCSMOO=.T ! .T to activate smoothing of temp & density
TAUSMNNC=0.050 ! smoothing convolution time (tau) for densities
TAUSMTNC=0.050 ! smoothing convolution time (tau) for temperatures
XSMN_NC=0.050 ! radial smoothing half width for densities (r/a units)
XSMT_NC=0.050 ! radial smoothing half width for temperatures (r/a units)
NLSAW_NC=.T ! .T to reinitialize time convolutions at each sawtooth
!
! FUSION PRODUCTS MODEL
NPTCLF=50000 !Monte Carlo particles for fusion ion species
NLFHE4=.T. !.T for DT->HE4
PLFHE4=5.0E3 !min energy (watts) for full NPTCLF
NLFHE3=.T !.T for DD->He3
PLFHE3=1.0E4 !energy (watts)
NLFST=.T !.T for DD->T
PLFST=1.0E4 !energy (watts)
NLBOUT=.T !.F for fast ion loss file
NLFATOM=.T ! atomic effects on/of fusion products ions (CX loss)
!
! MAGNETICS MODELING
NLQLIM0=.T
QLIM0=9. ! q0-limit as max(QLIM0,2*qmin),
NLSPIZ=.F !* .T TO USE SPITZER INSTEAD OF N.C. RESISTIVITY
NLMDIF=.T ! .T TO SOLVE POLOIDAL FIELD DIFFUSION EQN.
NLBOOT=.T !* =.T TO INCLUDE BOOTSTRAP CURRENTS
NMCURB=3 ! beam current model 0=none 1=classical 3=neoclassical
NLI2PB=.T !* COMPARE LI/2+BETA MEASURED DATA
NLVSUR=.F
NLALAM=.F
NLBDIA=.T !* =.T to read in diamagnetic beta
NLBPDA=.F
NLDFLX=.T !* =.T to read in diamagnetic flux
NEFLD=3
VLPMIN=0.1
VLPMML=0.05
!XPEFLD=2.0
XZFMIN=0.1
!QEFLD=4.0
!RQEFLD=0.0
NLQDATA=.T
NLPCUR=.T ! to match current profile using edge current
NQMODA=1,4,1 ! 4 for QPR input
NQMODB=1,1,1
TQMODA=2.0,8.0
TAUQMOD=0.1
!
! NEUTRALS MODEL - SEE ALSO PTCL BALANCE
NSOMOD=1 ! SELECT 1 FOR FRANTIC, 8 FOR SORCE8
MOD0ED=1 ! =1 TO SET T0(EDGE)=TIEDGE , =3 T0(EDGE) = TE(EDGE)
MODIEDG=3 ! =1 TO SET TI(EDGE)=TIEDGE , =3 TI(EDGE) = TE(EDGE)
TI0FRC=.033333 ! T0(A)=TI0FRC*TI(0) IF MOD0ED=2
E0IN=10., 2*3., 10., 2*3. ! T0(A) FOR COLD SCES, HOT SCES IF MOD0ED=1
FH0ESC=.3 ! FRAC. OF ESC. NEUTRALS REFLECTED TO PLASMA "HOT"
!
! EFLUX MODEL
NLEFLX=.F !* .T FOR NEUTRAL EFLUX CALCULATION
!
!================================================================
! ICH Model
!================================================================
!
NLICRF=.F
NICRF=1 !ICRF model switch (1=new SPRUCE; 5=old SPRUCE)
DTICRF=0.050 !HOW OFTEN TO CALL ICRF HEATING PACKAGE
XZMINI=1.0 !* MINORITY CHARGE STATE
AMINI= 1.0 !* ATOMIC WT OF MINORITY
FRMINI=0.058 !* MINORITY CONCENTRATION nmin/ne (if not set by a ufile)
!FRACMINI=0.10,0.90 !* MINORITY relative fractions
TAUMIN=0.2 ! Minority ptcl conf. time (not known)
!
! namelist for TORIC
NMDTORIC=63 ! N of poloidal modes: Npol=2^n-1, n is selected, Npol is calculated
RFARTR=5.0 ! distance from antenna to Faraday shield, cm
ANTLCTR=1.6 ! effective antenna propagation constant
NFLRTR=1 ! ion FLR contributions, =1 included, =0 ignored
! NFLRETR=1 ! electron FLR contribution, NOT USED IN ROBERT's run
! FLRFACTR=1.0 ! not used in Robert's run
NBPOLTR=1 ! poloidal magn. field, =1 included, =0 ignored
NQTORTR=1 ! toroidal broadenning of plasma dispersion
NCOLLTR=0 ! collisional contribution (=0 in RB namelist)
ENHCOLTR=1.0 ! electron collision enhancement factor
! ALFVNTR(20) ad hoc collisional broadenning of Alfven and ion-ion resonance
ALFVNTR(1)=0.0 ! =1 included, =0 ignored
ALFVNTR(2)=0.1 ! enhancement factor
ALFVNTR(3)=3.0 ! value of ABS((n//^2-S)/R) below which damping is added
ALFVNTR(4)=5.0 ! value of abs(w/(k//*v_te)) below which damping is calculated
!
NICHPSI=128 !* Number of radial grid points default=128
NICHCHI=64 !* Number of poloidal grid points default=64
ERFMAX=10.E6 !* Max grid energy (eV) default=10MeV
!
MSYM_NPHI=1 ! 1 for symmetric n_|| spectra
NUM_NPHI =1 ! Number of n_||
NNPHI(1,1)=+27 ! sets n_|| values
WNPHI(1,1)=1.0 ! weights for n_|| values
!
NICHA=4 !* Number of ICH antennas must agree with RFLIST
NGEOANT=13 !* 1 to specify antennae geom w RMNICHA, RMJICHA circles
RGEOANT=380.06,382.87,385.18,386.98,388.27,389.04,389.30,389.04,388.27,
RGEOANT(10)=386.98,385.18,382.87,380.06
YGEOANT=89.912,80.439,70.833,61.119,51.325,41.476,31.600,21.724,11.875,
YGEOANT(10)=2.0808,-7.6329,-17.239,-26.712
RMJICHA(1)=4*200.6 !* Major radius of each antenna (cm)
RMNICHA(1)=4*188.7 !* Minor radius of antenna (<0 for high field ant.)
THICHA(1)=4*36.0 !* Poloidal extent of each antenna (degrees)
WIDICHA(1)=4*16.0 !* Toroidal width of each antenna element (cm)
SEPICHA(1)=4*45.6 !* Toroidal Separation of antenna elements (cm)
!
PHICHA(1,1)=-90.0,+90.0,-90.0,+90.0
PHICHA(1,2)=-90.0,+90.0,-90.0,+90.0
PHICHA(1,3)=-90.0,+90.0,-90.0,+90.0
PHICHA(1,4)= 0.0, 0.0,-90.0,+90.0
VVRMOM(1)=292.8,134.5,10.4
VVZMOM(1)=0.0,209.6,-16.2
TCRFON= 3.90696716
TCRFOF= 9.86857605
!FI_OUTTIM(1)=5.
!FI_OUTTIM(2)=8.
!
NLLH=.F ! .T for LH model !TIMLSOUT=6.50,7.50,8.50,9.50 ! times for cd.dat files NYXINV=100 ! No of Y pts in 2D map NCUPLRLH=6,6 ! 6 for 'JET' NDOBRAM=1 ! 1 do Brambilla calc w JEStevens code (1) PHASEDLH(1)=0.0, ! phasing in deg between waveguids if NGRPS<1 NANTLH=1, ! Only 1 antenna modeled so far NGRPSLH=1, ! num of spectrun peaks: num of waveguids if DO>0 (3) !TOTPWRLH=1.9E+06, ! LH power (WATTS) used by TRANSP !TLHON=8.00, ! LH start time !TLHOFF=10.5, ! LH stop time DTLH=0.01 ! Frequency of calls to LH calculation (sec) ! $inpval TBT 1/30/92 from IGNAT$:[TRANSP]INPUT.lhh NVLH=399, ! number of v_el/c bins points btw -1 and +1 (199) NSLICELH=401, ! num of n_par slices used in Brambilla calc (301) WGTITRLH=0.2, NGTYPELH=3, ! 1 for linear grid, 2 or 3 for expon grid, 3 is smoother NFREQLH=50, ! steps btwn preparing ray plot in 3d (100) NSTEPLH=20000 ! max steps in following each ray (500) NPSILH=80, ! numr of psi shells for power/current dep (10) NZONESLH=2000, ! number of shell crossings tracked (100) NSMOOLH=5, ! num vel bins for smoothing D_ql (9) NSMWLH=3, ! characteristic width of fun covering NSMOO vs bins (3) FGHZLH=3.7, ! frequency in GHz (3.7 for JET) HLH=0.010, ! step length in m along pathlength (0.005 for PBXM) PARMINLH=1.5, ! left edge of n_parallel peak PARMAXLH=2.3, ! right edge of n parallel peak CENTERLH(1)=1.8, !center of n_parallel peak WIDTHSLH(1)=0.25, !array (size NGRPS) giving width of n_par peak (1.) POWERSLH(1)=1.00, ! relative LH power (usually 1.0) !irememfe=0, ! default=0 NPLFLGLH( 1)=0,0,0,0,0,0,0,0,0,0, ! plot controls NPLFLGLH(11)=0,0,0,0,0,0,0,0,0,0, NPLFLGLH(21)=0,0,0,0,0, NPRPFGLH( 1)=0,0,0,0,0,0,0,0,0,0, NPRPFGLH(11)=0,0,0,0,0,0,0,0,0,0, NPRPFGLH(21)=0,0,0,0,0,0,0,0,0,0, NDORFDLH=1, !deposition grafs NDORFHLH=0, NDOPASLH=0, NDOSCRLH=1, ! screen output NDODBGLH=0, ! debug? NRAMPULH=100, ! number of steps to ramp up power NFLATLH=101, ! number of times in each step NDIAGLH=10, NRAYSLH=15, ! num of rays in peak (10) ! ! ***************************************** ! Pellet data ! ***************************************** ! NPEL=0 ! Number of pellets TPEL=4.4999,4.7506,4.9996 ! List of pellet popping times PELVEL=1.289E+5,1.169E+5,1.374E+5 ! List of pellet velocities PELRAD=0.27,0.27,0.4 ! List of pellet radii APEL=2.0,2.0,2.0 ! List of atomic weights of pellet material ! NMPELA=1 ! 1=std model +ablation, 2=ablation only LPLPRI=3 ! ablation printout detail, larger n => more detail PLTEDG=1.0 ! edge region Te for pellet ablation PLDEDG=1.E11 ! edge region ne for pellet ablation PLRSTA=700.0 ! R start of pellet PLYSTA=0.0 ! Y start of pellet PLTHEA=0.0 ! Poloidal tilt angle of pellet gun (degrees) PLPHIA=0.00 ! Toroidal tilt angle of pellet gun (degrees) ! ! DATA HANDLING ! ! NSHOT = 77694 JETtransp Fri Apr 23 11:08:36 TGRID1=0.020 ! UPF F(T) TIME GRID DT TGRID2=0.020 ! UPF F(X,T) TIME GRID DT PRFAC=0.1 !* FRACTION OF CX EFLUX TOT POWER TO SUBTRACT FROM BOLO. ! NRITER=1 NSYTER=0 ! Set to 1.0 for slice and stack if NRITER=3 NRINER=1 NSYNER=0 ! Set to 1.0 for slice and stack if NRINER=3 NSYECF=0 NRIECF=1 NRITI2=1 ! Ti(R) map 1= input on major radius to right of axis only NRIOMG=1 NRIVP2=1 NRIZF2=1 NSYZF2=0 NRIQPR=-6 NRIBOL=-6 FNEMIN=1.0E8 ! ! MISCELLANEOUS NLEBAL=.T ! .T FOR ELECTRON ENERGY BALANCE CALCULATION NLTIBL=.T ! .T FOR TI BALANCE CALCULATION NLPBAL=.T ! .T FOR PTCL BALANCE CALCULATION ! ! ***************************************** ! Diagnostic simulation parameters ! ***************************************** ! ! Interferometer and Faraday rotation array ! geometry of KG4 pulse #49382 ! R=1.892,2.703,3.042,3.744,3.000,3.000,3.000,3.000 ! Z=0.000,0.000,0.000,0.000,-0.614,-0.376,-0.147,0.190 ! A=1.5708,1.5694,-1.5677,1.5708,0.3901,0.1681,-0.0389,-0.2437 ! Y.Baranov converted to TRANSP variables 14 Jan 2000 NLDA=8 PHLDA=8*0.00 RLDA=189.2,270.3,304.2,374.4,400.0,400.0,400.0,400.0 YLDA=+300.0,+300.0,+300.0,+300.0,-20.28,-20.63,-18.59,-5.86 THLDA=-90.00,-90.00,-90.00,-90.00,-22.35,-9.63,+2.228,+13.960 ! ! VB array ! Geometry from P. Morgan, 9 Sept 87 ! geometry correcte 2 August 1990 by PMS as in run 4078 !NVBA = 15 !PHVBA = 15*0. !THVBA = 66.0,69.4,72.0,74.4,81.1,79.8,84.6,86.8,85.1,92.1,90.9, !THVBA(12) = 98.3,103.2,90.0,0.0 !RVBA = 324.0,324.0,324.0,324.0,288.0,324.0,288.0,288.0, !RVBA(9) = 324.0,288.0,324.0,288.0,263.0,310.0,450.0 !YVBA = 14*-367.0,-29.0 ! ! Geometry for 1991 diagnostic 10 sept 91 ! VBINT(10) corresponds to S3-AD16 used for HAL/EFZ ! !NVBA = 10 !PHVBA = 10*0. !THVBA = -67.6,-69.9,-72.8,-75.8,-84.4,-86.3,-94.9,-98.4,-107.8,-90.0 !RVBA = 6*331.49,3*279.76,311.0 !YVBA = 10*351.65 ! ! Geometry for 1996 diagnostic / 22 Aug 1996 from P.Morgan ! VBINT(1) corresponds to KS3/BASH used for KS3/ZEFH ! VBINT(2) corresponds to KS3/BASV used for KS3/ZEFV ! NVBA = 2 PHVBA = 2*0. THVBA = -2.00522,90.0 RVBA = 590.08,311.0 YVBA = 2.15,-300. ! ! Collimated Neutron Array ! Geometry from G. Sadler 15 Sept 87 ! only good for '86 and winter '87 runs NNTA = 19 PHNTA = 19*0. THNTA = 18.62,14.68,10.6,6.41,2.14,-2.14,-6.41,-10.6,-14.68,-18.62, THNTA(11) = -77.58,-80.62,-83.72,-86.85,-90.0,-93.14,-96.28, THNTA(18) = -99.38,-102.42 RNTA = 10*605.5, 5*302.15, 4*302.15 ! focal point YNTA corrected 12/07/88 YNTA = 10*0.,9*344.2 ! ! CX detector array ! ! Crystal Spectrometer ! Geometry from E. Kallne 9 Sept 87 NXCA = 1 RXCA = 4.0 YXCA = 0. THXCA = 0. PHXCA = 27.065 ! ********************************************* ! TRDAT NAMELIST ! ********************************************* ! $TRDATA ! INPUTDIR='transp_data/77694/61097C03/data' NLRNNE = .FALSE. ! .T TO RENORM NE(R) TO LINE AVG NE in TRDAT LFIXUP=2 XRCTER=0.05 XRCNER=0.05 XRCECF=-1.0 XRCTI2=-0.1 !XRCKI2=0.5 !XRCKE2=0.5 !XRCVP2=-0.1 XRCOMG=-0.1 XRCZF2=-0.1 !TPELDA=11.900,12.100!list of bad data times around each pellet ! The TPELDA times must be at least 2*max(tgrid1,trgrid2) away ! from the pellet time TPEL !TPELDA(1,1)=4.4 ! TIME0 = 40.0 !JET plasma's don't begin at TIME0=0.0! ! PREQPR='EFIT' EXTQPR='Q_0030 PPF:2:tranppf:0030:$SDD:Q_PROFILE(PSI_NORM) :' PRECUR='MG2' EXTCUR='XIP_0063 PPF:1:JETPPF:0063:$SDD:PLASMA CURRENT :' PREDFL='MAGN' EXTDFL='FLX_0003 PPF:1:JETPPF:0003:$SDD:Diamagnetic flux :' PREVSF='MG3' EXTVSF='VPA_0063 PPF:1:JETPPF:0063:$SDD:VOLTAGE PARALLEL TO B :' PRETER='HRTS' EXTTER='TE_0045 PPF:2:JETPPF:0045:$SDD:TEMPERATURE :' PRENER='HRTS' EXTNER='NE_0045 PPF:2:JETPPF:0045:$SDD:DENSITY :' PREL2B='MG2' EXTL2B='XLM_0063 PPF:1:JETPPF:0063:$SDD:1/2 LI + BETAI :' PREMRY='MSCR' EXTMRY='MRY PPF:0:tranppf:0029:$SDD: :' PRERBZ='MAGN' EXTRBZ='ITOR_0003 PPF:1:JETPPF:0003:$SDD:Torroidal field current :' PREBOL='BOLO' EXTBOL='TOBU_0056 PPF:1:JETPPF:0056:$SDD:Tot. bulk rad. (vert.) :' PRENTX='TIN' EXTNTX='RNT_0024 PPF:1:JETPPF:0024:$SDD:TOTAL NEUTRON RATE :' PRENB2='PNBI' EXTNB2='NBI PPF:0:tranppf:0028:$SDD: :' !PREZF2='CXFM' !EXTZF2='ZFCX' PREZEF='KS3' EXTZEF='ZEFV_0037 PPF:1:JETPPF:0037:$SDD:<Zeff> vertical (AD16) :' PRETI2='CXFM' EXTTI2='TI_0008 PPF:2:cxsbatch:0008:$TDD:Ion Temperature :' PREBDI='MG2' EXTBDI='BET_0063 PPF:1:JETPPF:0063:$SDD:BETAI (DIAMAGNETIC) :' PREOMG='CXFM' EXTOMG='ANGF_0008 PPF:2:cxsbatch:0008:$TDD:Toroidal rotation :' PRERCY='EDG7' EXTRCY='FLW_0021 PPF:1:JETPPF:0021:$SDD:Total wall D-alp ph flux:' PREGFD='GASM' EXTGFD='MAJR_0004 PPF:1:JETPPF:0004:$SDD:D2 E L RATE :' !PRERFP='PICH' !EXTRFP='ICRH' !PRERFF='FREQ' !EXTRFF='ICRH' !PRESAW='A' !sawtooth correlation function !EXTSAW='SAW' ! ! $END ! ***** Beginning of the file 'JETTRN.TAPE4A.NBFIXD.DAT'. ***** ! ***** End of the file 'JETTRN.TAPE4A.NBFIXD.DAT' ***** ! JET VACUUM VESSEL DESCRIPTION-- TREAT AS LIMITER FOR ! BEAM DEPOSITION AND EFLUX ! ***** End of the file 'jet2000.dat' ***** ! ***** Beginning of the file jet2000.dat. ***** ! %NBLIST.FOR-- BEAM NAMELIST DATA WRITTEN TO END OF FILE ! JET VACUUM VESSEL DESCRIPTION-- TREAT AS LIMITER FOR ! BEAM DEPOSITION AND EFLUX ! ***** End of the file 'jet2000.dat' ***** NPTCLS = 50000 ! no. of MC ptcls ! TrFile (voits, 14-Jan-2009 18:38:07) SC_VSF = -1.0 SC_RBZ = 768.0 SC_RCY = 10.0 ! TrFile (voits, 15-Jan-2009 16:53:17) ! TrFile (voits, 16-Jan-2009 11:24:28) ! TrFile (voits, 16-Jan-2009 17:20:33) ! TrFile (voits, 16-Jan-2009 17:59:45) ! TrFile (voits, 20-Jan-2009 16:06:36) ! TrFile (voits, 27-Jan-2009 13:19:45) ! TrFile (ij, 20-Mar-2009 17:11:39) ! TrFile (ij, 24-Mar-2009 10:45:23) ! TrFile (ij, 26-Mar-2009 11:19:50) ! TrFile (ij, 16-Apr-2009 14:05:16) ! TrFile (ij, 21-Apr-2009 14:15:14) ! TrFile (ij, 22-Apr-2009 11:17:44) ! TrFile (voits, 13-May-2009 16:49:56) ! TrFile (voits, 14-May-2009 10:09:55) ! TrFile (voits, 14-May-2009 15:46:29) ! TrFile (voits, 15-May-2009 11:22:04) !NBLIST.FOR ! beam constants, configuration = JET2000 NLBEAM = .T DTBEAM = 0.025 ! beam code timestep NDEP0 = 500 ! no. of deposition tracks GOOMAX = 5000.0 ! max goosing factor DTN = 1.8E-6 ! orbit timestep control REDGE = 17.0 ! aperture dimensions XZEDGE = 19.0 NBSHAP = 1 ! 1=rectangular source grid FOCLR = 1.E3 ! focal lengths FOCLZ = 1.4E3 DIVR = 1.22E-2 ! divergences DIVZ = 1.22E-2 BMWIDR = 9.0 ! source grid dimensions BMWIDZ = 22.5 DN0OUT = 5.E11 ! "external" neutral density NCIRLM = 2 ! circle limiters CRLMR1 = 480.0, 189.3 CRLMY1 = 28.0, 31.6 CRLMRD = 300.6, 198.9 NLINLM = 3 ! line limiters ALNLMR = 0.0,0.0,421.0 ALNLMY = 215.0,-215.0,0.0 ALNLMT = 0.0,0.0,90.0 ! BEAMLINE P1U %NBLIST.FOR GENERATED OUTPUT NLCO(1) = .T EINJA(1) = 0.000000E+00 PINJA(1) = 0.000000E+00 TBONA(1) = 8.007500E+00 TBOFFA(1) = 1.600150E+01 ABEAMA(1) = 2.000000E+00 XZBEAMA(1) = 1.000000E+00 PDELTA(1) = 1.000000E-03 RTCENA(1) = 1.848000E+02 XLBAPA(1) = 8.300000E+02 XLBTNA(1) = 1.264600E+03 FFULLA(1) = 0.000000E+00 FHALFA(1) = 0.000000E+00 XYBAPA(1) = 2.600000E+01 XYBSCA(1) = 4.590000E+01 XBZETA(1) = 3.297900E+02 ! BEAMLINE P2 %NBLIST.FOR GENERATED OUTPUT NLCO(2) = .T EINJA(2) = 0.000000E+00 PINJA(2) = 0.000000E+00 TBONA(2) = 8.007500E+00 TBOFFA(2) = 1.600150E+01 ABEAMA(2) = 2.000000E+00 XZBEAMA(2) = 1.000000E+00 PDELTA(2) = 1.000000E-03 RTCENA(2) = 1.848000E+02 XLBAPA(2) = 8.205000E+02 XLBTNA(2) = 1.255200E+03 FFULLA(2) = 0.000000E+00 FHALFA(2) = 0.000000E+00 XYBAPA(2) = 1.500000E+01 XYBSCA(2) = 1.461000E+02 XBZETA(2) = 3.297900E+02 ! BEAMLINE P3 %NBLIST.FOR GENERATED OUTPUT NLCO(3) = .T EINJA(3) = 0.000000E+00 PINJA(3) = 0.000000E+00 TBONA(3) = 8.007500E+00 TBOFFA(3) = 1.600150E+01 ABEAMA(3) = 2.000000E+00 XZBEAMA(3) = 1.000000E+00 PDELTA(3) = 1.000000E-03 RTCENA(3) = 1.310000E+02 XLBAPA(3) = 8.205000E+02 XLBTNA(3) = 1.274300E+03 FFULLA(3) = 0.000000E+00 FHALFA(3) = 0.000000E+00 XYBAPA(3) = 1.500000E+01 XYBSCA(3) = 1.461000E+02 XBZETA(3) = 3.253000E+02 ! BEAMLINE P4 %NBLIST.FOR GENERATED OUTPUT NLCO(4) = .T EINJA(4) = 0.000000E+00 PINJA(4) = 0.000000E+00 TBONA(4) = 8.007500E+00 TBOFFA(4) = 1.600150E+01 ABEAMA(4) = 2.000000E+00 XZBEAMA(4) = 1.000000E+00 PDELTA(4) = 1.000000E-03 RTCENA(4) = 1.310000E+02 XLBAPA(4) = 8.300000E+02 XLBTNA(4) = 1.283800E+03 FFULLA(4) = 0.000000E+00 FHALFA(4) = 0.000000E+00 XYBAPA(4) = 6.200000E+00 XYBSCA(4) = 4.840000E+01 XBZETA(4) = 3.253000E+02 ! BEAMLINE P5 %NBLIST.FOR GENERATED OUTPUT NLCO(5) = .T EINJA(5) = 0.000000E+00 PINJA(5) = 0.000000E+00 TBONA(5) = 8.008499E+00 TBOFFA(5) = 1.600150E+01 ABEAMA(5) = 2.000000E+00 XZBEAMA(5) = 1.000000E+00 PDELTA(5) = 1.000000E-03 RTCENA(5) = 1.310000E+02 XLBAPA(5) = 8.299000E+02 XLBTNA(5) = 1.283700E+03 FFULLA(5) = 0.000000E+00 FHALFA(5) = 0.000000E+00 XYBAPA(5) = -6.000000E+00 XYBSCA(5) = -4.840000E+01 XBZETA(5) = 3.253000E+02 ! BEAMLINE P6U %NBLIST.FOR GENERATED OUTPUT NLCO(6) = .T EINJA(6) = 0.000000E+00 PINJA(6) = 0.000000E+00 TBONA(6) = 8.007500E+00 TBOFFA(6) = 1.600150E+01 ABEAMA(6) = 2.000000E+00 XZBEAMA(6) = 1.000000E+00 PDELTA(6) = 1.000000E-03 RTCENA(6) = 1.310000E+02 XLBAPA(6) = 8.190000E+02 XLBTNA(6) = 1.272700E+03 FFULLA(6) = 0.000000E+00 FHALFA(6) = 0.000000E+00 XYBAPA(6) = -3.000000E+00 XYBSCA(6) = -1.474000E+02 XBZETA(6) = 3.253000E+02 ! BEAMLINE P7U %NBLIST.FOR GENERATED OUTPUT NLCO(7) = .T EINJA(7) = 0.000000E+00 PINJA(7) = 0.000000E+00 TBONA(7) = 8.007500E+00 TBOFFA(7) = 1.600150E+01 ABEAMA(7) = 2.000000E+00 XZBEAMA(7) = 1.000000E+00 PDELTA(7) = 1.000000E-03 RTCENA(7) = 1.848000E+02 XLBAPA(7) = 8.190000E+02 XLBTNA(7) = 1.253600E+03 FFULLA(7) = 0.000000E+00 FHALFA(7) = 0.000000E+00 XYBAPA(7) = -3.000000E+00 XYBSCA(7) = -1.474000E+02 XBZETA(7) = 3.297900E+02 ! BEAMLINE P8 %NBLIST.FOR GENERATED OUTPUT NLCO(8) = .T EINJA(8) = 0.000000E+00 PINJA(8) = 0.000000E+00 TBONA(8) = 8.007500E+00 TBOFFA(8) = 1.600150E+01 ABEAMA(8) = 2.000000E+00 XZBEAMA(8) = 1.000000E+00 PDELTA(8) = 1.000000E-03 RTCENA(8) = 1.848000E+02 XLBAPA(8) = 8.299000E+02 XLBTNA(8) = 1.264600E+03 FFULLA(8) = 0.000000E+00 FHALFA(8) = 0.000000E+00 XYBAPA(8) = -6.000000E+00 XYBSCA(8) = -4.840000E+01 XBZETA(8) = 3.297900E+02 ! BEAMLINE Q3 %NBLIST.FOR GENERATED OUTPUT NLCO(9) = .T EINJA(9) = 0.000000E+00 PINJA(9) = 0.000000E+00 TBONA(9) = 8.007500E+00 TBOFFA(9) = 1.600050E+01 ABEAMA(9) = 2.000000E+00 XZBEAMA(9) = 1.000000E+00 PDELTA(9) = 1.000000E-03 RTCENA(9) = 1.310000E+02 XLBAPA(9) = 8.205000E+02 XLBTNA(9) = 1.274300E+03 FFULLA(9) = 0.000000E+00 FHALFA(9) = 0.000000E+00 XYBAPA(9) = 1.500000E+01 XYBSCA(9) = 1.461000E+02 XBZETA(9) = 1.453000E+02 ! BEAMLINE Q4 %NBLIST.FOR GENERATED OUTPUT NLCO(10) = .T EINJA(10) = 0.000000E+00 PINJA(10) = 0.000000E+00 TBONA(10) = 8.007500E+00 TBOFFA(10) = 1.600050E+01 ABEAMA(10) = 2.000000E+00 XZBEAMA(10) = 1.000000E+00 PDELTA(10) = 1.000000E-03 RTCENA(10) = 1.310000E+02 XLBAPA(10) = 8.300000E+02 XLBTNA(10) = 1.283800E+03 FFULLA(10) = 0.000000E+00 FHALFA(10) = 0.000000E+00 XYBAPA(10) = 6.200000E+00 XYBSCA(10) = 4.840000E+01 XBZETA(10) = 1.453000E+02 ! BEAMLINE Q5 %NBLIST.FOR GENERATED OUTPUT NLCO(11) = .T EINJA(11) = 0.000000E+00 PINJA(11) = 0.000000E+00 TBONA(11) = 8.007500E+00 TBOFFA(11) = 1.600050E+01 ABEAMA(11) = 2.000000E+00 XZBEAMA(11) = 1.000000E+00 PDELTA(11) = 1.000000E-03 RTCENA(11) = 1.310000E+02 XLBAPA(11) = 8.299000E+02 XLBTNA(11) = 1.283700E+03 FFULLA(11) = 0.000000E+00 FHALFA(11) = 0.000000E+00 XYBAPA(11) = -6.000000E+00 XYBSCA(11) = -4.840000E+01 XBZETA(11) = 1.453000E+02 ! BEAMLINE Q6 %NBLIST.FOR GENERATED OUTPUT NLCO(12) = .T EINJA(12) = 0.000000E+00 PINJA(12) = 0.000000E+00 TBONA(12) = 8.007500E+00 TBOFFA(12) = 1.600050E+01 ABEAMA(12) = 2.000000E+00 XZBEAMA(12) = 1.000000E+00 PDELTA(12) = 1.000000E-03 RTCENA(12) = 1.310000E+02 XLBAPA(12) = 8.190000E+02 XLBTNA(12) = 1.272700E+03 FFULLA(12) = 0.000000E+00 FHALFA(12) = 0.000000E+00 XYBAPA(12) = 0.000000E+00 XYBSCA(12) = 0.000000E+00 XBZETA(12) = 1.453000E+02 ! BEAMLINE Q7U %NBLIST.FOR GENERATED OUTPUT NLCO(13) = .T EINJA(13) = 0.000000E+00 PINJA(13) = 0.000000E+00 TBONA(13) = 8.007500E+00 TBOFFA(13) = 1.600050E+01 ABEAMA(13) = 2.000000E+00 XZBEAMA(13) = 1.000000E+00 PDELTA(13) = 1.000000E-03 RTCENA(13) = 1.848000E+02 XLBAPA(13) = 8.190000E+02 XLBTNA(13) = 1.253600E+03 FFULLA(13) = 0.000000E+00 FHALFA(13) = 0.000000E+00 XYBAPA(13) = -3.000000E+00 XYBSCA(13) = -1.474000E+02 XBZETA(13) = 1.497900E+02 ! BEAMLINE Q8 %NBLIST.FOR GENERATED OUTPUT NLCO(14) = .T EINJA(14) = 0.000000E+00 PINJA(14) = 0.000000E+00 TBONA(14) = 8.007500E+00 TBOFFA(14) = 1.600050E+01 ABEAMA(14) = 2.000000E+00 XZBEAMA(14) = 1.000000E+00 PDELTA(14) = 1.000000E-03 RTCENA(14) = 1.848000E+02 XLBAPA(14) = 8.299000E+02 XLBTNA(14) = 1.264600E+03 FFULLA(14) = 0.000000E+00 FHALFA(14) = 0.000000E+00 XYBAPA(14) = -6.000000E+00 XYBSCA(14) = -4.840000E+01 XBZETA(14) = 1.497900E+02 NBEAM= 14 ! NBLIST.FOR: NO. OF BEAMS ! NBLIST.FOR: 20% .GT. MAX(INJECTION ENERGY): EBDMAX= 1.409762E+05 NLBDAT=.T ! ** BEAM DATA UFILE FLAG ** ! TrFile (voits, 15-May-2009 11:47:49) ! TrFile (voits, 18-May-2009 17:09:04) ! TrFile (voits, 19-May-2009 17:18:56) ! TrFile (voits, 21-May-2009 18:23:20) |