fixed script to generate boussinesq results and figure

Author: Daniel Ruprecht

examples/boussinesq_2d_imex/plotgmrescounter.py

@@ -25,17 +25,17 @@
   rcParams['figure.figsize'] = 5.0, 2.5
   fig = plt.figure()
 
-  plt.plot(xx[:,5], udirk[2,:,5], '--', color='r', markersize=fs-2, label='DIRK', dashes=(3,3))
-  plt.plot(xx[:,5], uend[2,:,5], '-', color='b', label='SDC')
-  plt.plot(xx[:,5], uimex[2,:,5], '--', color='g', markersize=fs-2, label='RK-IMEX', dashes=(3,3))
+  plt.plot(xx[:,5], udirk[2,:,5], '--', color='g', markersize=fs-2, label='DIRK(4)', dashes=(3,3))
+  plt.plot(xx[:,5], uend[2,:,5], '-', color='b', label='SDC(4)')
+  plt.plot(xx[:,5], uimex[2,:,5], '--', color='r', markersize=fs-2, label='RK-IMEX(4)', dashes=(3,3))
   #plt.plot(xx[:,5], utrap[2,:,5], '--', color='k', markersize=fs-2, label='Trap', dashes=(3,3))
   plt.legend(loc='lower left', fontsize=fs, prop={'size':fs})
   plt.yticks(fontsize=fs)
   plt.xticks(fontsize=fs)
   plt.xlabel('x [km]', fontsize=fs, labelpad=0)
   plt.ylabel('Bouyancy', fontsize=fs, labelpad=1)
   #plt.show()
-  filename = 'sdc-fwsw-boussinesq.pdf'
+  filename = 'boussinesq.pdf'
   plt.savefig(filename, bbox_inches='tight')
   call(["pdfcrop", filename, filename])
 

examples/boussinesq_2d_imex/rungmrescounter.py

@@ -44,7 +44,7 @@
     sparams = {}
     sparams['maxiter'] = 4
 
-    dirk_order = 4
+    dirk_order = sparams['maxiter']
 
     # setup parameters "in time"
     t0     = 0