Now for non-DMDA results one can plot a single component using "-dof …

…N -c k" where k=0,1,...,N.
bueler · Oct 13, 2021 · 5f6aca9 · 5f6aca9
1 parent f863ec7
commit 5f6aca9
Showing 1 changed file with 6 additions and 3 deletions.
diff --git a/c/ch5/plotTS.py b/c/ch5/plotTS.py
@@ -103,13 +103,16 @@
             plt.pause(0.1)
     print('.')
 else:
-    for k in range(dims[0]):
-        plt.plot(t,U[k],label='y[%d]' % k)
+    if args.dof > 1:
+        print('plotting only component %d ...' % args.c)
+        plt.plot(t,U[args.c],label='y[%d]' % args.c)
+    else:
+        for k in range(dims[0]):
+            plt.plot(t,U[k],label='y[%d]' % k)
     plt.xlabel('t')
     plt.legend()
     if args.filename:
         print('writing file %s' % args.filename)
         plt.savefig(args.filename)
     else:
         plt.show()
-