- Timestamp:
- Jan 24, 2012, 6:35:35 PM (13 years ago)
- File:
-
- 1 edited
-
trunk/UTIL/PYTHON/README.PP (modified) (9 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/UTIL/PYTHON/README.PP
r489 r504 12 12 ************************************** 13 13 14 ************************************ *****************************15 MAPMODE 1... MAPPING MODE...SIMPLE EXAMPLES on a SAMPLE GCM FILE16 ************************************ *****************************14 ************************************ 15 SIMPLE EXAMPLES on a SAMPLE GCM FILE 16 ************************************ 17 17 Goal: The simplest, most minimal example. Mapping topography. 18 18 pp.py -f diagfired.nc … … 51 51 pp.py -f diagfired.nc -v tauice -m 0.2 -M 0.9 -H -b vishires -p ortho 52 52 53 Goal: I want the same map, but projected on the sphere and centered on longitude 100° 54 pp.py -f diagfired.nc -v tauice -m 0.2 -M 0.9 -H -b vishires -p ortho --blon -100 55 53 56 Goal: I want the same map, but projected with north polar stereographic view 54 57 pp.py -f diagfired.nc -v tauice -m 0.2 -M 0.9 -H -b vishires -p npstere … … 67 70 68 71 Goal: I want to plot results for two different times in the file next to one another 69 pp.py -f diagfi .nc -v tsurf --time 4 --time 772 pp.py -f diagfired.nc -v tsurf --time 0.5 --time 0.9 70 73 71 74 Goal: I want to plot averaged results in the file from one time to another time 72 pp.py -f diagfi .nc -v tsurf --time 4,775 pp.py -f diagfired.nc -v tsurf --time 0.5,0.9 73 76 74 77 Goal: I want to plot the minimum/maximum value over times stored in the file 75 pp.py -f diagfi.nc -v tsurf --redope mint 76 pp.py -f diagfi.nc -v tsurf --redope maxt 78 pp.py -f diagfired.nc -v tsurf --redope mint 79 pp.py -f diagfired.nc -v tsurf --redope maxt 80 81 Goal: I want to plot a section of temperature at longitude 0 82 pp.py -f diagfired.nc --var temp --lon 0 --time 0.5 83 84 Goal: I want to plot a section of zonally-averaged temperature 85 pp.py -f diagfired.nc -v temp --lon -180,180 --time 0.5 86 87 Goal: I want to plot a section of zonally-averaged temperature with contours of zonally-averaged zonal wind 88 pp.py -f diagfired.nc -v temp -w u --lon -180,180 --time 0.5 89 90 Goal: I want to plot a section of zonally-averaged temperature and zonally-averaged zonal wind 91 pp.py -f diagfired.nc -v temp,u --lon -180,180 --time 0.5 77 92 78 93 Goal: I want to plot a globally-averaged 1D temperature profile 79 pp.py -f diagfi .nc -v temp --time 4--lat -90,90 --lon -180,18094 pp.py -f diagfired.nc -v temp --time 0.5 --lat -90,90 --lon -180,180 80 95 81 96 Goal: I want to overplot few globally-averaged 1D temperature profiles at different times 82 pp.py -f diagfi.nc -v temp --time 4 --time 7 --lat -90,90 --lon -180,180 83 84 [only mesoscale for the moment] 97 pp.py -f diagfired.nc -v temp --time 0.5 --time 0.9 --lat -90,90 --lon -180,180 98 99 Goal: I want to overplot lat=0,lon=0 1D temperature profiles at different times 100 pp.py -f diagfired.nc -v temp --time 0.5 --time 0.9 --lat 0 --lon 0 --ymax 20 --xmin 180. 101 102 ******************** [specific mesoscale] 103 85 104 Goal: I want to plot results for two different LOCAL times in the file next to one another 86 pp.py -f wrfout**** -v TSURF --time 4 -- time 7 --axtime lt 105 pp.py -f wrfout_d01_2024-05-03_01:00:00 -v TSURF --time 4 --time 7 --axtime lt 106 107 Goal: I want to plot a section of temperature at longitude -120 between AMR altitudes 0 km and 30 km with 50 levels 108 pp.py -f wrfout_d01_2024-05-03_01:00:00 -v tk --time 4 --axtime lt --lon -120 -i 3 -l 0,30,50 109 110 Goal: I want to plot a section of temperature at longitude -120 between ALS altitudes 0 km and 30 km with 50 levels 111 pp.py -f wrfout_d01_2024-05-03_01:00:00 -v tk --time 4 --axtime lt --lon -120 -i 4 -l 0,30,50 87 112 88 113 *********************************************************************************** … … 152 177 relaxed towards prescribed GCM fields. 153 178 179 154 180 *********************************************************************************** 155 181 Simple 2D plot: Zonal mean. … … 158 184 Goal: 159 185 160 Plot the zonal mean temperature from a netcdf fiel representing one month.161 162 Command: 163 164 gcm.py -f POLAR_NIGHT_RUN/diagfi16.nc --var temp --lon 180,-180 --time 0,65186 Plot the zonal mean temperature from a netcdf field representing one month. 187 188 Command: 189 190 pp.py -f POLAR_NIGHT_RUN/diagfi16.nc --var temp --lon 180,-180 --time 0,65 165 191 166 192 Note: … … 180 206 Command: 181 207 182 gcm.py -f POLAR_NIGHT_RUN/diagfi16.nc --var temp --lon -180,180 --time 0,65 -i 4208 pp.py -f POLAR_NIGHT_RUN/diagfi16.nc --var temp --lon -180,180 --time 0,65 -i 4 183 209 184 210 Note: … … 199 225 Command: 200 226 201 gcm.py -f POLAR_NIGHT_RUN/stats16.nc --var temp --lon -180,180 --time 1 -i 4 --fref POLAR_NIGHT_REF/stats16.nc --operation - --mope -2 --Mope 2 --title "Polar temperatures with new parametrizations" --titleref "Reference run"227 pp.py -f POLAR_NIGHT_RUN/stats16.nc --var temp --lon -180,180 --time 1 -i 4 --fref POLAR_NIGHT_REF/stats16.nc --operation - --mope -2 --Mope 2 --title "Polar temperatures with new parametrizations" --titleref "Reference run" 202 228 203 229 Note: … … 207 233 One can combine this command with projections and means, so that for example, to compare co2 depletion at the south pole: 208 234 209 gcm.py -f POLAR_NIGHT_RUN/start16.nc --var co2 --vert 0,150 --proj spstere --time 1 -i 4 --fref POLAR_NIGHT_REF/start16.nc --operation - --mope -0.5 --Mope 0.5 --title "Polar co2 with new parametrizations" --titleref "Reference run"235 pp.py -f POLAR_NIGHT_RUN/start16.nc --var co2 --vert 0,150 --proj spstere --time 1 -i 4 --fref POLAR_NIGHT_REF/start16.nc --operation - --mope -0.5 --Mope 0.5 --title "Polar co2 with new parametrizations" --titleref "Reference run" 210 236 211 237 *********************************************************************************** … … 219 245 Command: 220 246 221 gcm.py -f TES.MappedClimatology.nadir.MY25.nc --var T_nadir_day --lat -76. --time 90 --ymin 500 --ymax 1 -m 128 -M 148 -H247 pp.py -f TES.MappedClimatology.nadir.MY25.nc --var T_nadir_day --lat -76. --time 90 --ymin 500 --ymax 1 -m 128 -M 148 -H 222 248 223 249 Note:
Note: See TracChangeset
for help on using the changeset viewer.
