now printing cfl numbers prior to simulation

Author: Daniel Ruprecht

examples/boussinesq_2d_imex/HookClass.py

@@ -43,7 +43,7 @@ def dump_step(self,status):
 
         fs = 18
 
-        if True:
+        if False:
           yplot = self.level.uend.values
           xx    = self.level.prob.xx
           zz    = self.level.prob.zz

examples/boussinesq_2d_imex/playground.py

@@ -28,7 +28,7 @@
 
     # This comes as read-in for the level class
     lparams = {}
-    lparams['restol'] = 1E-12
+    lparams['restol'] = 1E-6
 
     swparams = {}
     swparams['collocation_class'] = collclass.CollGaussLobatto
@@ -41,7 +41,7 @@
     # setup parameters "in time"
     t0     = 0
     Tend   = 3000
-    Nsteps = 500
+    Nsteps =  250
     dt = Tend/float(Nsteps)
 
     # This comes as read-in for the problem class
@@ -52,10 +52,10 @@
     pparams['Nfreq']    = 0.01
     pparams['x_bounds'] = [(-150.0, 150.0)]
     pparams['z_bounds'] = [(   0.0,  10.0)]
-    pparams['order']    = [0, 0] # [fine_level, coarse_level]
-    pparams['gmres_maxiter'] = [50, 50]
+    pparams['order']    = [0] # [fine_level, coarse_level]
+    pparams['gmres_maxiter'] = [50]
     pparams['gmres_restart'] = 20
-    pparams['gmres_tol']     = 1e-14
+    pparams['gmres_tol']     = 1e-6
 
     # This comes as read-in for the transfer operations
     tparams = {}
@@ -80,6 +80,13 @@
     # get initial values on finest level
     P = MS[0].levels[0].prob
     uinit = P.u_exact(t0)
+    
+    cfl_advection = pparams['u_adv']*dt/P.h[0]
+    cfl_acoustic_hor = pparams['c_s']*dt/P.h[0]
+    cfl_acoustic_ver = pparams['c_s']*dt/P.h[1]
+    print ("CFL number of advection: %4.2f" % cfl_advection)
+    print ("CFL number of acoustics (horizontal): %4.2f" % cfl_acoustic_hor)
+    print ("CFL number of acoustics (vertical):   %4.2f" % cfl_acoustic_ver)
 
     # call main function to get things done...
     uend,stats = mp.run_pfasst(MS,u0=uinit,t0=t0,dt=dt,Tend=Tend)