Changeset 2742 in lmdz_wrf for trunk

Timestamp:

Nov 1, 2019, 7:07:41 PM (6 years ago)

Author:

lfita

Message:

Adding:

• `DefineMap?': Function to provide the evnrionment for the projection of the maps

Location:

trunk/tools

Files:

• OBSstations.csv (modified) (1 diff)
• drawing_tools.py (modified) (3 diffs)

trunk/tools/OBSstations.csv

@@ -54 +54 @@
 VILLA REYNOLDS AERO,87448,33.733333,33,44,0.00,65.383333,65,23,0.00,486,SAN LUIS,Argentina,Villa Reynolds Airpot,,
 VILLAGUAY AERO,87385,31.850000,31,51,0.00,59.083333,59,05,0.00,43,ENTRE RIOS,Argentina,Villaguay Airpot,,
+

trunk/tools/drawing_tools.py

@@ -165 +165 @@
 # check_colorBar: Check if the given colorbar exists in matplotlib
 # colorbar_vals: Function to provide the colorbar values for a figure
+# DefineMap: Function to provide the evnrionment for the projection of the maps
 # draw_MerPara: Function to add the parallels and meridians into a figure with a 
 #   panel of plots
@@ -3002 +3003 @@
     return axisv, newaxist, newaxisL
 
+def DefineMap(proj,res,projvals):
+    """ Function to provide the evnrionment for the projection of the maps
+      proj: name of the projection
+        'cyl': cylindric
+        'lcc': lamber-conformal
+        'npstere': north polar stereographic
+        'spstere': south polar stereographic
+      res: resolution to use
+      projvals: values for the given projection
+    """
+    fname = 'DefineMap'
+
+    availproj = ['cyl', 'lcc', 'npstere', 'spstere']
+
+
+    if proj == 'cyl':
+        Nvals = 4
+        if len(projvals) != Nvals:
+            expvals = ['llcrnrlon', 'llcrnrlat', 'urcrnrlon','urcrnrlat']
+            print errrormsg
+            print '  ' + fname + ": projection '" + proj + "' reequires ", Nvals,    \
+              "values, ", len(projvals), 'provided!!'
+            print '    expected values', expvals
+            quit(-1)
+        defmap=Basemap(projection=proj, llcrnrlon=projvals[0], llcrnrlat=projvals[1],\
+         urcrnrlon=projvals[2], urcrnrlat=projvals[3], resolution=res)
+
+    elif proj == 'lcc':
+        Nvals = 6
+        if len(projvals) != Nvals:
+            expvals = ['lat_0', 'lon_0', 'llcrnrlon', 'llcrnrlat', 'urcrnrlon',      \
+              'urcrnrlat']
+            print errrormsg
+            print '  ' + fname + ": projection '" + proj + "' reequires ", Nvals,    \
+              "values, ", len(projvals), 'provided!!'
+            print '    expected values', expvals
+            quit(-1)
+            m = Basemap(projection=proj, lat_0=projvals[0], lon_0=projvals[1],       \
+              llcrnrlon=projvals[2], llcrnrlat=projvals[3], urcrnrlon=projvals[4],   \
+              urcrnrlat=projvals[5], resolution=res)
+
+    elif proj == 'npstere':
+        Nvals = 2
+        if len(projvals) != Nvals:
+            expvals = ['boundinglat', 'lon_0']
+            print errrormsg
+            print '  ' + fname + ": projection '" + proj + "' reequires ", Nvals,    \
+              "values, ", len(projvals), 'provided!!'
+            print '    expected values', expvals
+            quit(-1)
+        defmap = Basemap(projection=proj, boundinglat=projvals[0], lon_0=projvals[1],\
+          resolution=res)
+
+    elif proj == 'spstere':
+        Nvals = 2
+        if len(projvals) != Nvals:
+            expvals = ['boundinglat', 'lon_0']
+            print errrormsg
+            print '  ' + fname + ": projection '" + proj + "' reequires ", Nvals,    \
+              "values, ", len(projvals), 'provided!!'
+            print '    expected values', expvals
+            quit(-1)
+        defmap = Basemap(projection=proj, boundinglat=projvals[0], lon_0=projvals[1],\
+          resolution=res)
+    else:
+        print errormsg
+        print '  ' + fname + ": projection '" + + "' not available !!"
+        print '    avaialble ones:', availproj
+        quit(-1)
+
+    return defmap
+
 def figureMap_size(lonlatb, figwidth=8, Ncol=1, Nrow=1, titpercen=1.1, dpi=100):
     """ Function to determine the real size of a figure adjusted to the ratio of the 
@@ -4322 +4395 @@
 
         if map_proj == 'cyl':
-            m = Basemap(projection=map_proj, llcrnrlon=nlon, llcrnrlat=nlat,         \
-              urcrnrlon=xlon, urcrnrlat= xlat, resolution=map_res)
+            m = DefineMap(map_proj, map_res, [nlon,nlat,xlon,xlat])
         elif map_proj == 'lcc':
-            m = Basemap(projection=map_proj, lat_0=lat2, lon_0=lon2, llcrnrlon=nlon, \
-              llcrnrlat=nlat, urcrnrlon=xlon, urcrnrlat= xlat, resolution=map_res)
-        else:
-            print errormsg
-            print '  ' + fname + ": projection '" + map_proj + "' does not exist!!"
-            print '    existing ones: cyl, lcc'
-            quit(-1)
-
+            m = DefineMap(map_proj, map_res, [lat2, lon2, nlon, nlat, xlon, xlat])
+        elif map_proj == 'npstere':
+            m = DefineMap(map_proj, map_res, [lat2, lon2])
+        elif map_proj == 'spstere':
+            m = DefineMap(map_proj, map_res, [xlat, -65.])
         if len(olon[:].shape) == 1:
             lons, lats = np.meshgrid(olon[:], olat[:])