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/')
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')
# 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')
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')
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')
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) )