awage · awage · May 29, 2025 · May 29, 2025 · May 30, 2025 · Jul 13, 2026
diff --git a/plots/4d_sprott_memristor_g=0.03_m=11_r=5.8_res=1200_α=1_β=0.05_γ=0.5.png b/plots/4d_sprott_memristor_g=0.03_m=11_r=5.8_res=1200_α=1_β=0.05_γ=0.5.png
diff --git a/plots/alfven_res=1200_ν=0.0175.png b/plots/alfven_res=1200_ν=0.0175.png
diff --git a/plots/bairstow_map_a=0.8_res=1200.png b/plots/bairstow_map_a=0.8_res=1200.png
diff --git a/plots/basin_sommerer_res=1200.png b/plots/basin_sommerer_res=1200.png
diff --git a/plots/basins_CO2_laser_f=0.005_m=0.02_res=1200.png b/plots/basins_CO2_laser_f=0.005_m=0.02_res=1200.png
diff --git a/plots/basins_economic_res=1500_Φ=0.085.png b/plots/basins_economic_res=1500_Φ=0.085.png
diff --git a/plots/basins_kur_ni=14_res=300.png b/plots/basins_kur_ni=14_res=300.png
diff --git a/plots/basins_lasers_f=70200.0_m=0.8_res=1200.png b/plots/basins_lasers_f=70200.0_m=0.8_res=1200.png
diff --git a/plots/basins_multi_rand_net_N=100_g=2.5_res=1200_s=4.5.png b/plots/basins_multi_rand_net_N=100_g=2.5_res=1200_s=4.5.png
diff --git a/plots/basins_riddle_ott_res=1200.png b/plots/basins_riddle_ott_res=1200.png
diff --git a/plots/basins_tentacles_N=40_res=1200.png b/plots/basins_tentacles_N=40_res=1200.png
diff --git a/plots/basins_tlns_res=1200_θ=1.0.png b/plots/basins_tlns_res=1200_θ=1.0.png
diff --git a/plots/bell_yoke_Tmax=150_lr=-0.03_res=700.png b/plots/bell_yoke_Tmax=150_lr=-0.03_res=700.png
diff --git a/plots/biorhythm_ks=1.99_res=1200.png b/plots/biorhythm_ks=1.99_res=1200.png
diff --git a/plots/black_holes_res=1200_Δpϕ=0.03.png b/plots/black_holes_res=1200_Δpϕ=0.03.png
diff --git a/plots/bogdanov_k=1.2_res=2000_ε=0.0125_μ=-0.1.png b/plots/bogdanov_k=1.2_res=2000_ε=0.0125_μ=-0.1.png
diff --git a/plots/carpet_res=1200.png b/plots/carpet_res=1200.png
diff --git a/plots/cbic_map_j=0.66_res=1200_μ=2.9.png b/plots/cbic_map_j=0.66_res=1200_μ=2.9.png
diff --git a/plots/chen_S=0.77_res=1200.png b/plots/chen_S=0.77_res=1200.png
diff --git a/plots/chimera_states_res=1200_β=0.025_μ=0.55_ν=0.45.png b/plots/chimera_states_res=1200_β=0.025_μ=0.55_ν=0.45.png
diff --git a/plots/cold_atoms_res=1200_vx=0.1_x=-500.png b/plots/cold_atoms_res=1200_vx=0.1_x=-500.png
diff --git a/plots/comptetion_model_res=1200.png b/plots/comptetion_model_res=1200.png
diff --git a/plots/coupled_HR_res=1200.png b/plots/coupled_HR_res=1200.png
diff --git a/plots/coupled_lorenz_res=1200_α=10_β=24.8_γ=2.67_ε=1.1.png b/plots/coupled_lorenz_res=1200_α=10_β=24.8_γ=2.67_ε=1.1.png
diff --git a/plots/coupled_rulkov_res=1200_Δ=0.0_α=4.1_γ1=-1.75_γ2=-1.75_σ=0.035.png b/plots/coupled_rulkov_res=1200_Δ=0.0_α=4.1_γ1=-1.75_γ2=-1.75_σ=0.035.png
diff --git a/plots/cournot2d_res=1200_μ1=3.55_μ2=3.55.png b/plots/cournot2d_res=1200_μ1=3.55_μ2=3.55.png
diff --git a/plots/cpashwin_res=1200_α=0.7_ε=0.5_ν=1.28.png b/plots/cpashwin_res=1200_α=0.7_ε=0.5_ν=1.28.png
diff --git a/plots/cplog_a=3.6_res=1200_ε=-1.0.png b/plots/cplog_a=3.6_res=1200_ε=-1.0.png
diff --git a/plots/cyc_comp_r=3.35_res=1200.png b/plots/cyc_comp_r=3.35_res=1200.png
diff --git a/plots/dsnm_a=0.55_b=0.45_res=2000.png b/plots/dsnm_a=0.55_b=0.45_res=2000.png
diff --git a/plots/duffing_0.05_F=0.098_d=0.05_res=1200_ω=1.15.png b/plots/duffing_0.05_F=0.098_d=0.05_res=1200_ω=1.15.png
diff --git a/plots/earth_reversal_res=1200_Γ=0.9_μ=0.119_ν=0.1.png b/plots/earth_reversal_res=1200_Γ=0.9_μ=0.119_ν=0.1.png
diff --git a/plots/eckhardt_Re=425.0_res=1200.png b/plots/eckhardt_Re=425.0_res=1200.png
diff --git a/plots/feudel_a=0.0015_r=3.83_res=1200.png b/plots/feudel_a=0.0015_r=3.83_res=1200.png
diff --git a/plots/gear_rattle2_res=1200.png b/plots/gear_rattle2_res=1200.png
diff --git a/plots/gear_rattle_res=1200_β=0.6_δ=0.6_ε=0.1.png b/plots/gear_rattle_res=1200_β=0.6_δ=0.6_ε=0.1.png
diff --git a/plots/grebogi_res=1200.png b/plots/grebogi_res=1200.png
diff --git a/plots/gyrostat_res=1200.png b/plots/gyrostat_res=1200.png
diff --git a/plots/henon_j=0.9_res=1200_μ=1.08.png b/plots/henon_j=0.9_res=1200_μ=1.08.png
diff --git a/plots/hh_E=0.25_res=1200.png b/plots/hh_E=0.25_res=1200.png
diff --git a/plots/hidden_chua_m0=-0.177_m1=-1.15_res=1200_α=8.46_β=12.1_γ=0.00516.png b/plots/hidden_chua_m0=-0.177_m1=-1.15_res=1200_α=8.46_β=12.1_γ=0.00516.png
diff --git a/plots/hopfield_neuron_b1=3.0_b2=3.0_res=1200_α=-0.1_β=0.1_δ=1.5.png b/plots/hopfield_neuron_b1=3.0_b2=3.0_res=1200_α=-0.1_β=0.1_δ=1.5.png
diff --git a/plots/impacting_res=1200.png b/plots/impacting_res=1200.png
diff --git a/plots/josephson_idc=1.88_irf=10.2_res=1200_Ω=1_β=25.png b/plots/josephson_idc=1.88_irf=10.2_res=1200_Ω=1_β=25.png
diff --git a/plots/kapitaniak_d1=-0.935_d2=-0.935_l=1.41_pp=-1.41_res=1200.png b/plots/kapitaniak_d1=-0.935_d2=-0.935_l=1.41_pp=-1.41_res=1200.png
diff --git a/plots/kicked_rotor_4d_res=1200_ρ=6.5.png b/plots/kicked_rotor_4d_res=1200_ρ=6.5.png
diff --git a/plots/kicked_rotor_f0=4.0_res=1200_ν=0.02.png b/plots/kicked_rotor_f0=4.0_res=1200_ν=0.02.png
diff --git a/plots/kicked_rotor_ott_f0=0.1_res=1200_ν=0.1.png b/plots/kicked_rotor_ott_f0=0.1_res=1200_ν=0.1.png
diff --git a/plots/lai_pcw_res=1200.png b/plots/lai_pcw_res=1200.png
diff --git a/plots/lebovitz_Re=307.0_res=1200.png b/plots/lebovitz_Re=307.0_res=1200.png
diff --git a/plots/li_sprott_a=13.0_b=0.55_res=1200.png b/plots/li_sprott_a=13.0_b=0.55_res=1200.png
diff --git a/plots/lorenz84_F=6.85_G=1.29_a=0.25_b=4.0_res=1200.png b/plots/lorenz84_F=6.85_G=1.29_a=0.25_b=4.0_res=1200.png
diff --git a/plots/lorenz_complexity_res=1200.png b/plots/lorenz_complexity_res=1200.png
diff --git a/plots/mag_pend_N=3_d=0.3_res=1200_α=0.2_γ=1_ω=0.5.png b/plots/mag_pend_N=3_d=0.3_res=1200_α=0.2_γ=1_ω=0.5.png
diff --git a/plots/matryoshka_res=1200.png b/plots/matryoshka_res=1200.png
diff --git a/plots/megastable_res=1200.png b/plots/megastable_res=1200.png
diff --git a/plots/mira_map_a=-0.42_b=-1.32_res=1200.png b/plots/mira_map_a=-0.42_b=-1.32_res=1200.png
diff --git a/plots/nash_equilibrium_res=1200_λ=0.5_μ=1.95.png b/plots/nash_equilibrium_res=1200_λ=0.5_μ=1.95.png
diff --git a/plots/neuron_synapse_B=2_c=1.8_g=1.7_res=1200.png b/plots/neuron_synapse_B=2_c=1.8_g=1.7_res=1200.png
diff --git a/plots/newton__res=1200.png b/plots/newton__res=1200.png
diff --git a/plots/open_disks_res=1200_ww=0.001.png b/plots/open_disks_res=1200_ww=0.001.png
diff --git a/plots/open_sinai_res=1200_Δ=0.8.png b/plots/open_sinai_res=1200_Δ=0.8.png
diff --git a/plots/parametric_pendulum_A_1=0.0_A_2=4.1_res=1200_ω=1.5.png b/plots/parametric_pendulum_A_1=0.0_A_2=4.1_res=1200_ω=1.5.png
diff --git a/plots/pendulum_F=1.36_d=0.2_res=1200_ω=0.5.png b/plots/pendulum_F=1.36_d=0.2_res=1200_ω=0.5.png
diff --git a/plots/pendulum_F=1.66_d=0.2_res=1200_ω=1.0.png b/plots/pendulum_F=1.66_d=0.2_res=1200_ω=1.0.png
diff --git a/plots/photonic_coupler_k=5_res=1200_α=2_β=1.5_ε=0.5.png b/plots/photonic_coupler_k=5_res=1200_α=2_β=1.5_ε=0.5.png
diff --git a/plots/poon_disks_res=1200_ww=0.1.png b/plots/poon_disks_res=1200_ww=0.1.png
diff --git a/plots/pred_prey_k=7.94_m=0.104_res=1200.png b/plots/pred_prey_k=7.94_m=0.104_res=1200.png
diff --git a/plots/ricker_gatto_res=2000.png b/plots/ricker_gatto_res=2000.png
diff --git a/plots/rikitake_res=1200_α=1.0_μ=0.5.png b/plots/rikitake_res=1200_α=1.0_μ=0.5.png
diff --git a/plots/shrimali_res=1200_α=3.25_β=0.01_ϵp=0.5.png b/plots/shrimali_res=1200_α=3.25_β=0.01_ϵp=0.5.png
diff --git a/plots/slim_fractal_res=1200_μ=0.2.png b/plots/slim_fractal_res=1200_μ=0.2.png
diff --git a/plots/split_ring_res_res=1200_μ=35_ω=1.03.png b/plots/split_ring_res_res=1200_μ=35_ω=1.03.png
diff --git a/plots/sporadical_res=1200.png b/plots/sporadical_res=1200.png
diff --git a/plots/thomas_b=0.166_res=1200.png b/plots/thomas_b=0.166_res=1200.png
diff --git a/plots/triangular_img_res=3000_λ=1.03.png b/plots/triangular_img_res=3000_λ=1.03.png
diff --git a/scripts/4d_kicked.jl b/scripts/4d_kicked.jl
@@ -56,4 +56,4 @@ L = W1*exp(λ1*T) + W2*exp(λ2*T)
 M = W1*(exp(λ1*T)-1)/λ1 + W2*(exp(λ2*T) -1)/λ2
 # res = 400
 params = @strdict res  M L  ρ
-print_fig(params, "kicked_rotor_4d", compute_kicked_rotor_4d; ylab = L"\theta_2", xlab = L"\theta_1", force = false)
+print_fig(params, "kicked_rotor_4d", compute_kicked_rotor_4d; ylab = L"\theta_2", xlab = L"\theta_1", force = false, cmap)
diff --git a/scripts/HR_extreme_events.jl b/scripts/HR_extreme_events.jl
@@ -55,7 +55,7 @@ end
 
 res = 1200
 params = @strdict res 
-cmap = ColorScheme([ RGB(0.9,0.2,0.1), RGB(1,1,1), RGB(0,0,0) ] )
+# cmap = ColorScheme([ RGB(0.9,0.2,0.1), RGB(1,1,1), RGB(0,0,0) ] )
 print_fig(params, "coupled_HR", compute_HR; force = false, xlab = L"z_1", ylab = L"z_2", cmap) 
 # att = get_att(params, "coupled_HR", compute_HR; force = false) 
 

diff --git a/scripts/LuDP_gene_switch.jl b/scripts/LuDP_gene_switch.jl
@@ -57,4 +57,4 @@ end
 
 res = 400
     params = @strdict res
-    print_fig(params, "ludp", compute_basins_ludp; ylab = L"y", xlab = L"x", force = false) 
+    print_fig(params, "ludp", compute_basins_ludp; ylab = L"y", xlab = L"x", force = false, cmap) 
diff --git a/scripts/TLNs_network.jl b/scripts/TLNs_network.jl
@@ -56,6 +56,6 @@ end
 let res = 1200
 θ = 1.; 
 params = @strdict res θ
-print_fig(params, "basins_tlns", compute_basins_tlns; force = false, xlab = L"x_6", ylab = L"x_7")
+print_fig(params, "basins_tlns", compute_basins_tlns; force = false, xlab = L"x_6", ylab = L"x_7", cmap)
 end
 
diff --git a/scripts/aguirre_billiard.jl b/scripts/aguirre_billiard.jl
@@ -68,4 +68,4 @@ end
 
 ww = 0.001; #res = 400
 params = @strdict ww res
-print_fig(params, "open_disks", compute_exit_three_disk; xlab = L"x_0", ylab = L"\theta_0")
+print_fig(params, "open_disks", compute_exit_three_disk; xlab = L"x_0", ylab = L"\theta_0", cmap)
diff --git a/scripts/alfven_complexity.jl b/scripts/alfven_complexity.jl
@@ -38,7 +38,7 @@ end
 let res = 1200
 ν = 0.01747; 
 params = @strdict ν res
-cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(1,0.26, 0.26), RGB(0.34,0.34,1), RGB(0.1,0.1,0.1) ] )
+# cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(1,0.26, 0.26), RGB(0.34,0.34,1), RGB(0.1,0.1,0.1) ] )
 print_fig(params, "alfven", compute_alfven; ylab = L"b_z", xlab = L"b_y", force = false, cmap) 
 end
 
diff --git a/scripts/bell_yoke2.jl b/scripts/bell_yoke2.jl
@@ -94,7 +94,7 @@ end
 
 res = 700; Tmax = 150; lr = -0.03
 params = @strdict res Tmax lr
-cmap = ColorScheme([RGB(1,1,1),  RGB(0.9,0.2,0.1), RGB(0.2,0.2,0.2)] )
+# cmap = ColorScheme([RGB(1,1,1),  RGB(0.9,0.2,0.1), RGB(0.2,0.2,0.2)] )
 print_fig(params, "bell_yoke", compute_bell_yoke; ylab= L"\phi_1", xlab= L"\phi_2", force = false, cmap)
 att = get_att(params, "bell_yoke", compute_bell_yoke)
 

diff --git a/scripts/black_holes.jl b/scripts/black_holes.jl
@@ -126,5 +126,5 @@ end
 
 Δpϕ = 0.03; #res = 500
 params = @strdict Δpϕ res
-    cmap = ColorScheme([RGB(1,1,1), RGB(0.99,0.4,0.4), RGB(0., .9, .7),  RGB(0.2,0.2,0.2)] )
+    # cmap = ColorScheme([RGB(1,1,1), RGB(0.99,0.4,0.4), RGB(0., .9, .7),  RGB(0.2,0.2,0.2)] )
 print_fig(params, "black_holes", compute_BH; xlab = L"\rho", ylab = L"z", force = false, cmap)
diff --git a/scripts/buck_converter.jl b/scripts/buck_converter.jl
@@ -49,7 +49,7 @@ end
 
 res = 1200; Vin = 30.10
 params = @strdict res Vin
-cmap = ColorScheme([RGB(1,1,1),  RGB(0.7,0.7,0.7), RGB(0.3,0.3,0.3)] )
+# cmap = ColorScheme([RGB(1,1,1),  RGB(0.7,0.7,0.7), RGB(0.3,0.3,0.3)] )
 print_fig(params, "buck", compute_buck; ylab= L"V", xlab= L"i", force = false, cmap)
 
 
diff --git a/scripts/carpet_oscillator.jl b/scripts/carpet_oscillator.jl
@@ -38,6 +38,6 @@ end
 
 res = 1200
 params = @strdict  res
-print_fig(params, "carpet", compute_carpet; force = false)
+print_fig(params, "carpet", compute_carpet; force = false, cmap)
 
 
diff --git a/scripts/chimera_reduced_model.jl b/scripts/chimera_reduced_model.jl
@@ -59,5 +59,5 @@ res = 1200;
 A = 0.1; β = 0.025
 μ = (A+1)/2; ν = 1 - μ
 params = @strdict res μ ν β 
-print_fig(params, "chimera_states", compute_chimera; force = false, xlab = L"s", ylab = L"\psi") 
-att = get_att(params, "chimera_states", compute_chimera) 
+print_fig(params, "chimera_states", compute_chimera; force = false, xlab = L"s", ylab = L"\psi", cmap) 
+# att = get_att(params, "chimera_states", compute_chimera) 
diff --git a/scripts/co2_modulated_laser.jl b/scripts/co2_modulated_laser.jl
@@ -56,7 +56,7 @@ end
 let res = 1200
     m = 0.02; f = 0.005
     params = @strdict res m f
-    cmap = ColorScheme([RGB(1,1,1),  RGB(0.9,0.2,0.1)] )
+    # cmap = ColorScheme([RGB(1,1,1),  RGB(0.9,0.2,0.1)] )
     print_fig(params, "basins_CO2_laser", compute_basins_CO2_lasers; force = false, xlab = L"\log x", ylab = L"y", cmap)
 end
 
diff --git a/scripts/cold_atoms.jl b/scripts/cold_atoms.jl
@@ -85,5 +85,5 @@ end
 
 x = -500; vx = 0.1; #res = 300
 params = @strdict x vx res
-cmap = ColorScheme([RGB(1,1,1), RGB(0.55,0.9,0.35), RGB(0.9,0.4,0.1),  RGB(0.50,0.24,1)] )
+# cmap = ColorScheme([RGB(1,1,1), RGB(0.55,0.9,0.35), RGB(0.9,0.4,0.1),  RGB(0.50,0.24,1)] )
 print_fig(params, "cold_atoms", compute_cold_atoms; ylab = L"v_y", xlab = L"y", cmap)
diff --git a/scripts/competition_ecology_plancton.jl b/scripts/competition_ecology_plancton.jl
@@ -60,5 +60,5 @@ end
 
 let res = 1200
 params = @strdict res
-print_fig(params, "comptetion_model", compute_competition_model; xlab = L"N_2", ylab = L"N_4", force = false)
+print_fig(params, "comptetion_model", compute_competition_model; xlab = L"N_2", ylab = L"N_4", force = false, cmap)
 end
diff --git a/scripts/coupled_hopfield_neuron.jl b/scripts/coupled_hopfield_neuron.jl
@@ -50,6 +50,6 @@ end
 res = 1200
 δ = 1.5; α = -0.1; b1 = 3.0; b2 = 3.0; β = -α
 params = @strdict res δ α b1 b2 β
-cmap = :jet
+# cmap = :jet
 print_fig(params, "hopfield_neuron", compute_hpf_cpld; force = false, xlab = L"x", ylab = L"y", cmap) 
 # att = get_att(params, "hopfield_neuron", compute_hpf_cpld; force = false) 
diff --git a/scripts/coupled_lorenz.jl b/scripts/coupled_lorenz.jl
@@ -44,7 +44,7 @@ end
 let res = 1200
 α = 10; β = 24.76; ε = 1.1; γ = 8/3; 
 params = @strdict res α β ε γ
-cmap = ColorScheme([RGB(1,1,1),  RGB(0.9,0.2,0.1)] )
+# cmap = ColorScheme([RGB(1,1,1),  RGB(0.9,0.2,0.1)] )
 print_fig(params, "coupled_lorenz", compute_lorenz; force = false, xlab = L"z_1", ylab = L"z_2", cmap) 
 # att = get_att(params, "coupled_lorenz", compute_lorenz; force = false) 
 # @show att 

diff --git a/scripts/coupled_rulkov_neurons.jl b/scripts/coupled_rulkov_neurons.jl
@@ -36,7 +36,7 @@ end
 res = 1200
 Δ = 0.0; γ1 = -1.75; γ2 = γ1 + Δ; α = 4.1 ; σ = 0.035
 params = @strdict res γ1 γ2  α Δ σ
-cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(0.34,0.34,1), RGB(1,0.46,0.46), RGB(0.1,0.1,0.1) ] )
+# cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(0.34,0.34,1), RGB(1,0.46,0.46), RGB(0.1,0.1,0.1) ] )
 print_fig(params, "coupled_rulkov", compute_cpld_rulkov; force = false, cmap) 
 
 
diff --git a/scripts/cuencas_hh.jl b/scripts/cuencas_hh.jl
@@ -73,7 +73,7 @@ res = 1200
 E = 0.25; #res = 800
 params = @strdict E res
 # cmap = ColorScheme([RGB(1,1,1), RGB(0.55,0.9,0.35), RGB(0.1,0.4,0.1), RGB(0.25,0.25, 0.25), RGB(0.50,0.24,1)] )
-cmap = ColorScheme([RGB(1,1,1), RGB(0.55,0.9,0.35), RGB(0.9,0.4,0.1),  RGB(0.50,0.24,1)] )
+# cmap = ColorScheme([RGB(1,1,1), RGB(0.55,0.9,0.35), RGB(0.9,0.4,0.1),  RGB(0.50,0.24,1)] )
 print_fig(params, "hh", comp_basins_hh; cmap)
 
 
diff --git a/scripts/cyclic_competition.jl b/scripts/cyclic_competition.jl
@@ -60,6 +60,6 @@ end
 let res = 1200
 r = 3.35; #res = 1500;
 params = @strdict r res
-cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(1,0.16, 0.16), RGB(0.14,0.14,1), RGB(0.8,0.8,0.1), RGB(0.1,0.1,0.1) ] )
+# cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(1,0.16, 0.16), RGB(0.14,0.14,1), RGB(0.8,0.8,0.1), RGB(0.1,0.1,0.1) ] )
 print_fig(params , "cyc_comp", compute_cyc_comp; xlab = L"x", ylab = L"y", cmap) 
 end
diff --git a/scripts/decroly_biorhythm.jl b/scripts/decroly_biorhythm.jl
@@ -45,7 +45,7 @@ res = 1200;
 # for ks in 2.0:0.05:2.5
 ks = 1.99
 params = @strdict  ks res
-cmap = ColorScheme([RGB(1,1,1),  RGB(0.9,0.2,0.1)] )
+# cmap = ColorScheme([RGB(1,1,1),  RGB(0.9,0.2,0.1)] )
 print_fig(params, "biorhythm", compute_biorhythm; ylab = L"y", xlab = L"x", force = false, cmap)
 # att = get_att(params, "biorhythm", compute_biorhythm)
 # end

diff --git a/scripts/dsnm.jl b/scripts/dsnm.jl
@@ -43,5 +43,5 @@ end
 
 a = 0.55; b = 0.45; #res = 1000
 params = @strdict a b res
-print_fig(params, "dsnm", compute_dsnm)
+print_fig(params, "dsnm", compute_dsnm; cmap)
 
diff --git a/scripts/duffing_basins.jl b/scripts/duffing_basins.jl
@@ -30,7 +30,7 @@ function compute_basins_duffing(di::Dict)
 end
 
 
-cmap = ColorScheme([RGB(0,0,0), RGB(1,1,1)] )
+# cmap = ColorScheme([RGB(0,0,0), RGB(1,1,1)] )
 
 # res = 600; 
 # d = 0.1; F=0.1; ω=0.1;  # smooth boundary
@@ -44,5 +44,5 @@ cmap = ColorScheme([RGB(0,0,0), RGB(1,1,1)] )
 let
     d = 0.05; F=0.098; ω=1.15;  # smooth boundary
     params = @strdict d F ω res
-    print_fig(params, string("duffing_",d), compute_basins_duffing; ylab = L"$\dot{x}$", xlab = L"x") 
+    print_fig(params, string("duffing_",d), compute_basins_duffing; ylab = L"$\dot{x}$", xlab = L"x", cmap) 
 end
diff --git a/scripts/duffing_parametrical_cos.jl b/scripts/duffing_parametrical_cos.jl
@@ -54,4 +54,4 @@ end
 
 
 a = b = 1; c = 0.2; ω=1.;  # smooth boundary
-print_fig(600, 600, a, b, c, ω, 600) 
+print_fig(600, 600, a, b, c, ω, 600; cmap) 
diff --git a/scripts/earth_field_reversal.jl b/scripts/earth_field_reversal.jl
@@ -42,6 +42,6 @@ end
 let
     μ = 0.1193; ν = 0.1; Γ = 0.9; res = 1200
     params = @strdict  res μ ν Γ 
-    cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(0.7,0.7,0.7), RGB(1,0,0)] )
+    # cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(0.7,0.7,0.7), RGB(1,0,0)] )
     print_fig(params, "earth_reversal", compute_earth_reversal; cmap, xlab = L"Q", ylab = L"D", force = false)
 end
diff --git a/scripts/eckhardt_9D.jl b/scripts/eckhardt_9D.jl
@@ -90,5 +90,5 @@ end
 
 Re = 425.; res = 1200
 params = @strdict res Re
-print_fig(params, "eckhardt", compute_E9D; ylab = L"$V_{att}-V_{att,sym}$", xlab = L"V_{att}+V_{att,sym}") 
+print_fig(params, "eckhardt", compute_E9D; ylab = L"$V_{att}-V_{att,sym}$", xlab = L"V_{att}+V_{att,sym}", cmap) 
 
diff --git a/scripts/economic_geographic_model.jl b/scripts/economic_geographic_model.jl
@@ -90,6 +90,6 @@ end
 let res = 1500
     Φ = 0.085
     params = @strdict res Φ
-    print_fig(params, "basins_economic", compute_basins_economic; force = false, xlab = L"\lambda_1", ylab = L"\lambda_2")
+    print_fig(params, "basins_economic", compute_basins_economic; force = false, xlab = L"\lambda_1", ylab = L"\lambda_2", cmap)
 end
 
diff --git a/scripts/gear_rattle.jl b/scripts/gear_rattle.jl
@@ -57,7 +57,7 @@ end
 
 β = 0.6; δ = 0.6; ε = 0.1; res = 1200
 params = @strdict β δ ε res
-print_fig(params, "gear_rattle", compute_gear_rattle; ylab= L"\dot{\theta}", xlab= L"\theta", force = false)
+print_fig(params, "gear_rattle", compute_gear_rattle; ylab= L"\dot{\theta}", xlab= L"\theta", force = false, cmap)
 # att = get_att(params, "gear_rattle", compute_gear_rattle)
 
 # function gear_rattle!(du,u, p, t)

diff --git a/scripts/gear_rattle_souza.jl b/scripts/gear_rattle_souza.jl
@@ -51,7 +51,7 @@ end
 
 res = 1200
 params = @strdict res
-cmap = ColorScheme([RGB(1,1,1),  RGB(0.9,0.2,0.1)] )
+# cmap = ColorScheme([RGB(1,1,1),  RGB(0.9,0.2,0.1)] )
 print_fig(params, "gear_rattle2", compute_gear_rattle2; ylab= L"\dot{s}", xlab= L"s", force = false, cmap)
 # att = get_att(params, "gear_rattle2", compute_gear_rattle2)
 

diff --git a/scripts/gyrostat.jl b/scripts/gyrostat.jl
@@ -40,7 +40,7 @@ end
 
 res = 1200
 params = @strdict res
-cmap = ColorScheme([RGB(1,1,1),  RGB(0.9,0.25,0.2)] )
+# cmap = ColorScheme([RGB(1,1,1),  RGB(0.9,0.25,0.2)] )
 print_fig(params, "gyrostat", compute_gyrostat; ylab = L"x_2", xlab = L"x_1", force = false, cmap)
-att = get_att(params, "gyrostat", compute_gyrostat)
+# att = get_att(params, "gyrostat", compute_gyrostat)
 
diff --git a/scripts/hidden_chua.jl b/scripts/hidden_chua.jl
@@ -51,6 +51,6 @@ end
 # α = 8.5; β = 14.28; γ = 0; m0 = -8/7; m1 = -5/7
 α = 8.4562218418; β = 12.0732335925; γ = 0.0051631393; m0 = -0.1767573476; m1 = -1.1467573476; # res = 500
 params = @strdict α β γ m0 m1 res
-cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(1,0.46, 0.46), RGB(0.34,0.34,1), RGB(0.1,0.1,0.1) ] )
+# cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(1,0.46, 0.46), RGB(0.34,0.34,1), RGB(0.1,0.1,0.1) ] )
 print_fig(params, "hidden_chua", compute_basins_chua; force = false, cmap)
 
diff --git a/scripts/josephson_junction.jl b/scripts/josephson_junction.jl
@@ -49,5 +49,5 @@ end
 
 β = 25; idc = 1.878; irf = 10.198; Ω = 1; #res = 700
 params = @strdict β idc irf Ω res
-cmap = ColorScheme([RGB(1,1,1),  RGB(0.9,0.4,0.1),  RGB(0.50,0.24,1)] )
+# cmap = ColorScheme([RGB(1,1,1),  RGB(0.9,0.4,0.1),  RGB(0.50,0.24,1)] )
 print_fig(params, "josephson", compute_basins_josephson; ylab = L"\dot{\phi}", xlab = L"\phi", force = false, cmap)
diff --git a/scripts/kicked_ott.jl b/scripts/kicked_ott.jl
@@ -64,5 +64,5 @@ end
 # res = 200
 f0 = 0.1; ν = 0.1
 params = @strdict f0 ν res
-cmap = ColorScheme([RGB(1,1,1), RGB(1,0,0), RGB(0.7,0.7,0.7), RGB(1,0,0)] )
+# cmap = ColorScheme([RGB(1,1,1), RGB(1,0,0), RGB(0.7,0.7,0.7), RGB(1,0,0)] )
 print_fig(params, "kicked_rotor_ott", compute_kicked_rotor_ott; xlab = L"\theta_1", ylab = L"\theta_2", force = false, cmap)
diff --git a/scripts/kicked_rotor.jl b/scripts/kicked_rotor.jl
@@ -40,5 +40,5 @@ end
 f0 = 4.
 ν = 0.02
 params = @strdict f0 ν res
-print_fig(params, "kicked_rotor", compute_kicked_rotor; ylab = L"y", xlab = L"x", cmap = :flag)
+print_fig(params, "kicked_rotor", compute_kicked_rotor; ylab = L"y", xlab = L"x", cmap)
 # print_fig(600,600, nothing, f0, ν, 400) 
diff --git a/scripts/kur_halekotte.jl b/scripts/kur_halekotte.jl
@@ -94,7 +94,7 @@ let res = 300
 # print_fig(params, "basins_kur", compute_kur_halekotte; xlab = L"\phi", ylab = L"\omega") 
 ni = 14
 params = @strdict ni res
-print_fig(params, "basins_kur", compute_kur_halekotte; xlab = L"\Delta\phi", ylab = L"\omega", force = false) 
+print_fig(params, "basins_kur", compute_kur_halekotte; xlab = L"\Delta\phi", ylab = L"\omega", force = false, cmap) 
 # ni = 58
 # params = @strdict ni res
 # print_fig(params, "basins_kur", compute_kur_halekotte; xlab = L"\phi", ylab = L"\omega") 

diff --git a/scripts/kuramoto_basins_with_tentacles.jl b/scripts/kuramoto_basins_with_tentacles.jl
@@ -86,7 +86,7 @@ end
 let res = 1200
     N = 40
     params = @strdict res N
-    print_fig(params, "basins_tentacles", compute_basin_tentacle; force = false, xlab = L"\alpha_1", ylab = L"\alpha_2")
+    print_fig(params, "basins_tentacles", compute_basin_tentacle; force = false, xlab = L"\alpha_1", ylab = L"\alpha_2", cmap)
 end
 
 # dd = compute_basin_tentacle_feats(40, 100)

diff --git a/scripts/lebovitz_mariotti.jl b/scripts/lebovitz_mariotti.jl
@@ -75,6 +75,6 @@ end
 # res = 450; 
 Re = 307.
 params = @strdict res Re
-cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(0.1,0.1,0.1), RGB(1,0.46, 0.46), RGB(0.34,0.34,1)] )
+# cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(0.1,0.1,0.1), RGB(1,0.46, 0.46), RGB(0.34,0.34,1)] )
 print_fig(params, "lebovitz", compute_LM; cmap)
 # att = get_att(params, "lebovitz", compute_LM)
diff --git a/scripts/li_sprott.jl b/scripts/li_sprott.jl
@@ -44,7 +44,7 @@ end
 
 a = 13.; b = 0.55; #res = 500
 params = @strdict res a b
-cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(0.7,0.7,0.7), RGB(1,0,0)] )
+# cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(0.7,0.7,0.7), RGB(1,0,0)] )
 print_fig(params, "li_sprott", compute_li_sprott; cmap) 
 
 
diff --git a/scripts/lorenz84.jl b/scripts/lorenz84.jl
@@ -39,5 +39,5 @@ end
 let
     F=6.846; G=1.287; a=0.25; b=4.; res = 1200
     params = @strdict F G a b res
-    print_fig(params, "lorenz84", compute_lorenz84; ylab = L"y", xlab = L"x", force = false)
+    print_fig(params, "lorenz84", compute_lorenz84; ylab = L"y", xlab = L"x", force = false, cmap)
 end
diff --git a/scripts/lorenz_computational_chaos.jl b/scripts/lorenz_computational_chaos.jl
@@ -50,7 +50,7 @@ end
 
 res = 1200
 params = @strdict res
-cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(1,0.46, 0.46), RGB(0.34,0.34,1), RGB(0.1,0.1,0.1) ] )
+# cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(1,0.46, 0.46), RGB(0.34,0.34,1), RGB(0.1,0.1,0.1) ] )
 print_fig(params, "lorenz_complexity", compute_lorenz_comp;  force = false, cmap)
 att = get_att(params, "lorenz_complexity", compute_lorenz_comp)
 
diff --git a/scripts/magnetic_pendulum.jl b/scripts/magnetic_pendulum.jl
@@ -72,4 +72,4 @@ end
 
 γ=1; d=0.3; α=0.2; ω=0.5; N=3; #res = 300
 params = @strdict γ d α ω N res
-print_fig(params, "mag_pend", compute_mag_pend; ylab = L"y", xlab = L"x", force = false)
+print_fig(params, "mag_pend", compute_mag_pend; ylab = L"y", xlab = L"x", force = false, cmap)
diff --git a/scripts/map_ashwin.jl b/scripts/map_ashwin.jl
@@ -36,7 +36,7 @@ end
 
 # a = 3.57480493875920; ε = -0.2
 ν = 1.28; α = 0.7; ε = 0.5; 
-cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(0.7,0.7,0.7), RGB(1,0,0)] )
+# cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(0.7,0.7,0.7), RGB(1,0,0)] )
 params = @strdict res α ε ν
 print_fig(params, "cpashwin", compute_ashwin; cmap) 
 
diff --git a/scripts/map_bairstow.jl b/scripts/map_bairstow.jl
@@ -37,6 +37,6 @@ res = 1200
 a = 0.8
 # a = 0.15
 params = @strdict res a 
-cmap = ColorScheme([RGB(1,1,1),  RGB(1,1,1),  RGB(1,1,1),  RGB(0.9,0.1,0.1)] )
+# cmap = ColorScheme([RGB(1,1,1),  RGB(1,1,1),  RGB(1,1,1),  RGB(0.9,0.1,0.1)] )
 print_fig(params, "bairstow_map", compute_bairstow; xlab = L"u", ylab = L"v", force = false, cmap)
 # att = get_att(params, "bairstow_map", compute_bairstow)
diff --git a/scripts/map_bischi_cournot.jl b/scripts/map_bischi_cournot.jl
@@ -41,6 +41,6 @@ res = 1200
 μ1 = 3.55; μ2 = 3.55
 params = @strdict res μ1 μ2
 # cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(0.34,0.34,1), RGB(1,0.46,0.46), RGB(0.1,0.1,0.1) ] )
-print_fig(params, "cournot2d", compute_cournot_2d; force = false) 
+print_fig(params, "cournot2d", compute_cournot_2d; force = false, cmap) 
 # att = get_att(params, "cournot2d", compute_cplog) 
 
diff --git a/scripts/map_bogdanov.jl b/scripts/map_bogdanov.jl
@@ -41,5 +41,5 @@ k = 1.2
 ε = 0.0125 
 res = 2000
 params = @strdict res μ k ε
-cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(0.7,0.7,0.7), RGB(1,0,0)] )
+# cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(0.7,0.7,0.7), RGB(1,0,0)] )
 print_fig(params, "bogdanov", compute_bogdanov; force = false,  cmap) 
diff --git a/scripts/map_cbic.jl b/scripts/map_cbic.jl
@@ -40,5 +40,5 @@ end
 j = 0.66
 # res = 1000
 params = @strdict res μ j
-cmap = ColorScheme([RGB(1,1,1), RGB(1,0,0), RGB(0.8,0.8,0.8), RGB(0,0,0)] )
+# cmap = ColorScheme([RGB(1,1,1), RGB(1,0,0), RGB(0.8,0.8,0.8), RGB(0,0,0)] )
 print_fig(params, "cbic_map", compute_cbic_map; cmap) 
diff --git a/scripts/map_cpled_logstc.jl b/scripts/map_cpled_logstc.jl
@@ -35,6 +35,6 @@ res = 1200
 # a = 3.57480493875920; ε = -0.2
 a = 3.6; ε = -1.
 params = @strdict res a ε
-cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(0.34,0.34,1), RGB(1,0.46,0.46), RGB(0.1,0.1,0.1) ] )
+# cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(0.34,0.34,1), RGB(1,0.46,0.46), RGB(0.1,0.1,0.1) ] )
 print_fig(params, "cplog", compute_cplog; force = false, cmap) 
 
diff --git a/scripts/map_feudel.jl b/scripts/map_feudel.jl
@@ -41,4 +41,4 @@ a = 0.0015
 # a = 0.0024
 # res = 400
 params = @strdict res a r
-print_fig(params, "feudel", compute_feudel; xlab = L"x", ylab = L"\theta", force = false)
+print_fig(params, "feudel", compute_feudel; xlab = L"x", ylab = L"\theta", force = false, cmap)
diff --git a/scripts/map_grebogi.jl b/scripts/map_grebogi.jl
@@ -32,5 +32,5 @@ end
 
 
 params = @strdict res 
-cmap = ColorScheme([RGB(1,1,1),  RGB(0.9,0.2,0.1)] )
+# cmap = ColorScheme([RGB(1,1,1),  RGB(0.9,0.2,0.1)] )
 print_fig(params, "grebogi", compute_grebogi; ylab = L"x", xlab = L"\theta", cmap)
diff --git a/scripts/map_henon.jl b/scripts/map_henon.jl
@@ -31,5 +31,5 @@ end
 j = 0.9
 # res = 1000
 params = @strdict res μ j
-cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(1,0.46, 0.46), RGB(0.34,0.34,1), RGB(0.1,0.1,0.1) ] )
+# cmap = ColorScheme([RGB(1,1,1), RGB(0,1,0), RGB(1,0.46, 0.46), RGB(0.34,0.34,1), RGB(0.1,0.1,0.1) ] )
 print_fig(params, "henon", compute_henon; cmap)
diff --git a/scripts/map_holt_mcpeek.jl b/scripts/map_holt_mcpeek.jl
@@ -40,5 +40,5 @@ end
 
 res = 1200; 
 params = @strdict res 
-cmap = ColorScheme([RGB(1,1,1),  RGB(0.9,0.2,0.1)] )
+# cmap = ColorScheme([RGB(1,1,1),  RGB(0.9,0.2,0.1)] )
 print_fig(params, "cazelle", compute_cazelle; xlab = L"x", ylab = L"y", force = false, cmap)
Original file line number	Diff line number	Diff line change
Expand Up		@@ -54,4 +54,4 @@ end


		a = b = 1; c = 0.2; ω=1.; # smooth boundary
		print_fig(600, 600, a, b, c, ω, 600)
		print_fig(600, 600, a, b, c, ω, 600; cmap)