From 43a0eab66c8dc719bc9f12a48b68bdd802daafc6 Mon Sep 17 00:00:00 2001 From: John Krasting Date: Wed, 3 May 2017 08:02:51 -0400 Subject: [PATCH] Veritcal levels in TS_drift.py were reversed - Corrected bug where the vertical levels were flipped in the analysis script. --- tools/analysis/TS_drift.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/tools/analysis/TS_drift.py b/tools/analysis/TS_drift.py index 068bd5e080..87561ebca1 100755 --- a/tools/analysis/TS_drift.py +++ b/tools/analysis/TS_drift.py @@ -32,9 +32,9 @@ def main(cmdLineArgs,stream=False): if 'so_xyave' not in rootGroupS.variables: raise Exception('Could not find "so_xyave" files "%s"'%(cmdLineArgs.infile)) if 'zt' in rootGroupT.variables.keys(): - zt = rootGroupT.variables['zt'][::-1] * -1 + zt = rootGroupT.variables['zt'][:] * -1 elif 'z_l' in rootGroupT.variables.keys(): - zt = rootGroupT.variables['z_l'][::-1] * -1 + zt = rootGroupT.variables['z_l'][:] * -1 timeT = rootGroupT.variables['time'] timeS = rootGroupS.variables['time'] timeT = numpy.array([int(x.year) for x in netCDF4.num2date(timeT[:],timeT.units,calendar=timeT.calendar)])