STEP_60_Generate_SWC_PREC_Datasets.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Step 60: Generate SWC PREC Datasets\n",
    "\n",
    "Purpose: combine SWC stations with CREAT3 projected data, historic extreme event data and rainfall \n",
    "distribution datasets.  Partition data by station into individual feature classes modelling the \n",
    "final text file outputs.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Step 60: Generate SWC PREC Datasets\n",
      "  Station Count          : 5144\n",
      "  Interpolated Grid Count: 6332\n",
      "  Historic Station Count : 11165\n",
      "  Rainfall Distribution Coverage Count: 19614\n"
     ]
    }
   ],
   "source": [
    "import arcpy;\n",
    "import os,sys;\n",
    "import datetime,importlib;\n",
    "\n",
    "print(\"Step 60: Generate SWC PREC Datasets\");\n",
    "\n",
    "import swc_resources;\n",
    "importlib.reload(swc_resources);\n",
    "rez = swc_resources.rez();\n",
    "\n",
    "near_max_historic_val  = 1;\n",
    "near_max_historic_unit = 'degree';\n",
    "near_max_rainfall_val  = 1;\n",
    "near_max_rainfall_unit = 'degree';\n",
    "\n",
    "stations       = rez['results'] + os.sep + 'SWC_Station_Universe';\n",
    "stations_cnt   = arcpy.GetCount_management(stations)[0];\n",
    "\n",
    "gridinterp     = rez['results'] + os.sep + 'CRWU_CREAT_Grid_Projections_DATA_INTERP';\n",
    "gridinterp_cnt = arcpy.GetCount_management(gridinterp)[0];\n",
    "\n",
    "hist      = rez['source'] + os.sep + 'CRWU_CREAT_Historic_Climate_Stations';\n",
    "hist_cnt  = arcpy.GetCount_management(hist)[0];\n",
    "\n",
    "rfcov     = rez['results'] + os.sep + 'RainfallDistribution';\n",
    "rfcov_cnt = arcpy.GetCount_management(rfcov)[0];\n",
    "\n",
    "print(\"  Station Count          : \" + str(stations_cnt));\n",
    "print(\"  Interpolated Grid Count: \" + str(gridinterp_cnt));\n",
    "print(\"  Historic Station Count : \" + str(hist_cnt));\n",
    "print(\"  Rainfall Distribution Coverage Count: \" + str(rfcov_cnt));\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    " ### 60.010: Calculate nearest historical station"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "  copying to new work stations\n",
      "  add tracking fields for joined data\n",
      "  calculating the nearest historic station\n",
      "  adding indexes\n",
      "  joining in historic storms information\n",
      "  tracking unmatched stations\n",
      "  Station Count: 5144\n",
      "Wall time: 15.7 s\n"
     ]
    }
   ],
   "source": [
    "%%time\n",
    "\n",
    "input60010  = rez['results'] + os.sep + 'SWC_Station_Universe';\n",
    "output60010 = rez['working']  + os.sep + 'D4EMLite_PREC_Details_WRK2';\n",
    "\n",
    "if arcpy.Exists(output60010):\n",
    "    arcpy.Delete_management(output60010);\n",
    "   \n",
    "print(\"  copying to new work stations\");\n",
    "arcpy.conversion.FeatureClassToFeatureClass(\n",
    "     in_features = input60010\n",
    "    ,out_path    = os.path.dirname(output60010)\n",
    "    ,out_name    = os.path.basename(output60010)\n",
    ");\n",
    "\n",
    "print(\"  add tracking fields for joined data\");\n",
    "arcpy.management.AddFields(\n",
    "     in_table = output60010\n",
    "    ,field_description = [\n",
    "         ['has_creat3_grid_data'     ,'TEXT','Has CREAT3 Grid Data', 1]\n",
    "        ,['is_interp_grid_data'      ,'TEXT','Is Interpolated Grid Data', 1]\n",
    "        ,['has_creat3_histstorm_data','TEXT','Has CREAT3 Historic Storm Data',1]\n",
    "        ,['has_rainfall_data'        ,'TEXT','Has Rainfall Data',1]\n",
    "    ]\n",
    ");\n",
    "\n",
    "print(\"  calculating the nearest historic station\");\n",
    "arcpy.analysis.Near(\n",
    "     in_features   = output60010\n",
    "    ,near_features = hist\n",
    "    ,location      = 'NO_LOCATION'\n",
    "    ,method        = 'PLANAR'\n",
    ");\n",
    "\n",
    "print(\"  adding indexes\");\n",
    "arcpy.management.AddIndex(\n",
    "     in_table        = output60010\n",
    "    ,fields          = 'NEAR_FID'\n",
    "    ,index_name      = 'NEAR_FID_IDX'\n",
    ");\n",
    "\n",
    "print(\"  joining in historic storms information\");\n",
    "flds = [\n",
    "     'IntensePrecip_Historic_5Year_24'\n",
    "    ,'IntensePrecip_Historic_10Year_2'\n",
    "    ,'IntensePrecip_Historic_15Year_2'\n",
    "    ,'IntensePrecip_Historic_30Year_2'\n",
    "    ,'IntensePrecip_Historic_50Year_2'\n",
    "    ,'IntensePrecip_Historic_100Year_'\n",
    "];\n",
    "\n",
    "arcpy.management.JoinField(\n",
    "     in_data         = output60010\n",
    "    ,in_field        = 'NEAR_FID'\n",
    "    ,join_table      = hist\n",
    "    ,join_field      = 'OBJECTID'\n",
    "    ,fields          = flds\n",
    ");\n",
    "\n",
    "flds = [\n",
    "     'IntensePrecip_Historic_5Year_24'\n",
    "    ,'IntensePrecip_Historic_10Year_2'\n",
    "    ,'IntensePrecip_Historic_15Year_2'\n",
    "    ,'IntensePrecip_Historic_30Year_2'\n",
    "    ,'IntensePrecip_Historic_50Year_2'\n",
    "    ,'IntensePrecip_Historic_100Year_'\n",
    "];\n",
    "\n",
    "for fld in flds:\n",
    "\n",
    "    arcpy.management.CalculateField(\n",
    "         in_table        = output60010\n",
    "        ,field           = fld\n",
    "        ,expression      = \"nIfNone(!\" + fld + \"!,!NEAR_DIST!)\"\n",
    "        ,expression_type = 'PYTHON3'\n",
    "        ,code_block      = \"\"\"\n",
    "def nIfNone(pval,pdist):\n",
    "    if pdist > \"\"\" + str(near_max_historic_val) + \"\"\" or pdist == -1:\n",
    "        return 0;\n",
    "    else:\n",
    "        return pval;\n",
    "\"\"\"\n",
    "    );\n",
    "\n",
    "print(\"  tracking unmatched stations\");\n",
    "arcpy.management.CalculateField(\n",
    "     in_table        = output60010\n",
    "    ,field           = 'has_creat3_histstorm_data'\n",
    "    ,expression      = \"nIfNone(!IntensePrecip_Historic_5Year_24!)\"\n",
    "    ,expression_type = 'PYTHON3'\n",
    "    ,code_block      = \"\"\"\n",
    "def nIfNone(pin):\n",
    "    if pin is None or pin <= 0:\n",
    "       return 'N';\n",
    "    else:\n",
    "       return 'Y';\n",
    "\"\"\"\n",
    ");\n",
    "\n",
    "cnt = arcpy.GetCount_management(output60010)[0];\n",
    "print(\"  Station Count: \" + str(cnt));\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 60.020: Calculate the containing CREAT3 grid with temperture and precipitation data for stations"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "  calculating the containing CREAT3 grid with precipitation data\n",
      "  updating tracking fields for grid source\n",
      "  cleaning up unneeded fields\n",
      "  zeroing out unmatched stations\n",
      "  Station Count: 5144\n",
      "Wall time: 3min 15s\n"
     ]
    }
   ],
   "source": [
    "%%time\n",
    "\n",
    "flds = [\n",
    "     'TEMP_2035_HOTDRY_JAN_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_FEB_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_MAR_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_APR_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_MAY_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_JUN_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_JUL_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_AUG_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_SEP_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_OCT_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_NOV_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_DEC_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_ANNUAL_CHANGEF'\n",
    "    \n",
    "    ,'TEMP_2035_CENTRAL_JAN_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_FEB_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_MAR_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_APR_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_MAY_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_JUN_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_JUL_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_AUG_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_SEP_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_OCT_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_NOV_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_DEC_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_ANNUAL_CHANGE'\n",
    "    \n",
    "    ,'TEMP_2035_WETWARM_JAN_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_FEB_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_MAR_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_APR_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_MAY_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_JUN_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_JUL_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_AUG_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_SEP_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_OCT_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_NOV_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_DEC_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_ANNUAL_CHANGE'\n",
    "    \n",
    "    ,'TEMP_2060_HOTDRY_JAN_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_FEB_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_MAR_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_APR_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_MAY_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_JUN_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_JUL_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_AUG_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_SEP_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_OCT_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_NOV_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_DEC_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_ANNUAL_CHANGEF'\n",
    "    \n",
    "    ,'TEMP_2060_CENTRAL_JAN_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_FEB_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_MAR_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_APR_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_MAY_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_JUN_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_JUL_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_AUG_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_SEP_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_OCT_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_NOV_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_DEC_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_ANNUAL_CHANGE'\n",
    "    \n",
    "    ,'TEMP_2060_WETWARM_JAN_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_FEB_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_MAR_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_APR_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_MAY_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_JUN_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_JUL_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_AUG_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_SEP_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_OCT_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_NOV_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_DEC_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_ANNUAL_CHANGE'\n",
    "    \n",
    "    ,'PRECIP_2035_HOTDRY_JAN_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_FEB_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_MAR_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_APR_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_MAY_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_JUN_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_JUL_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_AUG_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_SEP_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_OCT_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_NOV_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_DEC_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_ANNUAL_CHANG'\n",
    "\n",
    "    ,'PRECIP_2035_CENTRAL_JAN_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_FEB_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_MAR_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_APR_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_MAY_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_JUN_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_JUL_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_AUG_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_SEP_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_OCT_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_NOV_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_DEC_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_ANNUAL_CHAN'\n",
    "\n",
    "    ,'PRECIP_2035_WETWARM_JAN_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_FEB_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_MAR_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_APR_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_MAY_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_JUN_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_JUL_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_AUG_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_SEP_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_OCT_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_NOV_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_DEC_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_ANNUAL_CHAN'\n",
    "\n",
    "    ,'PRECIP_2060_HOTDRY_JAN_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_FEB_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_MAR_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_APR_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_MAY_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_JUN_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_JUL_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_AUG_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_SEP_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_OCT_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_NOV_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_DEC_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_ANNUAL_CHANG'\n",
    "\n",
    "    ,'PRECIP_2060_CENTRAL_JAN_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_FEB_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_MAR_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_APR_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_MAY_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_JUN_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_JUL_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_AUG_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_SEP_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_OCT_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_NOV_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_DEC_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_ANNUAL_CHAN'\n",
    "\n",
    "    ,'PRECIP_2060_WETWARM_JAN_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_FEB_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_MAR_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_APR_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_MAY_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_JUN_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_JUL_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_AUG_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_SEP_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_OCT_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_NOV_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_DEC_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_ANNUAL_CHAN'\n",
    "\n",
    "    ,'STORM5YR_2035_STORMY_ANNUAL_CHA'\n",
    "    ,'STORM10YR_2035_STORMY_ANNUAL_CH'\n",
    "    ,'STORM15YR_2035_STORMY_ANNUAL_CH'\n",
    "    ,'STORM30YR_2035_STORMY_ANNUAL_CH'\n",
    "    ,'STORM50YR_2035_STORMY_ANNUAL_CH'\n",
    "    ,'STORM100YR_2035_STORMY_ANNUAL_C'\n",
    "    \n",
    "    ,'STORM5YR_2060_STORMY_ANNUAL_CHA'\n",
    "    ,'STORM10YR_2060_STORMY_ANNUAL_CH'\n",
    "    ,'STORM15YR_2060_STORMY_ANNUAL_CH'\n",
    "    ,'STORM30YR_2060_STORMY_ANNUAL_CH'\n",
    "    ,'STORM50YR_2060_STORMY_ANNUAL_CH'\n",
    "    ,'STORM100YR_2060_STORMY_ANNUAL_C'\n",
    "    \n",
    "    ,'STORM5YR_2035_LESSSTORMY_ANNUAL'\n",
    "    ,'STORM10YR_2035_LESSSTORMY_ANNUA'\n",
    "    ,'STORM15YR_2035_LESSSTORMY_ANNUA'\n",
    "    ,'STORM30YR_2035_LESSSTORMY_ANNUA'\n",
    "    ,'STORM50YR_2035_LESSSTORMY_ANNUA'\n",
    "    ,'STORM100YR_2035_LESSSTORMY_ANNU'\n",
    "    \n",
    "    ,'STORM5YR_2060_LESSSTORMY_ANNUAL'\n",
    "    ,'STORM10YR_2060_LESSSTORMY_ANNUA'\n",
    "    ,'STORM15YR_2060_LESSSTORMY_ANNUA'\n",
    "    ,'STORM30YR_2060_LESSSTORMY_ANNUA'\n",
    "    ,'STORM50YR_2060_LESSSTORMY_ANNUA'\n",
    "    ,'STORM100YR_2060_LESSSTORMY_ANNU'\n",
    "];\n",
    "\n",
    "wrkstat2 = rez['working'] + os.sep + 'D4EMLite_PREC_Details_WRK2';\n",
    "wrkstat3 = rez['working'] + os.sep + 'D4EMLite_PREC_Details_WRK3';\n",
    "\n",
    "if arcpy.Exists(wrkstat3):\n",
    "    arcpy.Delete_management(wrkstat3);\n",
    "    \n",
    "print(\"  calculating the containing CREAT3 grid with precipitation data\");\n",
    "arcpy.analysis.SpatialJoin(\n",
    "     target_features   = wrkstat2\n",
    "    ,join_features     = gridinterp\n",
    "    ,out_feature_class = wrkstat3\n",
    "    ,join_operation    = 'JOIN_ONE_TO_ONE'\n",
    "    ,join_type         = 'KEEP_ALL'\n",
    "    ,match_option      = 'WITHIN'\n",
    ");\n",
    "\n",
    "print(\"  updating tracking fields for grid source\");\n",
    "arcpy.management.CalculateField(\n",
    "     in_table        = wrkstat3\n",
    "    ,field           = 'has_creat3_grid_data'\n",
    "    ,expression      = \"nIfNone(!PRECIP_2035_HOTDRY_JAN_CHANGEPC!)\"\n",
    "    ,expression_type = 'PYTHON3'\n",
    "    ,code_block      = \"\"\"\n",
    "def nIfNone(pin):\n",
    "    if pin is None:\n",
    "       return 'N';\n",
    "    else:\n",
    "       return 'Y';\n",
    "\"\"\"\n",
    ");\n",
    "\n",
    "arcpy.management.CalculateField(\n",
    "     in_table        = wrkstat3\n",
    "    ,field           = 'is_interp_grid_data'\n",
    "    ,expression      = \"!is_interpolated!\"\n",
    "    ,expression_type = 'PYTHON3'\n",
    ");\n",
    "\n",
    "print(\"  cleaning up unneeded fields\")\n",
    "arcpy.management.DeleteField(\n",
    "     in_table   = wrkstat3\n",
    "    ,drop_field = [\n",
    "         'join_count'\n",
    "        ,'TARGET_FID'\n",
    "        ,'CREAT_ID'\n",
    "        ,'GRIDCODE'\n",
    "        ,'HurrGridID'\n",
    "        ,'grid_over_land'\n",
    "        ,'Longitude_Centroid'\n",
    "        ,'Latitude_Centroid'\n",
    "        ,'Longitude_Min'\n",
    "        ,'Latitude_Min'\n",
    "        ,'Longitude_Max'\n",
    "        ,'Latitude_Max'\n",
    "        ,'is_interpolated'\n",
    "    ]\n",
    ");\n",
    "\n",
    "print(\"  zeroing out unmatched stations\");\n",
    "for fld in flds:\n",
    "   \n",
    "    arcpy.management.CalculateField(\n",
    "         in_table        = wrkstat3\n",
    "        ,field           = fld\n",
    "        ,expression      = \"zeroIfNone(!\" + fld + \"!)\"\n",
    "        ,expression_type = 'PYTHON3'\n",
    "        ,code_block      = \"\"\"\n",
    "def zeroIfNone(pin):\n",
    "    if pin is None:\n",
    "       return 0;\n",
    "    else:\n",
    "       return pin;\n",
    "\"\"\"\n",
    "    );\n",
    "    \n",
    "cnt = arcpy.GetCount_management(wrkstat3)[0];\n",
    "print(\"  Station Count: \" + str(cnt));\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 60.030: Determine Precipitation RF value for stations"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "  copying to new work stations\n",
      "  calculating the nearest rainfall coverage polygon within 1 degree(s)\n",
      "  joining in rainfall distribution information\n",
      "  tracking unmatched stations\n",
      "  cleaning up unneeded fields\n",
      "  Station Count: 5144\n",
      "Wall time: 11min 34s\n"
     ]
    }
   ],
   "source": [
    "%%time\n",
    "\n",
    "wrkstat3 = rez['working'] + os.sep + 'D4EMLite_PREC_Details_WRK3';\n",
    "wrkstat4 = rez['working'] + os.sep + 'D4EMLite_PREC_Details_WRK4';\n",
    "\n",
    "if arcpy.Exists(wrkstat4):\n",
    "    arcpy.Delete_management(wrkstat4);\n",
    "    \n",
    "print(\"  copying to new work stations\");\n",
    "arcpy.conversion.FeatureClassToFeatureClass(\n",
    "     in_features = wrkstat3\n",
    "    ,out_path    = os.path.dirname(wrkstat4)\n",
    "    ,out_name    = os.path.basename(wrkstat4)\n",
    ");\n",
    "\n",
    "print(\"  calculating the nearest rainfall coverage polygon \",end='');\n",
    "print(\"within \" + str(near_max_rainfall_val) + \" \" + near_max_rainfall_unit + \"(s)\");\n",
    "arcpy.analysis.Near(\n",
    "     in_features   = wrkstat4\n",
    "    ,near_features = rfcov\n",
    "    ,location      = 'NO_LOCATION'\n",
    "    ,method        = 'PLANAR'\n",
    ");\n",
    "\n",
    "print(\"  joining in rainfall distribution information\");\n",
    "arcpy.management.JoinField(\n",
    "     in_data         = wrkstat4\n",
    "    ,in_field        = 'NEAR_FID'\n",
    "    ,join_table      = rfcov\n",
    "    ,join_field      = 'OBJECTID'\n",
    "    ,fields          = ['rf_value']\n",
    ");\n",
    "\n",
    "arcpy.management.CalculateField(\n",
    "     in_table        = wrkstat4\n",
    "    ,field           = 'rf_value'\n",
    "    ,expression      = \"nIfNone(!rf_value!,!NEAR_DIST!)\"\n",
    "    ,expression_type = 'PYTHON3'\n",
    "    ,code_block      = \"\"\"\n",
    "def nIfNone(pval,pdist):\n",
    "    if pdist > \"\"\" + str(near_max_rainfall_val) + \"\"\" or pdist == -1:\n",
    "        return None;\n",
    "    else:\n",
    "        return pval;\n",
    "\"\"\"\n",
    ");\n",
    "\n",
    "print(\"  tracking unmatched stations\");\n",
    "z = arcpy.management.CalculateField(\n",
    "     in_table        = wrkstat4\n",
    "    ,field           = 'has_rainfall_data'\n",
    "    ,expression      = \"nIfNone(!rf_value!)\"\n",
    "    ,expression_type = 'PYTHON3'\n",
    "    ,code_block      = \"\"\"\n",
    "def nIfNone(pin):\n",
    "    if pin is None:\n",
    "       return 'N';\n",
    "    else:\n",
    "       return 'Y';\n",
    "\"\"\"\n",
    ");\n",
    "\n",
    "print(\"  cleaning up unneeded fields\")\n",
    "z = arcpy.management.DeleteField(\n",
    "     in_table   = wrkstat4\n",
    "    ,drop_field = [\n",
    "         'NEAR_FID'\n",
    "        ,'NEAR_DIST'\n",
    "    ]\n",
    ");\n",
    "\n",
    "cnt = arcpy.GetCount_management(wrkstat4)[0];\n",
    "print(\"  Station Count: \" + str(cnt));\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### TEMP and PRECIP Table Makers\n",
    "\n",
    "Functions to create empty TEMP and PRECIP tables"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "def createTEMP(fgdb,fc):\n",
    "    \n",
    "    if arcpy.Exists(fgdb + os.sep + fc):\n",
    "        arcpy.Delete_management(fgdb + os.sep + fc);\n",
    "    \n",
    "    arcpy.CreateFeatureclass_management(\n",
    "         out_path      = fgdb\n",
    "        ,out_name      = fc\n",
    "        ,geometry_type = \"POINT\"\n",
    "        ,has_m         = \"DISABLED\"\n",
    "        ,has_z         = \"DISABLED\"\n",
    "        ,spatial_reference = arcpy.SpatialReference(4269) \n",
    "    );\n",
    "    \n",
    "    arcpy.management.AddFields(\n",
    "         in_table = fgdb + os.sep + fc\n",
    "        ,field_description = [\n",
    "             ['StationId','TEXT'  ,'StationId',14]\n",
    "            ,['Jan'      ,'DOUBLE','Jan']\n",
    "            ,['Feb'      ,'DOUBLE','Feb']\n",
    "            ,['Mar'      ,'DOUBLE','Mar']\n",
    "            ,['Apr'      ,'DOUBLE','Apr']\n",
    "            ,['May'      ,'DOUBLE','May']\n",
    "            ,['Jun'      ,'DOUBLE','Jun']\n",
    "            ,['Jul'      ,'DOUBLE','Jul']\n",
    "            ,['Aug'      ,'DOUBLE','Aug']\n",
    "            ,['Sep'      ,'DOUBLE','Sep']\n",
    "            ,['Oct'      ,'DOUBLE','Oct']\n",
    "            ,['Nov'      ,'DOUBLE','Nov']\n",
    "            ,['Dec'      ,'DOUBLE','Dec']\n",
    "            ,['Ann'      ,'DOUBLE','Ann']\n",
    "         ]\n",
    "    );\n",
    "    \n",
    "def createPREC(fgdb,fc):\n",
    "    \n",
    "    if arcpy.Exists(fgdb + os.sep + fc):\n",
    "        arcpy.Delete_management(fgdb + os.sep + fc);\n",
    "    \n",
    "    arcpy.CreateFeatureclass_management(\n",
    "         out_path      = fgdb\n",
    "        ,out_name      = fc\n",
    "        ,geometry_type = \"POINT\"\n",
    "        ,has_m         = \"DISABLED\"\n",
    "        ,has_z         = \"DISABLED\"\n",
    "        ,spatial_reference = arcpy.SpatialReference(4269) \n",
    "    );\n",
    "    \n",
    "    arcpy.management.AddFields(\n",
    "         in_table = fgdb + os.sep + fc\n",
    "        ,field_description = [\n",
    "             ['StationId','TEXT'  ,'StationId',14]\n",
    "            ,['Jan'      ,'DOUBLE','Jan']\n",
    "            ,['Feb'      ,'DOUBLE','Feb']\n",
    "            ,['Mar'      ,'DOUBLE','Mar']\n",
    "            ,['Apr'      ,'DOUBLE','Apr']\n",
    "            ,['May'      ,'DOUBLE','May']\n",
    "            ,['Jun'      ,'DOUBLE','Jun']\n",
    "            ,['Jul'      ,'DOUBLE','Jul']\n",
    "            ,['Aug'      ,'DOUBLE','Aug']\n",
    "            ,['Sep'      ,'DOUBLE','Sep']\n",
    "            ,['Oct'      ,'DOUBLE','Oct']\n",
    "            ,['Nov'      ,'DOUBLE','Nov']\n",
    "            ,['Dec'      ,'DOUBLE','Dec']\n",
    "            ,['Ann'      ,'DOUBLE','Ann']\n",
    "         ]\n",
    "    );\n",
    "    \n",
    "def createGEVdepth(fgdb,fc):\n",
    "    \n",
    "    if arcpy.Exists(fgdb + os.sep + fc):\n",
    "        arcpy.Delete_management(fgdb + os.sep + fc);\n",
    "    \n",
    "    arcpy.CreateFeatureclass_management(\n",
    "         out_path      = fgdb\n",
    "        ,out_name      = fc\n",
    "        ,geometry_type = \"POINT\"\n",
    "        ,has_m         = \"DISABLED\"\n",
    "        ,has_z         = \"DISABLED\"\n",
    "        ,spatial_reference = arcpy.SpatialReference(4269) \n",
    "    );\n",
    "    \n",
    "    arcpy.management.AddFields(\n",
    "         in_table = fgdb + os.sep + fc\n",
    "        ,field_description = [\n",
    "             ['StationId'   ,'TEXT'  ,'StationId',14]\n",
    "            ,['z5yrann_mm'  ,'DOUBLE','5yrann_mm']\n",
    "            ,['z10yrann_mm' ,'DOUBLE','10yrann_mm']\n",
    "            ,['z15yrann_mm' ,'DOUBLE','15yrann_mm']\n",
    "            ,['z30yrann_mm' ,'DOUBLE','30yrann_mm']\n",
    "            ,['z50yrann_mm' ,'DOUBLE','50yrann_mm']\n",
    "            ,['z100yrann_mm','DOUBLE','100yrann_mm']\n",
    "         ]\n",
    "    );\n",
    "    "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 60.040: Create fresh tables for TEMP and PREC extracts "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "  creating fresh TEMP feature classes\n",
      "  creating fresh PREC feature classes\n",
      "  creating fresh GEVddepth feature classes\n",
      "  creating fresh PREC_SCS_Types feature class\n",
      "Wall time: 48.8 s\n"
     ]
    }
   ],
   "source": [
    "%%time\n",
    "\n",
    "print(\"  creating fresh TEMP feature classes\");\n",
    "createTEMP(rez['results'],'TEMP2035HotDry');\n",
    "\n",
    "createTEMP(rez['results'],'TEMP2035Central');\n",
    "\n",
    "createTEMP(rez['results'],'TEMP2035WetWarm');\n",
    "\n",
    "createTEMP(rez['results'],'TEMP2060HotDry');\n",
    "\n",
    "createTEMP(rez['results'],'TEMP2060Central');\n",
    "\n",
    "createTEMP(rez['results'],'TEMP2060WetWarm');\n",
    "\n",
    "print(\"  creating fresh PREC feature classes\");\n",
    "createPREC(rez['results'],'PREC2035HotDry');\n",
    "\n",
    "createPREC(rez['results'],'PREC2035Central');\n",
    "\n",
    "createPREC(rez['results'],'PREC2035WetWarm');\n",
    "\n",
    "createPREC(rez['results'],'PREC2060HotDry');\n",
    "\n",
    "createPREC(rez['results'],'PREC2060Central');\n",
    "\n",
    "createPREC(rez['results'],'PREC2060WetWarm');\n",
    "\n",
    "print(\"  creating fresh GEVddepth feature classes\");\n",
    "createGEVdepth(rez['results'],'GEVdepth2035Stormy');\n",
    "\n",
    "createGEVdepth(rez['results'],'GEVdepth2035LessStormy');\n",
    "\n",
    "createGEVdepth(rez['results'],'GEVdepth2060Stormy');\n",
    "\n",
    "createGEVdepth(rez['results'],'GEVdepth2060LessStormy');\n",
    "\n",
    "createGEVdepth(rez['results'],'GEVdepth_historical');\n",
    "\n",
    "print(\"  creating fresh PREC_SCS_Types feature class\");\n",
    "if arcpy.Exists(rez['results'] + os.sep + 'PREC_SCS_Types'):\n",
    "    arcpy.Delete_management(rez['results'] + os.sep + 'PREC_SCS_Types');\n",
    "        \n",
    "arcpy.CreateFeatureclass_management(\n",
    "     out_path      = rez['results']\n",
    "    ,out_name      = 'PREC_SCS_Types'\n",
    "    ,geometry_type = \"POINT\"\n",
    "    ,has_m         = \"DISABLED\"\n",
    "    ,has_z         = \"DISABLED\"\n",
    "    ,spatial_reference = arcpy.SpatialReference(4269) \n",
    ");\n",
    "\n",
    "z = arcpy.management.AddFields(\n",
    "     in_table = rez['results'] + os.sep + 'PREC_SCS_Types'\n",
    "    ,field_description = [\n",
    "         ['StationId'                 ,'TEXT','StationId',14]\n",
    "        ,['Rainfall_Distribution_Type','TEXT','Rainfall_Distribution_Type',255]\n",
    "     ]\n",
    ");\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 60.050: Write out datasets to individual feature classes"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "  truncate datasets for safety\n",
      "  opening an editor session to support multiple cursors\n",
      "  writing the station data to individual feature classes\n",
      "Wall time: 19.7 s\n"
     ]
    }
   ],
   "source": [
    "%%time\n",
    "\n",
    "wrkstat4 = rez['working'] + os.sep + 'D4EMLite_PREC_Details_WRK4';\n",
    "\n",
    "fldin = [\n",
    "     'StationId'\n",
    "    # 1\n",
    "    ,'TEMP_2035_HOTDRY_JAN_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_FEB_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_MAR_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_APR_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_MAY_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_JUN_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_JUL_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_AUG_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_SEP_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_OCT_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_NOV_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_DEC_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_ANNUAL_CHANGEF'\n",
    "    # 14\n",
    "    ,'TEMP_2035_CENTRAL_JAN_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_FEB_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_MAR_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_APR_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_MAY_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_JUN_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_JUL_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_AUG_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_SEP_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_OCT_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_NOV_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_DEC_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_ANNUAL_CHANGE'\n",
    "    # 27\n",
    "    ,'TEMP_2035_WETWARM_JAN_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_FEB_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_MAR_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_APR_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_MAY_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_JUN_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_JUL_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_AUG_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_SEP_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_OCT_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_NOV_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_DEC_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_ANNUAL_CHANGE'\n",
    "    # 40\n",
    "    ,'TEMP_2060_HOTDRY_JAN_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_FEB_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_MAR_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_APR_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_MAY_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_JUN_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_JUL_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_AUG_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_SEP_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_OCT_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_NOV_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_DEC_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_ANNUAL_CHANGEF'\n",
    "    # 53\n",
    "    ,'TEMP_2060_CENTRAL_JAN_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_FEB_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_MAR_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_APR_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_MAY_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_JUN_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_JUL_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_AUG_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_SEP_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_OCT_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_NOV_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_DEC_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_ANNUAL_CHANGE'\n",
    "    # 66\n",
    "    ,'TEMP_2060_WETWARM_JAN_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_FEB_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_MAR_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_APR_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_MAY_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_JUN_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_JUL_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_AUG_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_SEP_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_OCT_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_NOV_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_DEC_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_ANNUAL_CHANGE'\n",
    "    # 79 \n",
    "    ,'PRECIP_2035_HOTDRY_JAN_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_FEB_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_MAR_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_APR_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_MAY_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_JUN_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_JUL_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_AUG_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_SEP_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_OCT_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_NOV_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_DEC_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_ANNUAL_CHANG'\n",
    "    # 92\n",
    "    ,'PRECIP_2035_CENTRAL_JAN_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_FEB_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_MAR_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_APR_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_MAY_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_JUN_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_JUL_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_AUG_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_SEP_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_OCT_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_NOV_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_DEC_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_ANNUAL_CHAN'\n",
    "    # 105\n",
    "    ,'PRECIP_2035_WETWARM_JAN_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_FEB_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_MAR_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_APR_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_MAY_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_JUN_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_JUL_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_AUG_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_SEP_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_OCT_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_NOV_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_DEC_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_ANNUAL_CHAN'\n",
    "    # 118\n",
    "    ,'PRECIP_2060_HOTDRY_JAN_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_FEB_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_MAR_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_APR_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_MAY_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_JUN_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_JUL_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_AUG_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_SEP_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_OCT_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_NOV_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_DEC_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_ANNUAL_CHANG'\n",
    "    # 131\n",
    "    ,'PRECIP_2060_CENTRAL_JAN_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_FEB_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_MAR_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_APR_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_MAY_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_JUN_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_JUL_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_AUG_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_SEP_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_OCT_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_NOV_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_DEC_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_ANNUAL_CHAN'\n",
    "    # 144\n",
    "    ,'PRECIP_2060_WETWARM_JAN_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_FEB_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_MAR_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_APR_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_MAY_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_JUN_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_JUL_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_AUG_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_SEP_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_OCT_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_NOV_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_DEC_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_ANNUAL_CHAN'\n",
    "    # 157\n",
    "    ,'STORM5YR_2035_STORMY_ANNUAL_CHA'\n",
    "    ,'STORM10YR_2035_STORMY_ANNUAL_CH'\n",
    "    ,'STORM15YR_2035_STORMY_ANNUAL_CH'\n",
    "    ,'STORM30YR_2035_STORMY_ANNUAL_CH'\n",
    "    ,'STORM50YR_2035_STORMY_ANNUAL_CH'\n",
    "    ,'STORM100YR_2035_STORMY_ANNUAL_C'\n",
    "    # 163\n",
    "    ,'STORM5YR_2060_STORMY_ANNUAL_CHA'\n",
    "    ,'STORM10YR_2060_STORMY_ANNUAL_CH'\n",
    "    ,'STORM15YR_2060_STORMY_ANNUAL_CH'\n",
    "    ,'STORM30YR_2060_STORMY_ANNUAL_CH'\n",
    "    ,'STORM50YR_2060_STORMY_ANNUAL_CH'\n",
    "    ,'STORM100YR_2060_STORMY_ANNUAL_C'\n",
    "    # 169\n",
    "    ,'STORM5YR_2035_LESSSTORMY_ANNUAL'\n",
    "    ,'STORM10YR_2035_LESSSTORMY_ANNUA'\n",
    "    ,'STORM15YR_2035_LESSSTORMY_ANNUA'\n",
    "    ,'STORM30YR_2035_LESSSTORMY_ANNUA'\n",
    "    ,'STORM50YR_2035_LESSSTORMY_ANNUA'\n",
    "    ,'STORM100YR_2035_LESSSTORMY_ANNU'\n",
    "    # 175\n",
    "    ,'STORM5YR_2060_LESSSTORMY_ANNUAL'\n",
    "    ,'STORM10YR_2060_LESSSTORMY_ANNUA'\n",
    "    ,'STORM15YR_2060_LESSSTORMY_ANNUA'\n",
    "    ,'STORM30YR_2060_LESSSTORMY_ANNUA'\n",
    "    ,'STORM50YR_2060_LESSSTORMY_ANNUA'\n",
    "    ,'STORM100YR_2060_LESSSTORMY_ANNU'\n",
    "    # 181\n",
    "    ,'IntensePrecip_Historic_5Year_24'\n",
    "    ,'IntensePrecip_Historic_10Year_2'\n",
    "    ,'IntensePrecip_Historic_15Year_2'\n",
    "    ,'IntensePrecip_Historic_30Year_2'\n",
    "    ,'IntensePrecip_Historic_50Year_2'\n",
    "    ,'IntensePrecip_Historic_100Year_'\n",
    "    # 187\n",
    "    ,'rf_value'\n",
    "];\n",
    "\n",
    "fldout = [\n",
    "     'StationId'\n",
    "    ,'Jan'\n",
    "    ,'Feb'\n",
    "    ,'Mar'\n",
    "    ,'Apr'\n",
    "    ,'May'\n",
    "    ,'Jun'\n",
    "    ,'Jul'\n",
    "    ,'Aug'\n",
    "    ,'Sep'\n",
    "    ,'Oct'\n",
    "    ,'Nov'\n",
    "    ,'Dec'\n",
    "    ,'Ann'\n",
    "];\n",
    "\n",
    "fldgev = [\n",
    "     'StationId'\n",
    "    ,'z5yrann_mm'\n",
    "    ,'z10yrann_mm'\n",
    "    ,'z15yrann_mm'\n",
    "    ,'z30yrann_mm'\n",
    "    ,'z50yrann_mm'\n",
    "    ,'z100yrann_mm'\n",
    "];\n",
    "\n",
    "fldscs = [\n",
    "     'StationId'\n",
    "    ,'Rainfall_Distribution_Type'\n",
    "];\n",
    "\n",
    "print(\"  truncate datasets for safety\");\n",
    "arcpy.TruncateTable_management(rez['results'] + os.sep + 'TEMP2035HotDry');\n",
    "arcpy.TruncateTable_management(rez['results'] + os.sep + 'TEMP2035Central');\n",
    "arcpy.TruncateTable_management(rez['results'] + os.sep + 'TEMP2035WetWarm');\n",
    "arcpy.TruncateTable_management(rez['results'] + os.sep + 'TEMP2060HotDry');\n",
    "arcpy.TruncateTable_management(rez['results'] + os.sep + 'TEMP2060Central');\n",
    "arcpy.TruncateTable_management(rez['results'] + os.sep + 'TEMP2060WetWarm');\n",
    "\n",
    "arcpy.TruncateTable_management(rez['results'] + os.sep + 'PREC2035HotDry');\n",
    "arcpy.TruncateTable_management(rez['results'] + os.sep + 'PREC2035Central');\n",
    "arcpy.TruncateTable_management(rez['results'] + os.sep + 'PREC2035WetWarm');\n",
    "arcpy.TruncateTable_management(rez['results'] + os.sep + 'PREC2060HotDry');\n",
    "arcpy.TruncateTable_management(rez['results'] + os.sep + 'PREC2060Central');\n",
    "arcpy.TruncateTable_management(rez['results'] + os.sep + 'PREC2060WetWarm');\n",
    "\n",
    "arcpy.TruncateTable_management(rez['results'] + os.sep + 'GEVdepth2035Stormy');\n",
    "arcpy.TruncateTable_management(rez['results'] + os.sep + 'GEVdepth2035LessStormy');\n",
    "arcpy.TruncateTable_management(rez['results'] + os.sep + 'GEVdepth2060Stormy');\n",
    "arcpy.TruncateTable_management(rez['results'] + os.sep + 'GEVdepth2060LessStormy');\n",
    "arcpy.TruncateTable_management(rez['results'] + os.sep + 'GEVdepth_historical');\n",
    "\n",
    "arcpy.TruncateTable_management(rez['results'] + os.sep + 'PREC_SCS_Types');\n",
    "\n",
    "print(\"  opening an editor session to support multiple cursors\");\n",
    "with arcpy.da.Editor(rez['results']) as edit:\n",
    "      \n",
    "    curs_TEMP2035HotDry = arcpy.da.InsertCursor(\n",
    "         in_table    = rez['results'] + os.sep + 'TEMP2035HotDry'   \n",
    "        ,field_names = fldout\n",
    "    );\n",
    "    \n",
    "    curs_TEMP2035Central = arcpy.da.InsertCursor(\n",
    "         in_table    = rez['results'] + os.sep + 'TEMP2035Central'   \n",
    "        ,field_names = fldout\n",
    "    );\n",
    "        \n",
    "    curs_TEMP2035WetWarm = arcpy.da.InsertCursor(\n",
    "         in_table    = rez['results'] + os.sep + 'TEMP2035WetWarm'   \n",
    "        ,field_names = fldout\n",
    "    );\n",
    "    \n",
    "    curs_TEMP2060HotDry = arcpy.da.InsertCursor(\n",
    "         in_table    = rez['results'] + os.sep + 'TEMP2060HotDry'   \n",
    "        ,field_names = fldout\n",
    "    );\n",
    "    \n",
    "    curs_TEMP2060Central = arcpy.da.InsertCursor(\n",
    "         in_table    = rez['results'] + os.sep + 'TEMP2060Central'   \n",
    "        ,field_names = fldout\n",
    "    );\n",
    "        \n",
    "    curs_TEMP2060WetWarm = arcpy.da.InsertCursor(\n",
    "         in_table    = rez['results'] + os.sep + 'TEMP2060WetWarm'   \n",
    "        ,field_names = fldout\n",
    "    );\n",
    "    \n",
    "    curs_PREC2035HotDry = arcpy.da.InsertCursor(\n",
    "         in_table    = rez['results'] + os.sep + 'PREC2035HotDry'   \n",
    "        ,field_names = fldout\n",
    "    );\n",
    "    \n",
    "    curs_PREC2035Central = arcpy.da.InsertCursor(\n",
    "         in_table    = rez['results'] + os.sep + 'PREC2035Central'   \n",
    "        ,field_names = fldout\n",
    "    );\n",
    "        \n",
    "    curs_PREC2035WetWarm = arcpy.da.InsertCursor(\n",
    "         in_table    = rez['results'] + os.sep + 'PREC2035WetWarm'   \n",
    "        ,field_names = fldout\n",
    "    );\n",
    "    \n",
    "    curs_PREC2060HotDry = arcpy.da.InsertCursor(\n",
    "         in_table    = rez['results'] + os.sep + 'PREC2060HotDry'   \n",
    "        ,field_names = fldout\n",
    "    );\n",
    "    \n",
    "    curs_PREC2060Central = arcpy.da.InsertCursor(\n",
    "         in_table    = rez['results'] + os.sep + 'PREC2060Central'   \n",
    "        ,field_names = fldout\n",
    "    );\n",
    "        \n",
    "    curs_PREC2060WetWarm = arcpy.da.InsertCursor(\n",
    "         in_table    = rez['results'] + os.sep + 'PREC2060WetWarm'   \n",
    "        ,field_names = fldout\n",
    "    );\n",
    "    \n",
    "    curs_GEVdepth2035Stormy = arcpy.da.InsertCursor(\n",
    "         in_table    = rez['results'] + os.sep + 'GEVdepth2035Stormy'   \n",
    "        ,field_names = fldgev\n",
    "    );\n",
    "        \n",
    "    curs_GEVdepth2035LessStormy = arcpy.da.InsertCursor(\n",
    "         in_table    = rez['results'] + os.sep + 'GEVdepth2035LessStormy'   \n",
    "        ,field_names = fldgev\n",
    "    );\n",
    "    \n",
    "    curs_GEVdepth2060Stormy = arcpy.da.InsertCursor(\n",
    "         in_table    = rez['results'] + os.sep + 'GEVdepth2060Stormy'   \n",
    "        ,field_names = fldgev\n",
    "    );\n",
    "        \n",
    "    curs_GEVdepth2060LessStormy = arcpy.da.InsertCursor(\n",
    "         in_table    = rez['results'] + os.sep + 'GEVdepth2060LessStormy'   \n",
    "        ,field_names = fldgev\n",
    "    );\n",
    "    \n",
    "    curs_GEVdepth_historical = arcpy.da.InsertCursor(\n",
    "         in_table    = rez['results'] + os.sep + 'GEVdepth_historical'   \n",
    "        ,field_names = fldgev\n",
    "    );\n",
    "    \n",
    "    curs_PREC_SCS_Types = arcpy.da.InsertCursor(\n",
    "         in_table    = rez['results'] + os.sep + 'PREC_SCS_Types'   \n",
    "        ,field_names = fldscs\n",
    "    );\n",
    "    \n",
    "    print(\"  writing the station data to individual feature classes\");\n",
    "    with arcpy.da.SearchCursor(\n",
    "         in_table    = wrkstat4\n",
    "        ,field_names = fldin\n",
    "    ) as incur:\n",
    "    \n",
    "        for row in incur:\n",
    "            \n",
    "            # TEMP 2035 HOTDRY\n",
    "            curs_TEMP2035HotDry.insertRow((\n",
    "                 row[0]\n",
    "                ,round(row[1]  ,2)\n",
    "                ,round(row[2]  ,2)\n",
    "                ,round(row[3]  ,2)\n",
    "                ,round(row[4]  ,2)\n",
    "                ,round(row[5]  ,2)\n",
    "                ,round(row[6]  ,2)\n",
    "                ,round(row[7]  ,2)\n",
    "                ,round(row[8]  ,2)\n",
    "                ,round(row[9]  ,2)\n",
    "                ,round(row[10] ,2)\n",
    "                ,round(row[11] ,2)\n",
    "                ,round(row[12] ,2)\n",
    "                ,round(row[13] ,2)\n",
    "            ));\n",
    "            \n",
    "            # TEMP 2035 CENTRAL\n",
    "            curs_TEMP2035Central.insertRow((\n",
    "                 row[0]\n",
    "                ,round(row[14] ,2)\n",
    "                ,round(row[15] ,2)\n",
    "                ,round(row[16] ,2)\n",
    "                ,round(row[17] ,2)\n",
    "                ,round(row[18] ,2)\n",
    "                ,round(row[19] ,2)\n",
    "                ,round(row[20] ,2)\n",
    "                ,round(row[21] ,2)\n",
    "                ,round(row[22] ,2)\n",
    "                ,round(row[23] ,2)\n",
    "                ,round(row[24] ,2)\n",
    "                ,round(row[25] ,2)\n",
    "                ,round(row[26] ,2)\n",
    "            ));\n",
    "            \n",
    "            # TEMP 2035 WETWARM\n",
    "            curs_TEMP2035WetWarm.insertRow((\n",
    "                 row[0]\n",
    "                ,round(row[27] ,2)\n",
    "                ,round(row[28] ,2)\n",
    "                ,round(row[29] ,2)\n",
    "                ,round(row[30] ,2)\n",
    "                ,round(row[31] ,2)\n",
    "                ,round(row[32] ,2)\n",
    "                ,round(row[33] ,2)\n",
    "                ,round(row[34] ,2)\n",
    "                ,round(row[35] ,2)\n",
    "                ,round(row[36] ,2)\n",
    "                ,round(row[37] ,2)\n",
    "                ,round(row[38] ,2)\n",
    "                ,round(row[39] ,2)\n",
    "            ));\n",
    "            \n",
    "            # TEMP 2060 HOTDRY\n",
    "            curs_TEMP2060HotDry.insertRow((\n",
    "                 row[0]\n",
    "                ,round(row[40] ,2)\n",
    "                ,round(row[41] ,2)\n",
    "                ,round(row[42] ,2)\n",
    "                ,round(row[43] ,2)\n",
    "                ,round(row[44] ,2)\n",
    "                ,round(row[45] ,2)\n",
    "                ,round(row[46] ,2)\n",
    "                ,round(row[47] ,2)\n",
    "                ,round(row[48] ,2)\n",
    "                ,round(row[49] ,2)\n",
    "                ,round(row[50] ,2)\n",
    "                ,round(row[51] ,2)\n",
    "                ,round(row[52] ,2)\n",
    "            ));\n",
    "            \n",
    "            # TEMP 2060 CENTRAL\n",
    "            curs_TEMP2060Central.insertRow((\n",
    "                 row[0]\n",
    "                ,round(row[53] ,2)\n",
    "                ,round(row[54] ,2)\n",
    "                ,round(row[55] ,2)\n",
    "                ,round(row[56] ,2)\n",
    "                ,round(row[57] ,2)\n",
    "                ,round(row[58] ,2)\n",
    "                ,round(row[59] ,2)\n",
    "                ,round(row[60] ,2)\n",
    "                ,round(row[61] ,2)\n",
    "                ,round(row[62] ,2)\n",
    "                ,round(row[63] ,2)\n",
    "                ,round(row[64] ,2)\n",
    "                ,round(row[65] ,2)\n",
    "            ));\n",
    "            \n",
    "            # TEMP 2060 WETWARM\n",
    "            curs_TEMP2060WetWarm.insertRow((\n",
    "                 row[0]\n",
    "                ,round(row[66] ,2)\n",
    "                ,round(row[67] ,2)\n",
    "                ,round(row[68] ,2)\n",
    "                ,round(row[69] ,2)\n",
    "                ,round(row[70] ,2)\n",
    "                ,round(row[71] ,2)\n",
    "                ,round(row[72] ,2)\n",
    "                ,round(row[73] ,2)\n",
    "                ,round(row[74] ,2)\n",
    "                ,round(row[75] ,2)\n",
    "                ,round(row[76] ,2)\n",
    "                ,round(row[77] ,2)\n",
    "                ,round(row[78] ,2)\n",
    "            ));\n",
    "            \n",
    "            # PRECIP 2035 HOTDRY\n",
    "            curs_PREC2035HotDry.insertRow((\n",
    "                 row[0]\n",
    "                ,round(row[79] ,2)\n",
    "                ,round(row[80] ,2)\n",
    "                ,round(row[81] ,2)\n",
    "                ,round(row[82] ,2)\n",
    "                ,round(row[83] ,2)\n",
    "                ,round(row[84] ,2)\n",
    "                ,round(row[85] ,2)\n",
    "                ,round(row[86] ,2)\n",
    "                ,round(row[87] ,2)\n",
    "                ,round(row[88] ,2)\n",
    "                ,round(row[89] ,2)\n",
    "                ,round(row[90] ,2)\n",
    "                ,round(row[91] ,2)\n",
    "            ));\n",
    "            \n",
    "            # PRECIP 2035 CENTRAL\n",
    "            curs_PREC2035Central.insertRow((\n",
    "                 row[0]\n",
    "                ,round(row[92] ,2)\n",
    "                ,round(row[93] ,2)\n",
    "                ,round(row[94] ,2)\n",
    "                ,round(row[95] ,2)\n",
    "                ,round(row[96] ,2)\n",
    "                ,round(row[97] ,2)\n",
    "                ,round(row[98] ,2)\n",
    "                ,round(row[99] ,2)\n",
    "                ,round(row[100],2)\n",
    "                ,round(row[101],2)\n",
    "                ,round(row[102],2)\n",
    "                ,round(row[103],2)\n",
    "                ,round(row[104],2)\n",
    "            ));\n",
    "            \n",
    "            # PRECIP 2035 WETWARM\n",
    "            curs_PREC2035WetWarm.insertRow((\n",
    "                 row[0]\n",
    "                ,round(row[105],2)\n",
    "                ,round(row[106],2)\n",
    "                ,round(row[107],2)\n",
    "                ,round(row[108],2)\n",
    "                ,round(row[109],2)\n",
    "                ,round(row[110],2)\n",
    "                ,round(row[111],2)\n",
    "                ,round(row[112],2)\n",
    "                ,round(row[113],2)\n",
    "                ,round(row[114],2)\n",
    "                ,round(row[115],2)\n",
    "                ,round(row[116],2)\n",
    "                ,round(row[117],2)\n",
    "            ));\n",
    "            \n",
    "            # PRECIP 2060 HOTDRY\n",
    "            curs_PREC2060HotDry.insertRow((\n",
    "                 row[0]\n",
    "                ,round(row[118],2)\n",
    "                ,round(row[119],2)\n",
    "                ,round(row[120],2)\n",
    "                ,round(row[121],2)\n",
    "                ,round(row[122],2)\n",
    "                ,round(row[123],2)\n",
    "                ,round(row[124],2)\n",
    "                ,round(row[125],2)\n",
    "                ,round(row[126],2)\n",
    "                ,round(row[127],2)\n",
    "                ,round(row[128],2)\n",
    "                ,round(row[129],2)\n",
    "                ,round(row[130],2)\n",
    "            ));\n",
    "            \n",
    "            # PRECIP 2060 CENTRAL\n",
    "            curs_PREC2060Central.insertRow((\n",
    "                 row[0]\n",
    "                ,round(row[131],2)\n",
    "                ,round(row[132],2)\n",
    "                ,round(row[133],2)\n",
    "                ,round(row[134],2)\n",
    "                ,round(row[135],2)\n",
    "                ,round(row[136],2)\n",
    "                ,round(row[137],2)\n",
    "                ,round(row[138],2)\n",
    "                ,round(row[139],2)\n",
    "                ,round(row[140],2)\n",
    "                ,round(row[141],2)\n",
    "                ,round(row[142],2)\n",
    "                ,round(row[143],2)\n",
    "            ));\n",
    "            \n",
    "            # PRECIP 2060 WETWARM\n",
    "            curs_PREC2060WetWarm.insertRow((\n",
    "                 row[0]\n",
    "                ,round(row[144],2)\n",
    "                ,round(row[145],2)\n",
    "                ,round(row[146],2)\n",
    "                ,round(row[147],2)\n",
    "                ,round(row[148],2)\n",
    "                ,round(row[149],2)\n",
    "                ,round(row[150],2)\n",
    "                ,round(row[151],2)\n",
    "                ,round(row[152],2)\n",
    "                ,round(row[153],2)\n",
    "                ,round(row[154],2)\n",
    "                ,round(row[155],2)\n",
    "                ,round(row[156],2)\n",
    "            ));\n",
    "            \n",
    "            IntensePrecip_Historic_5Year_mm   = row[181] * 25.4;\n",
    "            IntensePrecip_Historic_10Year_mm  = row[182] * 25.4;\n",
    "            IntensePrecip_Historic_15Year_mm  = row[183] * 25.4;\n",
    "            IntensePrecip_Historic_30Year_mm  = row[184] * 25.4;\n",
    "            IntensePrecip_Historic_50Year_mm  = row[185] * 25.4;\n",
    "            IntensePrecip_Historic_100Year_mm = row[186] * 25.4;\n",
    "            \n",
    "            # STORMY 2035\n",
    "            curs_GEVdepth2035Stormy.insertRow((\n",
    "                 row[0]\n",
    "                ,round(IntensePrecip_Historic_5Year_mm   * (1 + row[157] / 100),2)\n",
    "                ,round(IntensePrecip_Historic_10Year_mm  * (1 + row[158] / 100),2)\n",
    "                ,round(IntensePrecip_Historic_15Year_mm  * (1 + row[159] / 100),2)\n",
    "                ,round(IntensePrecip_Historic_30Year_mm  * (1 + row[160] / 100),2)\n",
    "                ,round(IntensePrecip_Historic_50Year_mm  * (1 + row[161] / 100),2)\n",
    "                ,round(IntensePrecip_Historic_100Year_mm * (1 + row[162] / 100),2)\n",
    "            ));\n",
    "            \n",
    "            # STORMY 2060\n",
    "            curs_GEVdepth2060Stormy.insertRow((\n",
    "                 row[0]\n",
    "                ,round(IntensePrecip_Historic_5Year_mm   * (1 + row[163] / 100),2)\n",
    "                ,round(IntensePrecip_Historic_10Year_mm  * (1 + row[164] / 100),2)\n",
    "                ,round(IntensePrecip_Historic_15Year_mm  * (1 + row[165] / 100),2)\n",
    "                ,round(IntensePrecip_Historic_30Year_mm  * (1 + row[166] / 100),2)\n",
    "                ,round(IntensePrecip_Historic_50Year_mm  * (1 + row[167] / 100),2)\n",
    "                ,round(IntensePrecip_Historic_100Year_mm * (1 + row[168] / 100),2)\n",
    "            ));\n",
    "            \n",
    "            # LESS STORMY 2035\n",
    "            curs_GEVdepth2035LessStormy.insertRow((\n",
    "                 row[0]\n",
    "                ,round(IntensePrecip_Historic_5Year_mm   * (1 + row[169] / 100),2)\n",
    "                ,round(IntensePrecip_Historic_10Year_mm  * (1 + row[170] / 100),2)\n",
    "                ,round(IntensePrecip_Historic_15Year_mm  * (1 + row[171] / 100),2)\n",
    "                ,round(IntensePrecip_Historic_30Year_mm  * (1 + row[172] / 100),2)\n",
    "                ,round(IntensePrecip_Historic_50Year_mm  * (1 + row[173] / 100),2)\n",
    "                ,round(IntensePrecip_Historic_100Year_mm * (1 + row[174] / 100),2)\n",
    "            ));\n",
    "            \n",
    "            # LESS STORMY 2060\n",
    "            curs_GEVdepth2060LessStormy.insertRow((\n",
    "                 row[0]\n",
    "                ,round(IntensePrecip_Historic_5Year_mm   * (1 + row[175] / 100),2)\n",
    "                ,round(IntensePrecip_Historic_10Year_mm  * (1 + row[176] / 100),2)\n",
    "                ,round(IntensePrecip_Historic_15Year_mm  * (1 + row[177] / 100),2)\n",
    "                ,round(IntensePrecip_Historic_30Year_mm  * (1 + row[178] / 100),2)\n",
    "                ,round(IntensePrecip_Historic_50Year_mm  * (1 + row[179] / 100),2)\n",
    "                ,round(IntensePrecip_Historic_100Year_mm * (1 + row[180] / 100),2)\n",
    "            ));\n",
    "            \n",
    "            # Intense Precipitation Historic in mm\n",
    "            curs_GEVdepth_historical.insertRow((\n",
    "                 row[0]\n",
    "                ,round(IntensePrecip_Historic_5Year_mm  ,2)\n",
    "                ,round(IntensePrecip_Historic_10Year_mm ,2)\n",
    "                ,round(IntensePrecip_Historic_15Year_mm ,2)\n",
    "                ,round(IntensePrecip_Historic_30Year_mm ,2)\n",
    "                ,round(IntensePrecip_Historic_50Year_mm ,2)\n",
    "                ,round(IntensePrecip_Historic_100Year_mm,2)\n",
    "            ));\n",
    "            \n",
    "            curs_PREC_SCS_Types.insertRow((\n",
    "                 row[0]\n",
    "                ,row[187]\n",
    "            ));\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 60.060: Finalize Results and cleanup unneeded fields"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "  copying to workstations results\n",
      "  cleaning up the fieldset\n",
      "Wall time: 1min\n"
     ]
    }
   ],
   "source": [
    "%%time\n",
    "\n",
    "wrkstat4 = rez['working'] + os.sep + 'D4EMLite_PREC_Details_WRK4';\n",
    "\n",
    "rezstat = rez['results'] + os.sep + 'D4EMLite_PREC_Details';\n",
    "\n",
    "if arcpy.Exists(rezstat):\n",
    "    arcpy.Delete_management(rezstat);\n",
    "   \n",
    "print(\"  copying to workstations results\");\n",
    "arcpy.conversion.FeatureClassToFeatureClass(\n",
    "     in_features = wrkstat4\n",
    "    ,out_path    = os.path.dirname(rezstat)\n",
    "    ,out_name    = os.path.basename(rezstat)\n",
    ");\n",
    "\n",
    "print(\"  cleaning up the fieldset\");\n",
    "flds_to_keep = [\n",
    "     'StationId'\n",
    "    ,'DataType'\n",
    "    ,'FileName'\n",
    "    ,'Id'\n",
    "    ,'Lat'\n",
    "    ,'Long'\n",
    "    ,'Scenario'\n",
    "    ,'Constituent'\n",
    "    ,'SDate'\n",
    "    ,'EDate'\n",
    "    ,'YrCount'\n",
    "    ,'Value'\n",
    "    ,'StaNam'\n",
    "    ,'has_creat3_grid_data'\n",
    "    ,'is_interp_grid_data'\n",
    "    ,'has_creat3_histstorm_data'\n",
    "    ,'has_rainfall_data'\n",
    "    ,'Shape'\n",
    "    ,'OBJECTID'\n",
    "];\n",
    "\n",
    "flds = [f.name for f in arcpy.ListFields(rezstat)]\n",
    "delete_flds = list(set(flds) - set(flds_to_keep));\n",
    "\n",
    "z = arcpy.management.DeleteField(\n",
    "     in_table   = rezstat\n",
    "    ,drop_field = delete_flds\n",
    ");\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 60.070: Review and QA\n",
    "\n",
    "QA Products:\n",
    "\n",
    "1. comparison of interpolated values to nearest values saved to **step60qa.txt**\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "  writing out QA report\n",
      "  D4EMLite_PREC_Details Count : 5144\n",
      "  TEMP2035HotDry Count        : 5144\n",
      "  TEMP2035Central Count       : 5144\n",
      "  TEMP2035WetWarm Count       : 5144\n",
      "  TEMP2060HotDry Count        : 5144\n",
      "  TEMP2060Central Count       : 5144\n",
      "  TEMP2060WetWarm Count       : 5144\n",
      "  PREC2035HotDry Count        : 5144\n",
      "  PREC2035Central Count       : 5144\n",
      "  PREC2035WetWarm Count       : 5144\n",
      "  PREC2060HotDry Count        : 5144\n",
      "  PREC2060Central Count       : 5144\n",
      "  PREC2060WetWarm Count       : 5144\n",
      "  GEVdepth2035Stormy Count    : 5144\n",
      "  GEVdepth2035LessStormy Count: 5144\n",
      "  GEVdepth2060Stormy Count    : 5144\n",
      "  GEVdepth2060LessStormy Count: 5144\n",
      "  GEVdepth_historical Count   : 5144\n",
      "  PREC_SCS_Types Count        : 5144\n",
      " \n",
      "Wall time: 3.68 s\n"
     ]
    }
   ],
   "source": [
    "%%time\n",
    "\n",
    "wrkstat4 = rez['working'] + os.sep + 'D4EMLite_PREC_Details_WRK4';\n",
    "\n",
    "qa_flds = [\n",
    "     # 0\n",
    "     'StationId'\n",
    "    ,'DataType'\n",
    "    ,'FileName'\n",
    "    ,'Lat'\n",
    "    ,'Long'\n",
    "    ,'Scenario'\n",
    "    ,'Constituent'\n",
    "    ,'SDate'\n",
    "    ,'EDate'\n",
    "    ,'YrCount'\n",
    "    # 10\n",
    "    ,'Value'\n",
    "    ,'StaNam'\n",
    "    ,'has_creat3_grid_data'\n",
    "    ,'is_interp_grid_data'\n",
    "    ,'has_creat3_histstorm_data'\n",
    "    ,'has_rainfall_data'\n",
    "    ,'IntensePrecip_Historic_5Year_24'\n",
    "    ,'IntensePrecip_Historic_10Year_2'\n",
    "    ,'IntensePrecip_Historic_15Year_2'\n",
    "    ,'IntensePrecip_Historic_30Year_2'\n",
    "    # 20\n",
    "    ,'IntensePrecip_Historic_50Year_2'\n",
    "    ,'IntensePrecip_Historic_100Year_'\n",
    "    ,'TEMP_2035_HOTDRY_JAN_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_FEB_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_MAR_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_APR_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_MAY_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_JUN_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_JUL_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_AUG_CHANGEF'\n",
    "    # 30\n",
    "    ,'TEMP_2035_HOTDRY_SEP_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_OCT_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_NOV_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_DEC_CHANGEF'\n",
    "    ,'TEMP_2035_HOTDRY_ANNUAL_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_JAN_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_FEB_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_MAR_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_APR_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_MAY_CHANGEF'\n",
    "    # 40\n",
    "    ,'TEMP_2035_CENTRAL_JUN_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_JUL_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_AUG_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_SEP_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_OCT_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_NOV_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_DEC_CHANGEF'\n",
    "    ,'TEMP_2035_CENTRAL_ANNUAL_CHANGE'\n",
    "    ,'TEMP_2035_WETWARM_JAN_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_FEB_CHANGEF'\n",
    "    # 50\n",
    "    ,'TEMP_2035_WETWARM_MAR_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_APR_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_MAY_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_JUN_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_JUL_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_AUG_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_SEP_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_OCT_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_NOV_CHANGEF'\n",
    "    ,'TEMP_2035_WETWARM_DEC_CHANGEF'\n",
    "    # 60\n",
    "    ,'TEMP_2035_WETWARM_ANNUAL_CHANGE'\n",
    "    ,'TEMP_2060_HOTDRY_JAN_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_FEB_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_MAR_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_APR_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_MAY_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_JUN_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_JUL_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_AUG_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_SEP_CHANGEF'\n",
    "    # 70\n",
    "    ,'TEMP_2060_HOTDRY_OCT_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_NOV_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_DEC_CHANGEF'\n",
    "    ,'TEMP_2060_HOTDRY_ANNUAL_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_JAN_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_FEB_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_MAR_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_APR_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_MAY_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_JUN_CHANGEF'\n",
    "    # 80\n",
    "    ,'TEMP_2060_CENTRAL_JUL_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_AUG_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_SEP_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_OCT_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_NOV_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_DEC_CHANGEF'\n",
    "    ,'TEMP_2060_CENTRAL_ANNUAL_CHANGE'\n",
    "    ,'TEMP_2060_WETWARM_JAN_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_FEB_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_MAR_CHANGEF'\n",
    "    # 90\n",
    "    ,'TEMP_2060_WETWARM_APR_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_MAY_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_JUN_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_JUL_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_AUG_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_SEP_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_OCT_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_NOV_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_DEC_CHANGEF'\n",
    "    ,'TEMP_2060_WETWARM_ANNUAL_CHANGE' \n",
    "    # 100\n",
    "    ,'PRECIP_2035_HOTDRY_JAN_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_FEB_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_MAR_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_APR_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_MAY_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_JUN_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_JUL_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_AUG_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_SEP_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_OCT_CHANGEPC'\n",
    "    # 110\n",
    "    ,'PRECIP_2035_HOTDRY_NOV_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_DEC_CHANGEPC'\n",
    "    ,'PRECIP_2035_HOTDRY_ANNUAL_CHANG'\n",
    "    ,'PRECIP_2035_CENTRAL_JAN_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_FEB_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_MAR_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_APR_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_MAY_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_JUN_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_JUL_CHANGEP'\n",
    "    # 120\n",
    "    ,'PRECIP_2035_CENTRAL_AUG_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_SEP_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_OCT_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_NOV_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_DEC_CHANGEP'\n",
    "    ,'PRECIP_2035_CENTRAL_ANNUAL_CHAN'\n",
    "    ,'PRECIP_2035_WETWARM_JAN_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_FEB_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_MAR_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_APR_CHANGEP'\n",
    "    # 130\n",
    "    ,'PRECIP_2035_WETWARM_MAY_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_JUN_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_JUL_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_AUG_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_SEP_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_OCT_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_NOV_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_DEC_CHANGEP'\n",
    "    ,'PRECIP_2035_WETWARM_ANNUAL_CHAN'\n",
    "    # 150\n",
    "    ,'PRECIP_2060_HOTDRY_JAN_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_FEB_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_MAR_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_APR_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_MAY_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_JUN_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_JUL_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_AUG_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_SEP_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_OCT_CHANGEPC'\n",
    "    # 160\n",
    "    ,'PRECIP_2060_HOTDRY_NOV_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_DEC_CHANGEPC'\n",
    "    ,'PRECIP_2060_HOTDRY_ANNUAL_CHANG'\n",
    "    ,'PRECIP_2060_CENTRAL_JAN_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_FEB_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_MAR_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_APR_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_MAY_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_JUN_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_JUL_CHANGEP'\n",
    "    # 170\n",
    "    ,'PRECIP_2060_CENTRAL_AUG_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_SEP_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_OCT_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_NOV_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_DEC_CHANGEP'\n",
    "    ,'PRECIP_2060_CENTRAL_ANNUAL_CHAN'\n",
    "    ,'PRECIP_2060_WETWARM_JAN_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_FEB_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_MAR_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_APR_CHANGEP'\n",
    "    # 180\n",
    "    ,'PRECIP_2060_WETWARM_MAY_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_JUN_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_JUL_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_AUG_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_SEP_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_OCT_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_NOV_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_DEC_CHANGEP'\n",
    "    ,'PRECIP_2060_WETWARM_ANNUAL_CHAN'\n",
    "    ,'STORM5YR_2035_STORMY_ANNUAL_CHA'\n",
    "    # 190\n",
    "    ,'STORM10YR_2035_STORMY_ANNUAL_CH'\n",
    "    ,'STORM15YR_2035_STORMY_ANNUAL_CH'\n",
    "    ,'STORM30YR_2035_STORMY_ANNUAL_CH'\n",
    "    ,'STORM50YR_2035_STORMY_ANNUAL_CH'\n",
    "    ,'STORM100YR_2035_STORMY_ANNUAL_C'\n",
    "    ,'STORM5YR_2060_STORMY_ANNUAL_CHA'\n",
    "    ,'STORM10YR_2060_STORMY_ANNUAL_CH'\n",
    "    ,'STORM15YR_2060_STORMY_ANNUAL_CH'\n",
    "    ,'STORM30YR_2060_STORMY_ANNUAL_CH'\n",
    "    ,'STORM50YR_2060_STORMY_ANNUAL_CH'\n",
    "    # 200\n",
    "    ,'STORM100YR_2060_STORMY_ANNUAL_C'\n",
    "    ,'STORM5YR_2035_LESSSTORMY_ANNUAL'\n",
    "    ,'STORM10YR_2035_LESSSTORMY_ANNUA'\n",
    "    ,'STORM15YR_2035_LESSSTORMY_ANNUA'\n",
    "    ,'STORM30YR_2035_LESSSTORMY_ANNUA'\n",
    "    ,'STORM50YR_2035_LESSSTORMY_ANNUA'\n",
    "    ,'STORM100YR_2035_LESSSTORMY_ANNU'\n",
    "    ,'STORM5YR_2060_LESSSTORMY_ANNUAL'\n",
    "    ,'STORM10YR_2060_LESSSTORMY_ANNUA'\n",
    "    ,'STORM15YR_2060_LESSSTORMY_ANNUA'\n",
    "    # 210\n",
    "    ,'STORM30YR_2060_LESSSTORMY_ANNUA'\n",
    "    ,'STORM50YR_2060_LESSSTORMY_ANNUA'\n",
    "    ,'STORM100YR_2060_LESSSTORMY_ANNU'\n",
    "    ,'IntensePrecip_Historic_5Year_24'\n",
    "    ,'IntensePrecip_Historic_10Year_2'\n",
    "    ,'IntensePrecip_Historic_15Year_2'\n",
    "    ,'IntensePrecip_Historic_30Year_2'\n",
    "    ,'IntensePrecip_Historic_50Year_2'\n",
    "    ,'IntensePrecip_Historic_100Year_'\n",
    "    # 219\n",
    "    ,'rf_value'\n",
    "]\n",
    "\n",
    "print(\"  writing out QA report\");\n",
    "with open(rez['qa'] + os.sep + 'step60qa.txt',\"w\",encoding=\"utf-8\") as out:\n",
    "    out.write(\"Step 60 QA Station Review\\n\");\n",
    "    out.write(datetime.datetime.now().isoformat() + \"\\n\");\n",
    "    \n",
    "    val = \"\";\n",
    "    for fld in qa_flds:\n",
    "        val += '\"' + fld + '\",';\n",
    "    out.write(val.rstrip(',') + \"\\n\");\n",
    "\n",
    "    with arcpy.da.SearchCursor(\n",
    "         in_table    = wrkstat4\n",
    "        ,field_names = qa_flds\n",
    "        ,sql_clause  = (None,'ORDER BY StationId')\n",
    "    ) as cur:\n",
    "        \n",
    "        for row in cur:\n",
    "            \n",
    "            val = \"\";\n",
    "            for item in row:\n",
    "                if item is None:\n",
    "                    val += \",\";\n",
    "                elif type(item) == int or type(item) == float:\n",
    "                    val += str(item) + \",\";\n",
    "                else:\n",
    "                    val += '\"' + item + '\",';\n",
    "                    \n",
    "            out.write(val.rstrip(',') + \"\\n\");\n",
    "\n",
    "D4EMLite_PREC_Details_cnt  = arcpy.GetCount_management(rez['results'] + os.sep + 'D4EMLite_PREC_Details')[0];\n",
    "\n",
    "TEMP2035HotDry_cnt         = arcpy.GetCount_management(rez['results'] + os.sep + 'TEMP2035HotDry')[0];\n",
    "TEMP2035Central_cnt        = arcpy.GetCount_management(rez['results'] + os.sep + 'TEMP2035Central')[0];\n",
    "TEMP2035WetWarm_cnt        = arcpy.GetCount_management(rez['results'] + os.sep + 'TEMP2035WetWarm')[0];\n",
    "\n",
    "TEMP2060HotDry_cnt         = arcpy.GetCount_management(rez['results'] + os.sep + 'TEMP2060HotDry')[0];\n",
    "TEMP2060Central_cnt        = arcpy.GetCount_management(rez['results'] + os.sep + 'TEMP2060Central')[0];\n",
    "TEMP2060WetWarm_cnt        = arcpy.GetCount_management(rez['results'] + os.sep + 'TEMP2060WetWarm')[0];\n",
    "\n",
    "PREC2035HotDry_cnt         = arcpy.GetCount_management(rez['results'] + os.sep + 'PREC2035HotDry')[0];\n",
    "PREC2035Central_cnt        = arcpy.GetCount_management(rez['results'] + os.sep + 'PREC2035Central')[0];\n",
    "PREC2035WetWarm_cnt        = arcpy.GetCount_management(rez['results'] + os.sep + 'PREC2035WetWarm')[0];\n",
    "\n",
    "PREC2060HotDry_cnt         = arcpy.GetCount_management(rez['results'] + os.sep + 'PREC2060HotDry')[0];\n",
    "PREC2060Central_cnt        = arcpy.GetCount_management(rez['results'] + os.sep + 'PREC2060Central')[0];\n",
    "PREC2060WetWarm_cnt        = arcpy.GetCount_management(rez['results'] + os.sep + 'PREC2060WetWarm')[0];\n",
    "\n",
    "GEVdepth2035Stormy_cnt     = arcpy.GetCount_management(rez['results'] + os.sep + 'GEVdepth2035Stormy')[0];\n",
    "GEVdepth2035LessStormy_cnt = arcpy.GetCount_management(rez['results'] + os.sep + 'GEVdepth2035LessStormy')[0];\n",
    "\n",
    "GEVdepth2060Stormy_cnt     = arcpy.GetCount_management(rez['results'] + os.sep + 'GEVdepth2060Stormy')[0];\n",
    "GEVdepth2060LessStormy_cnt = arcpy.GetCount_management(rez['results'] + os.sep + 'GEVdepth2060LessStormy')[0];\n",
    "GEVdepth_historical_cnt    = arcpy.GetCount_management(rez['results'] + os.sep + 'GEVdepth_historical')[0];\n",
    "\n",
    "PREC_SCS_Types_cnt         = arcpy.GetCount_management(rez['results'] + os.sep + 'PREC_SCS_Types')[0];\n",
    "\n",
    "print(\"  D4EMLite_PREC_Details Count : \" + str(D4EMLite_PREC_Details_cnt));\n",
    "\n",
    "print(\"  TEMP2035HotDry Count        : \" + str(TEMP2035HotDry_cnt));\n",
    "print(\"  TEMP2035Central Count       : \" + str(TEMP2035Central_cnt));\n",
    "print(\"  TEMP2035WetWarm Count       : \" + str(TEMP2035WetWarm_cnt));\n",
    "\n",
    "print(\"  TEMP2060HotDry Count        : \" + str(TEMP2060HotDry_cnt));\n",
    "print(\"  TEMP2060Central Count       : \" + str(TEMP2060Central_cnt));\n",
    "print(\"  TEMP2060WetWarm Count       : \" + str(TEMP2060WetWarm_cnt));\n",
    "\n",
    "print(\"  PREC2035HotDry Count        : \" + str(PREC2035HotDry_cnt));\n",
    "print(\"  PREC2035Central Count       : \" + str(PREC2035Central_cnt));\n",
    "print(\"  PREC2035WetWarm Count       : \" + str(PREC2035WetWarm_cnt));\n",
    "\n",
    "print(\"  PREC2060HotDry Count        : \" + str(PREC2060HotDry_cnt));\n",
    "print(\"  PREC2060Central Count       : \" + str(PREC2060Central_cnt));\n",
    "print(\"  PREC2060WetWarm Count       : \" + str(PREC2060WetWarm_cnt));\n",
    "\n",
    "print(\"  GEVdepth2035Stormy Count    : \" + str(GEVdepth2035Stormy_cnt));\n",
    "print(\"  GEVdepth2035LessStormy Count: \" + str(GEVdepth2035LessStormy_cnt));\n",
    "\n",
    "print(\"  GEVdepth2060Stormy Count    : \" + str(GEVdepth2060Stormy_cnt));\n",
    "print(\"  GEVdepth2060LessStormy Count: \" + str(GEVdepth2060LessStormy_cnt));\n",
    "\n",
    "print(\"  GEVdepth_historical Count   : \" + str(GEVdepth_historical_cnt));\n",
    "\n",
    "print(\"  PREC_SCS_Types Count        : \" + str(PREC_SCS_Types_cnt));\n",
    "print(\" \");\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.10"
  },
  "widgets": {
   "application/vnd.jupyter.widget-state+json": {
    "state": {
     "02244e374ff24848afa9f3d8747acb76": {
      "model_module": "@jupyter-widgets/base",
      "model_module_version": "1.1.0",
      "model_name": "LayoutModel",
      "state": {}
     },
     "0d6adac7a2df46cf9da30cfbaf520741": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "1.4.0",
      "model_name": "IntProgressModel",
      "state": {
       "layout": "IPY_MODEL_19cc71b7413e4f11bce11bea24a1a766",
       "max": 4,
       "style": "IPY_MODEL_34c04f6fabce4b668c698db59672ab88",
       "value": 4
      }
     },
     "19cc71b7413e4f11bce11bea24a1a766": {
      "model_module": "@jupyter-widgets/base",
      "model_module_version": "1.1.0",
      "model_name": "LayoutModel",
      "state": {}
     },
     "31789e4d94c3423c96fb7b514da4cb51": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "1.4.0",
      "model_name": "IntProgressModel",
      "state": {
       "layout": "IPY_MODEL_02244e374ff24848afa9f3d8747acb76",
       "max": 2,
       "style": "IPY_MODEL_68586d734a434d39b2653773a8a01441",
       "value": 2
      }
     },
     "34c04f6fabce4b668c698db59672ab88": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "1.4.0",
      "model_name": "ProgressStyleModel",
      "state": {
       "description_width": ""
      }
     },
     "370bb85d9d464a0bbeb6dd721d1ba113": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "1.4.0",
      "model_name": "IntProgressModel",
      "state": {
       "layout": "IPY_MODEL_b421f3cd575a483a8e6d489da015a584",
       "max": 17,
       "style": "IPY_MODEL_cf7de08ce2d34cb29cf4dfa6fb741d8f",
       "value": 17
      }
     },
     "54bf1e7fed054f9eb356d2cf8c5fb366": {
      "model_module": "@jupyter-widgets/base",
      "model_module_version": "1.1.0",
      "model_name": "LayoutModel",
      "state": {}
     },
     "68586d734a434d39b2653773a8a01441": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "1.4.0",
      "model_name": "ProgressStyleModel",
      "state": {
       "description_width": ""
      }
     },
     "7a67c834eaa64f0ab5d70ccdc6340037": {
      "model_module": "@jupyter-widgets/base",
      "model_module_version": "1.1.0",
      "model_name": "LayoutModel",
      "state": {}
     },
     "879f26f186c6497290c86210675c0bb6": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "1.4.0",
      "model_name": "ProgressStyleModel",
      "state": {
       "description_width": ""
      }
     },
     "8fdbeaa1869a48db9851c12ff597d8ae": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "1.4.0",
      "model_name": "ProgressStyleModel",
      "state": {
       "description_width": ""
      }
     },
     "b421f3cd575a483a8e6d489da015a584": {
      "model_module": "@jupyter-widgets/base",
      "model_module_version": "1.1.0",
      "model_name": "LayoutModel",
      "state": {}
     },
     "bbd0ca77f33c4c8ab1991442d20bfaa1": {
      "model_module": "@jupyter-widgets/base",
      "model_module_version": "1.1.0",
      "model_name": "LayoutModel",
      "state": {}
     },
     "c0f6e207d0194812b5842b0bb5df66c3": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "1.4.0",
      "model_name": "IntProgressModel",
      "state": {
       "layout": "IPY_MODEL_7a67c834eaa64f0ab5d70ccdc6340037",
       "max": 2,
       "style": "IPY_MODEL_e8b7c58ad376466d9ef5b51f4fa781b2",
       "value": 2
      }
     },
     "c4f9df5fa8fb4cb3ac4d41011cacbc70": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "1.4.0",
      "model_name": "IntProgressModel",
      "state": {
       "layout": "IPY_MODEL_54bf1e7fed054f9eb356d2cf8c5fb366",
       "max": 8,
       "style": "IPY_MODEL_879f26f186c6497290c86210675c0bb6",
       "value": 8
      }
     },
     "cdb5dd825db84ffc8e08a5762c7d0f3f": {
      "model_module": "@jupyter-widgets/base",
      "model_module_version": "1.1.0",
      "model_name": "LayoutModel",
      "state": {}
     },
     "cf7de08ce2d34cb29cf4dfa6fb741d8f": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "1.4.0",
      "model_name": "ProgressStyleModel",
      "state": {
       "description_width": ""
      }
     },
     "e3b72ef4f19e4d39bccc92f51a460d9c": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "1.4.0",
      "model_name": "IntProgressModel",
      "state": {
       "layout": "IPY_MODEL_bbd0ca77f33c4c8ab1991442d20bfaa1",
       "max": 4,
       "style": "IPY_MODEL_8fdbeaa1869a48db9851c12ff597d8ae",
       "value": 4
      }
     },
     "e8b7c58ad376466d9ef5b51f4fa781b2": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "1.4.0",
      "model_name": "ProgressStyleModel",
      "state": {
       "description_width": ""
      }
     },
     "f908859b1b0745d5a00664f2051af612": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "1.4.0",
      "model_name": "IntProgressModel",
      "state": {
       "layout": "IPY_MODEL_cdb5dd825db84ffc8e08a5762c7d0f3f",
       "max": 8,
       "style": "IPY_MODEL_fde806ce7f52424bba3888f93658d794",
       "value": 8
      }
     },
     "fde806ce7f52424bba3888f93658d794": {
      "model_module": "@jupyter-widgets/controls",
      "model_module_version": "1.4.0",
      "model_name": "ProgressStyleModel",
      "state": {
       "description_width": ""
      }
     }
    },
    "version_major": 2,
    "version_minor": 0
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}