Changeset 2150 in lmdz_wrf

Timestamp:

Sep 20, 2018, 4:46:39 PM (7 years ago)

Author:

lfita

Message:

Adding reading of values

File:

• trunk/tools/snd_DCAO_netCDF.py (modified) (3 diffs)

trunk/tools/snd_DCAO_netCDF.py

@@ -157 +157 @@
     return sttinf
 
-
 # Arguments
 ##
@@ -192 +191 @@
 stations = {}
 
+siglev = False
+tipolev = False
+fijolev = False
+indices = False
+Nsiglev = 0
+Ntipolev = 0
+Nfijolev = 0
+Nindices = 0
+
 # Recovering information from file
 for line in osnd:
@@ -200 +208 @@
         # Got trhought lines with data!
         if len(linevals) != 0:
-        
-            # Starting of record by findinf 'ESTACION:'
+            print linevals
+            print siglev, Nsiglev, tipolev, Ntipolev, fijolev, Nfijolev, indices, Nindices
+
+            # Starting of record by finding 'ESTACION:'
             if linevals[0] == 'ESTACION:':
                 statinf = station_inf(linevals)
-
+                siglev = False
+                tipolev = False
+                fijolev = False
+                indices = False
+
+            # date/time  of record by finding 'UTC'
             elif gen.searchInlist(linevals, 'UTC'):
                 statTinf = station_timeinf(linevals)
 
                 statinf = statinf + statTinf
-                print statinf
-
+
+            # significative levels
+            elif gen.searchInlist(linevals, 'SIGNIFICATIVOS'):
+                siglev = True
+                Nsiglev = 0
+
+            # tipo levels
+            elif gen.searchInlist(linevals, 'TIPO'):
+                siglev = False
+                # Removing headers
+                Nsiglev = Nsiglev - 2
+
+                tipolev = True
+                Ntipolev = 0
+
+            # fijo levels
+            elif gen.searchInlist(linevals, 'FIJOS'):
+                tipolev = False
+                # Removing headers
+                Ntipolev = Ntipolev - 2
+
+                fijolev = True
+                Nfijolev = 0
+
+            # indices
+            elif gen.searchInlist(linevals, 'TERMICO'):
+                fijolev = False
+                # Removing headers
+                Nfijolev = Nfijolev - 2
+
+                indices = True
+                Nindices = 0
+
+            # End record
+            elif linevals[0] == '1' and indices:
+                indices = False
+
+                stvals = statinf + [Nsiglev, Ntipolev, Nfijolev, Nindices]
+                if not stations.has_key(statinf[0]):
+                    stations[statinf[0]] = [stvals]
+                else:
+                    stv = stations[statinf[0]]
+                    stv.append(stvals)
+                    stations[statinf[0]] = stv
+
+                print stvals
+                quit(-1)
+
+            elif tipolev: 
+                Ntipolev = Ntipolev + 1
+
+            elif siglev: 
+                Nsiglev = Nsiglev + 1
+
+            elif fijolev: 
+                Nfijolev = Nfijolev + 1
+
+            elif indices: 
+                Nindices = Nindices + 1
+
+gen.printing_dictionary(stations)
 osnd.close()
-
-            
-            
-