source: trunk/UTIL/PYTHON/Intercomparison/Plots_MasterScript/bam.sh @ 812

Last change on this file since 812 was 764, checked in by acolaitis, 13 years ago

###################################################
# PYTHON / PLANETPLOT #
###################################################

# ------------------- XY plots ------------------ #

# Added a new category of plot to unidim, contour, etc... called "xy"
# - xy plots are plots that do not use time,vert,lat or lon as axis
# - variables to be plotted are stored in plot_x and plot_y, which is done in
# select_getfield. (there is no "what_I_plot" var for "wy" plots)
# - plot_x and plot_y are also subject to reduce field. A None value indicates
# the plot is not 'xy'
# - "xy" plots are used for a specific subset of plots : histograms, fourier
# transforms, hodographs
# - added option --analysis to perform certain kinds of analysis on the data
# (corresponding to xy plots).
# - One selects the fields he wants to plot (e.g. -v UV --lon 0 --lat 20) and
# chooses a kind of analysis :

--analysis fft # in the particular case given above, this mode corresponds

# to the mean of the ampitude spectrum of the vertical spatial
# fast fourier transform taken at each time index
# note that for now, fft are always done along the vertical
# axis. This could be made more flexible.
# not that the minimum wavelength you can plot depends on the
# vertical step of your simulation. THIS STEP HAS TO BE
# CONSTANT, hence you MUST use API or ZRECAST with constant
# spacing.

--analysis histo # histogram on the flattened data. If the user asks for --lon 0,20,

# the average is done before computing the histogram (contrary to the fft).
# However, if given a 2D array (with only --lon and --lat on a
# 4D field for example), the data is flattened before computation
# so that the result is still a 1D histogram.

--analysis histodensity # histogram with a kernel density estimate to a

# gaussian distribution, also giving the mean, variance,
# skewness and kurtosis. Other distributions are available in
# the scipy.stats package and could be implemented.

# - added variables "-v hodograph" and "-v hodograph_2". First one is a
# regular hodograph with "u" and "v" as axis, with labels of the local times
# (use --axtime lt). This is a "xy" plot (you must specify a vertical level
# as well, usually --vert 0 with an interpolation at 5m (-i 4 -l 0.005)).
# Second one is the variation with local time (for exemple) of the wind
# rotation (arctan(v/u)). This is not a "xy" plot but "unidim".
# For --ope plots in "xy" cases, only the operation plot is displayed.

# ------------------- Operations ------------------ #

# - For operations --ope -, the histogram (fourth plot) has been removed. To get it
# back, call --ope -_histo.
# - _only has been added to "+","-","-%" operations (eg "-_only")
# - For operations "-","+","-%","-_only","+_only","-%_only","-_histo", it is now
# possible to call multiple files (sill one variable) and compare each of them
# with the unique given reference file. Ex: -f file1,file2 --fref file3 --ope - will give 6 plots:

file1 file3 file1-file3
file2 file3 file2-file3

# In the case of "xy" plots, the multiple operation plots are regrouped on a
# single plot (using multiple lines). the title of this plot is not 'fig(2)-fig(1)'
# (default) but the argument of the --title command. Labels have to be given
# as following : --labels "dummy","dummy","file1-file3 label","dummy","dummy","file2-file3 label" (dummy can be anything. this is to be improved)
# To be able to run on multiple files and easily introduce the correct number and order of plots, these operations have been moved outside of the main loop
# on namefiles and variables. Operation of the type "add_var" or "cat" have
# been left inside the loop and are unchanged.

# ------------------- Localtime ------------------ #

# Changed the way localtime is computed. Reasons:
# - it was assuming one timestep per hour in computations mixing indices and
# actual times
# - to determine the starting date, it was using the name of the file (for Meso), instead of using the netcdf
# attribute (START_DATE)
# - it was using the computed mean longitude of the domain, which is not correct for
# hemispheric domain. => better to use the netcdf attribute CEN_LON
# - it was using this mean longitude even for profile plots at given
# longitude => better to get the local time at this given lon, especially
# important for large domains
# => new localtime() is in myplot (old one is commented). Interv is obsolete (but not removed yet). Case "Geo" has not been looked at.
# new localtime in myplot correctly account for starting date of the file in
# all cases
# accounts for local longitude of the plot
# accounts for files that do not have per hour outputs, but per timestep.
# specific cases can be added in myplot in localtime()

# ------------------ Misc ------------------ #

# - added option --xlog to get x logarithmic axis (--ylog already existed)
# - added the possibility to use 2 files of different gridding (-f
# file1,file2) although of the same type (meso for exemple). For that purpose,
# lon, lat, alt and vert arrays are now indexed with 'index_f', as for
# all_var. => all_lon, all_lat, all_lat, all_vert
# - added function teta_to_tk in myplot so that a call to pp.py on standard
# LMD mesoscale file with "-v tk" can be done without the need to call API.
# The temperature is computed from T and PTOT, knowing P0, T0 and R_CP.
# - bug correction in determineplot() that was causing wrong plot number and
# slices number when not calling averaged lon, lat, vert or times. (which is
# often !)
# - added new options to redope: --redope edge_x1, edge_x2, edge_y1, edge_y2
# which plots the boundaries of the domain. This is different compared to
# asking for a fixed longitude, because the domain boundary might not be at constant
# longitude (hemispheric domain for exemple). x1 is the western boundary, x2
# the eastern, y1 the southern and y2 the northern. x1 and x2 reduce the
# dimension along --lon, y1 and y2 reduce the dimension along --lat.
# - added control in windamplitude() to determine whether winds are staggered or
# unstaggered. This is usefull when dealing with non LMD_MMM files.
# - corrected a bug in reduce_field where the mean was computed on the wrong
# axis !!! (pretty serious bug)

# Exemple of plots you can do with these new options can be found in
$YOUR_SVN/trunk/UTIL/PYTHON/Intercomparison/Plots_MasterScript/bam.sh

# ------------------ API ------------------ #

# - changed maximum of levels from 299 to 1000 in API (interpolation on 1000 levels is
# usefull to get larger bandwidth in fourier transform)
# the following concerns users of MRAMS files.
# - API has not been modified for MRAMS files. Instead, a python script
# (ic.py) is run on MRAMS .ctl and .dat files, which automatically format those files
# to be API and pp.py compatible.
# - ic.py is in $YOUR_SVN/trunk/UTIL/PYTHON/Intercomparison/File_conversion

###################################################
# INTERCOMPARISON TOOLS: #
###################################################

#ic.py in
$YOUR_SVN/trunk/UTIL/PYTHON/Intercomparison/File_conversion
#CDO installer with import_binary in
$YOUR_SVN/trunk/UTIL/PYTHON/Intercomparison/CDO
#Plotting scripts in
$YOUR_SVN/trunk/UTIL/PYTHON/Intercomparison/Plots_MasterScript
#1D sensibility tool in
$YOUR_SVN/trunk/UTIL/PYTHON/Intercomparison/1D_sensibility_tool

# See README in each of these folders for details.
  • Property svn:executable set to *
File size: 32.6 KB
RevLine 
[764]1#!/bin/bash
2
3########## Constants
4mramsfiles=/home/local_home/colaitis/intercomp_test/converted_files
5lmdfiles=/planeto/acolmd/SwRI_files/lmd/HOLDEN_nests_mrams
6lat=-26. # landing site
7lon=-34.   # landing site
8#lon=-7.919998  #central longitude to test synchro
9#lon=0.
10filemrams='_2024-06-17_07:30:00'
11filemramsa='_2024-06-16_07:30:00'
12filemramsb='_2024-06-18_07:30:00'
13filelmd='_2024-06-17_06:00:00'
14filelmda='_2024-06-16_06:00:00'
15filelmdb='_2024-06-18_06:00:00'
16########## ---------
17if [ $1 == boundaries ]
18then
19echo 'Boundaries'
20# ------------------
21if [ $2 == temp ]
22then
23for i in 1 2 3
24do
25/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d0${i}$filemrams --div 20 --axtime lt --time 12 --redope edge_x1 -v tk --title "West-most T boundary MRAMS Nest $i at 12:00" -O mrams_lat1tbound_d0$i -S png -d --res 200 -i 4 -l 0,40,30 -c spectral -m 100 -M 250
26/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d0${i}$filemrams --div 20 --axtime lt --time 12 --redope edge_x2 -v tk --title "East-most T boundary MRAMS Nest $i at 12:00" -O mrams_lat2tbound_d0$i -S png -d --res 200 -i 4 -l 0,40,30 -N -c spectral -m 100 -M 250
27/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d0${i}$filemrams --div 20 --axtime lt --time 12 --redope edge_y1 -v tk --title "South-most T boundary MRAMS Nest $i at 12:00" -O mrams_lon1tbound_d0$i -S png -d --res 200 -i 4 -l 0,40,30 -N -c spectral -m 100 -M 250
28/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d0${i}$filemrams --div 20 --axtime lt --time 12 --redope edge_y2 -v tk --title "North-most T boundary MRAMS Nest $i at 12:00" -O mrams_lon2tbound_d0$i -S png -d --res 200 -i 4 -l 0,40,30 -N -c spectral -m 100 -M 250
29/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d0${i}$filelmd --div 20 --axtime lt --time 12 --redope edge_x1 -v tk --title "West-most T boundary LMD_MMM Nest $i at 12:00" -O lmd_lat1tbound_d0$i -S png -d --res 200 -i 4 -l 0,40,30 -c spectral -m 100 -M 250
30/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d0${i}$filelmd --div 20 --axtime lt --time 12 --redope edge_x2 -v tk --title "East-most T boundary LMD_MMM Nest $i at 12:00" -O lmd_lat2tbound_d0$i -S png -d --res 200 -i 4 -l 0,40,30 -N -c spectral -m 100 -M 250
31/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d0${i}$filelmd --div 20 --axtime lt --time 12 --redope edge_y1 -v tk --title "South-most T boundary LMD_MMM Nest $i at 12:00" -O lmd_lon1tbound_d0$i -S png -d --res 200 -i 4 -l 0,40,30 -N -c spectral -m 100 -M 250
32/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d0${i}$filelmd --div 20 --axtime lt --time 12 --redope edge_y2 -v tk --title "North-most T boundary LMD_MMM Nest $i at 12:00" -O lmd_lon2tbound_d0$i -S png -d --res 200 -i 4 -l 0,40,30 -N -c spectral -m 100 -M 250
33montage mrams_lat1tbound_d0${i}_200.png mrams_lat2tbound_d0${i}_200.png mrams_lon1tbound_d0${i}_200.png mrams_lon2tbound_d0${i}_200.png lmd_lat1tbound_d0${i}_200.png lmd_lat2tbound_d0${i}_200.png lmd_lon1tbound_d0${i}_200.png lmd_lon2tbound_d0${i}_200.png -mode concatenate -tile 4x2 bounding_d0${i}_temp.png
34rm -rf mrams_lat1tbound_d0${i}_200.png mrams_lat2tbound_d0${i}_200.png mrams_lon1tbound_d0${i}_200.png mrams_lon2tbound_d0${i}_200.png lmd_lat1tbound_d0${i}_200.png lmd_lat2tbound_d0${i}_200.png lmd_lon1tbound_d0${i}_200.png lmd_lon2tbound_d0${i}_200.png mrams_lat1tbound_d0${i}.sh mrams_lat2tbound_d0${i}.sh mrams_lon1tbound_d0${i}.sh mrams_lon2tbound_d0${i}.sh lmd_lat1tbound_d0${i}.sh lmd_lat2tbound_d0${i}.sh lmd_lon1tbound_d0${i}.sh lmd_lon2tbound_d0${i}.sh
35done
36fi
37# ------------------
38if [ $2 == uv ]
39then
40for i in 1 2 3
41do
42/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d0${i}$filemrams --div 20 --axtime lt --time 12 --redope edge_x1 -v uvmet --title "West-most UV boundary MRAMS Nest $i at 12:00" -O mrams_lat1uvbound_d0$i -S png -d --res 200 -i 4 -l 0,40,30 -c spectral -m 0 -M 50
43/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d0${i}$filemrams --div 20 --axtime lt --time 12 --redope edge_x2 -v uvmet --title "East-most UV boundary MRAMS Nest $i at 12:00" -O mrams_lat2uvbound_d0$i -S png -d --res 200 -i 4 -l 0,40,30 -N -c spectral -m 0 -M 50
44/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d0${i}$filemrams --div 20 --axtime lt --time 12 --redope edge_y1 -v uvmet --title "South-most UV boundary MRAMS Nest $i at 12:00" -O mrams_lon1uvbound_d0$i -S png -d --res 200 -i 4 -l 0,40,30 -N -c spectral -m 0 -M 50
45/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d0${i}$filemrams --div 20 --axtime lt --time 12 --redope edge_y2 -v uvmet --title "North-most UV boundary MRAMS Nest $i at 12:00" -O mrams_lon2uvbound_d0$i -S png -d --res 200 -i 4 -l 0,40,30 -N -c spectral -m 0 -M 50
46/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d0${i}$filelmd --div 20 --axtime lt --time 12 --redope edge_x1 -v uvmet --title "West-most UV boundary LMD_MMM Nest $i at 12:00" -O lmd_lat1uvbound_d0$i -S png -d --res 200 -i 4 -l 0,40,30 -c spectral -m 0 -M 50
47/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d0${i}$filelmd --div 20 --axtime lt --time 12 --redope edge_x2 -v uvmet --title "East-most UV boundary LMD_MMM Nest $i at 12:00" -O lmd_lat2uvbound_d0$i -S png -d --res 200 -i 4 -l 0,40,30 -N -c spectral -m 0 -M 50
48/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d0${i}$filelmd --div 20 --axtime lt --time 12 --redope edge_y1 -v uvmet --title "South-most UV boundary LMD_MMM Nest $i at 12:00" -O lmd_lon1uvbound_d0$i -S png -d --res 200 -i 4 -l 0,40,30 -N -c spectral -m 0 -M 50
49/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d0${i}$filelmd --div 20 --axtime lt --time 12 --redope edge_y2 -v uvmet --title "North-most UV boundary LMD_MMM Nest $i at 12:00" -O lmd_lon2uvbound_d0$i -S png -d --res 200 -i 4 -l 0,40,30 -N -c spectral -m 0 -M 50
50montage mrams_lat1uvbound_d0${i}_200.png mrams_lat2uvbound_d0${i}_200.png mrams_lon1uvbound_d0${i}_200.png mrams_lon2uvbound_d0${i}_200.png lmd_lat1uvbound_d0${i}_200.png lmd_lat2uvbound_d0${i}_200.png lmd_lon1uvbound_d0${i}_200.png lmd_lon2uvbound_d0${i}_200.png -mode concatenate -tile 4x2 bounding_d0${i}_uv.png
51rm -rf mrams_lat1uvbound_d0${i}_200.png mrams_lat2uvbound_d0${i}_200.png mrams_lon1uvbound_d0${i}_200.png mrams_lon2uvbound_d0${i}_200.png lmd_lat1uvbound_d0${i}_200.png lmd_lat2uvbound_d0${i}_200.png lmd_lon1uvbound_d0${i}_200.png lmd_lon2uvbound_d0${i}_200.png mrams_lat1uvbound_d0${i}.sh mrams_lat2uvbound_d0${i}.sh mrams_lon1uvbound_d0${i}.sh mrams_lon2uvbound_d0${i}.sh lmd_lat1uvbound_d0${i}.sh lmd_lat2uvbound_d0${i}.sh lmd_lon1uvbound_d0${i}.sh lmd_lon2uvbound_d0${i}.sh
52done
53fi
54# ------------------
55fi
56########## ---------
57if [ $1 == sync ]
58then
59echo 'Synchronization'
60for i in [1,2,3]
61do
62/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d0${i}$filemrams --div 20 --axtime lt --time 12 -v tsurf --title "Tsurf MRAMS Nest $i at 12:00" -O mrams_tsurf_d0$i -S png -d --res 200 ---proj moll --mark $lon,$lat
63/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d0${i}$filelmd --div 20 --axtime lt --time 12 -v tsurf --title "Tsurf LMD_MMM Nest $i at 12:00" -O lmd_tsurf_d0$i -S png -d --res 200 --proj moll --mark $lon,$lat
64done
65montage mrams_tsurf_d01_200.png mrams_tsurf_d02_200.png mrams_tsurf_d03_200.png lmd_tsurf_d01_200.png lmd_tsurf_d02_200.png lmd_tsurf_d03_200.png tsurf_sync.png
66rm -f mrams_tsurf_d01_200.png mrams_tsurf_d02_200.png mrams_tsurf_d03_200.png lmd_tsurf_d01_200.png lmd_tsurf_d02_200.png lmd_tsurf_d03_200.png mrams_tsurf_d01_200.sh mrams_tsurf_d02_200.sh mrams_tsurf_d03_200.sh lmd_tsurf_d01_200.sh lmd_tsurf_d02_200.sh lmd_tsurf_d03_200.sh
67
68/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams,$lmdfiles/wrfout_d03$filelmd  --lon $lon --lat $lat --axtime lt -v tsurf --ylabel "Surface Temperature" --xlabel "Local time (h)" --title "NEST 3 SYNCHRONIZATION" -O sync_tsurfnest3 -S png -d --res 200 --labels "MRAMS","LMD_MMM" 
69rm -f sync_tsurfnest3.sh
70
71fi
72########## ---------
73
74if [ $1 == nests ]
75then
76echo 'Safe Reference Landing Site'
77
78if [ $2 == power ]
79then
80echo 'Power spectrum'
81/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d01$filemrams,$mramsfiles/mramsout_d02$filemrams,$mramsfiles/mramsout_d03$filemrams --div 20  --lon $lon --lat $lat --axtime lt --time 8,20 0 -i 4 -l 0.1,10,999 -v uv --analysis fft --logy --logx --xlabel "Wavelength (m)" --ylabel "Spectrum amplitude" --title "MRAMS Nest $i Daytime Power spectrum of horizontal winds at landing site (0-10km AGL)" -O mrams_power_nests -S png -d --res 200 --xmin 10 --xmax 10000 --ymin 1 --ymax 10000 --labels "MRAMS Nest 1","MRAMS Nest 2","MRAMS Nest 3"
82/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d01$filelmd,$lmdfiles/wrfout_d02$filelmd,$lmdfiles/wrfout_d03$filelmd --div 20  --lon $lon --lat $lat --axtime lt --time 8,20 0 -i 4 -l 0.1,10,999 -v uv --analysis fft --logy --logx --xlabel "Wavelength (m)" --ylabel "Spectrum amplitude" --title "LMD_MMM Nest $i Daytime Power spectrum of horizontal winds at landing site (0-10km AGL)" -O lmd_power_nests -S png -d --res 200 --xmin 10 --xmax 10000 --ymin 1 --ymax 10000 --labels "LMD_MMM Nest 1","LMD_MMM Nest 2","LMD_MMM Nest 3"
83rm -f mrams_power_nests.sh lmd_power_nests.sh
84fi
85
86if [ $2 == histo ]
87then
88echo 'Wind histograms'
89/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d01$filemrams,$mramsfiles/mramsout_d02$filemrams,$mramsfiles/mramsout_d03$filemrams -v uv --analysis histo --lon $lon --lat $lat --div 20 --labels "MRAMS Nest 1","MRAMS Nest 2","MRAMS Nest 3" --title "MRAMS horizontal wind histogram at landing site for all time and altitudes" -d -S png -O mrams_histo2d_nests
90/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d01$filelmd,$lmdfiles/wrfout_d02$filelmd,$lmdfiles/wrfout_d03$filelmd -v uv --analysis histo --lon $lon --lat $lat --div 20 --labels "LMD_MMM Nest 1","LMD_MMM Nest 2","LMD_MMM Nest 3" --title "LMD_MMM horizontal wind histogram at landing site for all time and altitudes" -d -S png -O lmd_histo2d_nests
91rm -f mrams_histo2d_nests.sh lmd_histo2d_nests.sh
92fi
93fi
94
95if [ $1 == srlsn3 ]
96then
97echo 'Safe Reference Landing Site inner nest'
98
99if [ $2 == histo ]
100then
101
102echo 'WARNING'
103echo 'WARNING'
104echo 'WARNING'
105echo 'THOSE HIST ARE FOR ALL TIMES, CHECK THAT SYNCHRONIZATION IS OK OR SPECIFY COMMON LOCAL TIMES'
106echo 'THOSE HIST ARE FOR ALL TIMES, CHECK THAT SYNCHRONIZATION IS OK OR SPECIFY COMMON LOCAL TIMES'
107echo 'THOSE HIST ARE FOR ALL TIMES, CHECK THAT SYNCHRONIZATION IS OK OR SPECIFY COMMON LOCAL TIMES'
108
109/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams,$lmdfiles/wrfout_d03$filelmd -v uv --analysis histo --lon $lon --lat $lat --div 20 --labels "MRAMS","LMD_MMM" --title "Inner-Nest Horizontal wind histogram at landing site for all time and altitudes" -d -S png -O comp_histo2d_srlsn3 --ylabel "Density of probability" --xlabel "Horizontal wind speed (m/s)"
110/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams,$lmdfiles/wrfout_d03$filelmd -v uv --analysis histo --lon $lon --lat $lat --div 20 --labels "MRAMS","LMD_MMM" --title "Hor. wind histogram / 10m-600m AGL range: POWERED DESCENT" -d -S png -O comp_histo1d_10m_600m_AGL_srlsn3 --vert 0,10 -i 4 -l 0.01,0.6,20 --ylabel "Density of probability" --xlabel "Horizontal wind speed (m/s)"
111/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams,$lmdfiles/wrfout_d03$filelmd -v uv --analysis histo --lon $lon --lat $lat --div 20 --labels "MRAMS","LMD_MMM" --title "Hor. wind histogram / 1km-4km AML range: MAIN PARACHUTE DESCENT" -d -S png -O comp_histo1d_1km_4km_AML_srlsn3 --vert 0,20 -i 3 -l 1,4,20 --ylabel "Density of probability" --xlabel "Horizontal wind speed (m/s)"
112/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams,$lmdfiles/wrfout_d03$filelmd -v uv --analysis histo --lon $lon --lat $lat --div 20 --labels "MRAMS","LMD_MMM" --title "Hor. wind histogram / 4km-10km AML range: DROGUE AND PARACHUTE DEPLOYMENT" -d -S png -O comp_histo1d_4km_10km_AML_srlsn3 --vert 0,20 -i 3 -l 4,10,20 --ylabel "Density of probability" --xlabel "Horizontal wind speed (m/s)"
113/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams,$lmdfiles/wrfout_d03$filelmd -v uv --analysis histo --lon $lon --lat $lat --div 20 --labels "MRAMS","LMD_MMM" --title "Hor. wind histogram / 10km-40km AML range: DECELERATION PHASE" -d -S png -O comp_histo1d_10km_40km_AML_srlsn3 --vert 0,20 -i 3 -l 10,40,20 --ylabel "Density of probability" --xlabel "Horizontal wind speed (m/s)"
114montage comp_histo1d_10m_600m_AGL_srlsn3_200.png comp_histo1d_1km_4km_AML_srlsn3_200.png comp_histo1d_4km_10km_AML_srlsn3_200.png comp_histo1d_10km_40km_AML_srlsn3_200.png -mode concatenate -tile 2x2 comp_histo1d_decomposed_srlsn3.png
115rm -f comp_histo2d_srlsn3.sh comp_histo1d_10m_600m_AGL_srlsn3.sh comp_histo1d_1km_4km_AML_srlsn3.sh comp_histo1d_4km_10km_AML_srlsn3.sh comp_histo1d_10km_40km_AML_srlsn3.sh comp_histo1d_10m_600m_AGL_srlsn3_200.png comp_histo1d_1km_4km_AML_srlsn3_200.png comp_histo1d_4km_10km_AML_srlsn3_200.png comp_histo1d_10km_40km_AML_srlsn3_200.png
116
117#/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams,$lmdfiles/wrfout_d03$filelmd -v tsurf --lon $lon --lat $lat --div 20 --labels "MRAMS","LMD_MMM" --title "Inner-Nest Surface temperature histogram" -d -S png -O comp_histo1d_tsurf_srlsn3
118#rm -f comp_histo1d_tsurf_srlsn3.sh
119fi
120
121if [ $2 == power ]
122then
123/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams,$lmdfiles/wrfout_d03$filelmd  --lon $lon --lat $lat --axtime lt --time 8,20 0 -i 4 -l 0.1,10,999 -v uv --analysis fft --logy --logx --xlabel "Wavelength (m)" --ylabel "Spectrum amplitude" --title "Inner-Nest Daytime Power spectrum of horizontal winds at landing site (0-10km AGL)" -O comp_power_srlsn3 -S png -d --res 200 --xmin 10 --xmax 10000 --ymin 1 --ymax 10000 --labels "MRAMS","LMD_MMM" 
124rm -f comp_power_srlsn3.sh
125fi
126
127if [ $2 == series ]
128then
129
130/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams,$lmdfiles/wrfout_d03$filelmd -v UV --vert 0 --lon $lon --lat $lat --div 20 --labels "MRAMS","LMD_MMM" --title "Horizontal winds at 10m AGL (final drop and touchdown)" -d -S png -O comp_uv10m_serie_srlsn3 -i 4 -l 0.01 --xlabel "Local time (h)" --axtime lt --ylabel "Horizontal wind speed (m/s)"
131rm -f comp_uv10m_serie_srlsn3.sh
132
133/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams,$lmdfiles/wrfout_d03$filelmd -v tk --vert 0 --lon $lon --lat $lat --div 20 --labels "MRAMS","LMD_MMM" --title "Temperature at 10m AGL (final drop and touchdown)" -d -S png -O comp_temp10m_serie_srlsn3 -i 4 -l 0.01 --xlabel "Local time (h)" --axtime lt --ylabel "Temperature (K)"
134rm -f comp_temp10m_serie_srlsn3.sh
135
136/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams,$lmdfiles/wrfout_d03$filelmd -v UV --lon $lon --lat $lat --div 20 --title "MRAMS","LMD_MMM" -d -S png -O comp_uv_serie2d_srlsn3 -i 4 -l 0.01,40,40 --xlabel "Local time (h)" --axtime lt --ylabel "Altitude above surface (km)" --xmin 5 --xmax 24 -m 0 -M 75
137rm -f comp_uv_serie2d_srlsn3.sh
138
139/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams,$lmdfiles/wrfout_d03$filelmd -v tk --lon $lon --lat $lat --div 20 --title "MRAMS","LMD_MMM" -d -S png -O comp_temp_serie2d_srlsn3 -i 4 -l 0.01,40,40 --xlabel "Local time (h)" --axtime lt --ylabel "Altitude above surface (km)" -m 105 -M 250 -c spectral --xmin 5 --xmax 24
140rm -f comp_temp_serie2d_srlsn3.sh
141
142/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams,$lmdfiles/wrfout_d03$filelmd -v hodograph --vert 0 --lon $lon --lat $lat --div 20 --labels "MRAMS","LMD_MMM" --title "Hodograph of horizontal winds at 10m AGL (final drop and touchdown)" -d -S png -O comp_hodograph_srlsn3 -i 4 -l 0.01 --xlabel "U component (m/s)" --axtime lt --ylabel "V component (m/s)" --xmin -15 --xmax 15 --ymin -15 --ymax 15
143rm -f comp_hodograph_srlsn3.sh
144
145/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams,$lmdfiles/wrfout_d03$filelmd -v hodograph_2 --vert 0 --lon $lon --lat $lat --div 20 --labels "MRAMS","LMD_MMM" --title "Horiztonal winds rotation at 10m AGL (final drop and touchdown)" -d -S png -O comp_hodograph_2_srlsn3 -i 4 -l 0.01 --xlabel "Local Time (h)" --axtime lt --ylabel "Arbitrary wind rotation (degrees)"
146rm -f comp_hodograph_2_srlsn3.sh
147
148fi
149
150if [ $2 == profiles ]
151then
152
153/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams,$lmdfiles/wrfout_d03$filelmd -v UV --lon $lon --lat $lat --div 20 --labels "MRAMS","LMD_MMM" --title "Horizontal winds at LT 13:00" -d -S png -O comp_uv_profile_lt13_srlsn3 -i 4 -l 0.01,40,40 --xlabel "Horizontal wind speed (m/s)" --axtime lt --ylabel "Altitude above surface (km)" --axtime lt --time 13
154rm -f comp_uv_profile_lt13_srlsn3.sh
155
156/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams,$lmdfiles/wrfout_d03$filelmd -v tk --lon $lon --lat $lat --div 20 --labels "MRAMS","LMD_MMM" --title "Temperature serie at LT 13:00" -d -S png -O comp_temp_profile_lt13_srlsn3 -i 4 -l 0.01,40,40 --xlabel "Temperature (K)" --axtime lt --ylabel "Altitude above surface (km)" --axtime lt --time 13
157rm -f comp_temp_profile_lt13_srlsn3.sh
158
159fi
160
161if [ $2 == surface ]
162then
163
164/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams,$lmdfiles/wrfout_d03$filelmd -v PSFC --axtime lt --title "MRAMS Inner Nest surface pressure","LMD_MMM Inner Nest surface pressure" -d -S png -O comp_pressure_srlsn3 --blat -20 --div 20 --res 200 --lon $lon --lat $lat --xlabel "Local Time (h)" --ylabel "Surface Pressure (Pa)"
165rm -f comp_pressure_srlsn3.sh
166/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams,$lmdfiles/wrfout_d03$filelmd -v CDU --axtime lt --title "MRAMS Inner Nest surface stress","LMD_MMM Inner Nest surface stress" -d -S png -O comp_surfstress_srlsn3 --blat -20 --div 20 --res 200 --lon $lon --lat $lat --xlabel "Local Time (h)" --ylabel "Surface stress (m/s)"
167rm -f comp_surfstress_srlsn3.sh
168/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams,$lmdfiles/wrfout_d03$filelmd -v HFX --axtime lt --title "MRAMS Inner Nest sensib. heat flux ","LMD_MMM Inner Nest sensib. heat flux" -d -S png -O comp_sensibflux_srlsn3 --blat -20 --div 20 --res 200 --lon $lon --lat $lat --xlabel "Local Time (h)" --ylabel "Sensible heat flux (K.m/s)"
169rm -f comp_sensibflux_srlsn3.sh
170
171fi
172
173fi
174
175########## ---------
176
177if [ $1 == maps ]
178then
179
180if [ $2 == winds ]
181then
182
183/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams -v UV --axtime lt --time 14 -i 4 -l 0.01 --title "MRAMS horizontal winds at 10m AGL" -d -S png -W --vert 0 -O mrams_winds_10mAGL_maps --proj spstere --blat -20 --div 20 --res 200 -s 5 --mark $lon,$lat
184/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams  -v UV --axtime lt --time 14 -i 4 -l 0.6 --title "MRAMS horizontal winds at 600m AGL" -d -S png -W --vert 0 -O mrams_winds_600mAGL_maps --proj spstere --blat -20 --div 20 --res 200 -s 5 --mark $lon,$lat
185/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams  -v UV --axtime lt --time 14 -i 3 -l 1 --title "MRAMS horizontal winds at 1km AML" -d -S png -W --vert 0 -O mrams_winds_1kmAML_maps --proj spstere --blat -20 --div 20 --res 200 -s 5 --mark $lon,$lat
186/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams  -v UV --axtime lt --time 14 -i 3 -l 4 --title "MRAMS horizontal winds at 4km AML" -d -S png -W --vert 0 -O mrams_winds_4kmAML_maps --proj spstere --blat -20 --div 20 --res 200 -s 5 --mark $lon,$lat
187/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams -v UV --axtime lt --time 14 -i 3 -l 10 --title "MRAMS horizontal winds at 10km AML" -d -S png -W --vert 0 -O mrams_winds_10kmAML_maps --proj spstere --blat -20 --div 20 --res 200 -s 5 --mark $lon,$lat
188/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams  -v UV --axtime lt --time 14 -i 3 -l 20 --title "MRAMS horizontal winds at 20km AML" -d -S png -W --vert 0 -O mrams_winds_20kmAML_maps --proj spstere --blat -20 --div 20 --res 200 -s 5 --mark $lon,$lat
189
190/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d03$filelmd -v UV --axtime lt --time 14 -i 4 -l 0.01 --title "LMD_MMM horizontal winds at 10m AGL" -d -S png -W --vert 0 -O lmd_winds_10mAGL_maps --proj spstere --blat -20 --div 20 --res 200 -s 5 --mark $lon,$lat
191/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d03$filelmd  -v UV --axtime lt --time 14 -i 4 -l 0.6 --title "LMD_MMM horizontal winds at 600m AGL" -d -S png -W --vert 0 -O lmd_winds_600mAGL_maps --proj spstere --blat -20 --div 20 --res 200 -s 5 --mark $lon,$lat
192/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d03$filelmd  -v UV --axtime lt --time 14 -i 3 -l 1 --title "LMD_MMM horizontal winds at 1km AML" -d -S png -W --vert 0 -O lmd_winds_1kmAML_maps --proj spstere --blat -20 --div 20 --res 200 -s 5 --mark $lon,$lat
193/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d03$filelmd  -v UV --axtime lt --time 14 -i 3 -l 4 --title "LMD_MMM horizontal winds at 4km AML" -d -S png -W --vert 0 -O lmd_winds_4kmAML_maps --proj spstere --blat -20 --div 20 --res 200 -s 5 --mark $lon,$lat
194/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d03$filelmd -v UV --axtime lt --time 14 -i 3 -l 10 --title "LMD_MMM horizontal winds at 10km AML" -d -S png -W --vert 0 -O lmd_winds_10kmAML_maps --proj spstere --blat -20 --div 20 --res 200 -s 5 --mark $lon,$lat
195/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d03$filelmd  -v UV --axtime lt --time 14 -i 3 -l 20 --title "LMD_MMM horizontal winds at 20km AML" -d -S png -W --vert 0 -O lmd_winds_20kmAML_maps --proj spstere --blat -20 --div 20 --res 200 -s 5 --mark $lon,$lat
196
197montage mrams_winds_10mAGL_maps_200.png lmd_winds_10mAGL_maps_200.png mrams_winds_600mAGL_maps_200.png lmd_winds_600mAGL_maps_200.png mrams_winds_1kmAML_maps_200.png lmd_winds_1kmAML_maps_200.png -mode concatenate -tile 2x3 comp_winds_maps_10m_1km.png
198montage mrams_winds_4kmAML_maps_200.png lmd_winds_4kmAML_maps_200.png mrams_winds_10kmAML_maps_200.png lmd_winds_10kmAML_maps_200.png mrams_winds_20kmAML_maps_200.png lmd_winds_20kmAML_maps_200.png -mode concatenate -tile 2x3 comp_winds_maps_4km_20km.png
199rm -f mrams_winds_10mAGL_maps.sh lmd_winds_10mAGL_maps_200.png mrams_winds_600mAGL_maps_200.png lmd_winds_600mAGL_maps_200.png mrams_winds_1kmAML_maps_200.png lmd_winds_1kmAML_maps_200.png mrams_winds_4kmAML_maps_200.png lmd_winds_4kmAML_maps_200.png mrams_winds_10kmAML_maps_200.png lmd_winds_10kmAML_maps_200.png mrams_winds_20kmAML_maps_200.png lmd_winds_20kmAML_maps_200.png mrams_winds_10mAGL_maps.sh lmd_winds_10mAGL_maps.sh mrams_winds_600mAGL_maps.sh lmd_winds_600mAGL_maps.sh mrams_winds_1kmAML_maps.sh lmd_winds_1kmAML_maps.sh mrams_winds_4kmAML_maps.sh lmd_winds_4kmAML_maps.sh mrams_winds_10kmAML_maps.sh lmd_winds_10kmAML_maps.sh mrams_winds_20kmAML_maps.sh lmd_winds_20kmAML_maps.sh
200
201fi
202
203if [ $2 == surface ]
204then
205#/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams,$lmdfiles/wrfout_d03$filelmd -v PSFC --axtime lt --time 14 --title "MRAMS Inner Nest surface pressure","LMD_MMM Inner Nest surface pressure" -d -S png -O comp_pressure_maps --proj spstere --blat -20 --div 20 --res 200 --mark $lon,$lat
206#rm -f comp_pressure_maps.sh
207#/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams,$lmdfiles/wrfout_d03$filelmd -v CDU --axtime lt --time 14 --title "MRAMS Inner Nest surface stress","LMD_MMM Inner Nest surface stress" -d -S png -O comp_surfstress_maps --proj spstere --blat -20 --div 20 --res 200 --mark $lon,$lat
208#rm -f comp_surfstress_maps.sh
209#/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemrams,$lmdfiles/wrfout_d03$filelmd -v HFX --axtime lt --time 14 --title "MRAMS Inner Nest sensib. heat flux ","LMD_MMM Inner Nest sensib. heat flux" -d -S png -O comp_sensibflux_maps --proj spstere --blat -20 --div 20 --res 200 --mark $lon,$lat
210#rm -f comp_sensibflux_maps.sh
211
212/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d03$filelmd -v HFX --axtime lt --time 14 --title "MRAMS Inner Nest sensib. heat flux ","LMD_MMM Inner Nest sensib. heat flux" --proj spstere --blat -20 --div 20 --mark $lon,$lat
213
214fi
215
216fi
217
218########## ---------
219
220if [ $1 == variability ]
221then
222
223if [ $2 == histo ]
224then
225/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemramsa,$mramsfiles/mramsout_d03$filemrams,$mramsfiles/mramsout_d03$filemramsb -v uv --analysis histo --lon $lon --lat $lat --div 20 --labels "MRAMS day-1","MRAMS day 0","MRAMS day+1" --title "MRAMS Inner-Nest Horizontal wind histogram at landing site for all time and altitudes" -d -S png -O mrams_histo2d_variability --ylabel "Density of probability" --xlabel "Horizontal wind speed (m/s)" --ymin 0 --ymax 0.05
226/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemramsa,$mramsfiles/mramsout_d03$filemrams,$mramsfiles/mramsout_d03$filemramsb -v uv --analysis histo --lon $lon --lat $lat --div 20 --labels "MRAMS day-1","MRAMS day 0","MRAMS day+1" --title "Hor. wind histogram / 10m-600m AGL range: POWERED DESCENT" -d -S png -O mrams_histo1d_10m_600m_AGL_variability --vert 0,10 -i 4 -l 0.01,0.6,20 --ylabel "Density of probability" --xlabel "Horizontal wind speed (m/s)"
227/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemramsa,$mramsfiles/mramsout_d03$filemrams,$mramsfiles/mramsout_d03$filemramsb -v uv --analysis histo --lon $lon --lat $lat --div 20 --labels "MRAMS day-1","MRAMS day 0","MRAMS day+1" --title "Hor. wind histogram / 1km-4km AML range: MAIN PARACHUTE DESCENT" -d -S png -O mrams_histo1d_1km_4km_AML_variability --vert 0,20 -i 3 -l 1,4,20 --ylabel "Density of probability" --xlabel "Horizontal wind speed (m/s)"
228/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemramsa,$mramsfiles/mramsout_d03$filemrams,$mramsfiles/mramsout_d03$filemramsb -v uv --analysis histo --lon $lon --lat $lat --div 20 --labels "MRAMS day-1","MRAMS day 0","MRAMS day+1" --title "Hor. wind histogram / 4km-10km AML range: DROGUE AND PARACHUTE DEPLOYMENT" -d -S png -O mrams_histo1d_4km_10km_AML_variability --vert 0,20 -i 3 -l 4,10,20 --ylabel "Density of probability" --xlabel "Horizontal wind speed (m/s)"
229/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemramsa,$mramsfiles/mramsout_d03$filemrams,$mramsfiles/mramsout_d03$filemramsb -v uv --analysis histo --lon $lon --lat $lat --div 20 --labels "MRAMS day-1","MRAMS day 0","MRAMS day+1" --title "Hor. wind histogram / 10km-40km AML range: DECELERATION PHASE" -d -S png -O mrams_histo1d_10km_40km_AML_variability --vert 0,20 -i 3 -l 10,40,20 --ylabel "Density of probability" --xlabel "Horizontal wind speed (m/s)"
230montage mrams_histo1d_10m_600m_AGL_variability_200.png mrams_histo1d_1km_4km_AML_variability_200.png mrams_histo1d_4km_10km_AML_variability_200.png mrams_histo1d_10km_40km_AML_variability_200.png -mode concatenate -tile 2x2 mrams_histo1d_decomposed_variability.png
231rm -f mrams_histo1d_10m_600m_AGL_variability_200.png mrams_histo1d_1km_4km_AML_variability_200.png mrams_histo1d_4km_10km_AML_variability_200.png mrams_histo1d_10km_40km_AML_variability_200.png mrams_histo1d_10m_600m_AGL_variability.sh mrams_histo1d_1km_4km_AML_variability.sh mrams_histo1d_4km_10km_AML_variability.sh mrams_histo1d_10km_40km_AML_variability.sh
232
233/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d03$filelmda,$lmdfiles/wrfout_d03$filelmd,$lmdfiles/wrfout_d03$filelmdb -v uv --analysis histo --lon $lon --lat $lat --div 20 --labels "LMD_MMM day-1","LMD_MMM day 0","LMD_MMM day+1" --title "LMD_MMM Inner-Nest Horizontal wind histogram at landing site for all time and altitudes" -d -S png -O lmd_histo2d_variability --ylabel "Density of probability" --xlabel "Horizontal wind speed (m/s)" --ymin 0 --ymax 0.05
234/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d03$filelmda,$lmdfiles/wrfout_d03$filelmd,$lmdfiles/wrfout_d03$filelmdb -v uv --analysis histo --lon $lon --lat $lat --div 20 --labels "LMD_MMM day-1","LMD_MMM day 0","LMD_MMM day+1" --title "Hor. wind histogram / 10m-600m AGL range: POWERED DESCENT" -d -S png -O lmd_histo1d_10m_600m_AGL_variability --vert 0,10 -i 4 -l 0.01,0.6,20 --ylabel "Density of probability" --xlabel "Horizontal wind speed (m/s)"
235/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d03$filelmda,$lmdfiles/wrfout_d03$filelmd,$lmdfiles/wrfout_d03$filelmdb -v uv --analysis histo --lon $lon --lat $lat --div 20 --labels "LMD_MMM day-1","LMD_MMM day 0","LMD_MMM day+1" --title "Hor. wind histogram / 1km-4km AML range: MAIN PARACHUTE DESCENT" -d -S png -O lmd_histo1d_1km_4km_AML_variability --vert 0,20 -i 3 -l 1,4,20 --ylabel "Density of probability" --xlabel "Horizontal wind speed (m/s)"
236/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d03$filelmda,$lmdfiles/wrfout_d03$filelmd,$lmdfiles/wrfout_d03$filelmdb -v uv --analysis histo --lon $lon --lat $lat --div 20 --labels "LMD_MMM day-1","LMD_MMM day 0","LMD_MMM day+1" --title "Hor. wind histogram / 4km-10km AML range: DROGUE AND PARACHUTE DEPLOYMENT" -d -S png -O lmd_histo1d_4km_10km_AML_variability --vert 0,20 -i 3 -l 4,10,20 --ylabel "Density of probability" --xlabel "Horizontal wind speed (m/s)"
237/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d03$filelmda,$lmdfiles/wrfout_d03$filelmd,$lmdfiles/wrfout_d03$filelmdb -v uv --analysis histo --lon $lon --lat $lat --div 20 --labels "LMD_MMM day-1","LMD_MMM day 0","LMD_MMM day+1" --title "Hor. wind histogram / 10km-40km AML range: DECELERATION PHASE" -d -S png -O lmd_histo1d_10km_40km_AML_variability --vert 0,20 -i 3 -l 10,40,20 --ylabel "Density of probability" --xlabel "Horizontal wind speed (m/s)"
238montage lmd_histo1d_10m_600m_AGL_variability_200.png lmd_histo1d_1km_4km_AML_variability_200.png lmd_histo1d_4km_10km_AML_variability_200.png lmd_histo1d_10km_40km_AML_variability_200.png -mode concatenate -tile 2x2 lmd_histo1d_decomposed_variability.png
239rm -f lmd_histo1d_10m_600m_AGL_variability_200.png lmd_histo1d_1km_4km_AML_variability_200.png lmd_histo1d_4km_10km_AML_variability_200.png lmd_histo1d_10km_40km_AML_variability_200.png lmd_histo1d_10m_600m_AGL_variability.sh lmd_histo1d_1km_4km_AML_variability.sh lmd_histo1d_4km_10km_AML_variability.sh lmd_histo1d_10km_40km_AML_variability.sh
240
241montage mrams_histo2d_variability_200.png lmd_histo2d_variability_200.png -mode concatenate -tile 2x1 comp_histo2d_variability.png
242rm -f mrams_histo2d_variability_200.png lmd_histo2d_variability_200.png mrams_histo2d_variability.sh lmd_histo2d_variability.sh
243fi
244
245if [ $2 == density ]
246then
247/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemramsa,$mramsfiles/mramsout_d03$filemramsb -v uv --lon $lon --lat $lat --div 60 --labels "dummy","dummy","MRAMS day-1/0","dummy","dummy","MRAMS day+1/0" --title "MRAMS day-to-day horizontal winds variability at Landing Site" -d -S png -O mrams_histodensity_uv_variability --ope - --fref $mramsfiles/mramsout_d03$filemrams  --analysis histodensity --xmin -15 --xmax 15 --ymin 0 --ymax 0.35
248/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d03$filelmda,$lmdfiles/wrfout_d03$filelmdb -v uv --lon $lon --lat $lat --div 20 --labels "dummy","dummy","LMD_MMM day-1/0","dummy","dummy","LMD_MMM day+1/0" --title "LMD_MMM day-to-day horizontal winds variability at Landing Site" -d -S png -O lmd_histodensity_uv_variability --ope - --fref $lmdfiles/wrfout_d03$filelmd  --analysis histodensity  --xmin -15 --xmax 15 --ymin 0 --ymax 0.35
249montage mrams_histodensity_uv_variability_200.png lmd_histodensity_uv_variability_200.png -mode concatenate -tile 2x1 comp_histodensity_uv_variability.png
250rm -f mrams_histodensity_uv_variability_200.png lmd_histodensity_uv_variability_200.png mrams_histodensity_uv_variability.sh lmd_histodensity_uv_variability.sh
251
252/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $mramsfiles/mramsout_d03$filemramsa,$mramsfiles/mramsout_d03$filemramsb -v tk --lon $lon --lat $lat --div 200 --labels "dummy","dummy","MRAMS day-1/0","dummy","dummy","MRAMS day+1/0" --title "MRAMS day-to-day temperature variability at Landing Site" -d -S png -O mrams_histodensity_tk_variability --ope - --fref $mramsfiles/mramsout_d03$filemrams --analysis histodensity --xmin -5 --xmax 5 --ymin 0 --ymax 0.8
253/san0/acolmd/SVN/trunk//UTIL/PYTHON/pp.py -f $lmdfiles/wrfout_d03$filelmda,$lmdfiles/wrfout_d03$filelmdb -v tk --lon $lon --lat $lat --div 40 --labels "dummy","dummy","LMD_MMM day-1/0","dummy","dummy","LMD_MMM day+1/0" --title "LMD_MMM day-to-day temperature variability at Landing Site" -d -S png -O lmd_histodensity_tk_variability --ope - --fref $lmdfiles/wrfout_d03$filelmd --analysis histodensity  --xmin -5 --xmax 5 --ymin 0 --ymax 0.8
254montage mrams_histodensity_tk_variability_200.png lmd_histodensity_tk_variability_200.png -mode concatenate -tile 2x1 comp_histodensity_tk_variability.png
255rm -f mrams_histodensity_tk_variability_200.png lmd_histodensity_tk_variability_200.png mrams_histodensity_tk_variability.sh lmd_histodensity_tk_variability.sh
256fi
257
258fi
Note: See TracBrowser for help on using the repository browser.