Skip to content

Commit

Permalink
adds 1.4 name changes in mass_balance script, fixes documentation bug
Browse files Browse the repository at this point in the history
  • Loading branch information
ecoon committed Mar 10, 2023
1 parent 9036198 commit a5b9864
Show file tree
Hide file tree
Showing 2 changed files with 34 additions and 9 deletions.
Original file line number Diff line number Diff line change
Expand Up @@ -21,7 +21,7 @@ converge.
.. _compressible-porosity-evaluator-spec
.. admonition:: compressible-porosity-evaluator-spec
* `"compressible porosity model parameters`" ``[compressible-porosity-model-spec-list]``
* `"compressible porosity model parameters`" ``[compressible-porosity-standard-model-spec-list]``
KEYS:
Expand Down
41 changes: 33 additions & 8 deletions tools/utils/mass_balance.py
Original file line number Diff line number Diff line change
Expand Up @@ -91,7 +91,27 @@
import ats_xdmf # from $ATS_SRC_DIR/tools/utils/


names_dev = {"ponded depth":"surface-ponded_depth.cell.0",
names_dev = {"ponded depth":"surface-ponded_depth",
"surface pressure":"surface-pressure",
"pressure":"pressure",
"rain":"surface-precipitation_rain",
"snow":"surface-precipitation_snow",
"water content":"water_content",
"surface water content":"surface-water_content",
"saturation":"saturation_liquid",
"saturation gas":"saturation_gas",
"saturation ice":"saturation_ice",
"cell volume":"cell_volume",
"surface cell volume":"surface-cell_volume",
"surface mass density":"surface-mass_density_liquid",
"exfiltration":"surface-surface_subsurface_flux",
"evaporation":"surface-evaporative_flux",
"transpiration":"transpiration",
"surface density": "surface-molar_density_liquid",
"density": "molar_density_liquid",
}

names_13 = {"ponded depth":"surface-ponded_depth.cell.0",
"surface pressure":"surface-pressure.cell.0",
"pressure":"pressure.cell.0",
"rain":"surface-precipitation_rain.cell.0",
Expand Down Expand Up @@ -137,14 +157,20 @@ class MassBalanceFromVis(object):
_names = names_dev

@classmethod
def set_names_dev(cls):
def set_names(cls, namestr):
"""Change this class to use dev variable names.
Usage:
MassBalanceFromVis.set_names_dev()
sim = MassBalanceFromVis(...)
"""
cls._names = names_dev
if namestr == '1.3':
cls._names = names_13
elif namestr == '0.86':
cls._names = names_086
else:
cls._names = names_dev


@classmethod
def set_names_086(cls):
Expand Down Expand Up @@ -491,8 +517,8 @@ def legend(axs):
parser = argparse.ArgumentParser(description="Plot mass balance of a set of runs.")
parser.add_argument('directories', metavar='dirs', type=str, nargs='+',
help="list directories to plot")
parser.add_argument('--dev', action="store_true", default=False,
help="use ats-dev variable names")
parser.add_argument('--ats-version', type=str, default='dev',
help="use variable names from a given ats version")
parser.add_argument('--names', type=str, default=None,
help="string containing a pythonic dictionary of variable names")
parser.add_argument('--symbol', '-s', default=None, type=str,
Expand All @@ -513,13 +539,12 @@ def legend(axs):
else:
names = None

if args.dev:
MassBalanceFromVis.set_names_dev()
MassBalanceFromVis.set_names(args.ats_version)

# loop and plot
for d in args.directories:
sim = MassBalanceFromVis(d, names=names)
color = colors.next()
color = next(colors)
plot(sim, axs, color, symbol=args.symbol, label=d)

decorate(axs)
Expand Down

0 comments on commit a5b9864

Please sign in to comment.