From 3e3f90a67d2cdc576052717f5d494337effccdbe Mon Sep 17 00:00:00 2001 From: Cory Martin Date: Tue, 16 Aug 2022 14:16:47 -0500 Subject: [PATCH] Fix combine obsspace bug --- src/gsi-ncdiag/combine_obsspace.py | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/gsi-ncdiag/combine_obsspace.py b/src/gsi-ncdiag/combine_obsspace.py index e5bb294ef..b743814bd 100755 --- a/src/gsi-ncdiag/combine_obsspace.py +++ b/src/gsi-ncdiag/combine_obsspace.py @@ -218,12 +218,12 @@ def combine_obsspace(FileList, OutFile, GeoDir): if v in obsspace.variables: _var = obsspace.Variable(v) tmpdata = np.array(_var.read_data()) + del _var else: tmpdata = np.full((obsspace.nlocs), iconv.get_default_fill_val(GeoVarTypes[v]), dtype=GeoVarTypes[v]) tmpgeodata.append(tmpdata) tmpgeoidx.append(np.ones_like(tmpdata).astype(int)*int(idx2)) - del _var del obsspace tmpgeodata = np.hstack(tmpgeodata) tmpgeoidx = np.hstack(tmpgeoidx) @@ -244,12 +244,12 @@ def combine_obsspace(FileList, OutFile, GeoDir): if v in obsspace.variables: _var = obsspace.Variable(v) tmpdata = np.array(_var.read_data()) + del _var else: tmpdata = np.full((obsspace.nlocs, nlevs), iconv.get_default_fill_val(GeoVarTypes[v]), dtype=GeoVarTypes[v]) tmpgeodata.append(tmpdata) tmpgeoidx.append(np.ones_like(tmpdata).astype(int)*int(idx2)) - del _var del obsspace tmpgeodata = np.vstack(tmpgeodata) tmpgeoidx = np.vstack(tmpgeoidx) @@ -268,12 +268,12 @@ def combine_obsspace(FileList, OutFile, GeoDir): if v in obsspace.variables: _var = obsspace.Variable(v) tmpdata = np.array(_var.read_data()) + del _var else: tmpdata = np.full((obsspace.nlocs, ninterfaces), iconv.get_default_fill_val(GeoVarTypes[v]), dtype=GeoVarTypes[v]) tmpgeodata.append(tmpdata) tmpgeoidx.append(np.ones_like(tmpdata).astype(int)*int(idx2)) - del _var del obsspace tmpgeodata = np.vstack(tmpgeodata) tmpgeoidx = np.vstack(tmpgeoidx)