diff options
| author | schnetter <> | 2004-02-09 13:53:00 +0000 |
|---|---|---|
| committer | schnetter <> | 2004-02-09 13:53:00 +0000 |
| commit | eecd0bc7c4957139c1d250e089820ff2dbefbcaf (patch) | |
| tree | 5db8dfd3cb2e62e94a4524eaa28ba94a22eccf82 /Carpet/CarpetIOASCII/src | |
| parent | bf275ebfb21841dd531eae41d74b4c76d2848ebd (diff) | |
Allow output for grid variables with just one grid point in certain
Allow output for grid variables with just one grid point in certain
directions. This grid point will get the lower boundary of the
coordinate system.
darcs-hash:20040209135302-07bb3-6c022e1c5e263535fa22fe6b8f368a9e9b1889ec.gz
Diffstat (limited to 'Carpet/CarpetIOASCII/src')
| -rw-r--r-- | Carpet/CarpetIOASCII/src/ioascii.cc | 12 |
1 files changed, 8 insertions, 4 deletions
diff --git a/Carpet/CarpetIOASCII/src/ioascii.cc b/Carpet/CarpetIOASCII/src/ioascii.cc index b27570141..8e4a92e3b 100644 --- a/Carpet/CarpetIOASCII/src/ioascii.cc +++ b/Carpet/CarpetIOASCII/src/ioascii.cc @@ -30,7 +30,7 @@ #include "ioascii.hh" extern "C" { - static const char* rcsid = "$Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetIOASCII/src/ioascii.cc,v 1.57 2004/02/07 16:21:11 schnetter Exp $"; + static const char* rcsid = "$Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetIOASCII/src/ioascii.cc,v 1.58 2004/02/09 14:53:02 schnetter Exp $"; CCTK_FILEVERSION(Carpet_CarpetIOASCII_ioascii_cc); } @@ -912,9 +912,13 @@ namespace CarpetIOASCII { os << " " << coord_time << " "; for (int d=0; d<D; ++d) { assert (gfext.upper()[d] - gfext.lower()[d] >= 0); - os << " " << (coord_lower[d] + (index[d] - gfext.lower()[d]) - * (coord_upper[d] - coord_lower[d]) - / (gfext.upper()[d] - gfext.lower()[d])); + if (gfext.upper()[d] - gfext.lower()[d] == 0) { + os << " " << coord_lower[d]; + } else { + os << " " << (coord_lower[d] + (index[d] - gfext.lower()[d]) + * (coord_upper[d] - coord_lower[d]) + / (gfext.upper()[d] - gfext.lower()[d])); + } } os << " "; switch (CCTK_VarTypeI(vi)) { |