r.watershed el=panama_srtm3_90m th=100000 accum=testaccumsrtm3_90m basin=testbasinsrtm3_90m r.mapcalc "teststreams_srtm3_90=if(abs(testaccumsrtm3_90m)>100,1,null())" r.watershed el=panama_srtm3_90m th=100000 accum=testaccumsrtm3_90ms basin=testbasinsrtm3_90ms r.mapcalc "teststreams_srtm3_90s=if(abs(testaccumsrtm3_90ms)>100,1,null())" r.watershed el=ifs_n10w080d th=100000 accum=testaccumifs_10m basin=testbasinifs_10m r.mapcalc "teststreams_ifs_10m=if(abs(testaccumifs_10m)>500,1,null())" r.thin teststreams_ifs_10m out=teststreams_ifs_10mt r.to.vect teststreams_ifs_10mt out=teststreams_ifs_10mt r.watershed el=ifs_n10w080d th=100000 accum=testaccumifs_30mnn basin=testbasinifs_30mnn r.mapcalc "teststreams_ifs_30mnn=if(abs(testaccumifs_30mnn)>100,1,null())" r.watershed el=panama_srtm3_90m th=100000 accum=testaccumsrtm3_30mnn basin=testbasinsrtm3_30mnn r.mapcalc "teststreams_srtm3_30mnn=if(abs(testaccumsrtm3_30mnn)>100,1,null())" r.terraflow -s ifsall_10m_patch filled=ifsall_10m_fill dir=ifsall_10m_dir swater=ifsall_10m_swatersh accum=ifsall_10m_accum tci=ifsall_10m_tci r.terraflow -s srtm90.8378 filled=srtm90.fill dir=srtm90.dir swater=srtm90.swatersh accum=srtm90.accum tci=srtm90.tci r.terraflow [-sq] elev=string filled=string direction=string swatershed=string accumulation=string tci=string [d8cut=value] [memory=value] [STREAM_DIR=string] [stats=string] r.mapcalc "srtmstream90maj=if(abs(srtm90.accum)>5000 && srtm90.8378 > 0,1,null())" r.mapcalc "srtmstream90med=if(abs(srtm90.accum)>1000 && srtm90.8378 > 0,1,null())" r.mapcalc "srtmstream90med2=if(abs(srtm90.accum)>2000 && srtm90.8378 > 0,1,null())" r.mapcalc "srtmstream90=if(abs(srtm90.accum)>100 && srtm90.8378 > 0,1,null())" r.mapcalc "srtmstream90med3=if(abs(accum90ctr) > 1000 && srtm90ctr > 0, 1, null())" r.mapcalc="if(bbasin.50m30K == 34 || bbasin.50m30K == 36 || bbasin.50m30K == 38 || bbasin.50m30K == 32, 1, null())" r.mapcalc "pacorabasin_90m=if(basinsctr90mask==122 || basinsctr90mask==120 || basinsctr90mask==114 || basinsctr90mask==118, 1, null())" r.mapcalc "pacorabasino_90m=if(basinsctr90mask==122 || basinsctr90mask==120 || basinsctr90mask==114 || basinsctr90mask==118 || basinsctr90mask==116 || basinsctr90mask==110, 1, null())" r.mapcalc "chagresbasin=if(basinsctr90mask==102 || basinsctr90mask==100 || basinsctr90mask==96 || basinsctr90mask==98, 1, null())" r.mapcalc "basinsmask=if(basinsctr90mask==102 || basinsctr90mask==100 || basinsctr90mask==96 || basinsctr90mask==98 ||basinsctr90mask==122 || basinsctr90mask==120 || basinsctr90mask==114 || basinsctr90mask==118 ||basinsctr90mask==116 || basinsctr90mask==110 || basinsctr90mask==92 || basinsctr90mask==34, 1, null())" r.mapcalc "basinsmask_inv=if(isnull(basinsmask), 1,null())" r.mapcalc "pequinibasin90m=if(basinsctr90mask==92, 1, null())" r.mapcalc "cuangobasin90m=if(basinsctr90mask==34, 1, null())" r.mapcalc "basinscanal=if(basin200mask==62 || basin200mask==60 || basin200mask==50 || basin200mask==44 ||basin200mask==46 || basin200mask==48 || basin200mask==52 || basin200mask==56 ||basin200mask==54 || basin200mask==366, 1, null())" r.mapcalc "boundaryctr=if(panama_srtm3_90m,1,null())" r.to.vect boundaryctr out=boundaryctr feature=area --o d.rast -o pacorabasin_90m d.rast -o cuangobasin90m d.rast -o pequinibasin90m d.rast -o chagresbasin d.rast -o pacorabasino_90m d.vect srtmstream90med3t d.vect PanamaMaster07_Feb05_GIS icon=basic/circle size=8 fcolor=black d.rast -o basinsmask_inv d.vect boundaryctr type=boundary width=2 d.barscale -mt d.vect PanamaMaster07_Feb05_GIS icon=basic/circle size=8 d.vect PanamaMaster07_Feb05_GIS icon=basic/marker size=16 width=2 d.vect sstreamctr90t d.vect srtmstream90medt d.vect sstreamctr90majt co=black width=2 d.rast pan_srtm3_300msh d.rast lake1_90m -o d.rast lake2 -o d.vect canal width=3 d.vect basinscanal type=boundary width=2 d.vect basin_cuango type=boundary d.vect basin_pacorabig type=boundary d.vect stream90m_canal d.vect stream90m_cuango d.vect stream90m_pacora v.overlay ain=srtmstream90majt bin=basinscanal out=stream90m_canal atype=line operator=and v.overlay ain=srtmstream90majt bin=basin_cuango out=stream90m_cuango atype=line operator=and v.overlay ain=srtmstream90majt bin=basin_pacorabig out=stream90m_pacora atype=line operator=and v.in.ascii feb05pts.txt r.mapcalc "lsfac=exp(fccum90ctr*90.,1.4)*exp(sin(slopectr_srtm3_90m),1.3)" r.mapcalc "lsfac2=exp(accum90ctr,1.4)*exp(sin(slope.30m),1.3)" r.average base=basinsctr90mask cover=slope30m_percc out=basins_slopeperc r.average base=basinsctr90mask cover=lsfac out=basins_lsfac r.average base=basinsctr90mask cover=lsfac2 out=basins_lsfac2 d.his i=srtm90ctr_shade h=basins_lsfac r.mapcalc "basins5ctrBW=if(isnull(pacorabasbigBW),0,1)+if(isnull(chagbasinBW),0,2) \ + if(isnull(pequinibasBW),0,3)+if(isnull(cuangobasBW),0,4)+if(isnull(diosbasinBW),0,5)" 0 white 1 170 170 170 2 110 110 110 3 210 210 210 4 150 150 150 5 130 130 130 551 20:53 d.rast -o pacorabasbigBW 552 20:53 d.rast -o pequinibasBW 553 20:54 d.rast -o chagbasinBW 554 20:54 d.rast -o cuangobasBW 555 20:55 history 556 20:55 r.colors diosbasinBW co=rules