From 4b664ce551cbdd948808aca89c3dcf733a7bcbd9 Mon Sep 17 00:00:00 2001
From: Florian Sprung <florian@ecpspc09.epfl.ch>
Date: Tue, 19 Feb 2019 18:54:44 +0100
Subject: [PATCH] changed the number of grid points that is used to save padded
 fields in physical space using hdf5. The modification makes hdf5 consistent
 with netcdf and does not affect reading (backward compatibility)

---
 channelflow/flowfield.cpp | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/channelflow/flowfield.cpp b/channelflow/flowfield.cpp
index bfd1118a..1065ad81 100644
--- a/channelflow/flowfield.cpp
+++ b/channelflow/flowfield.cpp
@@ -2735,8 +2735,8 @@ void FlowField::hdf5Save(const string& filebase) const {
     // transfer to nopadded grid and save that
     int Nx, Nz;
     if (this->padded()) {
-        Nx = (2 * Nx_) / 3;
-        Nz = (2 * Nz_) / 3;
+        Nx = 2 * (Nx_ / 3);
+        Nz = 2 * (Nz_ / 3);
 
     } else {
         Nx = Nx_;