diff --git a/scripts/ns2d.strat/coeff_diss.py b/scripts/ns2d.strat/coeff_diss.py index 62a57be42a6d161df6bbc66be1995dae98bbf311_c2NyaXB0cy9uczJkLnN0cmF0L2NvZWZmX2Rpc3MucHk=..b02d9dd8090e8153264c14c1def457cf5ac604c2_c2NyaXB0cy9uczJkLnN0cmF0L2NvZWZmX2Rpc3MucHk= 100644 --- a/scripts/ns2d.strat/coeff_diss.py +++ b/scripts/ns2d.strat/coeff_diss.py @@ -34,7 +34,7 @@ ### PARAMETERS ### # CONDITIONS -nb_wavenumbers_y = 14 +nb_wavenumbers_y = 16 threshold_ratio = 1e1 min_factor = 0.7 @@ -263,5 +263,5 @@ gamma_not = make_float_value_for_path(gamma) # Create directory in path -path_root = Path("/fsnet/project/meige/2015/15DELDUCA/DataSim/Coef_Diss") +path_root = "/fsnet/project/meige/2015/15DELDUCA/DataSim/Coef_Diss" name_directory = "Coef_Diss_gamma{}".format(gamma_not) @@ -267,5 +267,5 @@ name_directory = "Coef_Diss_gamma{}".format(gamma_not) -path = path_root / name_directory +path = os.path.join(path_root, name_directory) if mpi.rank == 0 and not os.path.exists(path): os.mkdir(path) @@ -278,6 +278,12 @@ PLOT_FIGURES = PLOT_FIGURES and mpi.rank == 0 +# Write temporal results in .txt file +path_file_write2 = os.path.join(path, "results_check_nu8.txt") +if mpi.rank == 0: + if os.path.exists(path_file_write2): + os.remove(path_file_write2) + if len(paths_sim) == 0: # Make FIRST SIMULATION params = make_parameters_simulation(gamma, F, sigma, nu_8, t_end=t_end, NO_SHEAR_MODES=NO_SHEAR_MODES) @@ -328,6 +334,17 @@ energies.append(energy) viscosities.append(sim.params.nu_8) +# Write results into temporal file +if mpi.rank == 0: + to_print = ("####\n" + "t = {:.4e} \n" + "E = {:.4e} \n" + "nu8 = {:.4e} \n" + "factor = {:.4e} \n").format( + time_total, energy, sim.params.nu_8, 1) + + write_to_file(path_file_write2, to_print, mode="w") + # Compute the data spectra energy budget kxE, kyE, dissE_kx, dissE_ky = load_mean_spect_energy_budg(sim, tmin=2., tmax=1000) idx_diss_max, idy_diss_max = _compute_ikdiss(dissE_kx, dissE_ky) @@ -427,7 +444,8 @@ viscosities.append(sim.params.nu_8) # Computes new index k_max_dissipation - kxE, kyE, dissE_kx, dissE_ky = load_mean_spect_energy_budg(sim, tmin=2., tmax=1000) + kxE, kyE, dissE_kx, dissE_ky = load_mean_spect_energy_budg( + sim, tmin=2., tmax=1000) idx_diss_max, idy_diss_max = _compute_ikdiss(dissE_kx, dissE_ky) idx_dealiasing, idy_dealiasing = _compute_ikmax(kxE, kyE) @@ -431,6 +449,16 @@ idx_diss_max, idy_diss_max = _compute_ikdiss(dissE_kx, dissE_ky) idx_dealiasing, idy_dealiasing = _compute_ikmax(kxE, kyE) + # Write results into temporal file + to_print = ("####\n" + "t = {:.4e} \n" + "E = {:.4e} \n" + "nu8 = {:.4e} \n" + "factor = {:.4e} \n").format( + time_total, energy, sim.params.nu_8, factor) + + write_to_file(path_file_write2, to_print, mode="a") + if PLOT_FIGURES: ax.plot(kyE, dissE_ky, label="nu8 = {:.2e}, diff = {}".format( params.nu_8, abs(idy_diss_max - idy_dealiasing)))