<p><b>mhoffman@lanl.gov</b> 2012-03-12 16:13:16 -0600 (Mon, 12 Mar 2012)</p><p>BRANCH COMMIT -- land_ice<br>
Cleaned up the dome test case: modified the namelist file to use land_ice_grid.nc and cleaned up the output visualization script.<br>
</p><hr noshade><pre><font color="gray">Modified: branches/land_ice_projects/test_cases/dome/namelist.input.land_ice_core
===================================================================
--- branches/land_ice_projects/test_cases/dome/namelist.input.land_ice_core        2012-03-12 21:57:46 UTC (rev 1622)
+++ branches/land_ice_projects/test_cases/dome/namelist.input.land_ice_core        2012-03-12 22:13:16 UTC (rev 1623)
@@ -14,7 +14,7 @@
 /
 
 &amp;io
-   config_input_name = 'grid.nc'
+   config_input_name = 'land_ice_grid.nc'
    config_output_name = 'output.nc'
    config_restart_name = 'restart.nc'
    config_output_interval = '00_00:01:00'

Modified: branches/land_ice_projects/test_cases/dome/visualize_dome.py
===================================================================
--- branches/land_ice_projects/test_cases/dome/visualize_dome.py        2012-03-12 21:57:46 UTC (rev 1622)
+++ branches/land_ice_projects/test_cases/dome/visualize_dome.py        2012-03-12 22:13:16 UTC (rev 1623)
@@ -1,19 +1,19 @@
 #!/usr/bin/python
-import sys, os, glob, shutil, numpy
-from netCDF import *
-from pylab import *
-
+# import numpy
+# from netCDF import *
+from Scientific.IO.NetCDF import *
+# from pylab import *
 from optparse import OptionParser
-
-import time
-import matplotlib
 import matplotlib.pyplot as plt
-from matplotlib.contour import QuadContourSet
+# from matplotlib.contour import QuadContourSet
+# import time
 
 
 parser = OptionParser()
 parser.add_option(&quot;-f&quot;, &quot;--file&quot;, dest=&quot;filename&quot;, help=&quot;file to visualize&quot;, metavar=&quot;FILE&quot;)
 parser.add_option(&quot;-t&quot;, &quot;--time&quot;, dest=&quot;time&quot;, help=&quot;time step to visualize (0 based)&quot;, metavar=&quot;TIME&quot;)
+parser.add_option(&quot;-s&quot;, &quot;--save&quot;, action=&quot;store_true&quot;, dest=&quot;saveimages&quot;, help=&quot;include this flag to save plots as files&quot;)
+parser.add_option(&quot;-n&quot;, &quot;--nodisp&quot;, action=&quot;store_true&quot;, dest=&quot;hidefigs&quot;, help=&quot;include this flag to not display plots (usually used with -s)&quot;)
 # parser.add_option(&quot;-v&quot;, &quot;--var&quot;, dest=&quot;variable&quot;, help=&quot;variable to visualize&quot;, metavar=&quot;VAR&quot;)
 # parser.add_option(&quot;--max&quot;, dest=&quot;maximum&quot;, help=&quot;maximum for color bar&quot;, metavar=&quot;MAX&quot;)
 # parser.add_option(&quot;--min&quot;, dest=&quot;minimum&quot;, help=&quot;minimum for color bar&quot;, metavar=&quot;MIN&quot;)
@@ -92,6 +92,9 @@
 plt.colorbar()
 plt.title('thickness at time ' + str(time_slice) )
 plt.draw()
+if options.saveimages:
+        print &quot;Saving figures to files.&quot;
+        plt.savefig('dome_thickness.png')
 
 fig = plt.figure(2)
 ax = fig.add_subplot(121, aspect='equal')
@@ -104,6 +107,8 @@
 plt.colorbar()
 plt.title('upper surface at time ' + str(time_slice) )
 plt.draw()
+if options.saveimages:
+        plt.savefig('dome_surfaces.png')
 
 
 fig = plt.figure(3)
@@ -115,7 +120,12 @@
     plt.colorbar()
     plt.title('temperature at level '+ str(templevel) + ' at time ' + str(time_slice) )
     plt.draw()
+if options.saveimages:
+        plt.savefig('dome_temperature.png')
 
+if options.hidefigs:
+     print &quot;Plot display disabled with -n argument.&quot;
+else:     
+     plt.show()
 
-plt.show()
 f.close()

</font>
</pre>