minor bug in plotting example in facit.py

fsciortino · Aug 21, 2023 · 2d87f44 · 2d87f44
1 parent 2d4963e
commit 2d87f44
Showing 1 changed file with 6 additions and 4 deletions.
diff --git a/aurora/facit.py b/aurora/facit.py
@@ -12,7 +12,7 @@ class FACIT:
     References:  Maget et al 2020 Plasma Phys. Control. Fusion 62 105001
                  Fajardo et al 2022 Plasma Phys. Control. Fusion 64 055017
                  Maget et al 2022 Plasma Phys. Control. Fusion 64 069501
-                 Fajardo et al 2023 Plasma Phys. Control. Fusion 65 
+                 Fajardo et al 2023 Plasma Phys. Control. Fusion 65 035021
           
     License statement:
     Software name : FACIT
@@ -1448,9 +1448,9 @@ def Jacobian(self, R, Z, r, theta):
     ax1.tick_params(which = 'both', direction = 'in', axis = 'both', top = True, right = True)
 
 
-    ax2.plot(rho, fct.Vz_CL, label = 'CL component')
-    ax2.plot(rho, fct.Vz_BP, label = 'BP component')
-    ax2.plot(rho, fct.Vz_PS, label = 'PS component')
+    ax2.plot(rho, fct.Vconv_CL, label = 'CL component')
+    ax2.plot(rho, fct.Vconv_BP, label = 'BP component')
+    ax2.plot(rho, fct.Vconv_PS, label = 'PS component')
     ax2.plot(rho, fct.Vconv, label = ' Total')
     ax2.tick_params(which = 'both', direction = 'in', axis = 'both', top = True, right = True)
 
@@ -1459,4 +1459,6 @@ def Jacobian(self, R, Z, r, theta):
     ax2.set_ylabel('Convective velocity [m/s]')
 
     fig.tight_layout()
+
+    plt.show()