Example #1
0
import traj_tools as trj
c1 = trj.loadme('c1.trj')

c1.draw_centered_vs_t('P', 'WCB', 'P600', plottype='Smooth', ylim = (100, 1000), sigma = 1, idtext = '061114e', savebase = '/usr/users/stephan.rasp/Dropbox/figures/Case1_20070720/')
c1.draw_centered_vs_t('var4', 'WCB', 'P600', plottype='Smooth', ylim = (-3, 3), sigma = 1, idtext = '061114f', savebase = '/usr/users/stephan.rasp/Dropbox/figures/Case1_20070720/')

Example #2
0
import traj_tools as trj
c0 = trj.loadme('c0.save')
# c0.create_filter('WCB_nonconv1', [('P600', 0, 2880), ('P400', 120, 2880)])
c0.draw_hist('P600', 'WCB_nonconv1', '/usr/users/stephan.rasp/Dropbox/figures/Case0_20121013/', idtext = '290814a')
c0.draw_trj_all([], 'WCB_nonconv1_1080', savebase = '/usr/users/stephan.rasp/Dropbox/figures/Case0_20121013/', idtext = '290814b')
c0.draw_trj_all([], 'WCB_nonconv1_1620', savebase = '/usr/users/stephan.rasp/Dropbox/figures/Case0_20121013/', idtext = '290814c')
c0.draw_trj_all([], 'WCB_nonconv1_1080', savebase = '/usr/users/stephan.rasp/Dropbox/figures/Case0_20121013/onlyasc_', idtext = '290814d', onlyasc='P100', linewidth = 1.5)
c0.draw_trj_all([], 'WCB_nonconv1_1620', savebase = '/usr/users/stephan.rasp/Dropbox/figures/Case0_20121013/onlyasc_', idtext = '290814e', onlyasc='P100', linewidth = 1.5)
c0.draw_scatter('P600', 'P100', factor2 = 6, filtername='WCB_nonconv1', savebase = '/usr/users/stephan.rasp/Dropbox/figures/Case0_20121013/scatter_P600_P100_WCB_nonconv1.png', idtext = '290814f')
c0.draw_scatter('P600', 'P100', factor2 = 6, filtername='WCB_nonconv1_1080', savebase = '/usr/users/stephan.rasp/Dropbox/figures/Case0_20121013/scatter_P600_P100_WCB_nonconv1_1080.png', idtext = '290814g')
c0.draw_scatter('P600', 'P100', factor2 = 6, filtername='WCB_nonconv1_1620', savebase = '/usr/users/stephan.rasp/Dropbox/figures/Case0_20121013/scatter_P600_P100_WCB_nonconv1_1620.png', idtext = '290814h')
Example #3
0
# Script for 030914
import traj_tools as trj


c0 = trj.loadme('c0.trj')
c1 = trj.loadme('c1.trj')

c0.new_asc(100)
c0.new_asc(400)
c0.new_asc(600)

c0.create_filter('WCB', [('P600', 0, 2880)])
c0.create_filter('WCB_NonConv', [('P600', 0, 2880), ('P400', 120, 2880)])
c0.create_filter('WCB_Conv', [('P600', 0, 2880), ('P400', 0, 120)])

# Case 0, P600vP100 Scatter
c0.draw_scatter('P600', 'P100', factor2 = 6, carray = 'P100', 
                filtername = 'WCB_NonConv', 
                savebase = '/usr/users/stephan.rasp/Dropbox/figures/Case0_20121013/', 
                idtext = '030914a')
c0.draw_scatter('P600', 'P100', factor2 = 6, carray = 'P100', 
                filtername = 'WCB_Conv', 
                savebase = '/usr/users/stephan.rasp/Dropbox/figures/Case0_20121013/', 
                idtext = '030914b')
c0.draw_scatter('P600', 'P100', factor2 = 6, carray = 'P100', 
                filtername = 'WCB', 
                savebase = '/usr/users/stephan.rasp/Dropbox/figures/Case0_20121013/', 
                idtext = '030914c')


Example #4
0
import traj_tools as trj
c0 = trj.loadme('c0_rot.trj')
for level in [700, 500, 300]:
    array = c0.draw_intersect_hor(level, idtext='301014a', filtername='WCB', savebase='/usr/users/stephan.rasp/Dropbox/figures/Case0_20121013/intersect/')
c0.create_filter('WCB_Conv', [('P600', 0, 2880), ('P400', 0, 120)])    
c1 = trj.loadme('c1.trj')
for level in [700, 500, 300]:
    array = c1.draw_intersect_hor(level, idtext='301014b', filtername='WCB', savebase='/usr/users/stephan.rasp/Dropbox/figures/Case1_20070720//intersect/')

levellist = [300, 400, 500, 600, 700, 800, 900]
vellist1 = []
vellist2 = []
vellist3 = []
vellist4 = []
for level in levellist:
    vellist1.append((c0.draw_intersect_hor(level, filtername='WCB')).mean())
    vellist2.append((c1.draw_intersect_hor(level, filtername='WCB')).mean())
    vellist3.append((c0.draw_intersect_hor(level, filtername='WCB_NonConv')).mean())
    vellist4.append((c0.draw_intersect_hor(level, filtername='WCB_Conv')).mean())
    plt.close('all')

plt.plot(vellist1, levellist, 'r', label = 'Case0 P600')
plt.plot(vellist2, levellist, 'g', label = 'Case1 P600')
plt.plot(vellist3, levellist, 'b', label = 'Case0 P600 NonConvective')
plt.plot(vellist4, levellist, 'y', label = 'Case0 P600 Convective')
plt.gca().invert_yaxis()
plt.legend(loc=4)

plt.text(0.94, 1.02, '301014c', transform = plt.gca().transAxes, 
             fontsize = 6)
plt.savefig('/usr/users/stephan.rasp/Dropbox/figures/vel_vs_p_301014c.png')
Example #5
0

import traj_tools as trj
c0 = trj.loadme('c0_rot.trj')
c1 = trj.TrjObj('/home/scratch/users/stephan.rasp/Case1_20070720/d4deout/')

c0.new_delta('THETA', mode = 'climb')
c1.new_delta('THETA', mode = 'climb')


c0.draw_hist('deltaTHETAclimb', idtext='151014a', 
             savebase = '/usr/users/stephan.rasp/Dropbox/figures/Case0_20121013/', 
             range = (0, 100))
c1.draw_hist('deltaTHETAclimb', idtext='151014b', 
             savebase = '/usr/users/stephan.rasp/Dropbox/figures/Case1_20070720/',
             range = (0, 100))

c0.new_delta('P', mode = 'climb_r')
c1.new_delta('P', mode = 'climb_r')


c0.draw_hist('deltaPclimb_r', idtext='151014c', 
             savebase = '/usr/users/stephan.rasp/Dropbox/figures/Case0_20121013/', 
             range = (-2500, 0))
c1.draw_hist('deltaPclimb_r', idtext='151014d', 
             savebase = '/usr/users/stephan.rasp/Dropbox/figures/Case1_20070720/', 
             range = (-2500, 0))

c0.new_delta('QV', mode = 'climb_r')
c1.new_delta('QV', mode = 'climb_r')
Example #6
0
import traj_tools as trj

c0 = trj.loadme("c0_rot.trj")
c1 = trj.TrjObj("/home/scratch/users/stephan.rasp/Case1_20070720/d4deout/")

c0.new_delta("QV", mode="mintomax_r")
c1.new_delta("QV", mode="mintomax_r")

c0.draw_hist(
    "deltaQVmintomax_r",
    savebase="/usr/users/stephan.rasp/Dropbox/figures/Case0_20121013/",
    idtext="221014a",
    range=(-0.025, 0),
)
c1.draw_hist(
    "deltaQVmintomax_r",
    savebase="/usr/users/stephan.rasp/Dropbox/figures/Case1_20070720/",
    idtext="221014b",
    range=(-0.025, 0),
)

c0.new_max_diff("z")
c1.new_max_diff("z")

c0.draw_hist(
    "z_max_diff", savebase="/usr/users/stephan.rasp/Dropbox/figures/Case0_20121013/", idtext="221014c", range=(0, 15)
)
c1.draw_hist(
    "z_max_diff", savebase="/usr/users/stephan.rasp/Dropbox/figures/Case1_20070720/", idtext="221014d", range=(0, 15)
)