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

Timestamp:

Nov 30, 2011, 11:59:14 AM (13 years ago)

Author:

aslmd

Message:

GRAPHICS: previous commit by acolmd is actually useful for all idealized mesoscale simulations. corrected pp.py so that vert interp with those files is possible. some purely cosmetic changes. changed typefile='les' to 'mesoideal'

File:

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

trunk/UTIL/PYTHON/myplot.py

@@ -37 +37 @@
 ## Author: AS, AC
 def whatkindfile (nc):
-    if 'controle' in nc.variables:     typefile = 'gcm'
-    elif 'phisinit' in nc.variables:   typefile = 'gcm'
-    elif '9999' in getattr(nc,'START_DATE') : typefile = 'les'
-    elif 'vert' in nc.variables:       typefile = 'mesoapi'
-    elif 'U' in nc.variables:          typefile = 'meso'
-    elif 'HGT_M' in nc.variables:      typefile = 'geo'
+    if 'controle' in nc.variables:             typefile = 'gcm'
+    elif 'phisinit' in nc.variables:           typefile = 'gcm'
+    elif '9999' in getattr(nc,'START_DATE') :  typefile = 'mesoideal'
+    elif 'vert' in nc.variables:               typefile = 'mesoapi'
+    elif 'U' in nc.variables:                  typefile = 'meso'
+    elif 'HGT_M' in nc.variables:              typefile = 'geo'
     #else:                            errormess("whatkindfile: typefile not supported.")
-    else:                            typefile = 'gcm' # for lslin-ed files from gcm
+    else:                                      typefile = 'gcm' # for lslin-ed files from gcm
     return typefile
 
@@ -395 +395 @@
     elif typefile in ['geo']:
         [lon2d,lat2d] = getcoord2d(nc,nlat='XLAT_M',nlon='XLONG_M')
-    elif typefile in ['les']:
+    elif typefile in ['mesoideal']:
         nx=getattr(nc,'WEST-EAST_GRID_DIMENSION')
         ny=getattr(nc,'SOUTH-NORTH_GRID_DIMENSION')
@@ -490 +490 @@
     typefile = whatkindfile(nc)
     ###
-    if typefile is 'mesoapi':    [uchar,vchar] = ['Um','Vm']
-    elif typefile is 'gcm':      [uchar,vchar] = ['u','v']
-    elif typefile in ['meso','les']:     [uchar,vchar] = ['U','V']
-    else:                        [uchar,vchar] = ['not found','not found']
+    if typefile is 'mesoapi':                  [uchar,vchar] = ['Um','Vm']
+    elif typefile is 'gcm':                    [uchar,vchar] = ['u','v']
+    elif typefile in ['meso','mesoideal']:     [uchar,vchar] = ['U','V']
+    else:                                      [uchar,vchar] = ['not found','not found']
     ###
     if typefile in ['meso']:     metwind = False ## geometrical (wrt grid)