==>runtrx start: date: Thu Dec 26 14:57:44 EST 2024 ( mccune036.pppl.gov ) argv = 2 iarg = 2 cmd_opt = pre runtrx: tok.yy = NSTX.06 ==========(runtrx)====================== date: Thu Dec 26 14:57:44 EST 2024 ( mccune036.pppl.gov ) args: 120982J15 pre ==========(runtrx)====================== ==========(runtrx)====================== ==========TRANSP preprocessing========== date: Thu Dec 26 14:57:44 EST 2024 ( mccune036.pppl.gov ) ==========check for 120982J15_mdsmark.sh=== ==========(runtrx)====================== ==========TRDAT execution=============== date: Thu Dec 26 14:57:44 EST 2024 ( mccune036.pppl.gov ) TRDAT TRANSP Version: 24.4.0 TRANSP DOI: 10.11578/dc.20180627.4 Build Date: 2024/09/23 %trmpi_init.f90: LOG_LEVEL env. var.: 1 %trmpi_init.f90: logfile_level: warn (mpi_share_env) process myid= 0 cwd: /local/tr_gavdeeva/transp_compute/NSTX/120982J15 CTOKNM : Cdir = /local/tr_gavdeeva/transp_compute/NSTX/120982J15/ % DIRECTORY IS /local/tr_gavdeeva/transp_compute/NSTX/120982J15/; TOKAMAK ID SET TO NSTX %trmpi_openlog: LOGFILE_LEVEL = warn %splitn_module: update block detected, t= 0.390000000000000 namelist elements assigned the same value more than once: XYBAPA NLBDAT XYBSCA namelist element value field(s): decimal point(s) inserted: DTMING DTMINT RMJCHK GOOCON CRLMR1 CRLMY1 CRLMRD XKFAC %NLIST: open namelist file120982J15TR.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_NSTX 1209821015 %mds_setup: MDSplus tree open status: 265389633 -----------> entering UF0GET. %UF0DRD: DIMENSIONALITY OF INPUT FILE = 2 % NO. OF BEAMS (NB2 DATA) = 2 %UF0DRD: DIMENSIONALITY OF INPUT FILE = 3 UF0GET: -> Btor is clockwise from BTOR_CCW scalar in MMX data. UF0GET: -> Itor is counter-clockwise from ITOR_CCW scalar in MMX data. MMX: no max exponent found (OK). %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 TIDXSW defaulted: value of 0.05 assigned. %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 *DDEFIN: SET UP PHYSICS DATA FOR ANALYSIS: TIME RANGE IS TINIT= 1.9000E-01, FTIME= 6.6000E-01 SECONDS. SETSCD - SET SCALAR DATA FOR UPF SETTIM - SET UP TIMING INFORMATION FOR UPF *UF1GET: PROCESSING "CUR": plasma current amps ------------------- *UF1GET: PROCESSING "RBZ": [ext. B field] * R Tesla*cm ----------- *UF1GET: PROCESSING "DFL": diamagnetic flux Webers --------------- %SETDFLT: "GAS" INPUT DATA FILE NOT SPECIFIED ==>DEFAULT VALUE SUBSTITUTED, GASFLN = 0.0000E+00 *UF1GET: PROCESSING "NTX": Neutron flux n/sec -------------------- *UF1GET: PROCESSING "PLF": Poloidal flux Wb/rad ------------------ *UF1GET: PROCESSING "TRF": Toroidal flux Webers ------------------ *UF1GET: PROCESSING "VSF": surface voltage volts ----------------- *UF3GET: PROCESSING "MMX": Equilibrium moments cm ---------------- %FIXUNS: CONVERTING UFILE DATA UNITS, M , TO cm >CONVERSION FACTOR: * 100.000 ---> Geteq: #moments NMOM = 16 Geteq: Interpolate equilibrium data for UPF *UF2GET: PROCESSING "BOL": Power radiated W/cm**3 ---------------- %CKSETRD: NRIBOL=3 DATA: PLASMA COVERAGE CHECK: GREATEST LOWER BOUND (*) (R-R0)/A = -8.8884E-01 LEAST UPPER BOUND (*) (R-R0)/A = 1.1704E+00 *EXPECTED MINIMUM COVERAGE TO INCLUDE RANGE (R-R0)/A = -9.0000E-01 TO 9.0000E-01 AT ALL TIMES IN THE ANALYSIS. *REQUIREMENT RELAXED BY NAMELIST: XRCBOL= 1.000E+00 RMJCHK=-1.000E+00 %SETRD0: X VECTOR SET FROM UFILE *2-SIDED* MAJOR RADIUS DATA *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 %SETDAT2-- X axis is normalized flux, not modified, npts= 41 *UF2GET: PROCESSING "NER": electron density n/cm**3 -------------- %CKSETRD: NRINER=3 DATA: PLASMA COVERAGE CHECK: GREATEST LOWER BOUND (*) (R-R0)/A = -8.9219E-01 LEAST UPPER BOUND (*) (R-R0)/A = 1.1094E+00 *EXPECTED MINIMUM COVERAGE TO INCLUDE RANGE (R-R0)/A = -9.0000E-01 TO 9.0000E-01 AT ALL TIMES IN THE ANALYSIS. *REQUIREMENT RELAXED BY NAMELIST: XRCNER= 1.000E+00 RMJCHK=-1.000E+00 %SETRD0: X VECTOR SET FROM UFILE *2-SIDED* MAJOR RADIUS DATA *UF2GET: PROCESSING "NIM": impurity density n/cm**3 -------------- %CHKTIM: DATA TIME RANGE 8.5571E-02 TO 6.4557E-01 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). %SETRD0: X VECTOR SET FROM UFILE MINOR RADIUS DATA *UF2GET: PROCESSING "PRS": MHD Pressure profile Pascals ---------- %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 --------------------------- %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 --------------- %CKSETRD: NRITER=3 DATA: PLASMA COVERAGE CHECK: GREATEST LOWER BOUND (*) (R-R0)/A = -8.9219E-01 LEAST UPPER BOUND (*) (R-R0)/A = 1.1094E+00 *EXPECTED MINIMUM COVERAGE TO INCLUDE RANGE (R-R0)/A = -9.0000E-01 TO 9.0000E-01 AT ALL TIMES IN THE ANALYSIS. *REQUIREMENT RELAXED BY NAMELIST: XRCTER= 1.000E+00 RMJCHK=-1.000E+00 %SETRD0: X VECTOR SET FROM UFILE *2-SIDED* MAJOR RADIUS DATA *UF2GET: PROCESSING "TI2": ion temperature eV -------------------- %CHKTIM: DATA TIME RANGE 8.5571E-02 TO 6.4557E-01 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). %SETRD0: X VECTOR SET FROM UFILE MINOR RADIUS DATA *UF2GET: PROCESSING "VP2": plasma rotation rad/sec --------------- %CHKTIM: DATA TIME RANGE 8.5571E-02 TO 6.4557E-01 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). %SETRD0: X VECTOR SET FROM UFILE MINOR RADIUS DATA *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): 8.2500E-02 7.2350E-01 %setnbi: TBONMIN= 0.150000000000000 --delays beam turn-on from TBON= 8.250000327825546E-002 seconds. Non-simple Zeff data: Zeff inferred from impurity density data. %SETDFLT: "DTG" INPUT DATA FILE NOT SPECIFIED ==>DEFAULT VALUE SUBSTITUTED, DTMAXG = 2.0000E-03 %setmisc -- time dependent geometry timestep enabled (ZOOM is off). %SETDFLT: "DTS" INPUT DATA FILE NOT SPECIFIED ==>DEFAULT VALUE SUBSTITUTED, DT_SOURCES = 2.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-03 %setmisc -- time dependent output timestep enabled (ZOOM is off). %TDB_EXECSYM: PRESYMMETRIZE THE "TER" DATA -------------------------------------------------- %trdat_bdy_check -- checking plasma boundary data trdat_bdy_check: OK ...min(r_curvature/(Rmax-Rmin)) is: 4.420043231110674E-002 ...at time t= 0.190000000000000 seconds. ------------> thermal ion ptcl balance summary: Hydrogen NDEFINE=0 NMODEL=1 RC=Tau(p) GF=ZERO Deuterium NDEFINE=0 NMODEL=1 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) ------------------------------------ In TRDAT: calling COPY_PTNML COPY_PTNML_SUB: namelist: pt_template= " /p/omfit/tmp/jlestz/TRANSP/120982J15_tmp/120982J15_input_pt.nml" COPY_PTNML_SUB: PT_SOLVER namelist file= /p/omfit/tmp/jlestz/TRANSP/120982J15_tmp/120982J15_input_pt.nml COPY_PTNML_SUB: zcmd= cp /p/omfit/tmp/jlestz/TRANSP/120982J15_tmp/120982J15_input_pt.nml ./120982J15_ pt.nml COPY_PTNML_SUB: file copy status (0 is normal): 0 In TRDAT: calling COPY_TGLFNML COPY_TGLF_SUB: namelist: tglf_template= " /p/omfit/tmp/jlestz/TRANSP/120982J15_tmp/120982J15_input_tglf.nml" COPY_TGLF_SUB: TGLF namelist file= /p/omfit/tmp/jlestz/TRANSP/120982J15_tmp/120982J15_input_tglf.nml COPY_TGLF_SUB: zcmd= cp /p/omfit/tmp/jlestz/TRANSP/120982J15_tmp/120982J15_input_tglf.nml ./120982J1 5_tglf.nml COPY_TGLF_SUB: file copy status (0 is normal): 0 ==========(runtrx TRDAT sts = 0)====== ==>runtrx: trdat execution successful will write to mdsplus later ==========LABEL execution=============== date: Thu Dec 26 14:58:49 EST 2024 ( mccune036.pppl.gov ) ==>runtrx: label execution successful ==========add TF.PLN to MDSplus ======== ==========(runtrx)====================== date: Thu Dec 26 14:58:50 EST 2024 ( mccune036.pppl.gov ) ************************************************************************ Program: mds_start Version: v0.3.2 Revision: 953b182 Date: 03/20/2024 ************************************************************************ mds_conopn: option = 4 1209821015 TRANSP_NSTX ...connecting to server: TRANSPGRID.PPPL.GOV ...tcl("EDIT TRANSP_NSTX/SHOT=1209821015") nid = 155 open 120982J15TF.PLN add_file: 6289 lines - 63 tcl("WRITE") stat = 1 ==========add PH.CDF to MDSplus ======== ==========(runtrx)====================== date: Thu Dec 26 14:58:51 EST 2024 ( mccune036.pppl.gov ) ==>runtrx: LINK_ON_HOST environment variable detected. ==========(runtrx)====================== ==========>runtrx normal exit<========== date: Thu Dec 26 14:58:51 EST 2024 ( mccune036.pppl.gov )