TRANSP Grid Analysis 98004Q02 JET trdat.log
==>runtrx_r9 start: date: Fri Apr 3 11:17:51 AM EDT 2026 ( transpgrid-srv-r9.pppl.gov )
argv = 2
iarg = 2
cmd_opt = pre
runtrx_r9: tok.yy = JET.20
==========(runtrx_r9)======================
date: Fri Apr 3 11:17:51 AM EDT 2026 ( transpgrid-srv-r9.pppl.gov )
args: 98004Q02 pre
==========(runtrx_r9)======================
==========(runtrx_r9)======================
==========TRANSP preprocessing==========
date: Fri Apr 3 11:17:51 AM EDT 2026 ( transpgrid-srv-r9.pppl.gov )
==========check for 98004Q02_mdsmark.sh===
==========(runtrx_r9)======================
==========TRDAT execution===============
date: Fri Apr 3 11:17:51 AM EDT 2026 ( transpgrid-srv-r9.pppl.gov )
TRDAT
TRANSP Version: 21.X
TRANSP DOI: 10.11578/dc.20180627.4
Build Date: 2021/--/--
%trmpi_init.f90: LOG_LEVEL env. var.: 1
%trmpi_init.f90: logfile_level: warn
CTOKNM : Cdir = /scratch/shared/tr_jlombard/transp_compute/JET/98004Q02/
% DIRECTORY IS /scratch/shared/tr_jlombard/transp_compute/JET/98004Q02/; TOKAMAK ID SET TO JET
%trmpi_openlog: LOGFILE_LEVEL = warn
trcore/pre_tok -I- JET specific defaults set
namelist elements assigned more than once WITH change of value:
NNPHI NUM_NPHI
WNPHI
namelist element value field(s): decimal point(s) inserted:
AIMP XZIMP
XZIMPS XZMINI
AMINI FRACMINI
DTMINT DTMING
XZRCYC
%NLIST: open namelist file98004Q02TR.ZDA
%trcom_static_box: loading static data
-----------> entering UF0GET.
NAMED FILE: ./98004Q02_tmp/OMF98004.NB2
%UF0DRD: DIMENSIONALITY OF INPUT FILE = 2
% NO. OF BEAMS (NB2 DATA) = 16
NAMED FILE: ./98004Q02_tmp/OMF98004.SIM1
%UF0DRD: DIMENSIONALITY OF INPUT FILE = 2
%UF0GET: ATOMIC WEIGHT NOT FOUND IN SIM01 USING A = 9.012
%UF0GET: ATOMIC NUMBER NOT FOUND IN SIM01 USING Z = 4.000
NAMED FILE: ./98004Q02_tmp/OMF98004.SIM2
%UF0DRD: DIMENSIONALITY OF INPUT FILE = 2
%UF0GET: ATOMIC WEIGHT NOT FOUND IN SIM02 USING A = 58.693
%UF0GET: ATOMIC NUMBER NOT FOUND IN SIM02 USING Z = 28.000
NAMED FILE: ./98004Q02_tmp/OMF98004.MMX
%UF0DRD: DIMENSIONALITY OF INPUT FILE = 3
UF0GET: -> Btor is clockwise from BTOR_CCW scalar in MMX data.
UF0GET: -> Itor is clockwise from ITOR_CCW scalar in MMX data.
UF0GET: MMX max exponent: 2
%UF0GET -- finished.
%DATCHK: full equilibrium input
DATCHK: CHECKING ARRAY DIMENSIONS, SWITCHES, ETC
%DATCHK: NYXINV incremented to be odd: 101 161
%DATCHK: no ECH/ECCD, NLECH=F,
%DATCHK: no Lower Hybrid, NLLH=F
%DATCKICH: TORIC-5 solver, with numerical MHD equilibrium geometry.
%DATCKICH: RGEOANT_A is in use => ignore NGEOANT=1
NGEOANT reset to match number of non-zero RGEOANT values.
*** NGEOANT reset to: 10
Warning! Result depends on setting RGEOANT_A!
%datckich for TORIC5: increased NICHPSI from 128 to 320
%DATCKA: ACfile times pre-screen...
%DATCKA -- CHECKING NON-DEFAULT SPECIFICATIONS FOR AC FILE I/O
GIVEN IN NAMELIST CHARACTER DATA ARRAYS "SELOUT" AND "SELAVG"
%ASCXFIN -- AC FILE OPERATION COMPLETED SUCCESSFULLY.
dummy trdatusub vsn trgdat_check called (OK).
dummy trdatusub vsn trgdat_psload called (OK).
*UTRCHK: CHECKING TRDAT NAMELIST SPECIFICATIONS.
%UTRCHK: full equilibrium input
*** OFFSET (TIME0) OF 4.0000E+01 SECONDS APPLIED TO ALL INPUT DATA ***
*DDEFIN: SET UP PHYSICS DATA FOR ANALYSIS:
TIME RANGE IS TINIT= 8.0000E+00, FTIME= 1.3000E+01 SECONDS.
SETSCD - SET SCALAR DATA FOR UPF
SETTIM - SET UP TIMING INFORMATION FOR UPF
*UF1GET: PROCESSING "CUR": plasma current amps -------------------
NAMED FILE: ./98004Q02_tmp/OMF98004.CUR
*UF1GET: PROCESSING "RBZ": [ext. B field] * R Tesla*cm -----------
NAMED FILE: ./98004Q02_tmp/OMF98004.RBZ
%FIXUNS: CONVERTING UFILE RBZ DATA UNITS, M.T , TO Tesla*cm
>CONVERSION FACTOR: * 100.000
*UF1GET: PROCESSING "DFL": diamagnetic flux Webers ---------------
NAMED FILE: ./98004Q02_tmp/MAGN98004.FLX
*UF1GET: PROCESSING "GAS": gas flow source #/sec -----------------
NAMED FILE: ./98004Q02_tmp/GASM98004.ELER
*UF1GET: PROCESSING "NTX": Neutron flux n/sec --------------------
NAMED FILE: ./98004Q02_tmp/TIN98004.RNT
*UF1GET: PROCESSING "PF0": Psi(axis) Wb/rad ----------------------
NAMED FILE: ./98004Q02_tmp/OMF98004.PF0
*UF1GET: PROCESSING "PLF": Poloidal flux Wb/rad ------------------
NAMED FILE: ./98004Q02_tmp/OMF98004.PLF
*UF1GET: PROCESSING "RCY": recycling source #/sec ----------------
NAMED FILE: ./98004Q02_tmp/EDG798004.FLW
*UF1GET: PROCESSING "TRF": Toroidal flux Webers ------------------
NAMED FILE: ./98004Q02_tmp/OMF98004.TRF
*UF1GET: PROCESSING "VSF": surface voltage volts -----------------
NAMED FILE: ./98004Q02_tmp/EFP98004.VPA
*UF1GET: PROCESSING "ZEF": Zeff -------------------------------
NAMED FILE: ./98004Q02_tmp/KS398004.ZEFH
%FIXUNS: Rescaling Ufile ZEF DATA by 1.04400 ( SC_ZEF )
*UF3GET: PROCESSING "MMX": Equilibrium moments cm ----------------
NAMED FILE: ./98004Q02_tmp/OMF98004.MMX
%FIXUNS: CONVERTING UFILE DATA UNITS, M , TO cm
>CONVERSION FACTOR: * 100.000
---> Geteq: #moments NMOM = 16
Geteq: Interpolate equilibrium data for UPF
SETLIM - SET UP Axisymmetric limiter data for UPF
*UF1GET: PROCESSING "LIM": Axisymmetric Limiter cm ---------------
NAMED FILE: ./98004Q02_tmp/OMF98004.LIM
%UF1GET: X axis is not the timebase, processing axis units cm
%FIXUNS: CONVERTING UFILE LIM TIMEBASE UNITS, M , TO cm
>CONVERSION FACTOR: * 100.000
%FIXUNS: CONVERTING UFILE LIM DATA UNITS, M , TO cm
>CONVERSION FACTOR: * 100.000
*UF2GET: PROCESSING "BOL": Power radiated W/cm**3 ----------------
NAMED FILE: ./98004Q02_tmp/BOLP98004.TOBP
_%UF2GET FALLBACK: 1D UFILE INPUT => TIME SERIES F(T) DATA
%units check required: fallback to time series.
*UF2GET: PROCESSING "GRB": MHD g (R*Bt) profile Tesla*cm ---------
NAMED FILE: ./98004Q02_tmp/OMF98004.GRB
%UF2GET -- INPUT UFILE AXES OUT OF ORDER, WILL SWAP:
1ST DIMENSION = " "; 2ND DIMENSION = "Seconds "
%FIXUNS: CONVERTING UFILE GRB DATA UNITS, TESLA*M , TO Tesla*cm
>CONVERSION FACTOR: * 100.000
%SETDAT2-- X axis is normalized flux, not modified, npts= 41
*UF2GET: PROCESSING "NER": electron density n/cm**3 --------------
NAMED FILE: ./98004Q02_tmp/NE198004.NE
%UF2GET -- INPUT UFILE AXES OUT OF ORDER, WILL SWAP:
1ST DIMENSION = "SQRTPSI "; 2ND DIMENSION = "SEC "
%FIXUNS: CONVERTING UFILE NER DATA UNITS, M**-3 , TO n/cm**3
>CONVERSION FACTOR: *0.100000E-05
%SETDAT2-- X axis is normalized flux, not modified, npts=131
*UF2GET: PROCESSING "OMG": plasma rotation rad/sec ---------------
NAMED FILE: ./98004Q02_tmp/AF198004.AF
%UF2GET -- INPUT UFILE AXES OUT OF ORDER, WILL SWAP:
1ST DIMENSION = "PSIN "; 2ND DIMENSION = "SEC "
%CHKTIM: DATA TIME RANGE 8.2742E+00 TO 1.3027E+01 SECONDS
DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME).
%SETDAT2-- X axis is normalized flux, not modified, npts=131
*UF2GET: PROCESSING "PRS": MHD Pressure profile Pascals ----------
NAMED FILE: ./98004Q02_tmp/OMF98004.PRS
%UF2GET -- INPUT UFILE AXES OUT OF ORDER, WILL SWAP:
1ST DIMENSION = " "; 2ND DIMENSION = "Seconds "
%SETDAT2-- X axis is normalized flux, not modified, npts= 41
*UF2GET: PROCESSING "QPR": q profile ---------------------------
NAMED FILE: ./98004Q02_tmp/OMF98004.QPR
%UF2GET -- INPUT UFILE AXES OUT OF ORDER, WILL SWAP:
1ST DIMENSION = " "; 2ND DIMENSION = "Seconds "
%SETDAT2-- X axis is normalized flux, not modified, npts= 41
*UF2GET: PROCESSING "TER": electron temperature eV ---------------
NAMED FILE: ./98004Q02_tmp/TE198004.TE
%UF2GET -- INPUT UFILE AXES OUT OF ORDER, WILL SWAP:
1ST DIMENSION = "SQRTPSI "; 2ND DIMENSION = "SEC "
%SETDAT2-- X axis is normalized flux, not modified, npts=131
*UF2GET: PROCESSING "TI2": ion temperature eV --------------------
NAMED FILE: ./98004Q02_tmp/TI198004.TI
%UF2GET -- INPUT UFILE AXES OUT OF ORDER, WILL SWAP:
1ST DIMENSION = "PSIN "; 2ND DIMENSION = "SEC "
%CHKTIM: DATA TIME RANGE 8.1241E+00 TO 1.3027E+01 SECONDS
DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME).
%SETDAT2-- X axis is normalized flux, not modified, npts=131
*UF3GET: PROCESSING "SIM01": impurity species n/cm**3 ------------
NAMED FILE: ./98004Q02_tmp/OMF98004.SIM1
_%UF3GET FALLBACK: 2D UFILE INPUT => Y-AXIS INVARIANT PROFILE
%UF3GET -- INPUT UFILE AXES OUT OF ORDER, WILL SWAP:
1ST DIMENSION = "SQRTPSI "; 2ND DIMENSION = "SEC "
3RD DIMENSION = " "
%FIXUNS: CONVERTING UFILE DATA UNITS, M**-3 , TO n/cm**3
>CONVERSION FACTOR: *0.100000E-05
%SETDAT2-- X axis is normalized flux, not modified, npts=131
*UF3GET: PROCESSING "SIM02": impurity species n/cm**3 ------------
NAMED FILE: ./98004Q02_tmp/OMF98004.SIM2
_%UF3GET FALLBACK: 2D UFILE INPUT => Y-AXIS INVARIANT PROFILE
%UF3GET -- INPUT UFILE AXES OUT OF ORDER, WILL SWAP:
1ST DIMENSION = "SQRTPSI "; 2ND DIMENSION = "SEC "
3RD DIMENSION = " "
%FIXUNS: CONVERTING UFILE DATA UNITS, M**-3 , TO n/cm**3
>CONVERSION FACTOR: *0.100000E-05
%SETDAT2-- X axis is normalized flux, not modified, npts=131
*DDEFIN: Set Plasma Heating parameters :- ==========================
-> entering setnbi (neutral beam powers, voltages, etc.)
SETNBI - SET UP NEUTRAL BEAM DATA FOR UPF
*UF2GET: PROCESSING "NB2": beam heating data -------------------
NAMED FILE: ./98004Q02_tmp/OMF98004.NB2
%SETNBI- BEAM FULL ENERGY FRACTIONS SET BY NBI FILE **
%SETNBI- BEAM HALF ENERGY FRACTIONS SET BY NBI FILE **
%setnbi: NBI on/off times taken from data.
%NB_ONOFF: NBI on/off times (s): 7.5564E+00 1.3084E+01
-> entering setrfp (ICRF powers)
SETRFP - SET UP ICRF POWER DATA FOR UPF
*UF2GET: PROCESSING "RFP": RF heating data ---------------------
NAMED FILE: ./98004Q02_tmp/OMF98004.RFP
%ICRF on/off times taken from data.
%ICRF_ONOFF: ICRF on/off times (s): 8.0709E+00 1.5434E+01
SETRFP - SET UP ICRF FREQUENCY DATA FOR UPF
*UF2GET: PROCESSING "RFF": RF frequency data Hz ------------------
NAMED FILE: ./98004Q02_tmp/OMF98004.RFF
Zeff vs. time only:
Zeff vs. time input data provided.
%SETDFLT: "DTG" INPUT DATA FILE NOT SPECIFIED
==>DEFAULT VALUE SUBSTITUTED, DTMAXG = 2.0000E-02
%setmisc -- time dependent geometry timestep enabled (ZOOM is off).
%SETDFLT: "DTS" INPUT DATA FILE NOT SPECIFIED
==>DEFAULT VALUE SUBSTITUTED, DT_SOURCES = 5.0000E-03
%setmisc -- time dependent source timestep enabled (ZOOM is off).
%SETDFLT: "DTX" INPUT DATA FILE NOT SPECIFIED
==>DEFAULT VALUE SUBSTITUTED, DT_OUTPUT = 2.0000E-02
%setmisc -- time dependent output timestep enabled (ZOOM is off).
--------------------------------------------------
%trdat_bdy_check -- checking plasma boundary data
and checking plasma boundary is inside vacuum vessel
trdat_bdy_check: OK
...min(r_curvature/(Rmax-Rmin)) is: 9.5978722396684951E-002
...at time t= 13.000000000000000 seconds.
------------> thermal ion ptcl balance summary:
Hydrogen (absent)
Deuterium NDEFINE=0 NMODEL=4 RC=NLRCYC GF=NLRCYC
Tritium (absent)
Helium-3 (absent)
Helium-4 (absent)
Lithium-6 (absent)
Recycling fractions RFRAC are specified.
------------------------------------
Fast ion summary:
Beam species: D
Fusion product species: He4
Minority RF species: H
------------------------------------
(mpi_share_env) process myid= 0 cwd: /scratch/shared/tr_jlombard/transp_compute/JET/98004Q02
==========(runtrx_r9 TRDAT sts = 0)======
==========LABEL execution===============
date: Fri Apr 3 11:17:53 AM EDT 2026 ( transpgrid-srv-r9.pppl.gov )
==>runtrx_r9: label execution successful
==>runtrx_r9: LINK_ON_HOST environment variable detected.
==========(runtrx_r9)======================
==========>runtrx_r9 normal exit<==========
date: Fri Apr 3 11:17:53 AM EDT 2026 ( transpgrid-srv-r9.pppl.gov )