def plot_heat_power_density(self): ''' Function for plotting the heat flux profile Returns ------- None. ''' plt.plot(self._s, self._q) plt.xlabel('$s$') plt.ylabel('$q(s)$') imageFile = getOption('imageFile') if imageFile: plt.savefig(imageFile) else: plt.show()
FOOTPRINT.fx_in_out = 5. FOOTPRINT.calculate_heat_flux_density("lfs") Q_PARALLEL = np.array(FOOTPRINT._q) X_AFTER_LCFS = np.array(FOOTPRINT.get_global_coordinates()) EQUILIBRIUM = {} plt.figure() for i in X_AFTER_LCFS: P_0 = [i, 0, 0] FIELD_LINE_DICT = FIELD_LINE.follow_field_in_plane(p_0=P_0, max_length=15.0) plt.plot(FIELD_LINE_DICT['R'], FIELD_LINE_DICT['Z']) EQUILIBRIUM[i] = FIELD_LINE_DICT # X_AFTER_LCFS = HESELDATA.hesel_params.xaxis[I_AFTER_LCFS] + MID_PLANE_LOC DIVERTOR_COORDS = np.array((np.array([0.375, 0.675]), np.array([-0.78, -0.885]))) # EQUILIBRIUM = 'eq_0002' HEAT_FLUX_AT_OMP = np.array((X_AFTER_LCFS, Q_PARALLEL)) MAP_DICT = project_field_lines(X_AFTER_LCFS, DIVERTOR_COORDS, FIESTA) R_DIV = np.array([MAP_DICT[i]["R_pos"] for i in X_AFTER_LCFS]) F_X = np.array([MAP_DICT[i]["f_x"] for i in X_AFTER_LCFS]) ANGLES = np.array([MAP_DICT[i]["alpha"] for i in X_AFTER_LCFS]) FIG = plt.figure() plt.figure() plt.plot(R_DIV, Q_PARALLEL*X_AFTER_LCFS/(R_DIV*F_X/np.cos(ANGLES))) imageFile = getOption('imageFile') if imageFile : plt.savefig(imageFile) else : plt.show()