ERA5 and file checks implemented

Author: rzinke

methods/transient/Transient_Requirement_Validation.ipynb

@@ -499,12 +499,20 @@
    "outputs": [],
    "source": [
     "if 'do_iono' in site_info.keys() and site_info.get('do_iono') != \"False\":\n",
-    "    iono_stack_file = f\"{mintpy_dir}/inputs/ionStack.h5\"             # Input ionosphere stack file  \n",
+    "    # Input ionosphere stack file  \n",
+    "    iono_stack_file = f\"{mintpy_dir}/inputs/ionStack.h5\"\n",
     "\n",
-    "    dirpath, filename = os.path.split(ifgs_file)\n",
-    "    name, ext = os.path.splitext(filename)\n",
-    "    output_ifgs = os.path.join(dirpath, f\"{name}_iono{ext}\")\n",
-    "    print('#'*10, 'Ionosphere Correction set to True', '#'*10)\n",
+    "    print(f\"Ionosphere Correction dataset: {iono_stack_file:s}\")\n",
+    "\n",
+    "    # Check it tropo correction file exists\n",
+    "    if os.path.exists(iono_stack_file):\n",
+    "        dirpath, filename = os.path.split(ifgs_file)\n",
+    "        name, ext = os.path.splitext(filename)\n",
+    "        output_ifgs = os.path.join(dirpath, f\"{name}_iono{ext}\")\n",
+    "        print('#'*10, 'Ionosphere Correction set to True', '#'*10)\n",
+    "\n",
+    "    else:\n",
+    "        site_info['do_iono'] = \"False\"\n",
     "\n",
     "else:\n",
     "    site_info['do_iono'] = \"False\"\n",
@@ -594,22 +602,28 @@
     "    if site_info.get(\"tropo_model\") != \"HRRR\":\n",
     "        # ERA5-based correction\n",
     "        tropo_source = \"ERA5\"\n",
-    "        tropo_cor_file = os.path.join(mintpy_dir, \"inputs\", \"ERA5.h5\")\n",
+    "        tropo_cor_file = os.path.join(mintpy_dir, \"inputs\", f\"{tropo_source}.h5\")\n",
     "\n",
     "    else:\n",
     "        # HRRR-based correction\n",
     "        tropo_source = \"HRRR_ARIA\"\n",
-    "        tropo_cor_file = os.path.join(mintpy_dir, \"inputs\", \"HRRR_ARIA.h5\")\n",
+    "        tropo_cor_file = os.path.join(mintpy_dir, \"inputs\", f\"{tropo_source}.h5\")\n",
     "\n",
     "    print(f\"Troposphere Correction dataset: {tropo_cor_file:s}\")\n",
     "\n",
-    "    dirpath, filename = os.path.split(ifgs_file)\n",
-    "    name, ext = os.path.splitext(filename)\n",
-    "    output_ifgs = os.path.join(dirpath, f\"{name}_{tropo_source}{ext}\")\n",
-    "    print('#'*10, 'Troposphere Correction set to True', '#'*10)\n",
+    "    # Check it tropo correction file exists\n",
+    "    if os.path.exists(tropo_cor_file):\n",
+    "        dirpath, filename = os.path.split(ifgs_file)\n",
+    "        name, ext = os.path.splitext(filename)\n",
+    "        output_ifgs = os.path.join(dirpath, f\"{name}_{tropo_source}{ext}\")\n",
+    "        print('#'*10, 'Troposphere Correction set to True', '#'*10)\n",
+    "    else:\n",
+    "        site_info['do_tropo'] = \"False\"\n",
     "\n",
     "else:\n",
     "    site_info['do_tropo'] = \"False\"\n",
+    "\n",
+    "if site_info['do_tropo'] == \"False\":\n",
     "    print(\"#\" * 10, \"Troposphere Correction set to False\", \"#\" * 10)"
    ]
   },