Changeset 1689 in lmdz_wrf

Timestamp:

Dec 5, 2017, 1:20:45 PM (7 years ago)

Author:

lfita

Message:

Adding 'None' for each of the possible files with the list of stations

File:

• trunk/tools/get_stations.py (modified) (1 diff)

trunk/tools/get_stations.py

@@ -280 +280 @@
 
 # surface stations 
-osfcstatf = open(sfcstatfile, 'r')
-
-sfclonvals = []
-sfclatvals = []
-sfclabvals = []
-sfchgtvals = []
-
-for line in osfcstatf:
-    if line[0:1] != comchar and len(line) > 1:
-        vals = line.replace('\n','').split(sepchar)
-        if not gen.searchInlist(missvalS, vals[sfcloncol]) and                       \
-          not gen.searchInlist(missvalS, vals[sfclatcol]):
-            stlon = np.float(vals[sfcloncol])
-            stlat = np.float(vals[sfclatcol])
-            if stlon >= nlon and stlon <= xlon and stlat >= nlat and stlat <= xlat:
-                sfclonvals.append(stlon)
-                sfclatvals.append(stlat)
-                sfclabvals.append(vals[sfclabcol])
-                if gen.searchInlist(missvalS, vals[sfchgtcol]):
-                    if vals[sfchgtcol] == '0': 
+if stcstatfile != 'None':
+    osfcstatf = open(sfcstatfile, 'r')
+
+    sfclonvals = []
+    sfclatvals = []
+    sfclabvals = []
+    sfchgtvals = []
+
+    for line in osfcstatf:
+        if line[0:1] != comchar and len(line) > 1:
+            vals = line.replace('\n','').split(sepchar)
+            if not gen.searchInlist(missvalS, vals[sfcloncol]) and                       \
+              not gen.searchInlist(missvalS, vals[sfclatcol]):
+                stlon = np.float(vals[sfcloncol])
+                stlat = np.float(vals[sfclatcol])
+                if stlon >= nlon and stlon <= xlon and stlat >= nlat and stlat <= xlat:
+                    sfclonvals.append(stlon)
+                    sfclatvals.append(stlat)
+                    sfclabvals.append(vals[sfclabcol])
+                    if gen.searchInlist(missvalS, vals[sfchgtcol]):
+                        if vals[sfchgtcol] == '0': 
+                            sfchgtvals.append(np.float(vals[sfchgtcol]))
+                        else: 
+                            sfchgtvals.append(gen.fillValueF)
+                    else:
                         sfchgtvals.append(np.float(vals[sfchgtcol]))
-                    else: 
-                        sfchgtvals.append(gen.fillValueF)
-                else:
-                    sfchgtvals.append(np.float(vals[sfchgtcol]))
-
-osfcstatf.close()
-
-Nsfcstats = len(sfclonvals)
-print '  Number of surface stations: ', Nsfcstats
+    
+    osfcstatf.close()
+    Nsfcstats = len(sfclonvals)
+    print '  Number of surface stations: ', Nsfcstats
+else:
+    Nsfcstats = 0
+    print '  No surface stations !!'
 
 # sounding stations 
-osndstatf = open(sndstatfile, 'r')
-
-sndlonvals = []
-sndlatvals = []
-sndlabvals = []
-sndhgtvals = []
-
-for line in osndstatf:
-    if line[0:1] != comchar and len(line) > 1:
-        vals = line.replace('\n','').split(sepchar)
-        if vals[sndloncol] != missvalS and vals[sndlatcol] != missvalS:
-            stlon = np.float(vals[sndloncol])
-            stlat = np.float(vals[sndlatcol])
-            if stlon >= nlon and stlon <= xlon and stlat >= nlat and stlat <= xlat:
-                sndlonvals.append(stlon)
-                sndlatvals.append(stlat)
-                sndlabvals.append(vals[sndlabcol])
-                sndhgtvals.append(np.float(vals[sndhgtcol]))
-
-osndstatf.close()
-Nsndstats = len(sndlonvals)
-print '  Number of sounding stations: ', Nsndstats
+if stcstatfile != 'None':
+    osndstatf = open(sndstatfile, 'r')
+
+    sndlonvals = []
+    sndlatvals = []
+    sndlabvals = []
+    sndhgtvals = []
+
+    for line in osndstatf:
+        if line[0:1] != comchar and len(line) > 1:
+            vals = line.replace('\n','').split(sepchar)
+            if vals[sndloncol] != missvalS and vals[sndlatcol] != missvalS:
+                stlon = np.float(vals[sndloncol])
+                stlat = np.float(vals[sndlatcol])
+                if stlon >= nlon and stlon <= xlon and stlat >= nlat and stlat <= xlat:
+                    sndlonvals.append(stlon)
+                    sndlatvals.append(stlat)
+                    sndlabvals.append(vals[sndlabcol])
+                    sndhgtvals.append(np.float(vals[sndhgtcol]))
+
+    osndstatf.close()
+    Nsndstats = len(sndlonvals)
+    print '  Number of sounding stations: ', Nsndstats
+else:
+    Nsndstats = 0
+    print '  No sounding stations !!'
+
 
 if not os.path.isfile(simfilen):