07/08/2024, yujin.zeng@nasa.gov

This directory contains the code for running the river routing model introduced in the paper.

The script river_routing_model_CAL_CAL5%_UNCAL.f90 is for running experiments CAL, CAL5% and UNCAL.

The script river_routing_model_LIN.f90 is for running experiment LIN.

The input dir contains all the needed input files for the model.

The runoff forcings (OL7000) can be downloaded through: https://gmao.gsfc.nasa.gov/gmaoftp/yzeng3/river_routing_paper/data/OL7000

To run the model:

1. Set the path of input (var inputdir in the *.f90).
2. Set the path of run off forcing (var forcingdir in the *.f90).
3. Set the path of output (var outputdir in the *.f90).
4. Set the experiment name (CAL, CAL5%, UNCAL, LIN) in the *.f90.
5. Set the start date (var step_start in the *.f90) and end date (var step_end in the *.f90) in days. Day 1 is for Jan-1-1990. Should be within 9,221 (Mar-31-2015) and 11,596 (Sep-30-2021) as the limit of OL7000 runoff forcing data.
6. Build and run the model.

The outputs are in *.txt format with a daily frequency. The date is included in the name of each file. There are three variables:
1. Main river discharge [m3/s] in the Pfaf_Qr_*.txt
2. Main river storage [kg] in the Pfaf_Wr_*.txt
3. Local stream storage [kg] in the Pfaf_Ws_*.txt

Each *.txt file contains a list of 291,284 numbers corresponding to the variable for the catchment indexed from 1 to 291,284. To convert the list to a map, users need to use the catchment distribution map at 1-min: CatchIndex in SRTM_PfafData.nc (https://gmao.gsfc.nasa.gov/gmaoftp/yzeng3/river_routing_paper/SRTM_PfafData.nc). 

Please refer the paper or contact yujin.zeng@nasa.gov for more information.