|
52 | 52 | # generate Euler ensemble simulations for |
53 | 53 | # section 5c1 with assimilation of all |
54 | 54 | # fields |
55 | | -gen_5c1_euler_full = True |
| 55 | +gen_5c1_euler_full = False |
56 | 56 | # generate rising bubble ensemble |
57 | 57 | # simulations for section 5c2 with |
58 | 58 | # assimilation of the momenta fields |
59 | 59 | # along with the truth / obs |
60 | 60 | gen_5b_5c2_rb = False |
61 | 61 |
|
62 | | -gen_loc_errors_tc = False |
63 | | -gen_loc_errors_rb = False |
| 62 | +gen_loc_errors_tv = False |
| 63 | +gen_loc_errors_rb = True |
64 | 64 |
|
65 | 65 | # Otherwise, if gen_all = True, generate all |
66 | 66 | # results |
|
623 | 623 |
|
624 | 624 |
|
625 | 625 |
|
626 | | -if gen_loc_errors_tc or gen_all: |
| 626 | +if gen_loc_errors_tv or gen_all: |
627 | 627 | # Set ensemble with 10 members |
628 | 628 | rp.N = 10 |
629 | 629 | # Euler travelling vortex experiment |
|
770 | 770 | # Euler travelling vortex experiment |
771 | 771 | rp.tc = 'rb' |
772 | 772 | ud = { |
773 | | - 'aux' : 'wda_CFLfixed_63', |
| 773 | + 'aux' : 'wda_CFLfixed_179_79', |
774 | 774 | # Do blending for initial time-step |
775 | 775 | 'initial_blending' : True |
776 | 776 | } |
|
781 | 781 | 'obs_attrs' : ['rhou', 'rhov'], |
782 | 782 | # Path to the generated observation |
783 | 783 | 'obs_path' : path_to_obs + 'output_rising_bubble/output_rising_bubble_ensemble=1_160_80_1.0_truth_CFLfixed_ib-0.h5', |
784 | | - 'loc_setter' : (63,63) |
| 784 | + 'loc_setter' : (179,79) |
785 | 785 | } |
786 | 786 |
|
787 | 787 | # run simulation |
788 | 788 | rp.ud = json.dumps(ud) |
789 | 789 | rp.dap = json.dumps(dap) |
790 | | - # rp.queue_run() |
| 790 | + rp.queue_run() |
791 | 791 |
|
792 | 792 | ud = { |
793 | | - 'aux' : 'wda_CFLfixed_63', |
| 793 | + 'aux' : 'wda_CFLfixed_179_79', |
794 | 794 | # Do blending for initial time-step |
795 | 795 | 'initial_blending' : True, |
796 | 796 | 'continuous_blending' : True |
|
799 | 799 | # run simulation |
800 | 800 | rp.ud = json.dumps(ud) |
801 | 801 | rp.dap = json.dumps(dap) |
802 | | - # rp.queue_run() |
| 802 | + rp.queue_run() |
803 | 803 |
|
804 | 804 |
|
805 | 805 | ud = { |
806 | | - 'aux' : 'wda_CFLfixed_31', |
| 806 | + 'aux' : 'wda_CFLfixed_79', |
807 | 807 | # Do blending for initial time-step |
808 | 808 | 'initial_blending' : True |
809 | 809 | } |
|
816 | 816 | 'obs_attrs' : ['rhou', 'rhov'], |
817 | 817 | # Path to the generated observation |
818 | 818 | 'obs_path' : path_to_obs + 'output_rising_bubble/output_rising_bubble_ensemble=1_160_80_1.0_truth_CFLfixed_ib-0.h5', |
819 | | - 'loc_setter' : (31,31) |
| 819 | + 'loc_setter' : (79,79) |
820 | 820 | } |
821 | 821 |
|
822 | 822 | # run simulation |
|
825 | 825 | rp.queue_run() |
826 | 826 |
|
827 | 827 | ud = { |
828 | | - 'aux' : 'wda_CFLfixed_31', |
| 828 | + 'aux' : 'wda_CFLfixed_79', |
829 | 829 | # Do blending for initial time-step |
830 | 830 | 'initial_blending' : True, |
831 | 831 | 'continuous_blending' : True |
|
838 | 838 |
|
839 | 839 |
|
840 | 840 | ud = { |
841 | | - 'aux' : 'wda_CFLfixed_21', |
| 841 | + 'aux' : 'wda_CFLfixed_41', |
842 | 842 | # Do blending for initial time-step |
843 | 843 | 'initial_blending' : True |
844 | 844 | } |
|
850 | 850 | 'obs_attrs' : ['rhou', 'rhov'], |
851 | 851 | # Path to the generated observation |
852 | 852 | 'obs_path' : path_to_obs + 'output_rising_bubble/output_rising_bubble_ensemble=1_160_80_1.0_truth_CFLfixed_ib-0.h5', |
853 | | - 'loc_setter' : (21,21) |
| 853 | + 'loc_setter' : (41,41) |
854 | 854 | } |
855 | 855 |
|
856 | 856 | # run simulation |
|
859 | 859 | rp.queue_run() |
860 | 860 |
|
861 | 861 | ud = { |
862 | | - 'aux' : 'wda_CFLfixed_21', |
| 862 | + 'aux' : 'wda_CFLfixed_41', |
863 | 863 | # Do blending for initial time-step |
864 | 864 | 'initial_blending' : True, |
865 | 865 | 'continuous_blending' : True |
|
919 | 919 | gen_5c1_euler_momenta == False and \ |
920 | 920 | gen_5c1_euler_full == False and \ |
921 | 921 | gen_5b_5c2_rb == False and \ |
922 | | -gen_loc_errors_tc == False and \ |
| 922 | +gen_loc_errors_tv == False and \ |
923 | 923 | gen_loc_errors_rb == False and \ |
924 | 924 | gen_all == False: |
925 | 925 | print("no results generated, all switches were set to False") |
|
0 commit comments