Raytracing-based Atmospheric Delay Estimation - RADE

Getting set up to use the code

1. Create an account on the Copernicus servers here: https://cds.climate.copernicus.eu/user
2. Confirm your email, etc.
3. Install the public API key and client (https://cds.climate.copernicus.eu/api-how-to): a. Copy the URL and API key from the webpage into a file in your home directory name ~/.cdsapirc b. Install the CDS API using pip: pip install cdsapi.
4. You must accept the license for each product you wish to download: https://cds.climate.copernicus.eu/cdsapp/#!/terms/licence-to-use-copernicus-products
5. See the test_cdsapi.py script for details of the API. You can test that you can connect to the servers by running the test suite (described below).

TODO:

Create a test script which makes the plot and does a sanity check

############################## ######## Installation ######## ##############################

Installing Geometry Module:

1- create a directory somewhre that you want to build the module mkdir build cd build

2- setup the following variables:

setenv GEOMETRY_DIR /home/fattahi/tools/ray_tracing/raytracing_current/Geometry setenv GEOMETRY_LIB_DIR /home/fattahi/tools/ray_tracing/build

append-path PYTHONPATH /home/fattahi/tools/ray_tracing/build append-path LD_LIBRARY_PATH /home/fattahi/tools/ray_tracing/build

3- inside your build directory run the following command:

cmake . $src_dir/Geometry/cpp/ make python3 $src_dir/Geometry/cython/setup.py build_ext -b /home/fattahi/tools/ray_tracing/build

4- Now you should be able to call the geometry module from python e.g.:

import Geo2rdr

5- Look at an example in the test directory

############################## ######### Testing ########### ##############################

To test out the code, cd into the main directory of the repository, which should contain the test/ folder and the Geometry/ folder. Run the following:

*Step 1: Select the test cd test vi/nano/emacs test_scenario.py (whichever editor you like to use)

Then choose one of the scenarios to test run. scenario_0 uses ERA-I and is a single point, scenario_5 uses WRF and is a small region, etc. You can then run

cd .. python -m unittest

and the code will run that scenario in debug mode, which should produce a pdf plot of the weather model at 500 and 15000 meters, along with the weather model in a pickled object called weatherObj.dat. If the test completes successfully there will be a set of files in the test/ folder containing the estimated delay in meters in ENVI-format files. You can run

gdal2isce_xml.py -i mdx.py

to convert the files to ISCE-readable format and view them.