Fix example ChenShen98

LLNL · May 22, 2018 · a19f3de · a19f3de
1 parent 61b241b
commit a19f3de
Show file tree

Hide file tree

Showing 3 changed files with 16 additions and 4 deletions.
diff --git a/examples/ChenShen98/README b/examples/ChenShen98/README
@@ -37,16 +37,17 @@
 #generate initial conditions: grain of radius 100 cells in 256 x 256 square
 python ../../utils/make_sphere.py -x 256 -y 256 -z 1 -r 100 --qlen 1 initial256.nc
 
-set exe = /usr/gapps/phasefield/bin/ampe2d
+set exe = ampe2d
 
-srun -ppdebug -N1 -n16 $exe gg2d_256.input
+srun -ppdebug -n16 $exe gg2d_256.input
 
-python dtvst.py gg2d_256.log > 256.dat
+python dtvst.py gg2d_256.log > dtvst.dat
 
+gnuplot dtvst.plt
 
 #increase resolution: 512 x 512 mesh
 python ../../utils/make_sphere.py -x 512 -y 512 -z 1 -r 200 --qlen 1 initial512.nc
 
 srun -ppdebug -N1 -n16 $exe gg2d_512.input
 
-python dtvst.py gg2d_512.log > 512.dat
+python dtvst.py gg2d_512.log > dtvst.dat
diff --git a/examples/ChenShen98/dtvst.plt b/examples/ChenShen98/dtvst.plt
@@ -0,0 +1,8 @@
+set terminal png
+set pointsize 1
+set output "dtvst.png"
+
+set logscale y
+
+plot "dtvst.dat" lt 1 pt 5 with linespoints
+
diff --git a/examples/ChenShen98/gg2d_256.input b/examples/ChenShen98/gg2d_256.input
@@ -95,6 +95,8 @@ ModelParameters {
 
    free_energy_liquid = 0.
    free_energy_solid  = 0.
+
+   molar_volume = 1. // arbitrary value > 0.
 }
 
 ScalarDiagnostics {
@@ -104,6 +106,7 @@ ScalarDiagnostics {
 
 Integrator {
    atol = 3.e-5
+   Preconditioner{}
 }
 
 InitialConditions {