TRANSP Grid Analysis 36759D01 KSTR trdat.log
==>runtrx_r9 start: date: Fri Apr 3 11:21:46 AM EDT 2026 ( flux-login2.pppl.gov )
argv = 2
iarg = 2
cmd_opt = pre
runtrx_r9: tok.yy = KSTR.25
==========(runtrx_r9)======================
date: Fri Apr 3 11:21:46 AM EDT 2026 ( flux-login2.pppl.gov )
args: 36759D01 pre
==========(runtrx_r9)======================
==========(runtrx_r9)======================
==========TRANSP preprocessing==========
date: Fri Apr 3 11:21:46 AM EDT 2026 ( flux-login2.pppl.gov )
==========check for 36759D01_mdsmark.sh===
==========(runtrx_r9)======================
==========TRDAT execution===============
date: Fri Apr 3 11:21:46 AM EDT 2026 ( flux-login2.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_halkhawa/transp_compute/KSTR/36759D01/
% DIRECTORY IS /scratch/shared/tr_halkhawa/transp_compute/KSTR/36759D01/; TOKAMAK ID SET TO KSTR
%trmpi_openlog: LOGFILE_LEVEL = warn
%splitn_module: update block detected, t= 4.5099999999999998
namelist element value field(s): decimal point(s) inserted:
DTN FOCLRA
DTMINT DTMING
FOCLZA
%NLIST: open namelist file36759D01TR.ZDA
%trcom_static_box: loading static data
%mds_setup: setting MDSplus server to: TRANSPGRID.PPPL.GOV
%mds_setup: server set OK.
%mds_setup: MDSplus opening tree: TRANSP_KSTR 367590401
%mds_setup: MDSplus tree open status: 265389633
-----------> entering UF0GET.
%UF0DRD: DIMENSIONALITY OF INPUT FILE = 2
% NO. OF BEAMS (NB2 DATA) = 5
%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 Lower Hybrid, NLLH=F
%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
%UTRCZEF - ZF2 DATA SPECIFIED - PLEASE SET NLZFI2=.T TO USE.
*DDEFIN: SET UP PHYSICS DATA FOR ANALYSIS:
TIME RANGE IS TINIT= 4.5000E+00, FTIME= 9.0100E+00 SECONDS.
SETSCD - SET SCALAR DATA FOR UPF
SETTIM - SET UP TIMING INFORMATION FOR UPF
*UF1GET: PROCESSING "CUR": plasma current amps -------------------
%CHKTIM: DATA TIME RANGE 1.5000E+00 TO 8.9900E+00 SECONDS
DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME).
*UF1GET: PROCESSING "RBZ": [ext. B field] * R Tesla*cm -----------
%FIXUNS: CONVERTING UFILE RBZ DATA UNITS, M.T , TO Tesla*cm
>CONVERSION FACTOR: * 100.000
%CHKTIM: DATA TIME RANGE 1.5000E+00 TO 8.9900E+00 SECONDS
DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME).
*UF1GET: PROCESSING "BDI": beta diamagnetic --------------------
%CHKTIM: DATA TIME RANGE 1.5000E+00 TO 8.9900E+00 SECONDS
DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME).
*UF1GET: PROCESSING "DFL": diamagnetic flux Webers ---------------
%CHKTIM: DATA TIME RANGE 1.5000E+00 TO 8.9900E+00 SECONDS
DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME).
%SETDFLT: "GAS" INPUT DATA FILE NOT SPECIFIED
==>DEFAULT VALUE SUBSTITUTED, GASFLN = 0.0000E+00
*UF1GET: PROCESSING "L2B": li/2 + beta(pol) --------------------
%CHKTIM: DATA TIME RANGE 1.5000E+00 TO 8.9900E+00 SECONDS
DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME).
*UF1GET: PROCESSING "PF0": Psi(axis) Wb/rad ----------------------
%CHKTIM: DATA TIME RANGE 1.5000E+00 TO 8.9900E+00 SECONDS
DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME).
*UF1GET: PROCESSING "PLF": Poloidal flux Wb/rad ------------------
%CHKTIM: DATA TIME RANGE 1.5000E+00 TO 8.9900E+00 SECONDS
DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME).
*UF1GET: PROCESSING "TPI": particle confinement seconds ----------
*UF1GET: PROCESSING "TRF": Toroidal flux Webers ------------------
%CHKTIM: DATA TIME RANGE 1.5000E+00 TO 8.9900E+00 SECONDS
DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME).
*UF1GET: PROCESSING "VSF": surface voltage volts -----------------
%CHKTIM: DATA TIME RANGE 1.5000E+00 TO 8.9900E+00 SECONDS
DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME).
*UF3GET: PROCESSING "MMX": Equilibrium moments cm ----------------
%FIXUNS: CONVERTING UFILE DATA UNITS, M , TO cm
>CONVERSION FACTOR: * 100.000
%CHKTIM: DATA TIME RANGE 1.5000E+00 TO 8.9900E+00 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 ---------------
%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 ----------------
%UF2GET -- INPUT UFILE AXES OUT OF ORDER, WILL SWAP:
1ST DIMENSION = " "; 2ND DIMENSION = "seconds "
%SETDAT2-- X axis is normalized flux, not modified, npts=201
*UF2GET: PROCESSING "GRB": MHD g (R*Bt) profile Tesla*cm ---------
%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 1.5000E+00 TO 8.9900E+00 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 --------------
%UF2GET -- INPUT UFILE AXES OUT OF ORDER, WILL SWAP:
1ST DIMENSION = " "; 2ND DIMENSION = "seconds "
%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=201
*UF2GET: PROCESSING "NIM": impurity density n/cm**3 --------------
%UF2GET -- INPUT UFILE AXES OUT OF ORDER, WILL SWAP:
1ST DIMENSION = " "; 2ND DIMENSION = "seconds "
%FIXUNS: CONVERTING UFILE NIM DATA UNITS, M**-3 , TO n/cm**3
>CONVERSION FACTOR: *0.100000E-05
%SETDAT2-- X axis is normalized flux, not modified, npts=201
*UF2GET: PROCESSING "OMG": plasma rotation rad/sec ---------------
%UF2GET -- INPUT UFILE AXES OUT OF ORDER, WILL SWAP:
1ST DIMENSION = " "; 2ND DIMENSION = "seconds "
%SETDAT2-- X axis is normalized flux, not modified, npts=201
*UF2GET: PROCESSING "PRS": MHD Pressure profile Pascals ----------
%UF2GET -- INPUT UFILE AXES OUT OF ORDER, WILL SWAP:
1ST DIMENSION = " "; 2ND DIMENSION = "seconds "
%CHKTIM: DATA TIME RANGE 1.5000E+00 TO 8.9900E+00 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 ---------------------------
%UF2GET -- INPUT UFILE AXES OUT OF ORDER, WILL SWAP:
1ST DIMENSION = " "; 2ND DIMENSION = "seconds "
%CHKTIM: DATA TIME RANGE 1.5000E+00 TO 8.9900E+00 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 ---------------
%UF2GET -- INPUT UFILE AXES OUT OF ORDER, WILL SWAP:
1ST DIMENSION = " "; 2ND DIMENSION = "seconds "
%FIXUNS: CONVERTING UFILE TER DATA UNITS, KEV , TO eV
>CONVERSION FACTOR: * 1000.00
%SETDAT2-- X axis is normalized flux, not modified, npts=201
*UF2GET: PROCESSING "TI2": ion temperature eV --------------------
%UF2GET -- INPUT UFILE AXES OUT OF ORDER, WILL SWAP:
(mpi_share_env) process myid= 0 cwd: /scratch/shared/tr_halkhawa/transp_compute/KSTR/36759D01
In TRDAT: calling copy_toraynml
1ST DIMENSION = " "; 2ND DIMENSION = "seconds "
%FIXUNS: CONVERTING UFILE TI2 DATA UNITS, KEV , TO eV
>CONVERSION FACTOR: * 1000.00
%SETDAT2-- X axis is normalized flux, not modified, npts=201
*UF2GET: PROCESSING "ZF2": Zeff profile ------------------------
%UF2GET -- INPUT UFILE AXES OUT OF ORDER, WILL SWAP:
1ST DIMENSION = " "; 2ND DIMENSION = "seconds "
%SETDAT2-- X axis is normalized flux, not modified, npts=201
*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 -------------------
%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): 1.2901E+00 9.2099E+00
-> entering setecp (ECH/ECCD powers)
SETECP - SET UP ECRH POWER DATA FOR UPF
*UF2GET: PROCESSING "ECP": ECRF antenna powers watts -------------
%ECH/ECCD on/off times taken from data.
%EC_ONOFF: ECH on/off times (s): 2.7000E+00 9.2099E+00
SETECP - read polar/poloidal angles vs. time
*UF2GET: PROCESSING "ECA": ECRF poloidal aim degrees -------------
SETECP - read azimuthal/toroidal angles vs. time
*UF2GET: PROCESSING "ECB": ECRF toroidal aim degrees -------------
Non-simple Zeff data:
Zeff inferred from impurity density data.
%SETDFLT: "DTG" INPUT DATA FILE NOT SPECIFIED
==>DEFAULT VALUE SUBSTITUTED, DTMAXG = 1.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
trdat_bdy_check: OK
...min(r_curvature/(Rmax-Rmin)) is: 5.0079398599491759E-002
...at time t= 4.5699999999999985 seconds.
------------> thermal ion ptcl balance summary:
Hydrogen (absent)
Deuterium NDEFINE=0 NMODEL=4 RC=Tau(p) GF=ZERO
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: (None)
Minority RF species: (None)
------------------------------------
copy_toraynml: template fetch via "mds_get_inf" subprocess...
COPY_TORAYNML: $UFILES_ROOT/bin/mds_get_inf s TRANSPGRID.PPPL.GOV t TRANSP_KSTR r 36759D01 l .NML:TORAY q
************************************************************************
Program: mds_get_inf
Version: v0.6.0
Revision: 673c23d
Date: 02/11/2026
************************************************************************
mds_get_inf retrievs a text node from a TRANSP run
and writes it into a file <runid>.<node_label>
or uses TRANSP conventions:
name_list -> <runid>TR.DAT
contents -> <runid>TF.PLN
mds_shot = 367590401
mds_conopn: option = 3 367590401 TRANSP_KSTR
...connecting to server: TRANSPGRID.PPPL.GOV
mds_shot = 367590401
.NML:TORAY does not exist
STOP 1
COPY_TORAYNML: failed to extract .NML:TORAY
(now look for template file).
COPY_TORAYNML: namelist: TORAY_TEMPLATE= "/p/ptransp/halkhawa/template_new/input_toray.nml"
COPY_TORAYNML: Toray template file=
/p/ptransp/halkhawa/template_new/input_toray.nml
COPY_TORAYNML: zcmd= cp /p/ptransp/halkhawa/template_new/input_toray.nml ./36759D01_toray.in
In TRDAT: calling copy_gafitnml
COPY_TORAYNML: file copy status (0 is normal): 0
COPY_GAFITNML: namelist: GAFIT_TEMPLATE= ""
COPY_GAFITNML: GAFIT_TEMPLATE is blank; Default used
COPY_GAFITNML: Gafit template file=
$CODESYSDIR/tables/rf_nml/gafit.in
COPY_GAFITNML: zcmd= cp $CODESYSDIR/tables/rf_nml/gafit.in ./36759D01_gafit.in
In TRDAT: calling COPY_PTNML
COPY_GAFITNML: file copy status (0 is normal): 0
COPY_PTNML_SUB: namelist: pt_template= "/p/ptransp/halkhawa/template_new/pt_inp.nml"
COPY_PTNML_SUB: PT_SOLVER namelist file=
/p/ptransp/halkhawa/template_new/pt_inp.nml
COPY_PTNML_SUB: zcmd= cp /p/ptransp/halkhawa/template_new/pt_inp.nml ./36759D01_pt.nml
In TRDAT: calling COPY_MMMNML
COPY_PTNML_SUB: file copy status (0 is normal): 0
COPY_MMM_SUB: namelist: mmm_template= "/p/ptransp/halkhawa/template_new/mmm_ht.nml"
COPY_MMM_SUB: MMM namelist file=
/p/ptransp/halkhawa/template_new/mmm_ht.nml
COPY_MMM_SUB: zcmd= cp /p/ptransp/halkhawa/template_new/mmm_ht.nml ./36759D01_mmm.nml
COPY_MMM_SUB: file copy status (0 is normal): 0
==========(runtrx_r9 TRDAT sts = 0)======
==>runtrx_r9: trdat execution successful
will write to mdsplus later
==========LABEL execution===============
date: Fri Apr 3 11:21:52 AM EDT 2026 ( flux-login2.pppl.gov )
==>runtrx_r9: label execution successful
==========add TF.PLN to MDSplus ========
==========(runtrx_r9)======================
date: Fri Apr 3 11:21:53 AM EDT 2026 ( flux-login2.pppl.gov )
************************************************************************
Program: mds_start
Version: v0.6.0
Revision: 673c23d
Date: 02/11/2026
************************************************************************
mds_conopn: option = 4 367590401 TRANSP_KSTR
...connecting to server: TRANSPGRID.PPPL.GOV
...tcl("EDIT TRANSP_KSTR/SHOT=367590401")
nid = 1
open 36759D01TF.PLN
add_file: 6639 lines - 63
tcl("WRITE")
stat = 1
chmod: changing permissions of '/p/transpgrid/mds_trees/transp_kstr/transp_kstr_367590401.tree': Operation not permitted
==========add PH.CDF to MDSplus ========
==========(runtrx_r9)======================
date: Fri Apr 3 11:21:53 AM EDT 2026 ( flux-login2.pppl.gov )
==>runtrx_r9: LINK_ON_HOST environment variable detected.
==========(runtrx_r9)======================
==========>runtrx_r9 normal exit<==========
date: Fri Apr 3 11:21:53 AM EDT 2026 ( flux-login2.pppl.gov )