Changeset 566 in lmdz_wrf for trunk/tools/TS_ASCII_netCDF.py

Timestamp:

Jul 6, 2015, 2:38:18 PM (9 years ago)

Author:

lfita

Message:

Adding the use of `ts_header'

File:

• trunk/tools/TS_ASCII_netCDF.py (modified) (2 diffs)

trunk/tools/TS_ASCII_netCDF.py

@@ -13 +13 @@
 
 fillValue = 1.e20
+
+def ts_header(ln):
+    """ Function to get the values of the header of the *.TS files
+      line=ASCII lines with the header of the TS file
+      getting the line format from WRFV3.3 'EMCORE' in file 'share/wrf_timeseries.F'
+    """
+    fname = 'ts_header'
+
+    fmt=['A26', 'I2', 'I3', 'A6', 'A2', 'F7.3', 'A1', 'F8.3', 'A3', 'I4', 'A1', 'I4',\
+      'A3', 'F7.3', 'A1', 'F8.3', 'A2', 'F6.1', 'A7']
+
+    ncvar.values_fortran_fmt(ln,fmt)
+
+    return headervalues
 
 ####### ###### ##### #### ### ## #
@@ -87 +101 @@
     if iline == 0:
         newvar = objofile.createVariable('station','c')
-        valuespar=line.split('(')
+        valueshead = ts_header(line)
 
-        sec = ncvar.reduce_spaces(valuespar[0])
-        Nsec = len(sec)
-        ncvar.set_attribute(newvar, 'name', ncvar.numVector_String(sec[0:Nsec-3],' '))
-        ncvar.set_attribute(newvar, 'acronym',sec[Nsec-1])
+        ncvar.set_attribute(newvar, 'name', valueshead[0])
+        ncvar.set_attribute(newvar, 'acronym',valueshead[3].replace(' ',''))
 
-        sec = ncvar.reduce_spaces(valuespar[1])
-        Nsec = len(sec)
-        ncvar.set_attribute(newvar, 'real_lon',                                      \
-          np.float(sec[1].replace(',','').replace('(','').replace(')','')) )
-        ncvar.set_attribute(newvar, 'real_lat',                                      \
-          np.float(sec[0].replace(',','').replace('(','').replace(')','')) )
+        ncvar.set_attribute(newvar, 'real_lon', valueshead[5])
+        ncvar.set_attribute(newvar, 'real_lat', valueshead[7])
 
-        sec = ncvar.reduce_spaces(valuespar[2])
-        Nsec = len(sec)
-        ncvar.set_attribute(newvar, 'x_grid_point',                                  \
-          int(sec[0].replace(',','').replace('(','').replace(')','')) )
-        ncvar.set_attribute(newvar, 'y_grid_point',                                  \
-          int(sec[1].replace(',','').replace('(','').replace(')','')) )
+        ncvar.set_attribute(newvar, 'x_grid_point', valueshead[9])
+        ncvar.set_attribute(newvar, 'y_grid_point', valueshead[11])
 
-        sec = ncvar.reduce_spaces(valuespar[3])
-        Nsec = len(sec)
-        ncvar.set_attribute(newvar, 'model_lon',                                     \
-          np.float(sec[1].replace(',','').replace('(','').replace(')','')) )
-        ncvar.set_attribute(newvar, 'model_lat',                                     \
-          np.float(sec[0].replace(',','').replace('(','').replace(')','')) )
-        ncvar.set_attribute(newvar, 'model_height',                                  \
-          np.float(sec[2].replace(',','').replace('(','').replace(')','')) )
+        ncvar.set_attribute(newvar, 'model_lon', valueshead[13])
+        ncvar.set_attribute(newvar, 'model_lat', valueshead[15])
+        ncvar.set_attribute(newvar, 'model_height', valueshead[17])
         simstarttime = refdate
     else: