Changeset 2730 in lmdz_wrf

Timestamp:

Oct 25, 2019, 5:00:12 PM (5 years ago)

Author:

lfita

Message:

Adding:

• `list_listKs': Function to provide a new list transforming the values from a given list by a list of type of values
  • 'F-I' value in `typemod' as transform to integer, but passing first by transformation to float

File:

• trunk/tools/generic_tools.py (modified) (4 diffs)

trunk/tools/generic_tools.py

@@ -172 +172 @@
 # list_coincidences: Function to provide the coincidences between two lists
 # list_differences: Function to provide the differences between two lists
+# list_listKs: Function to provide a new list transforming the values from a given list by
+#    a list of type of values
 # list_norepeatcombos: Function to all possible (Num-1)-combinations of a Num values without repetitions
 # lstring_values: Function to provide a new list-string from a string which is a list of word separated by a character if some values are repeated they are not included
@@ -253 +255 @@
 # units_lunits: Function to provide LaTeX equivalences from a given units
 # vals_around: Function to provide the 3x3 values around a given j,i point
+# variables_values: Function to provide values to plot the different variables values from ASCII file
+#   'variables_values.dat'
 # xtrm_nx: Function to provide the extreme (a percentage above min,max) of a series of values
 # WRF_percenlevels: Function to compute eta-levels for WRF following a percentage distribution
@@ -5850 +5854 @@
         if typeval == 'int' or typeval == 'I':
             return int(value)
+        if typeval == 'float-int' or typeval == 'F-I':
+            return int(np.float(value))
         elif typeval == 'long':
             return long(value)
@@ -17215 +17221 @@
     return gamma
 
+def list_listKs(inlist, Klist):
+    """ Function to provide a new list transforming the values from a given list by
+      a list of type of values
+      inlist: input list
+      Klist: list of types of values to consecutively transform each value of the list
+        'S': String
+        'I': Integer
+        'F': Float
+        'D': Double float
+        'B': Bolean
+    >>> list_listKs(['2.2333', 'Hola', '231.21', '231.21'], ['F', 'S', 'I', 'D'])
+    [2.2333, 'Hola', 231, 231.21000000000001]
+    """
+    fname = 'list_listKs'
+
+    Nin = len(inlist)
+    NK = len(Klist)
+    if Nin != NK:
+        print errormsg
+        print '  ' + fname + ': number of values in list ', Nin, 'and ' +            \
+          'number of types to transform to', NK, 'do not coincide !!'
+        print '    passed values _______'
+        if Nin > NK:
+            for iv in range(NK):
+                print iv, ':', inlist[iv], Klist[iv]
+            print 'exceeded inlist', inlist[NK:Nin]
+        else:
+            for iv in range(Nin):
+                print iv, ':', inlist[iv], Klist[iv]
+            print 'exceeded Klist', Klist[Nin:NK]
+        quit(-1)
+
+    newlist = []
+    for iv in range(Nin): newlist.append(typemod(inlist[iv], Klist[iv]))
+
+    return newlist
+
+print list_listKs(['2.2333', 'Hola', '231.21', '231.21'], ['F', 'S', 'F-I', 'D'])
+
+quit()
+
+def stations_values(stcrit,stval):
+    """ Function to provide information from a given station from one of its values 
+      from ASCII file 'OBStations.csv'
+    variables_values([values])
+        [stcrit]: name of the criteria from OBStations.csv
+          one of: 'WMOid', 'station_name', 'longitude', 'latitude', 'lon_deg',         
+          'lon_min', 'lon_sec', 'lat_deg', 'lat_min', 'lat_sec', 'height', 'country', 
+          'prov', 'nice_name', 'add1', 'add2'
+          [WMOid]: WMO id value
+          [station_name]: station raw name 
+          [longitude]: longitude of the station in degrees
+          [latitude]: latitude of the stations in degrees
+          [lon_deg]: longitude degree of the station
+          [lon_min]: longitude minute of the station
+          [lon_sec]: longitude second of the station
+          [lat_deg]: latitude degree of the station
+          [lat_min]: latitude mniute of the station
+          [lat_sec]: latitude second of the station
+          [height]: height of the station
+          [country]: country of the station
+          [prov]: province/state/...
+          [nice_name]: Complete nice name of the station
+          [add1]: additional value 1
+          [add2]: additional value 2
+        [stval]: value of [stcrit] to be matched
+        return: [WMOid], [station_name], [longitude], [latitude], [lon_deg],         
+          [lon_min], [lon_sec], [lat_deg], [lat_min], [lat_sec], [height], [country], 
+          [prov], [nice_name], [add1], [add2]
+    >>> variables_values('WMOid','87582')
+
+    """
+    import subprocess as sub
+    fname='stations_values'
+
+    criteriavail = ['WMOid', 'station_name', 'longitude', 'latitude', 'lon_deg',     \
+          'lon_min', 'lon_sec', 'lat_deg', 'lat_min', 'lat_sec', 'height', 'country',\
+          'prov', 'nice_name', 'add1', 'add2']
+
+    if stcrit == 'h':
+        print fname + '_____________________________________________________________'
+        print stations_values.__doc__
+        quit()
+
+    folder = os.path.dirname(os.path.realpath(__file__))
+
+    infile = folder + '/OBSstations.csv'
+
+    if not os.path.isfile(infile):
+        print errormsg
+        print '  ' + fname + ": File '" + infile + "' does not exist !!"
+        quit(-1)
+
+    ifst = False
+
+    ncf = open(infile, 'r')
+
+    if not searchInlist(criteriavail,stcrit):
+        print errormsg
+        print '  ' + fname + ": criteria '" + stcrit + "' not available !!"
+        print '    available ones:', criteriavail
+        quit(-1)
+    icrit = criteriavail.index(stcrit)
+    stations = {}
+    Nst = 0
+    for line in ncf:
+        if line[0:1] != '#' and len(line) > 10:
+            values = line.replace('\n','').split(',')
+            if len(values) != 16:
+                print errormsg
+                print "problem in station: '", values[0],                            \
+                  ' it should have 16 values and it has', len(values)
+                quit(-1)
+
+            if values[icrit] == stval:                
+                stations[Nst] = [int(values[0]), values[1]]
+
+    print errormsg
+    print '  ' + fname + ": variable '" + varn + "' not defined !!!"
+    ncf.close()
+    quit(-1)
+
+    return 
+
 #quit()