Changeset 2015 in lmdz_wrf

Timestamp:

Aug 2, 2018, 5:56:36 PM (6 years ago)

Author:

lfita

Message:

Adding 'lonlatxtrms' into `draw_stations_map' in order to use lon,lat stations extremes as lonlatbox

File:

• trunk/tools/drawing.py (modified) (6 diffs)

trunk/tools/drawing.py

@@ -10195 +10195 @@
             * 'f', full
         SWNEbox: ',' list of the vertexs of the map [SWlon, NElon, SWlat, NElat]
+          'lonlatxtrms': to use the extremes of the lon,lat stations
         imgtit: title of the image ('!' for spaces)
         imgkind: kind of file output of the image (ps, pns, pdf, ...)
@@ -10211 +10212 @@
     stypevals = values.split(':')[0].split('@')
     mapvals = values.split(':')[1]
-    SWNEbox = gen.str_list_k(values.split(':')[2], ',', 'F')
+    SWNEboxS = values.split(':')[2]
     imgtit = values.split(':')[3].replace('!', ' ')
     imgkind = values.split(':')[4]
@@ -10236 +10237 @@
     ofile = open(filename, 'r')
     dicstations = {}
+    nlonst = 100000.
+    xlonst = -100000.
+    nlatst = 100000.
+    xlatst = -100000.
+
     for line in ofile:
         if line[0:1] != '#' and len(line) > 1: 
@@ -10244 +10250 @@
             llat = np.float(linev[3])
             lhgt = np.float(linev[4])
+            if nlonst > llon: nlonst=llon
+            if xlonst < llon: xlonst=llon
+            if nlatst > llat: nlatst=llat
+            if xlatst < llat: xlatst=llat
 
             if not gen.searchInlist(sttyps, ltyp):
@@ -10257 +10267 @@
     ofile.close()
 
+    if SWNEboxS == 'lonlatxtrms':
+        SWNEbox = [nlonst, xlonst, nlatst, xlatst]
+    else:
+        SWNEbox = gen.str_list_k(SWNEboxS, ',', 'F')
+
     Nd = 100
     nlon = SWNEbox[0]
@@ -10265 +10280 @@
     ddlon = (xlon - nlon)/Nd
     ddlat = (xlat - nlat)/Nd
+
+    if SWNEboxS == 'lonlatxtrms':
+        SWNEbox = [nlonst-ddlon, xlonst+ddlon, nlatst-ddlat, xlatst+ddlat]
+        nlon = SWNEbox[0]
+        xlon = SWNEbox[1]
+        nlat = SWNEbox[2]
+        xlat = SWNEbox[3]
 
     lon1D = np.arange(nlon, xlon+ddlon, ddlon)