run scripts used to generate results in the final production files and zenodo datasets

Author: ray-chew

RKLM_Python/queue_run_lamb_forcing.py

@@ -16,16 +16,16 @@
 omega = 7.292 * 1e-5
 # resol_x = [301]
 # resol_y = [120]
-resol_x = [301]
-resol_y = [120]
+resol_x = [151]
+resol_y = [60]
 resol_t = [1,2,4,8,10,12,14,16]
-resol_t = [10.0]
+# resol_t = [10.0]
 # resol_t = [10,12,14,16.0]
 # omegas = [0.0, 2.0 * omega * t_ref]
 # omegas = [2.0 * omega * t_ref]
 omegas = [0.0]
 
-tsteps = [1800, 900, 450, 360, 300, 258, 225]
+tsteps = [3600, 1800, 900, 450, 360, 300, 258, 225]
 tsteps = [360, 300, 258, 225]
 # tsteps = [21]
 
@@ -39,7 +39,7 @@
     ud['iny'] = y+1
 
     ud['ymax'] = 8.0
-    # ud['do_advection'] = False
+    ud['do_advection'] = True
     ud['rayleigh_forcing'] = False
 
     for t_idx, t in enumerate(resol_t):

RKLM_Python/queue_run_lamb_mixed.py

@@ -38,8 +38,9 @@
     ud['ymax'] = 8.0
     ud['rayleigh_forcing'] = False
     ud['mixed_run'] = True
-    # ud['do_advection'] = False
+    ud['do_advection'] = True
     ud['output_timesteps'] = False
+    ud['stepmax'] = 10001
 
 
     for t in resol_t:
@@ -54,7 +55,10 @@
             if om > 0:
                 ud['aux'] = 'test_run_S%i_a05' %t
             else:
-                ud['aux'] = 'test_run_S%i_mix' %t
+                if ud['mixed_run'] == True:
+                    ud['aux'] = 'test_run_S%i_mix' %t
+                else:
+                    ud['aux'] = 'test_run_S%i_adv' %t
 
             print(ud)
             # run simulation

RKLM_Python/queue_run_lamb_unstable.py

@@ -12,24 +12,25 @@
 
 t_ref = 100.0
 omega = 7.292 * 1e-5
-resol_x = [151,301,601]
-resol_y = [15,30,60]
+resol_x = [151,301]
+resol_y = [15,30]
 # resol_x = [601]
 # resol_y = [60]
 resol_x = [301]
 resol_y = [30]
 # resol_t = [10,12,14,16]
 # resol_t = [200,400,600,800,1000,1200,1400,1600]
-# resol_t = [1,2,4,8]
-resol_t = [10.0]
+resol_t = [1,2,4,8]
+# resol_t = [10.0]
 # omegas = [0.0, 2.0 * omega * t_ref]
 omegas = [2.0 * omega * t_ref]
 # omegas = [0.0]
 
 
 # tsteps = [3600, 1800, 900, 450, 360, 300, 258, 225]
-# tsteps = [3600, 1800, 900, 450]
-tsteps = [360]
+tsteps = [3600, 1800, 900, 450]
+# tsteps = [360, 300, 258, 225]
+# tsteps = [360]
 
 ud = {}
 dap = {
@@ -42,6 +43,7 @@
     ud['rayleigh_bdry_switch'] = True
     ud['rayleigh_forcing'] = True
     ud['do_advection'] = True
+    ud['trad_forcing'] = False
 
     for t_idx, t in enumerate(resol_t):
         ud['dtfixed0'] = t / t_ref
@@ -55,9 +57,15 @@
         for om in omegas:
             ud['coriolis_strength'] = [0.0, 0.0, om]
             if om > 0:
-                ud['aux'] = 'bdl_run_S%i_a05' %t
+                if ud['trad_forcing'] == True:
+                    ud['aux'] = 'bdl_run_S%i_a05_trad_forcing_file_test' %t    
+                else:
+                    ud['aux'] = 'bdl_run_S%i_a05' %t
             else:
-                ud['aux'] = 'bdl_run_S%i_noom_a05' %t
+                if ud['trad_forcing'] == True:
+                    ud['aux'] = 'bdl_run_S%i_noom_a05_trad_forcing_func' %t
+                else:
+                    ud['aux'] = 'bdl_run_S%i_noom_a05' %t
 
             print(ud)
             # run simulation