Add comments on 2d open channel example

examples/2d_open_channel_through_disorder/README.md

@@ -139,6 +139,15 @@ nz_low  = round(Int,(L_tot-L)/2/dx)
 nz_high = nz_low
 opts.nz_low = nz_low
 opts.nz_high = nz_high
+# opts.use_L0_threads = true would enhances the time performance, 
+# but marginally increases the memoery usage. 
+# the default of opts.use_L0_threads is true for 2D system.
+# here we turn opts.use_L0_threads off, because this L0-threads layer
+# in multithreads makes the numerical results slightly depends on the
+# number of threads with acceptable error. 
+# later we want to show that the result from mesti() and mesti2s()
+# should be exact the same later, so we turn it off to exclude the influence from
+# the this L0-threads layer.
 opts.use_L0_threads = false
 
 # for field-profile computations

examples/2d_open_channel_through_disorder/open_channel_through_disorder.ipynb

@@ -186,6 +186,15 @@
     "nz_high = nz_low\n",
     "opts.nz_low = nz_low\n",
     "opts.nz_high = nz_high\n",
+    "# opts.use_L0_threads = true would enhances the time performance, \n",
+    "# but marginally increases the memoery usage. \n",
+    "# the default of opts.use_L0_threads is true for 2D system.\n",
+    "# here we turn opts.use_L0_threads off, because this L0-threads layer\n",
+    "# in multithreads makes the numerical results slightly depends on the\n",
+    "# number of threads with acceptable error. \n",
+    "# later we want to show that the result from mesti() and mesti2s()\n",
+    "# should be exact the same later, so we turn it off to exclude the influence from\n",
+    "# the this L0-threads layer.\n",
     "opts.use_L0_threads = false\n",
     "\n",
     "# for field-profile computations\n",

examples/2d_open_channel_through_disorder/open_channel_through_disorder.jl

@@ -108,6 +108,15 @@ nz_low  = round(Int,(L_tot-L)/2/dx)
 nz_high = nz_low
 opts.nz_low = nz_low
 opts.nz_high = nz_high
+# opts.use_L0_threads = true would enhances the time performance, 
+# but marginally increases the memoery usage. 
+# the default of opts.use_L0_threads is true for 2D system.
+# here we turn opts.use_L0_threads off, because this L0-threads layer
+# in multithreads makes the numerical results slightly depends on the
+# number of threads with acceptable error. 
+# later we want to show that the result from mesti() and mesti2s()
+# should be exact the same later, so we turn it off to exclude the influence from
+# the this L0-threads layer.
 opts.use_L0_threads = false
 
 # compute the field-profiles through mesti()