Changeset c77630a in flex_extract.git for Source/Python/Classes/EcFlexpart.py

Timestamp:

Aug 13, 2019, 12:32:22 PM (5 years ago)

Author:

Anne Philipp <anne.philipp@…>

Branches:

master, ctbto, dev

Children:

0690a6c

Parents:

33a4ccc

Message:

BUGFIX: for python3-eccodes the file openings need to be in binary mode! added binary marker in filemode

File:

• Source/Python/Classes/EcFlexpart.py (modified) (10 diffs)

Source/Python/Classes/EcFlexpart.py

@@ -1094 +1094 @@
 
             print("outputfile = " + fnout)
-            f_handle = open(fnout, 'w')
-            h_handle = open(hnout, 'w')
-            g_handle = open(gnout, 'w')
+            f_handle = open(fnout, 'wb')
+            h_handle = open(hnout, 'wb')
+            g_handle = open(gnout, 'wb')
 
             # read message for message and store relevant data fields, where
@@ -1402 +1402 @@
             # write original time step to flux file as usual
             fluxfile = GribUtil(os.path.join(c.inputdir, fluxfilename))
-            fluxfile.set_keys(tmpfile, filemode='a', strict=True,
+            fluxfile.set_keys(tmpfile, filemode='ab', strict=True,
                               wherekeynames=['paramId'], wherekeyvalues=[142],
                               keynames=['perturbationNumber','date','time','stepRange','values'],
@@ -1408 +1408 @@
                                          date.hour*100, 0, lsp_new_np[inumb,:,it]],
                              )
-            fluxfile.set_keys(tmpfile, filemode='a', strict=True,
+            fluxfile.set_keys(tmpfile, filemode='ab', strict=True,
                               wherekeynames=['paramId'], wherekeyvalues=[143],
                               keynames=['perturbationNumber','date','time','stepRange','values'],
@@ -1416 +1416 @@
 
             # rr for first subgrid point is identified by step = 1
-            fluxfile.set_keys(tmpfile, filemode='a', strict=True,
+            fluxfile.set_keys(tmpfile, filemode='ab', strict=True,
                               wherekeynames=['paramId'], wherekeyvalues=[142],
                               keynames=['perturbationNumber','date','time','stepRange','values'],
@@ -1422 +1422 @@
                                          date.hour*100, '1', lsp_new_np[inumb,:,it+1]]
                               )
-            fluxfile.set_keys(tmpfile, filemode='a', strict=True,
+            fluxfile.set_keys(tmpfile, filemode='ab', strict=True,
                               wherekeynames=['paramId'], wherekeyvalues=[143],
                               keynames=['perturbationNumber','date','time','stepRange','values'],
@@ -1430 +1430 @@
 
             # rr for second subgrid point is identified by step = 2
-            fluxfile.set_keys(tmpfile, filemode='a', strict=True,
+            fluxfile.set_keys(tmpfile, filemode='ab', strict=True,
                               wherekeynames=['paramId'], wherekeyvalues=[142],
                               keynames=['perturbationNumber','date','time','stepRange','values'],
@@ -1436 +1436 @@
                                          date.hour*100, '2', lsp_new_np[inumb,:,it+2]]
                               )
-            fluxfile.set_keys(tmpfile, filemode='a', strict=True,
+            fluxfile.set_keys(tmpfile, filemode='ab', strict=True,
                               wherekeynames=['paramId'], wherekeyvalues=[143],
                               keynames=['perturbationNumber','date','time','stepRange','values'],
@@ -1466 +1466 @@
 
         gribfile.copy_dummy_msg(ifile, keynames=['paramId'],
-                      keyvalues=[142], filemode='w')
+                      keyvalues=[142], filemode='wb')
 
         gribfile.copy_dummy_msg(ifile, keynames=['paramId'],
-                      keyvalues=[143], filemode='a')
+                      keyvalues=[143], filemode='ab')
 
         return
@@ -1573 +1573 @@
                 fortfile = os.path.join(c.inputdir, 'fort.' + k)
                 silent_remove(fortfile)
-                fdict[k] = open(fortfile, 'w')
+                fdict[k] = open(fortfile, 'wb')
 #============================================================================================
             # create correct timestamp from the three time informations
@@ -1608 +1608 @@
             #    if 'olddate' not in locals() or cdate != olddate:
             #        fwrf = open(os.path.join(c.outputdir,
-            #                    'WRF' + cdate + '.' + ctime + '.000.grb2'), 'w')
+            #                    'WRF' + cdate + '.' + ctime + '.000.grb2'), 'wb')
             #        olddate = cdate[:]
 #============================================================================================