==>runtrx start: date: Fri Jun 28 09:10:47 EDT 2024 ( mccune008.pppl.gov ) argv = 2 iarg = 2 cmd_opt = pre runtrx: tok.yy = JET.21 ==========(runtrx)====================== date: Fri Jun 28 09:10:48 EDT 2024 ( mccune008.pppl.gov ) args: 104522M11 pre ==========(runtrx)====================== ==========(runtrx)====================== ==========TRANSP preprocessing========== date: Fri Jun 28 09:10:48 EDT 2024 ( mccune008.pppl.gov ) ==========check for 104522M11_mdsmark.sh=== ==========(runtrx)====================== ==========TRDAT execution=============== date: Fri Jun 28 09:10:48 EDT 2024 ( mccune008.pppl.gov ) TRDAT TRANSP Version: 24.0.1 TRANSP DOI: 10.11578/dc.20180627.4 Build Date: 2023/10/21 %trmpi_init.f90: LOG_LEVEL env. var.: 1 %trmpi_init.f90: logfile_level: warn (mpi_share_env) process myid= 0 cwd: /local/tr_mporadzi/transp_compute/JET/104522M11 CTOKNM : Cdir = /local/tr_mporadzi/transp_compute/JET/104522M11/ % DIRECTORY IS /local/tr_mporadzi/transp_compute/JET/104522M11/; 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: XZIMP XZIMPS XZMINI AMINI FRACMINI DTMINT DTMING XZRCYC %NLIST: open namelist file104522M11TR.ZDA %trcom_static_box: loading static data -----------> entering UF0GET. NAMED FILE: ./104522M11_tmp/OMF104522.NB2 %UF0DRD: DIMENSIONALITY OF INPUT FILE = 2 % NO. OF BEAMS (NB2 DATA) = 14 NAMED FILE: ./104522M11_tmp/OMF104522.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: ./104522M11_tmp/OMF104522.SIM2 %UF0DRD: DIMENSIONALITY OF INPUT FILE = 2 %UF0GET: ATOMIC WEIGHT NOT FOUND IN SIM02 USING A = 183.840 %UF0GET: ATOMIC NUMBER NOT FOUND IN SIM02 USING Z = 74.000 NAMED FILE: ./104522M11_tmp/OMF104522.SIM3 %UF0DRD: DIMENSIONALITY OF INPUT FILE = 2 %UF0GET: ATOMIC WEIGHT NOT FOUND IN SIM03 USING A = 58.693 %UF0GET: ATOMIC NUMBER NOT FOUND IN SIM03 USING Z = 28.000 NAMED FILE: ./104522M11_tmp/OMF104522.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: ndefine( 2 )=2; frac( 2 ) ignored as density profile is prescribed. %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= 7.0050E+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: ./104522M11_tmp/OMF104522.CUR %CHKTIM: DATA TIME RANGE 7.0364E+00 TO 1.3000E+01 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). *UF1GET: PROCESSING "RBZ": [ext. B field] * R Tesla*cm ----------- NAMED FILE: ./104522M11_tmp/OMF104522.RBZ %FIXUNS: CONVERTING UFILE RBZ DATA UNITS, M.T , TO Tesla*cm >CONVERSION FACTOR: * 100.000 %CHKTIM: DATA TIME RANGE 7.0364E+00 TO 1.3000E+01 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). *UF1GET: PROCESSING "DFL": diamagnetic flux Webers --------------- NAMED FILE: ./104522M11_tmp/MAGN104522.FLX *UF1GET: PROCESSING "GAS": gas flow source #/sec ----------------- NAMED FILE: ./104522M11_tmp/GASM104522.ELER ?FIXUNS: UFILE GAS DATA UNITS LABEL "e/s " NOT RECOGNIZED ?UF1GET: UFILE UNITS ERROR ./104522M11_tmp/GASM104522.ELER (GAS) ***** FATAL ERROR ***** UF1GET -- ERROR(S) ACQUIRING 1D INPUT DATA ?SETD1D: ERROR PROCESSING => ./104522M11_tmp/GASM104522.ELER (GAS) *UF1GET: PROCESSING "NTX": Neutron flux n/sec -------------------- NAMED FILE: ./104522M11_tmp/TIN104522.RNT *UF1GET: PROCESSING "PF0": Psi(axis) Wb/rad ---------------------- NAMED FILE: ./104522M11_tmp/OMF104522.PF0 %CHKTIM: DATA TIME RANGE 7.0364E+00 TO 1.3000E+01 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). *UF1GET: PROCESSING "PLF": Poloidal flux Wb/rad ------------------ NAMED FILE: ./104522M11_tmp/OMF104522.PLF %CHKTIM: DATA TIME RANGE 7.0364E+00 TO 1.3000E+01 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). *UF1GET: PROCESSING "RCY": recycling source #/sec ---------------- NAMED FILE: ./104522M11_tmp/EDG7104522.FLW *UF1GET: PROCESSING "TRF": Toroidal flux Webers ------------------ NAMED FILE: ./104522M11_tmp/OMF104522.TRF %CHKTIM: DATA TIME RANGE 7.0364E+00 TO 1.3000E+01 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). *UF1GET: PROCESSING "VSF": surface voltage volts ----------------- NAMED FILE: ./104522M11_tmp/EFP104522.VPA *UF1GET: PROCESSING "ZEF": Zeff ------------------------------- NAMED FILE: ./104522M11_tmp/KS3104522.ZEFV *UF3GET: PROCESSING "MMX": Equilibrium moments cm ---------------- NAMED FILE: ./104522M11_tmp/OMF104522.MMX %FIXUNS: CONVERTING UFILE DATA UNITS, M , TO cm >CONVERSION FACTOR: * 100.000 %CHKTIM: DATA TIME RANGE 7.0364E+00 TO 1.3000E+01 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). %Geteq: MMX file does not cover time range of run; equilibrium geometry will be extrapolated flat to fill the time range ---> 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: ./104522M11_tmp/OMF104522.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: ./104522M11_tmp/BOLP104522.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: ./104522M11_tmp/OMF104522.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 %CHKTIM: DATA TIME RANGE 7.0364E+00 TO 1.3000E+01 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). %SETDAT2-- X axis is normalized flux, not modified, npts= 41 *UF2GET: PROCESSING "NER": electron density n/cm**3 -------------- NAMED FILE: ./104522M11_tmp/PRFL104522.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 %CHKTIM: DATA TIME RANGE 7.1700E+00 TO 1.2990E+01 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). %SETDAT2-- X axis is normalized flux, not modified, npts=131 *UF2GET: PROCESSING "NID": D+ ion density n/cm**3 ---------------- NAMED FILE: ./104522M11_tmp/D104522.DENS %UF2GET -- INPUT UFILE AXES OUT OF ORDER, WILL SWAP: 1ST DIMENSION = "SQRTPSI "; 2ND DIMENSION = "SEC " %FIXUNS: CONVERTING UFILE NID DATA UNITS, M**-3 , TO n/cm**3 >CONVERSION FACTOR: *0.100000E-05 %CHKTIM: DATA TIME RANGE 7.1700E+00 TO 1.2990E+01 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). %SETDAT2-- X axis is normalized flux, not modified, npts=131 *UF2GET: PROCESSING "OMG": plasma rotation rad/sec --------------- NAMED FILE: ./104522M11_tmp/PRFL104522.AF %UF2GET -- INPUT UFILE AXES OUT OF ORDER, WILL SWAP: 1ST DIMENSION = "SQPSINT "; 2ND DIMENSION = "SEC " %CHKTIM: DATA TIME RANGE 7.1700E+00 TO 1.2990E+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: ./104522M11_tmp/OMF104522.PRS %UF2GET -- INPUT UFILE AXES OUT OF ORDER, WILL SWAP: 1ST DIMENSION = " "; 2ND DIMENSION = "Seconds " %CHKTIM: DATA TIME RANGE 7.0364E+00 TO 1.3000E+01 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). %SETDAT2-- X axis is normalized flux, not modified, npts= 41 *UF2GET: PROCESSING "QPR": q profile --------------------------- NAMED FILE: ./104522M11_tmp/OMF104522.QPR %UF2GET -- INPUT UFILE AXES OUT OF ORDER, WILL SWAP: 1ST DIMENSION = " "; 2ND DIMENSION = "Seconds " %CHKTIM: DATA TIME RANGE 7.0364E+00 TO 1.3000E+01 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). %SETDAT2-- X axis is normalized flux, not modified, npts= 41 *UF2GET: PROCESSING "TER": electron temperature eV --------------- NAMED FILE: ./104522M11_tmp/PRFL104522.TE %UF2GET -- INPUT UFILE AXES OUT OF ORDER, WILL SWAP: 1ST DIMENSION = "SQRTPSI "; 2ND DIMENSION = "SEC " %CHKTIM: DATA TIME RANGE 7.1700E+00 TO 1.2990E+01 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). %SETDAT2-- X axis is normalized flux, not modified, npts=131 *UF2GET: PROCESSING "TI2": ion temperature eV -------------------- NAMED FILE: ./104522M11_tmp/PRFL104522.TI %UF2GET -- INPUT UFILE AXES OUT OF ORDER, WILL SWAP: 1ST DIMENSION = "SQPSINT "; 2ND DIMENSION = "SEC " %CHKTIM: DATA TIME RANGE 7.1700E+00 TO 1.2990E+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: ./104522M11_tmp/OMF104522.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 %CHKTIM: DATA TIME RANGE 7.1700E+00 TO 1.2990E+01 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). %SETDAT2-- X axis is normalized flux, not modified, npts=131 *UF3GET: PROCESSING "SIM02": impurity species n/cm**3 ------------ NAMED FILE: ./104522M11_tmp/OMF104522.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 %CHKTIM: DATA TIME RANGE 7.1700E+00 TO 1.2990E+01 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). %SETDAT2-- X axis is normalized flux, not modified, npts=131 *UF3GET: PROCESSING "SIM03": impurity species n/cm**3 ------------ NAMED FILE: ./104522M11_tmp/OMF104522.SIM3 _%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 %CHKTIM: DATA TIME RANGE 7.1700E+00 TO 1.2990E+01 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). %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: ./104522M11_tmp/OMF104522.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.0064E+00 1.3172E+01 -> entering setrfp (ICRF powers) SETRFP - SET UP ICRF POWER DATA FOR UPF *UF2GET: PROCESSING "RFP": RF heating data --------------------- NAMED FILE: ./104522M11_tmp/OMF104522.RFP %ICRF on/off times taken from data. %ICRF_ONOFF: ICRF on/off times (s): 7.0261E+00 1.2997E+01 SETRFP - SET UP ICRF FREQUENCY DATA FOR UPF *UF2GET: PROCESSING "RFF": RF frequency data Hz ------------------ NAMED FILE: ./104522M11_tmp/OMF104522.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 = 1.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.980888851222844E-002 ...at time t= 11.3349999999999 seconds. ------------> thermal ion ptcl balance summary: Hydrogen (absent) Deuterium NDEFINE=2 RC=NLRCYC GF=NLRCYC Tritium NDEFINE=0 NMODEL=4 RC=NLRCYC GF=NLRCYC 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 ------------------------------------ ==========(runtrx TRDAT sts = 0)====== ==========LABEL execution=============== date: Fri Jun 28 09:10:49 EDT 2024 ( mccune008.pppl.gov ) ==>runtrx: label execution successful ==>runtrx: LINK_ON_HOST environment variable detected. ==========(runtrx)====================== ==========>runtrx normal exit<========== date: Fri Jun 28 09:10:50 EDT 2024 ( mccune008.pppl.gov )