dev_name = 'W7X' dev = pyfusion.getDevice(dev_name) shot = da['info']['shotdata']['shot'][0] chan = da['info']['channels'][0] echdata = dev.acq.getdata(shot,'W7X_TotECH') wech = np.where(echdata.signal > 100)[0] tech = echdata.timebase[wech[0]] t0_utc = int(tech * 1e9) + echdata.utc[0] dtprobe = (da['info']['params']['i_diag_utc'][0] - t0_utc)/1e9 tslice = da['t_mid'][s-step//2] + dtprobe dtaverage = da['t_mid'][s] - da['t_mid'][s-step] dtstr = '{dt:.2f}s avg'.format(dt=dtaverage) # axu.set_title('{fn}, time={t:.4f} {dt}'.format(fn=da.name, t=tslice, dt=dtstr)) tit = axu.get_title() cfg, cfg_dict = get_mag_config(shot) if cfg is None: ind = '' else: ind = 13 - 12*cfg[2]/.3904 if abs(ind - round(ind)) < .03: ind = int(round((ind))) ind = 'Ind {ind}'.format(ind=ind) axu.set_title(tit+'\nProgram Id {s}, time={t:.4f} {ind}'# {dt}' .format(s=shot, fn=da.name, t=tslice, dt=dtstr, ind=ind)) if cfg_dict is not None: plot_trim(axu, mag_dict=cfg_dict, color='g', aspect=1.2) if time_strip_h > 0:
tr=t_range, sh=shot_number)) fitsigs.append(fitsig1) sold, td, sig2, ebar2 = LCFS_plot(da, diag2, nrm=nrm, ax=axLCne, Tesmooth=fitsig1, **kwargs) sigs.append(sig2) tds.append(td) solds.append(sold) for fp in fp_list: fp.write(', Integral = {intg:.1f}\n'.format(intg=integrate(sig2, td))) # if diag1 in fits: params = fitgauss(sold, sig1, ebar1, ax=axLCne) maxTe = max(maxTe, np.nanmax(sig1)) maxne = max(maxne, np.nanmax(sig2)) #1/0 cfg, cfg_dict = get_mag_config(shot_number) if cfg_dict is not None: plot_trim(axLCTe, mag_dict=cfg_dict, aspect=2.6, size=0.06, color='g') ind = 13 - 12*cfg[2]/.3904 if abs(ind - round(ind)) < .03: ind = int(round((ind))) ind = 'Ind {ind}'.format(ind=ind) else: ind = '' tb = 'LP t_mid' if use_t_mid else 'into ECH' figLCFS.subplots_adjust(hspace=0, right=0.87) # figLCFS.suptitle('shot {s}, {avname} over time {fr}-{t}s {tb}: {ind}' title = str('shot {s}, {avname} filter over time {fr}-{t}s {tb}: {ind}' .format(s=[da['date'][0], da['progId'][0]], fr=t_range[0],t=t_range[1], tb=tb,
dev_name = 'W7X' dev = pyfusion.getDevice(dev_name) shot = da['info']['shotdata']['shot'][0] chan = da['info']['channels'][0] echdata = dev.acq.getdata(shot, 'W7X_TotECH') wech = np.where(echdata.signal > 100)[0] tech = echdata.timebase[wech[0]] t0_utc = int(tech * 1e9) + echdata.utc[0] dtprobe = (da['info']['params']['i_diag_utc'][0] - t0_utc) / 1e9 tslice = da['t_mid'][s - step // 2] + dtprobe dtaverage = da['t_mid'][s] - da['t_mid'][s - step] dtstr = '{dt:.2f}s avg'.format(dt=dtaverage) # axu.set_title('{fn}, time={t:.4f} {dt}'.format(fn=da.name, t=tslice, dt=dtstr)) tit = axu.get_title() cfg, cfg_dict = get_mag_config(shot) if cfg is None: ind = '' else: ind = 13 - 12 * cfg[2] / .3904 if abs(ind - round(ind)) < .03: ind = int(round((ind))) ind = 'Ind {ind}'.format(ind=ind) axu.set_title(tit + '\nProgram Id {s}, time={t:.4f} {ind}' # {dt}' .format(s=shot, fn=da.name, t=tslice, dt=dtstr, ind=ind)) if cfg_dict is not None: plot_trim(axu, mag_dict=cfg_dict, color='g', aspect=1.2) if time_strip_h > 0: