Small modifications to 1D input files to make them start from a small

x/r position to avoid huge fields in Radial case when r -> 0 due to
charge accumulation.

Added a Gnuplot script to plot both potential profiles extracted from
Gmsh.

runs/1D_Cathode/curve_potential.gp

@@ -0,0 +1,22 @@
+#Export the potential curves from Gmsh into a txt file.
+reset
+set terminal qt persist enhanced 1 size 600, 400
+
+set style line 1 pt 1 lc rgb "red"
+set style line 2 pt 2 lc rgb "blue"
+set style line 3 pt 4 lc rgb "#006400"
+set style line 4 pt 6 lc rgb "orange"
+set style line 5 pt 8 lc rgb "black"
+set style line 6 pt 10 lc rgb "#ADD8E6"
+
+set xlabel "x/r (mm)"
+
+set autoscale y
+set ylabel "Potential (V)"
+
+set key box at 90,-10
+
+set title "Effect of geometry"
+
+plot "Cartesian.dat" u ($5*1e3):($8) every 100 ls 1 t "Cartesian", \
+     "Radial.dat"    u ($5*1e3):($8) every 100 ls 2 t "Radial"

runs/1D_Cathode/inputCart.json

@@ -34,12 +34,12 @@
     {"name": "Infinite", "type": "dirichlet", "potential":    0.0, "physicalSurface": 2}
   ],
   "inject": [
-    {"name": "Cathode Electron", "species": "Electron", "flow":  9.0e-5, "units": "A", "v": 27500.0, "T": [2500.0, 2500.0, 2500.0],
+    {"name": "Cathode Electron", "species": "Electron", "flow":  1.0e-3, "units": "A", "v": 27500.0, "T": [2500.0, 2500.0, 2500.0],
                                  "velDist": ["Maxwellian", "Maxwellian", "Maxwellian"],     "n": [ 1, 0, 0], "physicalSurface": 1}
   ],
   "case": {
     "tau":  [1.0e-11],
-    "time": 1.0e-6,
+    "time": 1.0e-7,
     "pusher": ["1DCartCharged"],
     "EMSolver": "Electrostatic"
   },

runs/1D_Cathode/inputRad.json

@@ -34,12 +34,12 @@
     {"name": "Infinite", "type": "dirichlet", "potential":    0.0, "physicalSurface": 2}
   ],
   "inject": [
-    {"name": "Cathode Electron", "species": "Electron", "flow":  9.0e-5, "units": "A", "v": 27500.0, "T": [2500.0, 2500.0, 2500.0],
+    {"name": "Cathode Electron", "species": "Electron", "flow":  1.0e-3, "units": "A", "v": 27500.0, "T": [2500.0, 2500.0, 2500.0],
                                  "velDist": ["Maxwellian", "Maxwellian", "Maxwellian"],     "n": [ 1, 0, 0], "physicalSurface": 1}
   ],
   "case": {
     "tau":  [1.0e-11],
-    "time": 1.0e-6,
+    "time": 1.0e-7,
     "pusher": ["1DRadCharged"],
     "EMSolver": "Electrostatic"
   },

src/modules/moduleInput.f90

@@ -366,12 +366,12 @@ MODULE moduleInput
       !Gets the basename of the folder
       CALL config%get(object // '.folder', baseName, found)
       IF (found) THEN
-        folder = baseName
+        folder = baseName // '_'
 
       END IF
 
       !Compose the folder name
-      folder = folder // '_' // date_now(1:4) // '-' // date_now(5:6) // '-' // date_now(7:8) // '_' &
+      folder = folder // date_now(1:4) // '-' // date_now(5:6) // '-' // date_now(7:8) // '_' &
                       // time_now(1:2) // '.' // time_now(3:4) // '.' // time_now(5:6)
       
       !Creates the folder