<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 @@
/
&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("-f", "--file", dest="filename", help="file to visualize", metavar="FILE")
parser.add_option("-t", "--time", dest="time", help="time step to visualize (0 based)", metavar="TIME")
+parser.add_option("-s", "--save", action="store_true", dest="saveimages", help="include this flag to save plots as files")
+parser.add_option("-n", "--nodisp", action="store_true", dest="hidefigs", help="include this flag to not display plots (usually used with -s)")
# parser.add_option("-v", "--var", dest="variable", help="variable to visualize", metavar="VAR")
# parser.add_option("--max", dest="maximum", help="maximum for color bar", metavar="MAX")
# parser.add_option("--min", dest="minimum", help="minimum for color bar", metavar="MIN")
@@ -92,6 +92,9 @@
plt.colorbar()
plt.title('thickness at time ' + str(time_slice) )
plt.draw()
+if options.saveimages:
+ print "Saving figures to files."
+ 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 "Plot display disabled with -n argument."
+else:
+ plt.show()
-plt.show()
f.close()
</font>
</pre>