Prepare the maps required by TRIP2
| Attach:ic_elevtn.gif Δ |
Use the script >TRIP2/map/extract.sh which uses the NCO commands to extract flow direction and elevation first; then use >TRIP2/map/calc_rivseq.f to calculate river sequence based on the extracted flow direction map.
Example: Extract the Indochina region from the 1.0 degree global maps
[thanh@shui map]> more extract.sh
# extract indochina region from the Global map
ncks -d lon,93.,109.5 -d lat,8.,34. 1.0_deg_elevtn.nc ic_1.0_deg_elevtn.nc
ncks -d lon,93.,109.5 -d lat,8.,34. 1.0_deg_flwdir.nc ic_1.0_deg_flwdir.nc