From 3faddc48c762e55159025747ab72a5438d8e9f1a Mon Sep 17 00:00:00 2001 From: Chi Yan Date: Wed, 14 Jun 2023 12:06:11 -0400 Subject: [PATCH 1/6] Breuer_Case0_Benchmark --- B_Bench/B_Bench_init | Bin 0 -> 44 bytes B_Bench/main_input_0T | 90 ++++++++++++++++++++++++++++++++++++++++ B_Bench/main_input_100T | 90 ++++++++++++++++++++++++++++++++++++++++ B_Bench/main_input_80T | 90 ++++++++++++++++++++++++++++++++++++++++ 4 files changed, 270 insertions(+) create mode 100644 B_Bench/B_Bench_init create mode 100644 B_Bench/main_input_0T create mode 100644 B_Bench/main_input_100T create mode 100644 B_Bench/main_input_80T diff --git a/B_Bench/B_Bench_init b/B_Bench/B_Bench_init new file mode 100644 index 0000000000000000000000000000000000000000..893fcff3440c17c8a1822f5ff43c69419fbfa1f6 GIT binary patch literal 44 ecmcCwWMBXyFo4mlKpqPa&zcDWJMAH|3=9B8Nd(pa literal 0 HcmV?d00001 diff --git a/B_Bench/main_input_0T b/B_Bench/main_input_0T new file mode 100644 index 00000000..0aff72d9 --- /dev/null +++ b/B_Bench/main_input_0T @@ -0,0 +1,90 @@ +&problemsize_namelist + n_r = 48 + n_theta = 96 + nprow = 8 + npcol = 12 + aspect_ratio = 0.35d0 + shell_depth = 1.0d0 +/ +&numerical_controls_namelist +/ +&physical_controls_namelist + benchmark_mode = 0 + rotation = .True. + magnetism = .false. + viscous_heating = .false. + ohmic_heating = .false. + n_active_scalars = 1 +/ +&temporal_controls_namelist + max_time_step = 7.5e-5 + max_iterations = 100000 + checkpoint_interval = 10000 + cflmin = 0.4d0 + cflmax = 0.6d0 +/ +&io_controls_namelist +/ +&output_namelist + + +globalavg_values = 401,405,409 +globalavg_frequency = 1000 +globalavg_nrec = 10 + + +shellavg_values = 401,405,409,501,507,10001,10002 +shellavg_frequency = 1000 +shellavg_nrec = 10 + +equatorial_values = 1,501,10001 +equatorial_frequency = 1000 +equatorial_nrec = 10 + +shellslice_levels_nrm = 0.675 +shellslice_values = 1 +shellslice_frequency = 1000 +shellslice_nrec = 10 + + +azavg_values = 1,2,3,201,202,501 +azavg_frequency = 1000 +azavg_nrec = 10 + + +/ +&Boundary_Conditions_Namelist +no_slip_boundaries = .true. +strict_L_Conservation = .false. +T_Top = 0.0d0 +T_Bottom = 1.0d0 +fix_tvar_top = .true. +fix_tvar_bottom = .true. +chi_a_Top = 0.0d0 +chi_a_Bottom = 1.0d0 +fix_chivar_a_top = .true. +fix_chivar_a_bottom = .true. +/ +&Initial_Conditions_Namelist +init_type = 8 +!restart_iter = 0 +t_init_file = 'Bbench_t_init' +chi_a_init_file = 'Bbench_t_init' +/ +&Test_Namelist +/ +&Reference_Namelist + +Ekman_Number = 1.0d-3 +Rayleigh_Number = 0.0d0 +Prandtl_Number = 3.0d-1 +chi_a_Rayleigh_Number = -6.0d5 +chi_a_Prandtl_Number = 3.0d0 + +Magnetic_Prandtl_Number = 5.0d0 +reference_type = 1 +heating_type = 0 +gravity_power = 1.0d0 +/ +&Transport_Namelist +/ diff --git a/B_Bench/main_input_100T b/B_Bench/main_input_100T new file mode 100644 index 00000000..5f7a7de9 --- /dev/null +++ b/B_Bench/main_input_100T @@ -0,0 +1,90 @@ +&problemsize_namelist + n_r = 48 + n_theta = 96 + nprow = 8 + npcol = 12 + aspect_ratio = 0.35d0 + shell_depth = 1.0d0 +/ +&numerical_controls_namelist +/ +&physical_controls_namelist + benchmark_mode = 0 + rotation = .True. + magnetism = .false. + viscous_heating = .false. + ohmic_heating = .false. + n_active_scalars = 0 +/ +&temporal_controls_namelist + max_time_step = 7.5e-5 + max_iterations = 100000 + checkpoint_interval = 10000 + cflmin = 0.4d0 + cflmax = 0.6d0 +/ +&io_controls_namelist +/ +&output_namelist + + +globalavg_values = 401,405,409 +globalavg_frequency = 1000 +globalavg_nrec = 10 + + +shellavg_values = 401,405,409,501,507,10001,10002 +shellavg_frequency = 1000 +shellavg_nrec = 10 + +equatorial_values = 1,501,10001 +equatorial_frequency = 1000 +equatorial_nrec = 10 + +shellslice_levels_nrm = 0.675 +shellslice_values = 1 +shellslice_frequency = 1000 +shellslice_nrec = 10 + + +azavg_values = 1,2,3,201,202,501 +azavg_frequency = 1000 +azavg_nrec = 10 + + +/ +&Boundary_Conditions_Namelist +no_slip_boundaries = .true. +strict_L_Conservation = .false. +T_Top = 0.0d0 +T_Bottom = 1.0d0 +fix_tvar_top = .true. +fix_tvar_bottom = .true. +chi_a_Top = 0.0d0 +chi_a_Bottom = 1.0d0 +fix_chivar_a_top = .true. +fix_chivar_a_bottom = .true. +/ +&Initial_Conditions_Namelist +init_type = 8 +!restart_iter = 0 +t_init_file = 'bench_t_init' +chi_a_init_file = 'bench_t_init' +/ +&Test_Namelist +/ +&Reference_Namelist + +Ekman_Number = 1.0d-3 +Rayleigh_Number = 6.0d4 +Prandtl_Number = 3.0d-1 +chi_a_Rayleigh_Number = -0.0d5 +chi_a_Prandtl_Number = 3.0d0 + +Magnetic_Prandtl_Number = 5.0d0 +reference_type = 1 +heating_type = 0 +gravity_power = 1.0d0 +/ +&Transport_Namelist +/ diff --git a/B_Bench/main_input_80T b/B_Bench/main_input_80T new file mode 100644 index 00000000..2bf76185 --- /dev/null +++ b/B_Bench/main_input_80T @@ -0,0 +1,90 @@ +&problemsize_namelist + n_r = 48 + n_theta = 96 + nprow = 8 + npcol = 12 + aspect_ratio = 0.35d0 + shell_depth = 1.0d0 +/ +&numerical_controls_namelist +/ +&physical_controls_namelist + benchmark_mode = 0 + rotation = .True. + magnetism = .false. + viscous_heating = .false. + ohmic_heating = .false. + n_active_scalars = 1 +/ +&temporal_controls_namelist + max_time_step = 7.5e-5 + max_iterations = 100000 + checkpoint_interval = 10000 + cflmin = 0.4d0 + cflmax = 0.6d0 +/ +&io_controls_namelist +/ +&output_namelist + + +globalavg_values = 401,405,409 +globalavg_frequency = 1000 +globalavg_nrec = 10 + + +shellavg_values = 401,405,409,501,507,10001,10002 +shellavg_frequency = 1000 +shellavg_nrec = 10 + +equatorial_values = 1,501,10001 +equatorial_frequency = 1000 +equatorial_nrec = 10 + +shellslice_levels_nrm = 0.675 +shellslice_values = 1 +shellslice_frequency = 1000 +shellslice_nrec = 10 + + +azavg_values = 1,2,3,201,202,501 +azavg_frequency = 1000 +azavg_nrec = 10 + + +/ +&Boundary_Conditions_Namelist +no_slip_boundaries = .true. +strict_L_Conservation = .false. +T_Top = 0.0d0 +T_Bottom = 1.0d0 +fix_tvar_top = .true. +fix_tvar_bottom = .true. +chi_a_Top = 0.0d0 +chi_a_Bottom = 1.0d0 +fix_chivar_a_top = .true. +fix_chivar_a_bottom = .true. +/ +&Initial_Conditions_Namelist +init_type = 8 +!restart_iter = 0 +t_init_file = 'bench_t_init' +chi_a_init_file = 'bench_t_init' +/ +&Test_Namelist +/ +&Reference_Namelist + +Ekman_Number = 1.0d-3 +Rayleigh_Number = 4.8d4 +Prandtl_Number = 3.0d-1 +chi_a_Rayleigh_Number = -1.2d5 +chi_a_Prandtl_Number = 3.0d0 + +Magnetic_Prandtl_Number = 5.0d0 +reference_type = 1 +heating_type = 0 +gravity_power = 1.0d0 +/ +&Transport_Namelist +/ From 52d140b020374caa14d6a0d7af7c42ab72a4dcfc Mon Sep 17 00:00:00 2001 From: Chi Yan Date: Wed, 14 Jun 2023 18:42:27 -0400 Subject: [PATCH 2/6] Breuer_2010_case0_benchmark --- input_examples/b2010_case0_0T_input | 90 ++++++++++++++++++++++++++ input_examples/b2010_case0_100T_input | 90 ++++++++++++++++++++++++++ input_examples/b2010_case0_80T_input | 90 ++++++++++++++++++++++++++ input_examples/b2010_case0_initcond | Bin 0 -> 44 bytes 4 files changed, 270 insertions(+) create mode 100644 input_examples/b2010_case0_0T_input create mode 100644 input_examples/b2010_case0_100T_input create mode 100644 input_examples/b2010_case0_80T_input create mode 100644 input_examples/b2010_case0_initcond diff --git a/input_examples/b2010_case0_0T_input b/input_examples/b2010_case0_0T_input new file mode 100644 index 00000000..b5b6ba73 --- /dev/null +++ b/input_examples/b2010_case0_0T_input @@ -0,0 +1,90 @@ +&problemsize_namelist + n_r = 48 + n_theta = 96 + nprow = 8 + npcol = 12 + aspect_ratio = 0.35d0 + shell_depth = 1.0d0 +/ +&numerical_controls_namelist +/ +&physical_controls_namelist + benchmark_mode = 0 + rotation = .True. + magnetism = .false. + viscous_heating = .false. + ohmic_heating = .false. + n_active_scalars = 1 +/ +&temporal_controls_namelist + max_time_step = 7.5e-5 + max_iterations = 100000 + checkpoint_interval = 10000 + cflmin = 0.4d0 + cflmax = 0.6d0 +/ +&io_controls_namelist +/ +&output_namelist + + +globalavg_values = 401,405,409 +globalavg_frequency = 1000 +globalavg_nrec = 10 + + +shellavg_values = 401,405,409,501,507,10001,10002 +shellavg_frequency = 1000 +shellavg_nrec = 10 + +equatorial_values = 1,501,10001 +equatorial_frequency = 1000 +equatorial_nrec = 10 + +shellslice_levels_nrm = 0.675 +shellslice_values = 1 +shellslice_frequency = 1000 +shellslice_nrec = 10 + + +azavg_values = 1,2,3,201,202,501 +azavg_frequency = 1000 +azavg_nrec = 10 + + +/ +&Boundary_Conditions_Namelist +no_slip_boundaries = .true. +strict_L_Conservation = .false. +T_Top = 0.0d0 +T_Bottom = 1.0d0 +fix_tvar_top = .true. +fix_tvar_bottom = .true. +chi_a_Top = 0.0d0 +chi_a_Bottom = 1.0d0 +fix_chivar_a_top = .true. +fix_chivar_a_bottom = .true. +/ +&Initial_Conditions_Namelist +init_type = 8 +!restart_iter = 0 +t_init_file = 'b2010_case0_initcond' +chi_a_init_file = 'b2010_case0_initcond' +/ +&Test_Namelist +/ +&Reference_Namelist + +Ekman_Number = 1.0d-3 +Rayleigh_Number = 0.0d0 +Prandtl_Number = 3.0d-1 +chi_a_Rayleigh_Number = -6.0d5 +chi_a_Prandtl_Number = 3.0d0 + +Magnetic_Prandtl_Number = 5.0d0 +reference_type = 1 +heating_type = 0 +gravity_power = 1.0d0 +/ +&Transport_Namelist +/ diff --git a/input_examples/b2010_case0_100T_input b/input_examples/b2010_case0_100T_input new file mode 100644 index 00000000..942a9d98 --- /dev/null +++ b/input_examples/b2010_case0_100T_input @@ -0,0 +1,90 @@ +&problemsize_namelist + n_r = 48 + n_theta = 96 + nprow = 8 + npcol = 12 + aspect_ratio = 0.35d0 + shell_depth = 1.0d0 +/ +&numerical_controls_namelist +/ +&physical_controls_namelist + benchmark_mode = 0 + rotation = .True. + magnetism = .false. + viscous_heating = .false. + ohmic_heating = .false. + n_active_scalars = 0 +/ +&temporal_controls_namelist + max_time_step = 7.5e-5 + max_iterations = 100000 + checkpoint_interval = 10000 + cflmin = 0.4d0 + cflmax = 0.6d0 +/ +&io_controls_namelist +/ +&output_namelist + + +globalavg_values = 401,405,409 +globalavg_frequency = 1000 +globalavg_nrec = 10 + + +shellavg_values = 401,405,409,501,507,10001,10002 +shellavg_frequency = 1000 +shellavg_nrec = 10 + +equatorial_values = 1,501,10001 +equatorial_frequency = 1000 +equatorial_nrec = 10 + +shellslice_levels_nrm = 0.675 +shellslice_values = 1 +shellslice_frequency = 1000 +shellslice_nrec = 10 + + +azavg_values = 1,2,3,201,202,501 +azavg_frequency = 1000 +azavg_nrec = 10 + + +/ +&Boundary_Conditions_Namelist +no_slip_boundaries = .true. +strict_L_Conservation = .false. +T_Top = 0.0d0 +T_Bottom = 1.0d0 +fix_tvar_top = .true. +fix_tvar_bottom = .true. +chi_a_Top = 0.0d0 +chi_a_Bottom = 1.0d0 +fix_chivar_a_top = .true. +fix_chivar_a_bottom = .true. +/ +&Initial_Conditions_Namelist +init_type = 8 +!restart_iter = 0 +t_init_file = 'b2010_case0_initcond' +chi_a_init_file = 'b2010_case0_initcond' +/ +&Test_Namelist +/ +&Reference_Namelist + +Ekman_Number = 1.0d-3 +Rayleigh_Number = 6.0d4 +Prandtl_Number = 3.0d-1 +chi_a_Rayleigh_Number = -0.0d5 +chi_a_Prandtl_Number = 3.0d0 + +Magnetic_Prandtl_Number = 5.0d0 +reference_type = 1 +heating_type = 0 +gravity_power = 1.0d0 +/ +&Transport_Namelist +/ diff --git a/input_examples/b2010_case0_80T_input b/input_examples/b2010_case0_80T_input new file mode 100644 index 00000000..25065179 --- /dev/null +++ b/input_examples/b2010_case0_80T_input @@ -0,0 +1,90 @@ +&problemsize_namelist + n_r = 48 + n_theta = 96 + nprow = 8 + npcol = 12 + aspect_ratio = 0.35d0 + shell_depth = 1.0d0 +/ +&numerical_controls_namelist +/ +&physical_controls_namelist + benchmark_mode = 0 + rotation = .True. + magnetism = .false. + viscous_heating = .false. + ohmic_heating = .false. + n_active_scalars = 1 +/ +&temporal_controls_namelist + max_time_step = 7.5e-5 + max_iterations = 100000 + checkpoint_interval = 10000 + cflmin = 0.4d0 + cflmax = 0.6d0 +/ +&io_controls_namelist +/ +&output_namelist + + +globalavg_values = 401,405,409 +globalavg_frequency = 1000 +globalavg_nrec = 10 + + +shellavg_values = 401,405,409,501,507,10001,10002 +shellavg_frequency = 1000 +shellavg_nrec = 10 + +equatorial_values = 1,501,10001 +equatorial_frequency = 1000 +equatorial_nrec = 10 + +shellslice_levels_nrm = 0.675 +shellslice_values = 1 +shellslice_frequency = 1000 +shellslice_nrec = 10 + + +azavg_values = 1,2,3,201,202,501 +azavg_frequency = 1000 +azavg_nrec = 10 + + +/ +&Boundary_Conditions_Namelist +no_slip_boundaries = .true. +strict_L_Conservation = .false. +T_Top = 0.0d0 +T_Bottom = 1.0d0 +fix_tvar_top = .true. +fix_tvar_bottom = .true. +chi_a_Top = 0.0d0 +chi_a_Bottom = 1.0d0 +fix_chivar_a_top = .true. +fix_chivar_a_bottom = .true. +/ +&Initial_Conditions_Namelist +init_type = 8 +!restart_iter = 0 +t_init_file = 'b2010_case0_initcond' +chi_a_init_file = 'b2010_case0_initcond' +/ +&Test_Namelist +/ +&Reference_Namelist + +Ekman_Number = 1.0d-3 +Rayleigh_Number = 4.8d4 +Prandtl_Number = 3.0d-1 +chi_a_Rayleigh_Number = -1.2d5 +chi_a_Prandtl_Number = 3.0d0 + +Magnetic_Prandtl_Number = 5.0d0 +reference_type = 1 +heating_type = 0 +gravity_power = 1.0d0 +/ +&Transport_Namelist +/ diff --git a/input_examples/b2010_case0_initcond b/input_examples/b2010_case0_initcond new file mode 100644 index 0000000000000000000000000000000000000000..893fcff3440c17c8a1822f5ff43c69419fbfa1f6 GIT binary patch literal 44 ecmcCwWMBXyFo4mlKpqPa&zcDWJMAH|3=9B8Nd(pa literal 0 HcmV?d00001 From 5230b6e17dd01bb86fa357a53331b8fc1b48913c Mon Sep 17 00:00:00 2001 From: Chi Yan Date: Wed, 14 Jun 2023 18:57:05 -0400 Subject: [PATCH 3/6] Breuer_2010_Case0_Benchmark input lists --- input_examples/b2010_case0_0T_input | 1 + input_examples/b2010_case0_100T_input | 1 + input_examples/b2010_case0_80T_input | 1 + input_examples/b2010_case0_InitCond | Bin 0 -> 44 bytes 4 files changed, 3 insertions(+) create mode 100644 input_examples/b2010_case0_InitCond diff --git a/input_examples/b2010_case0_0T_input b/input_examples/b2010_case0_0T_input index b5b6ba73..1bd2cad0 100644 --- a/input_examples/b2010_case0_0T_input +++ b/input_examples/b2010_case0_0T_input @@ -1,3 +1,4 @@ +! The Breuer_2010 Benchmark for thermal-chemical buoyancy systems. &problemsize_namelist n_r = 48 n_theta = 96 diff --git a/input_examples/b2010_case0_100T_input b/input_examples/b2010_case0_100T_input index 942a9d98..e1e9079a 100644 --- a/input_examples/b2010_case0_100T_input +++ b/input_examples/b2010_case0_100T_input @@ -1,3 +1,4 @@ +! The Breuer_2010 Benchmark for thermal-chemical buoyancy systems. &problemsize_namelist n_r = 48 n_theta = 96 diff --git a/input_examples/b2010_case0_80T_input b/input_examples/b2010_case0_80T_input index 25065179..b53e4d23 100644 --- a/input_examples/b2010_case0_80T_input +++ b/input_examples/b2010_case0_80T_input @@ -1,3 +1,4 @@ +! The Breuer_2010 Benchmark for thermal-chemical buoyancy systems. &problemsize_namelist n_r = 48 n_theta = 96 diff --git a/input_examples/b2010_case0_InitCond b/input_examples/b2010_case0_InitCond new file mode 100644 index 0000000000000000000000000000000000000000..893fcff3440c17c8a1822f5ff43c69419fbfa1f6 GIT binary patch literal 44 ecmcCwWMBXyFo4mlKpqPa&zcDWJMAH|3=9B8Nd(pa literal 0 HcmV?d00001 From efeec1596ece5a2ad1eb96d6b39d2fd74e8ef817 Mon Sep 17 00:00:00 2001 From: Chi Yan Date: Thu, 15 Jun 2023 12:13:22 -0400 Subject: [PATCH 4/6] Breuer Bench clean unuseful files --- B_Bench/B_Bench_init | Bin 44 -> 0 bytes B_Bench/main_input_0T | 90 ---------------------------- B_Bench/main_input_100T | 90 ---------------------------- B_Bench/main_input_80T | 90 ---------------------------- input_examples/b2010_case0_initcond | Bin 44 -> 0 bytes 5 files changed, 270 deletions(-) delete mode 100644 B_Bench/B_Bench_init delete mode 100644 B_Bench/main_input_0T delete mode 100644 B_Bench/main_input_100T delete mode 100644 B_Bench/main_input_80T delete mode 100644 input_examples/b2010_case0_initcond diff --git a/B_Bench/B_Bench_init b/B_Bench/B_Bench_init deleted file mode 100644 index 893fcff3440c17c8a1822f5ff43c69419fbfa1f6..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 44 ecmcCwWMBXyFo4mlKpqPa&zcDWJMAH|3=9B8Nd(pa diff --git a/B_Bench/main_input_0T b/B_Bench/main_input_0T deleted file mode 100644 index 0aff72d9..00000000 --- a/B_Bench/main_input_0T +++ /dev/null @@ -1,90 +0,0 @@ -&problemsize_namelist - n_r = 48 - n_theta = 96 - nprow = 8 - npcol = 12 - aspect_ratio = 0.35d0 - shell_depth = 1.0d0 -/ -&numerical_controls_namelist -/ -&physical_controls_namelist - benchmark_mode = 0 - rotation = .True. - magnetism = .false. - viscous_heating = .false. - ohmic_heating = .false. - n_active_scalars = 1 -/ -&temporal_controls_namelist - max_time_step = 7.5e-5 - max_iterations = 100000 - checkpoint_interval = 10000 - cflmin = 0.4d0 - cflmax = 0.6d0 -/ -&io_controls_namelist -/ -&output_namelist - - -globalavg_values = 401,405,409 -globalavg_frequency = 1000 -globalavg_nrec = 10 - - -shellavg_values = 401,405,409,501,507,10001,10002 -shellavg_frequency = 1000 -shellavg_nrec = 10 - -equatorial_values = 1,501,10001 -equatorial_frequency = 1000 -equatorial_nrec = 10 - -shellslice_levels_nrm = 0.675 -shellslice_values = 1 -shellslice_frequency = 1000 -shellslice_nrec = 10 - - -azavg_values = 1,2,3,201,202,501 -azavg_frequency = 1000 -azavg_nrec = 10 - - -/ -&Boundary_Conditions_Namelist -no_slip_boundaries = .true. -strict_L_Conservation = .false. -T_Top = 0.0d0 -T_Bottom = 1.0d0 -fix_tvar_top = .true. -fix_tvar_bottom = .true. -chi_a_Top = 0.0d0 -chi_a_Bottom = 1.0d0 -fix_chivar_a_top = .true. -fix_chivar_a_bottom = .true. -/ -&Initial_Conditions_Namelist -init_type = 8 -!restart_iter = 0 -t_init_file = 'Bbench_t_init' -chi_a_init_file = 'Bbench_t_init' -/ -&Test_Namelist -/ -&Reference_Namelist - -Ekman_Number = 1.0d-3 -Rayleigh_Number = 0.0d0 -Prandtl_Number = 3.0d-1 -chi_a_Rayleigh_Number = -6.0d5 -chi_a_Prandtl_Number = 3.0d0 - -Magnetic_Prandtl_Number = 5.0d0 -reference_type = 1 -heating_type = 0 -gravity_power = 1.0d0 -/ -&Transport_Namelist -/ diff --git a/B_Bench/main_input_100T b/B_Bench/main_input_100T deleted file mode 100644 index 5f7a7de9..00000000 --- a/B_Bench/main_input_100T +++ /dev/null @@ -1,90 +0,0 @@ -&problemsize_namelist - n_r = 48 - n_theta = 96 - nprow = 8 - npcol = 12 - aspect_ratio = 0.35d0 - shell_depth = 1.0d0 -/ -&numerical_controls_namelist -/ -&physical_controls_namelist - benchmark_mode = 0 - rotation = .True. - magnetism = .false. - viscous_heating = .false. - ohmic_heating = .false. - n_active_scalars = 0 -/ -&temporal_controls_namelist - max_time_step = 7.5e-5 - max_iterations = 100000 - checkpoint_interval = 10000 - cflmin = 0.4d0 - cflmax = 0.6d0 -/ -&io_controls_namelist -/ -&output_namelist - - -globalavg_values = 401,405,409 -globalavg_frequency = 1000 -globalavg_nrec = 10 - - -shellavg_values = 401,405,409,501,507,10001,10002 -shellavg_frequency = 1000 -shellavg_nrec = 10 - -equatorial_values = 1,501,10001 -equatorial_frequency = 1000 -equatorial_nrec = 10 - -shellslice_levels_nrm = 0.675 -shellslice_values = 1 -shellslice_frequency = 1000 -shellslice_nrec = 10 - - -azavg_values = 1,2,3,201,202,501 -azavg_frequency = 1000 -azavg_nrec = 10 - - -/ -&Boundary_Conditions_Namelist -no_slip_boundaries = .true. -strict_L_Conservation = .false. -T_Top = 0.0d0 -T_Bottom = 1.0d0 -fix_tvar_top = .true. -fix_tvar_bottom = .true. -chi_a_Top = 0.0d0 -chi_a_Bottom = 1.0d0 -fix_chivar_a_top = .true. -fix_chivar_a_bottom = .true. -/ -&Initial_Conditions_Namelist -init_type = 8 -!restart_iter = 0 -t_init_file = 'bench_t_init' -chi_a_init_file = 'bench_t_init' -/ -&Test_Namelist -/ -&Reference_Namelist - -Ekman_Number = 1.0d-3 -Rayleigh_Number = 6.0d4 -Prandtl_Number = 3.0d-1 -chi_a_Rayleigh_Number = -0.0d5 -chi_a_Prandtl_Number = 3.0d0 - -Magnetic_Prandtl_Number = 5.0d0 -reference_type = 1 -heating_type = 0 -gravity_power = 1.0d0 -/ -&Transport_Namelist -/ diff --git a/B_Bench/main_input_80T b/B_Bench/main_input_80T deleted file mode 100644 index 2bf76185..00000000 --- a/B_Bench/main_input_80T +++ /dev/null @@ -1,90 +0,0 @@ -&problemsize_namelist - n_r = 48 - n_theta = 96 - nprow = 8 - npcol = 12 - aspect_ratio = 0.35d0 - shell_depth = 1.0d0 -/ -&numerical_controls_namelist -/ -&physical_controls_namelist - benchmark_mode = 0 - rotation = .True. - magnetism = .false. - viscous_heating = .false. - ohmic_heating = .false. - n_active_scalars = 1 -/ -&temporal_controls_namelist - max_time_step = 7.5e-5 - max_iterations = 100000 - checkpoint_interval = 10000 - cflmin = 0.4d0 - cflmax = 0.6d0 -/ -&io_controls_namelist -/ -&output_namelist - - -globalavg_values = 401,405,409 -globalavg_frequency = 1000 -globalavg_nrec = 10 - - -shellavg_values = 401,405,409,501,507,10001,10002 -shellavg_frequency = 1000 -shellavg_nrec = 10 - -equatorial_values = 1,501,10001 -equatorial_frequency = 1000 -equatorial_nrec = 10 - -shellslice_levels_nrm = 0.675 -shellslice_values = 1 -shellslice_frequency = 1000 -shellslice_nrec = 10 - - -azavg_values = 1,2,3,201,202,501 -azavg_frequency = 1000 -azavg_nrec = 10 - - -/ -&Boundary_Conditions_Namelist -no_slip_boundaries = .true. -strict_L_Conservation = .false. -T_Top = 0.0d0 -T_Bottom = 1.0d0 -fix_tvar_top = .true. -fix_tvar_bottom = .true. -chi_a_Top = 0.0d0 -chi_a_Bottom = 1.0d0 -fix_chivar_a_top = .true. -fix_chivar_a_bottom = .true. -/ -&Initial_Conditions_Namelist -init_type = 8 -!restart_iter = 0 -t_init_file = 'bench_t_init' -chi_a_init_file = 'bench_t_init' -/ -&Test_Namelist -/ -&Reference_Namelist - -Ekman_Number = 1.0d-3 -Rayleigh_Number = 4.8d4 -Prandtl_Number = 3.0d-1 -chi_a_Rayleigh_Number = -1.2d5 -chi_a_Prandtl_Number = 3.0d0 - -Magnetic_Prandtl_Number = 5.0d0 -reference_type = 1 -heating_type = 0 -gravity_power = 1.0d0 -/ -&Transport_Namelist -/ diff --git a/input_examples/b2010_case0_initcond b/input_examples/b2010_case0_initcond deleted file mode 100644 index 893fcff3440c17c8a1822f5ff43c69419fbfa1f6..0000000000000000000000000000000000000000 GIT binary patch literal 0 HcmV?d00001 literal 44 ecmcCwWMBXyFo4mlKpqPa&zcDWJMAH|3=9B8Nd(pa From cf76b0adeb6846533473eb6d9a198c9944bc025d Mon Sep 17 00:00:00 2001 From: Chi Yan Date: Thu, 15 Jun 2023 13:35:03 -0400 Subject: [PATCH 5/6] Description of B2010 benchmark in documentation --- doc/source/User_Guide/getting_started.rst | 31 ++++++++++++++++++++++- doc/source/User_Guide/references.bib | 12 +++++++++ 2 files changed, 42 insertions(+), 1 deletion(-) diff --git a/doc/source/User_Guide/getting_started.rst b/doc/source/User_Guide/getting_started.rst index f1231d8e..ad096174 100644 --- a/doc/source/User_Guide/getting_started.rst +++ b/doc/source/User_Guide/getting_started.rst @@ -565,6 +565,8 @@ in the directory: Rayleigh/input_examples/ +-----------------------+-----------------+--------------------------------+--------------------------------+ | Jones et al. 2011 | Steady MHD | j2011_steady_MHD_minimal | benchmark_mode = 4 | +-----------------------+-----------------+--------------------------------+--------------------------------+ +| Breuer et al. 2010 | Case 0 | b2010_case0_*T_input | | ++-----------------------+-----------------+--------------------------------+--------------------------------+ Standard benchmarks that generate minimal output files are discussed in the next four benchmarks: @@ -573,7 +575,7 @@ benchmarks: * :ref:`cookbookCase1Minimal` * :ref:`cookbookHydroAnelastic` * :ref:`cookbookMhdAnelastic` - +* :ref:`cookbookDDCBreuer` .. _cookbookCase0Minimal: Boussinesq non-MHD Benchmark: c2001_case0_minimal @@ -698,3 +700,30 @@ differing from the Jones et al. (2011) anelastic hydro benchmark quicksave_interval = 25000 num_quicksaves = 2 / + +.. _cookbookDDCBreuer: + +Steady Thermal-Chemical Boussinesq Convection Benchmark: b2010_case0_*T_input +~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ + +This is a Boussinesq convection benchmark described in Breuer et al. (2010) :cite:`Breuer2010` +in a dual buoyancy system that allows both thermal and chemical buoyancy sources. +The case 0 contains three input lists that describes varying contributions of +thermal vs chemical Rayleigh numbers whereas the total Rayleigh number stays the same. +This benchmark is specified for Rayleigh with input file b2010_case0_*T_input. +Below is an example for 80% thermal and 20% chemical convection scene for the +relevant Fortran namelists: + +:: + + &problemsize_namelist + n_r = 128 + n_theta = 192 + nprow = 32 + npcol = 16 + &Reference_Namelist + Ekman_Number = 1.0d-3 + Rayleigh_Number = 4.8d4 + Prandtl_Number = 3.0d-1 + chi_a_Rayleigh_Number = -1.2d5 + chi_a_Prandtl_Number = 3.0d0 diff --git a/doc/source/User_Guide/references.bib b/doc/source/User_Guide/references.bib index c3e4ee47..f2a3844b 100644 --- a/doc/source/User_Guide/references.bib +++ b/doc/source/User_Guide/references.bib @@ -81,3 +81,15 @@ @article{JONES2011120 keywords = "Magnetic fields, Jovian planets, Interiors", abstract = "Benchmark solutions for fully nonlinear anelastic compressible convection and dynamo action in a rotating spherical shell are proposed. Three benchmarks are specified. The first is a purely hydrodynamic case, which is steady in a uniformly drifting frame. The second is a self-excited saturated dynamo solution, also steady in a drifting frame. The third is again a self-excited dynamo but is unsteady in time, and it has a higher Rayleigh number than the steady dynamo benchmark. Four independent codes have been tested against these benchmarks, and very satisfactory agreement has been found. This provides an accurate reference standard against which new anelastic codes can be tested." } + +@article{Breuer2010, +title = "Thermochemically driven convection in a rotating spherical shell", +journal = "Geophysical Journal International", +volume = "183", +pages = "150-162", +year = "2010", +doi = "http://doi.org/10.1111/j.1365-246X.2010.04722.x", +url = "https://academic.oup.com/gji/article/183/1/150/592285", +author = "M. Breuer and A. Manglik and J. Wicht and T. Trumper and H. Harder and U. Hansen1", +keywords = "Numerical solutions, Dynamo: theories and simulations, Planetary interiors", +} \ No newline at end of file From c34ccfe6160c2ae57350288aa95aee8c617076c0 Mon Sep 17 00:00:00 2001 From: Chi Yan Date: Thu, 15 Jun 2023 13:48:30 -0400 Subject: [PATCH 6/6] minor fix to the reference file in Breuer bench --- doc/source/User_Guide/references.bib | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/doc/source/User_Guide/references.bib b/doc/source/User_Guide/references.bib index f2a3844b..ed9b321d 100644 --- a/doc/source/User_Guide/references.bib +++ b/doc/source/User_Guide/references.bib @@ -92,4 +92,4 @@ @article{Breuer2010 url = "https://academic.oup.com/gji/article/183/1/150/592285", author = "M. Breuer and A. Manglik and J. Wicht and T. Trumper and H. Harder and U. Hansen1", keywords = "Numerical solutions, Dynamo: theories and simulations, Planetary interiors", -} \ No newline at end of file +}