def test_axis_segfault():
from cdms2 import MV2
# Contributed by Lawson Hanson
month = 'January'
year = '2012'
miss_value = -9999.0
data = ['34.006348', '28.314002', '29.269668', '33.698551', '34.177242']
Rad_global_month = MV2.zeros([len(data)], MV2.float)
time_day = MV2.zeros([len(data)], MV2.float)
tim = cdms.createAxis(time_day)
tim.designateTime()
tim.id = "time"
tim.units = 'days since 2012-01-01'
for i in range(len(time_day)):
Rad_global_month[i] = data[i]
tim[i] = float(i)
# Create a temporary file
fd, f = tempfile.mkstemp('.nc')
os.close(fd)
try:
out = cdms.open(f, 'w')
rad_total_month = cdms.createVariable(Rad_global_month,
id='accum_swfcdown_' + month,
axis=[tim],
typecode='f')
rad_total_month.setAxisList([tim])
print rad_total_month, tim
print len(rad_total_month), len(tim)
out.write(rad_total_month)
finally:
out.close()
os.remove(f)
def test_axis_segfault():
from cdms2 import MV2
# Contributed by Lawson Hanson
month = 'January'
year = '2012'
miss_value = -9999.0
data = ['34.006348', '28.314002', '29.269668', '33.698551', '34.177242']
Rad_global_month = MV2.zeros([len(data)],MV2.float)
time_day = MV2.zeros([len(data)],MV2.float)
tim = cdms.createAxis(time_day)
tim.designateTime()
tim.id = "time"
tim.units = 'days since 2012-01-01'
for i in range(len(time_day)):
Rad_global_month[i] = data[i]
tim[i] = float(i)
# Create a temporary file
fd, f = tempfile.mkstemp('.nc')
os.close(fd)
try:
out = cdms.open(f,'w')
rad_total_month = cdms.createVariable(Rad_global_month,id = 'accum_swfcdown_'+month,axis=[tim],typecode='f')
rad_total_month.setAxisList([tim])
print rad_total_month,tim
print len(rad_total_month),len(tim)
out.write(rad_total_month)
finally:
out.close()
os.remove(f)
x.Graph[1].Set[0].line.linewidth = 2
# Let's set the colors
x.Graph[0].Set[0].line.color = 'red'
x.Graph[0].Set[1].line.color = 4 # color 4 is blue
# Let's set the second's set style to dash
x.Graph[1].Set[0].line.linestyle = 'dot'
x.Graph[0].Set[1].line.linestyle = 'dash'
x.Graph[0].Set[0].line.linestyle = 2 # 2 is dash
# Now let's set the name for each set, to be used in the legend
x.Graph[0].Set[0].legend = 'Year 1'
x.Graph[0].Set[1].legend = 'Year 2'
x.Graph[1].Set[0].legend = 'Year 2-1'
# ok one last trick, in order to have the 0 line plotted
# let's create a dummy array with zeros in it:
zero = MV.zeros((2), MV.Float)
zero = MV.asVariable(zero)
zero.setdimattribute(0, 'values', [-1., 1.])
# add a dataset and assign it graph 1 black color
x.add_set(1, 'black')
# Finally just for fun let's place a big red "Sample" accross it
x.add_string(0.5, 0.5, 'Sample', color='red', char_size=9, rot=55, just=14)
# And plot these babies:
x.plot([z1, z2, z3,
zero]) # You MUST pass a list of slab, even if only one slab
# or you can plot then 1 by 1
# first let's clear
x('kill G0.S0')
x('kill G0.S1')
x.Graph[1].Set[0].line.linewidth = 2
# Let's set the colors
x.Graph[0].Set[0].line.color = 'red'
x.Graph[0].Set[1].line.color = 4 # color 4 is blue
# Let's set the second's set style to dash
x.Graph[1].Set[0].line.linestyle = 'dot'
x.Graph[0].Set[1].line.linestyle = 'dash'
x.Graph[0].Set[0].line.linestyle = 2 # 2 is dash
# Now let's set the name for each set, to be used in the legend
x.Graph[0].Set[0].legend = 'Year 1'
x.Graph[0].Set[1].legend = 'Year 2'
x.Graph[1].Set[0].legend = 'Year 2-1'
# ok one last trick, in order to have the 0 line plotted
# let's create a dummy array with zeros in it:
zero = MV.zeros((2), MV.Float)
zero = MV.asVariable(zero)
zero.setdimattribute(0, 'values', [-1., 1.])
# add a dataset and assign it graph 1 black color
x.add_set(1, 'black')
# Finally just for fun let's place a big red "Sample" accross it
x.add_string(0.5, 0.5, 'Sample', color='red', char_size=9, rot=55, just=14)
# And plot these babies:
# You MUST pass a list of slab, even if only one slab
x.plot([z1, z2, z3, zero])
# or you can plot then 1 by 1
# first let's clear
x('kill G0.S0')
x('kill G0.S1')
