Precip distribution

pcmdi_metrics/precip_distribution/README.md

@@ -0,0 +1,23 @@
+# Precip distribution metrics
+
+Reference: Ahn, M.-S., P. A. Ullrich, P. J. Gleckler, J. Lee,  A. C. Ordonez, A. G. Pendergrass, and C. Jakob, 2022: Framework for Benchmarking Simulated Precipitation Distributions at Regional Scales. Geoscientific Model Development (in prep)
+
+## Driver code:
+- `precip_distribution_driver.py`
+
+## Parameter codes:
+- `param/`
+  - `precip_distribution_params_IMERG.py`
+  - `precip_distribution_params_TRMM.py`
+  - `precip_distribution_params_CMORPH.py`
+  - `precip_distribution_params_GPCP.py`
+  - `precip_distribution_params_PERSIANN.py`
+  - `precip_distribution_params_ERA5.py`
+  - `precip_distribution_params_cmip5.py`
+  - `precip_distribution_params_cmip6.py`
+
+## Run scripts:
+- `scripts_pcmdi/`
+  - `run_obs.bash`
+  - `run_parallel.wait.bash`
+

pcmdi_metrics/precip_distribution/lib/__init__.py

@@ -0,0 +1,19 @@
+from .argparse_functions import AddParserArgument  # noqa
+from .lib_precip_distribution import (  # noqa
+    precip_distribution_frq_amt,
+    precip_distribution_cum,
+    Regrid,
+    getDailyCalendarMonth,
+    CalcBinStructure,
+    MakeDists,
+    CalcRainMetrics,
+    CalcMetricsDomain,
+    CalcMetricsDomain3Clust,
+    CalcMetricsDomainAR6,
+    CalcPscore,
+    CalcP10P90,
+    oneyear,
+    MedDomain,
+    MedDomain3Clust,
+    MedDomainAR6,
+)

pcmdi_metrics/precip_distribution/lib/argparse_functions.py

@@ -0,0 +1,91 @@
+def AddParserArgument(P):
+    P.add_argument("--mip",
+                   type=str,
+                   dest='mip',
+                   default=None,
+                   help="cmip5, cmip6 or other mip")
+    P.add_argument("--exp",
+                   type=str,
+                   dest='exp',
+                   default=None,
+                   help="amip, cmip or others")
+    P.add_argument("--mod",
+                   type=str,
+                   dest='mod',
+                   default=None,
+                   help="model")
+    P.add_argument("--var",
+                   type=str,
+                   dest='var',
+                   default=None,
+                   help="pr or other variable")
+    P.add_argument("--frq",
+                   type=str,
+                   dest='frq',
+                   default=None,
+                   help="day, 3hr or other frequency")
+    P.add_argument("--modpath",
+                   type=str,
+                   dest='modpath',
+                   default=None,
+                   help="data directory path")
+    P.add_argument("--results_dir",
+                   type=str,
+                   dest='results_dir',
+                   default=None,
+                   help="results directory path")
+    P.add_argument("--case_id",
+                   type=str,
+                   dest='case_id',
+                   default=None,
+                   help="case_id with date")
+    P.add_argument("--prd",
+                   type=int,
+                   dest='prd',
+                   nargs='+',
+                   default=None,
+                   help="start- and end-year for analysis (e.g., 1985 2004)")
+    P.add_argument("--fac",
+                   type=str,
+                   dest='fac',
+                   default=None,
+                   help="factor to make unit of [mm/day]")
+    P.add_argument("--res",
+                   type=int,
+                   dest='res',
+                   nargs='+',
+                   default=None,
+                   help="list of target horizontal resolution [degree] for interporation (lon, lat)")
+    P.add_argument("--ref",
+                   type=str,
+                   dest='ref',
+                   default=None,
+                   help="reference data")
+    P.add_argument("--ref_dir",
+                   type=str,
+                   dest='ref_dir',
+                   default=None,
+                   help="reference directory path")
+    P.add_argument("--exp",
+                   type=str,
+                   dest='exp',
+                   default=None,
+                   help="e.g., historical or amip")
+    P.add_argument("--ver",
+                   type=str,
+                   dest='ver',
+                   default=None,
+                   help="version")
+    P.add_argument("--cmec",
+                   dest="cmec",
+                   default=False,
+                   action="store_true",
+                   help="Use to save CMEC format metrics JSON")
+    P.add_argument("--no_cmec",
+                   dest="cmec",
+                   default=False,
+                   action="store_false",
+                   help="Do not save CMEC format metrics JSON")
+    P.set_defaults(cmec=False)
+
+    return P