version=pmwiki-2.2.0-beta67 ordered=1 urlencoded=1 agent=Mozilla/5.0 (Windows; U; Windows NT 6.0; en-US; rv:1.8.1.16) Gecko/20080702 Firefox/2.0.0.16 author= charset=UTF-8 csum= ctime=1217988584 host=203.160.1.71 name=TRIP2.ExtractARegionOfInterest rev=1 targets=TRIP2.PrepareTheMapsRequiredByTRIP2,TRIP2.FlowDirectionMap,TRIP2.RiverSequenceMap,TRIP2.ElevationMap,TRIP2.MaskTheUpstreamRegionOfOneGivenPoint text=[[TRIP2/Prepare the maps required by TRIP2]]%0a(:table align='left' width=100%25 cellpadding=5 cellspacing=0:)%0a(:cell:)%0a* [[TRIP2/Flow direction map]]%0a* [[TRIP2/River sequence map]]%0a* [[TRIP2/Elevation map]]%0a* Extract a region of interest%0a* [[TRIP2/Mask the upstream region of one given point]]%0a(:cell:) Attach:ic_elevtn.gif %0a(:tableend:)%0a%0a%0a%0aUse 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.%0a%0aExample: Extract the Indochina region from the 1.0 degree global maps%0a%0a->'''[thanh@shui map]> more extract.sh''' %0a-># extract indochina region from the Global map%0a->ncks -d lon,93.,109.5 -d lat,8.,34. 1.0_deg_elevtn.nc ic_1.0_deg_elevtn.nc%0a->ncks -d lon,93.,109.5 -d lat,8.,34. 1.0_deg_flwdir.nc ic_1.0_deg_flwdir.nc time=1217988584 author:1217988584= diff:1217988584:1217988584:=1,21d0%0a%3c [[TRIP2/Prepare the maps required by TRIP2]]%0a%3c (:table align='left' width=100%25 cellpadding=5 cellspacing=0:)%0a%3c (:cell:)%0a%3c * [[TRIP2/Flow direction map]]%0a%3c * [[TRIP2/River sequence map]]%0a%3c * [[TRIP2/Elevation map]]%0a%3c * Extract a region of interest%0a%3c * [[TRIP2/Mask the upstream region of one given point]]%0a%3c (:cell:) Attach:ic_elevtn.gif %0a%3c (:tableend:)%0a%3c %0a%3c %0a%3c %0a%3c 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.%0a%3c %0a%3c Example: Extract the Indochina region from the 1.0 degree global maps%0a%3c %0a%3c ->'''[thanh@shui map]> more extract.sh''' %0a%3c -># extract indochina region from the Global map%0a%3c ->ncks -d lon,93.,109.5 -d lat,8.,34. 1.0_deg_elevtn.nc ic_1.0_deg_elevtn.nc%0a%3c ->ncks -d lon,93.,109.5 -d lat,8.,34. 1.0_deg_flwdir.nc ic_1.0_deg_flwdir.nc%0a\ No newline at end of file%0a host:1217988584=203.160.1.71