Changeset 475 for trunk/UTIL/PYTHON/myplot.py

Timestamp:

Dec 15, 2011, 4:05:49 PM (13 years ago)

Author:

aslmd

Message:

UTIL: Python graphics: 1. corrected locations lat/lon for meso. 2. fixed yaxis settings which did not work in some 1D and 2D cases. 3. moved 1D stuff in imov loop, so the loop is now consistent and generic. Also updated a bit farm_tour to avoid displaying e.g. grep as running job.

File:

• trunk/UTIL/PYTHON/myplot.py (modified) (3 diffs)

trunk/UTIL/PYTHON/myplot.py

@@ -414 +414 @@
     if typefile in ['mesoapi','meso']:
         [lon2d,lat2d] = getcoord2d(nc)
-        lon2d = dumpbdy(lon2d,6)
-        lat2d = dumpbdy(lat2d,6)
     elif typefile in ['gcm']: 
         [lon2d,lat2d] = getcoord2d(nc,nlat="latitude",nlon="longitude",is1d=True)
@@ -958 +956 @@
 
 ## Author: AS
-def bidimfind(lon2d,lat2d,vlon,vlat):
+def bidimfind(lon2d,lat2d,vlon,vlat,file=None):
    import numpy as np
+   import matplotlib.pyplot as mpl
    if vlat is None:    array = (lon2d - vlon)**2
    elif vlon is None:  array = (lat2d - vlat)**2
@@ -965 +964 @@
    idy,idx = np.unravel_index( np.argmin(array), lon2d.shape )
    if vlon is not None:
-       #print lon2d[idy,idx],vlon
-       if (np.abs(lon2d[idy,idx]-vlon)) > 5: errormess("longitude not found ",printvar=lon2d)
+      if (np.abs(lon2d[idy,idx]-vlon)) > 5: errormess("longitude not found ",printvar=lon2d)
    if vlat is not None:
-       #print lat2d[idy,idx],vlat
-       if (np.abs(lat2d[idy,idx]-vlat)) > 5: errormess("latitude not found ",printvar=lat2d)
-   return idx,idy
+      if (np.abs(lat2d[idy,idx]-vlat)) > 5: errormess("latitude not found ",printvar=lat2d)
+   if file is not None:
+      print idx,idy,lon2d[idy,idx],vlon
+      print idx,idy,lat2d[idy,idx],vlat
+      var = file.variables["HGT"][:,:,:]
+      mpl.contourf(var[0,:,:],30,cmap = mpl.get_cmap(name="Greys_r") ) ; mpl.axis('off') ; mpl.plot(idx,idy,'bo')
+      mpl.show()
+   return idy,idx
 
 ## Author: TN