Changeset 2e62398 in flex_extract.git for source/python/classes/GribTools.py


Ignore:
Timestamp:
Nov 16, 2018, 8:31:15 PM (5 years ago)
Author:
Anne Philipp <anne.philipp@…>
Branches:
master, ctbto, dev
Children:
403cbf1
Parents:
adb0b63
Message:

moved from grib_api to eccodes

File:
1 edited

Legend:

Unmodified
Added
Removed
  • source/python/classes/GribTools.py

    r274f9ef r2e62398  
    4444# ------------------------------------------------------------------------------
    4545import os
    46 from gribapi import grib_new_from_file, grib_is_defined, grib_get, \
    47      grib_release, grib_set, grib_write, grib_index_read, \
    48      grib_index_new_from_file, grib_index_add_file,  \
    49      grib_index_write
    50 
    51 
    52 # from eccodes import (codes_grib_new_from_file, codes_is_defined, codes_get,
    53                      # codes_release, codes_set, codes_write, codes_index_read,
    54                      # codes_index_new_from_file, codes_index_add_file,
    55                      # codes_index_write)
     46
     47from eccodes import (codes_grib_new_from_file, codes_is_defined, codes_get,
     48                     codes_release, codes_set, codes_write, codes_index_read,
     49                     codes_index_new_from_file, codes_index_add_file,
     50                     codes_index_write)
    5651
    5752# ------------------------------------------------------------------------------
     
    109104
    110105        while 1:
    111             gid = grib_new_from_file(fileid)
     106            gid = codes_new_from_file(fileid)
    112107
    113108            if gid is None:
     
    121116            i = 0
    122117            for wherekey in wherekeynames:
    123                 if not grib_is_defined(gid, wherekey):
     118                if not codes_is_defined(gid, wherekey):
    124119                    raise Exception("where key was not defined")
    125120
    126121                select = (select and (str(wherekeyvalues[i]) ==
    127                                       str(grib_get(gid, wherekey))))
     122                                      str(codes_get(gid, wherekey))))
    128123                i += 1
    129124
     
    131126                llist = []
    132127                for key in keynames:
    133                     llist.extend([str(grib_get(gid, key))])
     128                    llist.extend([str(codes_get(gid, key))])
    134129                return_list.append(llist)
    135130
    136             grib_release(gid)
     131            codes_release(gid)
    137132
    138133        fileid.close()
     
    182177
    183178        while 1:
    184             gid = grib_new_from_file(fin)
     179            gid = codes_new_from_file(fin)
    185180
    186181            if gid is None:
     
    193188            i = 0
    194189            for wherekey in wherekeynames:
    195                 if not grib_is_defined(gid, wherekey):
     190                if not codes_is_defined(gid, wherekey):
    196191                    raise Exception("where Key was not defined")
    197192
    198193                select = (select and (str(wherekeyvalues[i]) ==
    199                                       str(grib_get(gid, wherekey))))
     194                                      str(codes_get(gid, wherekey))))
    200195                i += 1
    201196
     
    203198                i = 0
    204199                for key in keynames:
    205                     grib_set(gid, key, keyvalues[i])
     200                    codes_set(gid, key, keyvalues[i])
    206201                    i += 1
    207202
    208             grib_write(gid, fout)
    209 
    210             grib_release(gid)
     203            codes_write(gid, fout)
     204
     205            codes_release(gid)
    211206
    212207        fin.close()
     
    250245
    251246        while 1:
    252             gid = grib_new_from_file(fin)
     247            gid = codes_new_from_file(fin)
    253248
    254249            if gid is None:
     
    261256            i = 0
    262257            for key in keynames:
    263                 if not grib_is_defined(gid, key):
     258                if not codes_is_defined(gid, key):
    264259                    raise Exception("Key was not defined")
    265260
    266261                if selectWhere:
    267262                    select = (select and (str(keyvalues[i]) ==
    268                                           str(grib_get(gid, key))))
     263                                          str(codes_get(gid, key))))
    269264                else:
    270265                    select = (select and (str(keyvalues[i]) !=
    271                                           str(grib_get(gid, key))))
     266                                          str(codes_get(gid, key))))
    272267                i += 1
    273268
    274269            if select:
    275                 grib_write(gid, fout)
    276 
    277             grib_release(gid)
     270                codes_write(gid, fout)
     271
     272            codes_release(gid)
    278273
    279274        fin.close()
     
    306301
    307302        if os.path.exists(index_file):
    308             iid = grib_index_read(index_file)
     303            iid = codes_index_read(index_file)
    309304            print("Use existing index file: %s " % (index_file))
    310305        else:
     
    312307                print("Inputfile: %s " % (filename))
    313308                if iid is None:
    314                     iid = grib_index_new_from_file(filename, index_keys)
     309                    iid = codes_index_new_from_file(filename, index_keys)
    315310                else:
    316                     grib_index_add_file(iid, filename)
     311                    codes_index_add_file(iid, filename)
    317312
    318313            if iid is not None:
    319                 grib_index_write(iid, index_file)
     314                codes_index_write(iid, index_file)
    320315
    321316        print('... index done')
Note: See TracChangeset for help on using the changeset viewer.
hosted by ZAMG