==>runtrx start: date: Tue Jun 25 12:20:58 EDT 2024 ( transpgrid3.pppl.gov ) argv = 2 iarg = 2 cmd_opt = pre runtrx: tok.yy = D3D.17 ==========(runtrx)====================== date: Tue Jun 25 12:20:58 EDT 2024 ( transpgrid3.pppl.gov ) args: 170339K01 pre ==========(runtrx)====================== ==========(runtrx)====================== ==========TRANSP preprocessing========== date: Tue Jun 25 12:20:58 EDT 2024 ( transpgrid3.pppl.gov ) ==========check for 170339K01_mdsmark.sh=== ==========(runtrx)====================== ==========TRDAT execution=============== date: Tue Jun 25 12:20:58 EDT 2024 ( transpgrid3.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_kshabbir/transp_compute/D3D/170339K01 CTOKNM : Cdir = /local/tr_kshabbir/transp_compute/D3D/170339K01/ % DIRECTORY IS /local/tr_kshabbir/transp_compute/D3D/170339K01/; TOKAMAK ID SET TO D3D %trmpi_openlog: LOGFILE_LEVEL = warn namelist element value field(s): decimal point(s) inserted: ECB_CURV_HOR ECB_CURV_VERT GOOMAX FOCLRA FOCLZA PDELTA %NLIST: open namelist file170339K01TR.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_D3D 1703391101 %mds_setup: MDSplus tree open status: 265389633 -----------> entering UF0GET. %UF0DRD: DIMENSIONALITY OF INPUT FILE = 2 % NO. OF BEAMS (NB2 DATA) = 14 %UF0DRD: DIMENSIONALITY OF INPUT FILE = 3 UF0GET: -> Btor is counter-clockwise from BTOR_CCW scalar in MMX data. UF0GET: -> Itor is counter-clockwise from ITOR_CCW scalar in MMX data. UF0GET: MMX max exponent: 2 %UF0GET -- finished. %datchk: reset MRSTRT to 0 %datchk: NLFBM set to .TRUE., no option. %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 ***** FATAL ERROR ***** ?UTRCZEF - ZEFF - NLZFI2=.T BUT NO ZF2 DATA SPECIFIED %UTRCHK warning: NO TOROIDAL ROTATION DATA & NLVPHI=.TRUE. _UTRCHK- ROTATION DATA TRIGRAPHS: VPH VP2 OMG ***** FATAL ERROR ***** ?UTRCHK - NO NER FILE SPECIFIED - NE(R) DATA REQUIRED!! ***** FATAL ERROR ***** ?UTRCHK - NO TE(R) FILE - TER OR ECF - TE(R) DATA REQUIRED! ***** FATAL ERROR ***** ?UTRCHK - NLTI2 SET AND NO TI2 FILE SPECIFIED ?? UTRCHK -- FATAL NAMELIST ERRORS DETECTED ! *DDEFIN: SET UP PHYSICS DATA FOR ANALYSIS: TIME RANGE IS TINIT= 2.1000E-01, FTIME= 5.0900E+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 4.0500E-01 TO 4.7550E+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 4.0500E-01 TO 4.7550E+00 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). *UF1GET: PROCESSING "BDI": beta diamagnetic -------------------- %CHKTIM: DATA TIME RANGE 4.0500E-01 TO 4.7550E+00 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). *UF1GET: PROCESSING "DFL": diamagnetic flux Webers --------------- %CHKTIM: DATA TIME RANGE 4.0500E-01 TO 4.7550E+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 4.0500E-01 TO 4.7550E+00 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). *UF1GET: PROCESSING "NTX": Neutron flux n/sec -------------------- *UF1GET: PROCESSING "PF0": Psi(axis) Wb/rad ---------------------- %CHKTIM: DATA TIME RANGE 4.3500E-01 TO 4.7550E+00 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). *UF1GET: PROCESSING "PLF": Poloidal flux Wb/rad ------------------ %CHKTIM: DATA TIME RANGE 4.3500E-01 TO 4.7550E+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 4.3500E-01 TO 4.7550E+00 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). *UF1GET: PROCESSING "VSF": surface voltage volts ----------------- %CHKTIM: DATA TIME RANGE 4.0500E-01 TO 4.7550E+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 4.3500E-01 TO 4.7550E+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 "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 4.3500E-01 TO 4.7550E+00 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). %SETDAT2-- X axis is normalized flux, not modified, npts= 41 *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 4.3500E-01 TO 4.7550E+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 4.3500E-01 TO 4.7550E+00 SECONDS DOES NOT COVER THE ENTIRE ANALYSIS TIME RANGE (TINIT,FTIME). %SETDAT2-- X axis is normalized flux, not modified, npts= 41 *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): 3.9990E-01 5.2899E+00 Zeff vs. (x,t) only: Zeff vs. (x,t) input data provided. %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 = 2.0000E-02 %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). -------------------------------------------------- %trdat_bdy_check -- checking plasma boundary data trdat_bdy_check: OK ...min(r_curvature/(Rmax-Rmin)) is: 4.823835904506524E-002 ...at time t= 1.73000000000000 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) ------------------------------------ In TRDAT: calling COPY_PTNML COPY_PTNML_SUB: namelist: pt_template= " /p/ptransp/trafiq/transp/work/D3D/pt_inp.nml" COPY_PTNML_SUB: PT_SOLVER namelist file= /p/ptransp/trafiq/transp/work/D3D/pt_inp.nml COPY_PTNML_SUB: zcmd= cp /p/ptransp/trafiq/transp/work/D3D/pt_inp.nml ./170339K01_pt.nml COPY_PTNML_SUB: file copy status (0 is normal): 0 In TRDAT: calling COPY_MMMNML COPY_MMM_SUB: namelist: mmm_template= " /p/ptransp/trafiq/transp/work/D3D/mmm_ht.nml" COPY_MMM_SUB: MMM namelist file= /p/ptransp/trafiq/transp/work/D3D/mmm_ht.nml COPY_MMM_SUB: zcmd= cp /p/ptransp/trafiq/transp/work/D3D/mmm_ht.nml ./170339K01_mmm.nml COPY_MMM_SUB: file copy status (0 is normal): 0 TRDAT ERROR: There were errors reading the data USER TYPEAHEAD BUFFER CONTENTS: w q %ABORTR: TA = 0.00000E+00 ** TRANSP ABORTR SUBROUTINE CALLED ** %bad_exit: generic f77 error exit call (errset status 999) ERRSET called, istat= 999 program exiting... ==========(runtrx TRDAT sts = 1)====== ==========add BAD PH.CDF to MDSplus ======== ==========(runtrx)====================== date: Tue Jun 25 12:21:10 EDT 2024 ( transpgrid3.pppl.gov ) ==========(runtrx)====================== date: Tue Jun 25 12:21:20 EDT 2024 ( transpgrid3.pppl.gov ) ==========>runtrx error exit status: 1