Skip to content

Commit 89c9122

Browse files
updates for cdo > 2.1 (#217)
* updates for cdo > 2.1 * Update whats_new.rst
1 parent 91cdda2 commit 89c9122

File tree

2 files changed

+11
-8
lines changed

2 files changed

+11
-8
lines changed

docs/whats_new.rst

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -14,6 +14,7 @@ UNRELEASED
1414
Internal Changes
1515
~~~~~~~~~~~~~~~~
1616

17+
- Updates for ``cdo > 2.1`` (:pull:`217`).
1718
- Removed ``nblink`` dependency (:pull:`196`).
1819

1920
v0.7.1 (2 August 2024)

pyremo/preproc/era5.py

Lines changed: 10 additions & 8 deletions
Original file line numberDiff line numberDiff line change
@@ -189,7 +189,7 @@ def _get_files(self, date):
189189
)
190190

191191
def _seldate(self, filename, date):
192-
return f"--seldate,{date} {filename}"
192+
return f"-seldate,{date} {filename}"
193193
# return self.cdo.seldate(date, input=filename)
194194

195195
def _seldates(self, filenames, date):
@@ -252,11 +252,11 @@ def _to_regular(self, filename, gridtype=None, setname="", table="ecmwf"):
252252
gridtype = self._gridtype(filename)
253253
# options = f"-f nc4 -t {table}"
254254
if setname:
255-
setname = f"--setname,{setname}" # {filename}"
255+
setname = f"-setname,{setname}" # {filename}"
256256
if gridtype == "gaussian_reduced":
257-
gaussian = "--setgridtype,regular"
257+
gaussian = "-setgridtype,regular"
258258
elif gridtype == "spectral":
259-
gaussian = "--sp2gpl"
259+
gaussian = "-sp2gpl"
260260
elif gridtype == "gaussian":
261261
gaussian = ""
262262
else:
@@ -268,7 +268,7 @@ def _to_regular(self, filename, gridtype=None, setname="", table="ecmwf"):
268268

269269
def _compute_wind(self, vort, div):
270270
"""compute wind from vorticity and divergence"""
271-
return f"--chname,u,ua,v,va --dv2uvl --merge {vort} {div}"
271+
return f"-chname,u,ua,v,va -dv2uvl -merge [ {vort} {div} ]"
272272

273273
def gfile(self, date, path=None, expid=None, filename=None):
274274
"""Create an ERA5 gfile dataset.
@@ -314,10 +314,12 @@ def gfile(self, date, path=None, expid=None, filename=None):
314314
print("computing wind...")
315315
wind = self._compute_wind(seldates["svo"], seldates["sd"])
316316

317-
merge = f"--setgrid,{self.gridfile} --merge " + " ".join(
318-
list(regulars.values()) + [wind]
317+
merge = (
318+
f"-setgrid,{self.gridfile} -merge [ "
319+
+ " ".join(list(regulars.values()) + [wind])
320+
+ " ]"
319321
)
320-
call = f"cdo {self.options} invertlev --invertlat {merge} {filename}"
322+
call = f"cdo {self.options} invertlev -invertlat {merge} {filename}"
321323
print(f"execute: {call}")
322324

323325
subprocess.run(

0 commit comments

Comments
 (0)