remove some references to defunct fxns in README, documentation

Author: ehinman

README.md

@@ -6,7 +6,7 @@
 
 ## Latest Announcements
 
-**02/24/2026** The `get_gwlevels`, `get_discharge_measurements` in the `nwis` module are deprecated and have been replaced with the `get_field_measurements` function in the `waterdata` module. The `get_pmcodes` function in the `nwis` module has been replaced with the `get_reference_table(collection='parameter_code)` function.
+**02/24/2026** The `get_gwlevels`, `get_discharge_measurements` functions in the `nwis` module are defunct and have been replaced with the `get_field_measurements` function in the `waterdata` module. The `get_pmcodes` function in the `nwis` module has been replaced with the `get_reference_table(collection='parameter_code)` function. Finally, the `get_water_use` function in the `nwis` module is defunct with no current replacement.
 
 :mega: **01/16/2026:** `dataretrieval` now features the `waterdata` module,
 which provides access to USGS's modernized [Water Data
@@ -138,7 +138,7 @@ To log messages to a file, you can specify a filename in the
 logging.basicConfig(filename='waterdata.log', level=logging.INFO)
 ```
 
-### NWIS Legacy Services (Deprecated but still functional)
+### Legacy NWIS Services (Deprecated but still functional)
 
 The `nwis` module accesses legacy NWIS Water Services:
 

demos/R Python Vignette equivalents.ipynb

@@ -62,9 +62,7 @@
     "\n",
     "# sample data Nitrate:\n",
     "parameterCd = \"00618\"\n",
-    "samples_data, md = waterdata.get_samples(monitoringLocationIdentifier=f\"USGS-{siteNumber}\", usgsPCode=parameterCd, activityStartDateLower=\"1980-01-01\", activityStartDateUpper=\"2010-01-01\")\n",
-    "\n",
-    "pCode, md = nwis.get_pmcodes(parameterCd=parameterCd)"
+    "samples_data, md = waterdata.get_samples(monitoringLocationIdentifier=f\"USGS-{siteNumber}\", usgsPCode=parameterCd, activityStartDateLower=\"1980-01-01\", activityStartDateUpper=\"2010-01-01\")\n"
    ]
   },
   {
@@ -98,25 +96,6 @@
     "dailyDataAvailable, md = nwis.get_dv(sites=siteNumbers, statCd=\"00003\")"
    ]
   },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {
-    "pycharm": {
-     "is_executing": false
-    }
-   },
-   "outputs": [],
-   "source": [
-    "'''\n",
-    "# Using defaults:\n",
-    "parameterCd <- \"00618\" \n",
-    "parameterINFO <- readNWISpCode(parameterCd)\n",
-    "'''\n",
-    "\n",
-    "pCode, md = nwis.get_pmcodes(parameterCd=\"00618\")"
-   ]
-  },
   {
    "cell_type": "code",
    "execution_count": null,
@@ -211,20 +190,6 @@
     "                                      activityStartDateLower=startDate, activityStartDateUpper=endDate)"
    ]
   },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "'''\n",
-    "siteNumber <- \"434400121275801\"\n",
-    "groundWater <- readNWISgwl(siteNumber)\n",
-    "'''\n",
-    "siteNumber = \"434400121275801\"\n",
-    "groundWater, md = nwis.get_gwlevels(sites=siteNumber)"
-   ]
-  },
   {
    "cell_type": "code",
    "execution_count": null,
@@ -252,35 +217,6 @@
     "ratings_data, md = nwis.get_ratings(site='01594440', file_type=\"base\")"
    ]
   },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "'''surfaceData <- readNWISmeas(siteNumber)'''\n",
-    "siteNumber = '01594440'\n",
-    "surface_data, md = nwis.get_discharge_measurements(sites=siteNumber)"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "'''\n",
-    "allegheny <- readNWISuse(stateCd = \"Pennsylvania\",\n",
-    "                         countyCd = \"Allegheny\")\n",
-    "national <- readNWISuse(stateCd = NULL, \n",
-    "                        countyCd = NULL, \n",
-    "                        transform = TRUE)\n",
-    "'''\n",
-    "allegheny, md = nwis.get_water_use(state=\"PA\", counties=\"003\")\n",
-    "\n",
-    "national, md = nwis.get_water_use()"
-   ]
-  },
   {
    "cell_type": "code",
    "execution_count": null,
@@ -295,22 +231,6 @@
     "discharge_stats, md = nwis.get_stats(sites='02319394', parameterCd=\"00060\", statReportType='annual', statTypeCd='all')"
    ]
   },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "# '''\n",
-    "# specificCond <- readWQPqw('WIDNR_WQX-10032762',\n",
-    "#                           'Specific conductance',\n",
-    "#                           '2011-05-01','2011-09-30')\n",
-    "# '''\n",
-    "# specific_cond, md = wqp.get_results(siteid='WIDNR_WQX-10032762',\n",
-    "#                           characteristicName = 'Specific conductance',\n",
-    "#                           startDateLo='2011-05-01', startDateHi='2011-09-30')"
-   ]
-  },
   {
    "cell_type": "code",
    "execution_count": null,