|
Last change
on this file since 796 was
349,
checked in by aslmd, 13 years ago
|
|
PYTHON GRAPHICS: A REFERENCE VERSION FOR BOTH GCM AND MESOSCALE USE!
Merged recent changes by TN for better GCM handling + 1D/2D plots with AS mesoscale version [winds.py]
Now
- there is a common function named planetoplot in planetoplot.py which supports both GCM and mesoscale plots
this will be a common reference for future development
- mesoscale users might want to use the script meso.py
- GCM users might want to use the script gcm.py
Todo
- a few capabilities that were in winds.py are broken in meso.py, this will be fixed soon.
the 'old' reference script winds.py is still here anyway!
- there is still a bug with 'ortho' projection with GCM fields. still searching...
Misc
- A few modifications were done in myplot.py and mymath.py
- Added a simple stupid little.py script illustrating very simple use of Python graphics
- Added the script doing statistical analysis for dust devils
|
-
Property svn:executable set to
*
|
|
File size:
349 bytes
|
| Line | |
|---|
| 1 | #! /usr/bin/env python |
|---|
| 2 | |
|---|
| 3 | import matplotlib.pyplot as plt |
|---|
| 4 | import numpy as np |
|---|
| 5 | |
|---|
| 6 | vvv = np.array([ 8.9, 7.3, 3.9, 2.0]) |
|---|
| 7 | rrr = np.array([ 2.0, 5.0, 8.0,11.0]) |
|---|
| 8 | |
|---|
| 9 | logv = np.log(vvv) |
|---|
| 10 | logr = np.log(rrr) |
|---|
| 11 | |
|---|
| 12 | vvv1 = 1./rrr |
|---|
| 13 | vvv2 = 1./np.sqrt(rrr) |
|---|
| 14 | |
|---|
| 15 | logv1 = 3. + np.log(vvv1) |
|---|
| 16 | logv2 = 3. + np.log(vvv2) |
|---|
| 17 | |
|---|
| 18 | plt.plot(logv,logr,'bo',logv1,logr,'r',logv2,logr,'y') |
|---|
| 19 | plt.show() |
|---|
Note: See
TracBrowser
for help on using the repository browser.
